当前位置：文档库 › 河海大学dsp实验3

河海大学dsp实验3

实验三用FFT 对信号进行频谱分析

班级：2012级通信一班

一、实验目的

1、加深对离散信号DFT 的理解；

2、掌握FFT 算法的流程及其MATLAB 实现；

3、利用FFT 对典型信号进行频谱分析；

4、结合理论知识，对频谱分析中出现的有关现象进行理论分析。

二、实验原理

混叠：采样序列的频谱是被采样信号频谱的周期延拓，当采样频率不满足奈奎斯特采样定理的时候，就会发生混叠，使得刺痒后的序列信号的频谱不能真实的反映原采样信号的频谱。

泄露：根据理论分析，一个时间的信号其频带宽度为无限，一个时间无限的信号其频带宽度则为有限。因此对一个时间有限的信号，应用DFT 进行分析，频谱混叠难以避免。对一个时间无限的信号虽然频带有限，但在实际运算中，时间总是取有限值，在将信号截断的过程中，出现了分散的扩展谱线的现象，称之为频谱泄露或功率泄露。

栅栏效应：DFT 是对单位圆上Z 变换的均匀采样，所以它不可能将频谱视为一个连续函数，就在一定意义上看，用DFT 来观察频谱就好象通过一个栅栏来观看一个景象一样，只能在离散点上看到真实的频谱，这样就有可能发生一些频谱的峰点和谷点被“尖桩的栅栏”所挡住，不能被我们观察到。

实验中用到的函数：y=exp(x): 以e 为底的指数。 conj(x): 取x 的共轭，即改变x 的虚部符号。 real(x): 取复数x 的实部。

rand(1,N): 生成在0和1之间均匀分布的随机序列，长度为N 。 randn(1,N): 生成正态分布(高斯分布)的随机序列，长度为N 。

X=fft(x,N): 计算序列x 的N 点FFT 。如果x 的长度小于N ，则在x 后面补零；如果x 的长度大于N ，则对x 进行截取；如果不指定参数N ，则以x 的实际长度作为FFT 的点数。

x=ifft(X,N): 计算序列X 的N 点IFFT 。

Y=fftshift(X): 将序列X 分成左右两部分并交换位置。

三、实验内容

1、对高斯序列x a (n )，令p =8， q 分别等于

2、4和8，观察q 值的改变对高斯序列时

域特性和幅频特性的影响，并给出理论解释：

()

()n p q a e

n x n --??

≤≤=???2

0150

其他

2、对衰减正弦序列x b (n )，a =0.1，f =0.0625，观察其时域和幅频特性，检查谱峰位置是否正确；改变f ，使f 分别等于0.4375和0.5625，比较这两种情况下频谱形状和谱峰位置，

t/T

x (n )

时域特性 p=8,q=2

X (k )

频域特性 p=8,q=2

t/T

x (n )

时域特性 p=8,q=4

X (k )

频域特性 p=8,q=4

t/T

x (n )

时域特性 p=8,q=8

X (k )

频域特性 p=8,q=8

并给出理论解释：

sin()()an b e fn n x n π-?≤≤=?

2015

0其他 3、观察三角波序列x c (n )和反三角波序列x d (n )的时域特性和幅频特性，用8点FFT 分析幅频特性，观察两者的序列形状和频谱曲线有什么异同？为什么？在x c (n )和x d (n )的末尾补零，用32点FFT 分析幅频特性，观察幅频特性发生了什么变化？这些变化说明了什么？

()c n n x n n n ≤≤??

=-≤≤???

038470其他

()d n n x n n n -≤≤??

=-≤≤???

4034470其他

4、对连续单频周期信号sin()a f t π2，按采样频率s a f f =8进行采样，截取长度分别选N =16和N ＝20，比较其幅频特性，并给出理论解释；当截取长度选N ＝164时，与N ＝20比较，频谱发生了哪些变化？为什么？

5、一个连续信号含两个频率分量，经采样得：

()sin(.)cos((.)) , ,,,x n n f n n N ππ=?+?+?=???-2012520125011

令N =16，Δf 分别为1/16和1/64，观察其频谱，两个谱峰是否可以区分开？为什么？当N =128时，Δf 不变，结果有何不同？为什么？

四、实验结果与分析

1.实验结果：

定理的时候，就会发生混叠。x a (n )中的参数p 为高斯序列的峰值位置，q 则表示高斯序列峰的尖锐度，q 值越大，时域图中图象越平缓，序列变化越慢，其幅频特性图中高频分量越少，频谱越窄，在同样的采样频率下越不容易产生混叠。 2.实验结果：

t/T

x (n )

时域特性 a=-0.1,f=0.0625051015

X (k )

频域特性 a=-0.1,f=0.0625

t/T

x (n )

时域特性 a=-0.1,f=0.4375k

X (k )

频域特性 a=-0.1,f=0.4375

t/T

x (n )

时域特性 a=-0.1,f=0.5625k

X (k )

频域特性 a=-0.1,f=0.5625

结论：由图可知：a =0.1，f =0.0625时谱峰位置不正确，此时采样频率f 太小，导致谱峰出现的位置不正确。频谱主瓣较宽，呈现出主瓣中间较为平缓，两侧较高的现象。

当f 分别等于0.4375和0.5625时，函数的幅频特性完全相同，谱峰位置也完全相同，且没有产生明显的混叠和泄露，随采样频率的增加频谱主瓣越来越窄。

3.实验结果：

t/T

x (n )

三角波序列时域特性 k

X (k )

三角波序列频域特性（8点FFT ）

t/T

x (n )

反三角波序列时域特性

X (k )

反三角波序列频域特性（8点FFT ）

X (k )

三角波序列频域特性（32点FFT ）

X (k )

反三角波序列频域特性（32点FFT ）

理论分析：

由图知：当N=8时，正反三角波的幅频特性相同，这是因为进行FFT 运算时，要先对序列进行周期延拓，之后再取主值部分，而两个函数在周期延拓后是相同的，只差一个相位，因此得到的幅频特性图形也是相同的。

当N=32时，两个函数幅频特性不同。这是因为在补零后再，二者的周期延拓就不再相同。栅栏效应，当N=8时，一些谱线被挡住。通过在原序列的末端补零，增加采样的点数，增加了谱线的数量，弱化了栅栏效应。而且反三角波补零后在n=8处会有一个突变，因此它的幅频特性曲线中频谱分量会增多。 4．实验结果：

k X (k )

频域特性 N=16

X (k )

频域特性 N=20

k X (k )

频域特性 N=164

理论分析：①比较N=16与N=20时FFT 的幅频特性，可以发现N=20时频谱出现泄漏，n=16时得到了单一谱线的频谱。这是由于时域中对信号的非整周期截断产生的频谱泄漏。N=20时截取了两个半周期，N=16时截取了两个周期，因此可以看到截取长度N 选得不合理会出现泄漏；

②比较N=20与N=164的幅频特性，可以看出虽然二者都产生了频谱泄漏，但N=164时谱峰要尖锐很多，谱峰两端的小谱峰更多更密集。这是因为在出现频谱泄漏后，随着随着截短长度的增加，截短后信号的频谱越接近原信号频谱的理论值，泄露误差越小。 5．实验结果：

X (k )

幅频特性 N=16，Δf=1/16

X (k )

幅频特性 N=16，Δf=1/64

X (k )

幅频特性 N=128，Δf=1/16

X (k )

幅频特性 N=128，Δf=1/64

理论分析：

观察实验所得频谱，当N=16，Δf=1/16时，两个谱峰可以区分开，这是因为此时的分辨率足够将两谱线分开；当N=16，Δf=1/64时，两个谱峰不能区分开，这是因为此时截取长度短，两谱线靠的太近，1/N的分辨率不能分开他们，出现了栅栏效应。而且由于时域窗太小会有泄漏现象。

当N=128时，两谱线都能分清，这是因为截取长度的提高导致了分辨率的提高，Δf≥1/N时可以分辨，不会发生栅栏效应。而且时域窗增大了，泄漏现象也有所改善。

五、实验小结

通过本次实验，我加深了对离散信号DFT的理解；掌握了FFT算法的流程及其MATLAB 实现；学会了利用FFT对典型信号进行频谱分析，对频谱分析中出现的有关现象及其原因有了进一步的理解，对今后的数字信号处理学习有很大的帮助。

1.程序：

n=0:1:15;

xa1=exp(-((n-8).^2)/2);

subplot(3,2,1);

plot(n,xa1,'-*');

xlabel('t/T');

ylabel('x(n)');

title('时域特性p=8,q=2');

xk1=abs(fft(xa1));

subplot(3,2,2);

stem(n,xk1)

xlabel('k');

ylabel('X(k)');

title('频域特性p=8,q=2');

xa2=exp(-((n-8).^2)/4);

subplot(3,2,3);

plot(n,xa2,'-*');

xlabel('t/T');

ylabel('x(n)');

title('时域特性p=8,q=4');

xk2=abs(fft(xa2));

subplot(3,2,4);

stem(n,xk2)

xlabel('k');

ylabel('X(k)');

title('频域特性p=8,q=4');

xa3=exp(-((n-8).^2)/8);

subplot(3,2,5);

plot(n,xa3,'-*');

xlabel('t/T');

ylabel('x(n)');

title('时域特性p=8,q=8');

xk3=abs(fft(xa3));

subplot(3,2,6);

stem(n,xk3)

xlabel('k');

ylabel('X(k)');

title('频域特性p=8,q=8');

2.程序：

n=0:1:15;

xa1=exp(-0.1*n).*sin(2*pi*0.0625*n); subplot(3,2,1);

plot(n,xa1,'-*');

xlabel('t/T');

ylabel('x(n)');

title('时域特性a=-0.1,f=0.0625');

xk1=abs(fft(xa1));

subplot(3,2,2);

stem(n,xk1)

xlabel('k');

ylabel('X(k)');

title('频域特性a=-0.1,f=0.0625');

xa2=exp(-0.1.*n).*sin(2*pi*0.4375.*n); subplot(3,2,3);

plot(n,xa2,'-*');

xlabel('t/T');

ylabel('x(n)');

title('时域特性a=-0.1,f=0.4375');

xk2=abs(fft(xa2));

subplot(3,2,4);

stem(n,xk2)

xlabel('k');

ylabel('X(k)');

title('频域特性a=-0.1,f=0.4375');

xa3=exp(-0.1.*n).*sin(2*pi*0.5625.*n); subplot(3,2,5);

plot(n,xa3,'-*');

xlabel('t/T');

ylabel('x(n)');

title('时域特性a=-0.1,f=0.5625');

xk3=abs(fft(xa3));

subplot(3,2,6);

stem(n,xk3)

xlabel('k');

ylabel('X(k)');

title('频域特性a=-0.1,f=0.5625');

3.程序：

n=0:1:7;

xa1=[0 1 2 3 4 3 2 1];

subplot(3,2,1);

plot(n,xa1,'-*');

xlabel('t/T');

ylabel('x(n)');

title('三角波序列时域特性');

xk1=abs(fft(xa1));

subplot(3,2,3);

stem(n,xk1)

xlabel('k');

ylabel('X(k)');

title('三角波序列频域特性（8点FFT）');

xa2=[4 3 2 1 0 1 2 3];

subplot(3,2,2);

plot(n,xa2,'-*');

xlabel('t/T');

ylabel('x(n)');

title('反三角波序列时域特性');

xk2=abs(fft(xa2));

subplot(3,2,4);

stem(n,xk2)

xlabel('k');

ylabel('X(k)');

title('反三角波序列频域特性（8点FFT）'); N=0:1:31;

xk3=abs(fft(xa1,32));

subplot(3,2,5);

stem(N,xk3)

xlabel('k');

ylabel('X(k)');

title('反三角波序列频域特性（32点FFT）'); xk4=abs(fft(xa2,32));

subplot(3,2,6);

stem(N,xk4)

xlabel('k');

ylabel('X(k)');

title('反三角波序列频域特性（32点FFT）'); 4．程序：

n1=0:1:15;

xa1=sin(2*pi*n1/8);

xk1=abs(fft(xa1));

subplot(3,1,1);

stem(n1,xk1)

xlabel('k');

ylabel('X(k)');

title('频域特性N=16');

n2=0:1:19;

xa2=sin(2*pi*n2/8);

xk2=abs(fft(xa2));

subplot(3,1,2);

stem(n2,xk2)

xlabel('k');

ylabel('X(k)');

title('频域特性N=20');

n3=0:1:163;

xa3=sin(2*pi*n3/8);

xk3=abs(fft(xa3));

subplot(3,1,3);

stem(n3,xk3)

xlabel('k');

ylabel('X(k)');

title('频域特性N=164');

5.程序：

n=0:1:15;

xa1=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1/16)*n); xk1=abs(fft(xa1));

subplot(2,1,1);

stem(n,xk1)

xlabel('k');

ylabel('X(k)');

title('幅频特性N=16，Δf=1/16');

xa2=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1/64)*n); xk2=abs(fft(xa2));

subplot(2,1,2);

stem(n,xk2)

xlabel('k');

ylabel('X(k)');

title('幅频特性N=16，Δf=1/64');

n=0:1:128;

xa3=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1/16)*n); xk3=abs(fft(xa3));

figure();

subplot(2,1,1);

stem(n,xk3)

xlabel('k');

ylabel('X(k)');

title('幅频特性N=128 ，Δf=1/16');

xa4=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1/64)*n); xk4=abs(fft(xa4));

subplot(2,1,2);

stem(n,xk4)

xlabel('k');

ylabel('X(k)');

title('幅频特性N=128，Δf=1/64');

南昌大学DSP实验报告

实验报告实验课程：DSP原理及应用学生姓名：学号：专业班级： 2012年 5月 25日

目录实验一定点除法运算实验二FIR滤波器实验三FFT算法实验四卷积计算实验五数码管显示实验六语音录放

实验一定点除法运算一、实验目的 1、熟悉C54指令系统，掌握常用汇编指令，学会设计程序和算法的技巧。 2、学习用指令实现除法运算。二、实验设备计算机；DSP 硬件仿真器；DSP 实验开发平台。三、实验原理由内置的硬件模块支持，数字信号处理器可以高速的完成加法和乘法运算。但TMS320 系列DSP不提供除法指令，为实现除法运算，需要编写除法子程序来实现。二进制除法是乘法的逆运算。乘法包括一系列的移位和加法，而除法可分解为一系列的减法和移位。本实验要求编写一个16 位的定点除法子程序。 1．除法运算的过程设累加器为8 位，且除法运算为10 除以3，除的过程包括与除数有关的除数逐步移位，然后进行减法运算，若所得商为正，则在商中置1，否则该位商为0 例如：4 位除法示例：(1)数的最低有效位对齐被除数的最高有效位00001010 - 00011000 11110010 (2)由于减法结果为负，丢弃减法结果，将被除数左移一位再减00010100 - 00011000 11111000 (3)结果仍为负，丢弃减法结果，将被除数左移一位再减00101000 - 00011000 00010000 (4)结果为正，将减法结果左移一位后把商置1，做最后一次减00100001 - 00011000 00001001 (5)结果为正，将减法结果左移一位加1 得最后结果，高4 位是余数，低4 位商：00010011 2．除法运算的实现为了尽量提高除法运算的效率，’C54x 系列提供了条件减指令SUBC 来完成除法操作。四、实验步骤 1．用Simulator 方式启动Code Composer。 2 ．执行Project New 建立新的项目，输入chuf作为项目的名称，将程序定位在D:\ti\myprojects\chuf目录。 3．执行File New Source File 建立新的程序文件，为创建新的程序文件命名为chuf.asm 并保存；执行Project Add Files to Project，把chuf.asm 加入项目中。4．执行File New Source File 建立新的文件并保存为chuf.cmd；执行Project Add Files to Project，把chuf.cmd 加入项目中。 5．编辑chuf.asm 加入如下内容： ;*** 编制计算除法运算的程序段。其中|被除数|<|除数|，商为小数*** .title "chuf.asm" .mmregs .def start,_c_int00

实验三 FIR滤波器的DSP实验报告

实验三FIR滤波器的DSP实验报告一、实验目的 1. 了解MATLAB的FIR滤波器设计方法与编程； 2. 掌握FIR滤波器算法基本原理和C语言的编程方法；二、实验内容 FIR滤波器设计流程为，根据给出的滤波器设计（通带带宽、阻带衰减、相位要求）要求，在Matlab环境下仿真得到滤波器系数，得到了滤波器系数后在CCS下编程设计DSP 的FIR滤波器。 1.使用Matlab中的FDAtool设计FIR滤波器设计实现fir低通滤波器，阶数为20，采样频率fs=8000Hz，截止频率fc=1500Hz，窗口设置为汉明窗Hamming。经过的信号频率是1000HZ和2000HZ的混叠波。 2.参阅FIR基本理论及C语言编程方法，研读、分析实验指导书中的代码； 3. 双击，启动CCS的配置程序选项，选择“C5502 Simulator”； 4. 启动CCS，打开实验工程文件，再编译并装载程序； 5. 仿真波形 A．输入信号波形 B．输入信号频 C输出信号波形

D 输出信号频谱从输入信号和输出信号的频谱对比中可以看出，输出信号有1KHZ 和2KHZ 的混叠信号。经过FIR 低通录波器后，滤掉了2KHZ 的频谱。只剩下1KHZ 的频谱。 4.修改代码：滤波器设计要求：线性相位低通滤波器，模拟信号通带为0~1kHz ，阻带衰减在5kHz 处为-60dB 。采样率为20kHz 。采用凯塞窗口设计此低通滤波器，设计过程如下：通带边界频率为1KHZ ，阻带边界频率是5KHZ 有60dB 的衰减，可得： rad rad s p 22052,102012ππωππ ω==== 过渡带πωωω5 2 = -=?p s ，A=60dB 设计凯塞窗口的长度N 和参数β分别为： 29285.28 =?-≈ ω A N 65326.5)7.8(1102.0=-=a β 使用MA TLAB 设计此滤波器，过程如下： 1）在Matlab 的Start 菜单中选择Toolboxes -> Filter Design -> Filter Design & Analysis Tools(fdatool)，或者在命令行中输入fdatool 来启动滤波器设计分析器。

河海大学2005-2006学年第二学期水力学期末试卷(A卷)

2005-2006学年第二学期水力学期末试卷(A卷) （水文、给排水、环工、资环04级、资环03级）班级学号姓名成绩一．是非题（每小题1分，共10分） 1．水流在边壁处的流速为零，因此该处的流速梯度为零。 ( ) 2．水力坡度就是单位长度流程上的水头损失。 ( ) 3．在恒定紊流中时均流速不随时间变化。 ( ) 4．不论平面在静止液体内如何放置，其静水总压力的作用点永远在平面形心之下。（）5．同一种管径和粗糙度的管道，雷诺数不同时，可以在管中形成紊流光滑区、粗糙区或过渡粗糙区。（）6．恒定总流的总水头线沿流程下降,而测压管水头线沿程可升可降。 ( ) 7．只要下游水位高于堰顶，一定是堰的淹没出流。 ( ) 8．棱柱形明渠中形成M3 型水面曲线时，其断面单位能量E s 沿程增大。( ) 9．如果两个液流中作用于相应点上的相同性质的力的大小成一定比例，则两个液流就是动力相似的。（）10．渐变无压渗流中任意过水断面各点的渗流流速相等，且等于断面平均流速。 ( ) 二．单项选择题（每小题2分，共20分） 1．已知水流某处的绝对压强值为2.94N/cm 2，当地大气压强为10 米水柱高，则该处真空度是（） (1) 3ｍ (2) 7ｍ (3) 2.94ｍ (4) 没有真空 2．在恒定流中 ( ) (1)流线一定互相平行 (2)断面平均流速必定沿程不变 (3)不同瞬时流线有可能相交 (4)同一点处不同时刻的动水压强相等 3．紊流过渡粗糙区的沿程水头损失系数λ ( ) (1) 只与雷诺数有关 (2) 只与相对粗糙度有关 (3) 只与绝对粗糙度有关 (4) 与相对粗糙度和雷诺数有关 4．按普朗特动量传递理论，紊流的断面流速分布规律符合() (1)对数分布(2)椭圆分布(3)抛物线分布(4)直线分布。 5．陡坡渠道上如发生非均匀流 ( ) (1) 一定是急流 (2) 一定是缓流 (3) 可以是缓流或急流 (4) 一定是渐变流 6．矩形断面渠道的均匀流，若断面单位能量与水深之比E s/h =2。则水流为（） (1)急流 (2)缓流 (3)临界流 (4) 急变流共5页第1页

DSP实验3

实验三快速傅里叶变换 FFT 一、实验目的 1．理解DFT 算法，并能用MATLAB 实现DFT 。 2. 加深对FFT 的理解，体会DFT 和FFT 之间的关系。 3．熟悉应用FFT 实现两个序列的线性卷积的方法。二、实验原理 N 点序列x(n) 的DFT 和IDFT 定义：若将DFT 变换的定义写成矩阵形式，则得到 X=A ﹒x ，其中DFT 变换矩阵A 为 ?????????? ????=---2)1(111...1......... ......11...11N N N N N N N W W W W A 可以用函数U=fft(u,N)和u=ifft(U,N)计算N 点序列的DFT 正、反变换。三、实验内容（一）离散傅里叶变换（DFT ） 1.用MATLAB 求N=16的有限序列)4/sin()8/sin()(ππn n n x +=的DFT 结果，并画出结果图。参考程序如下： N=16; n=0:1:N-1; %时域采样 xn=sin(n*pi/8)+sin(n*pi/4); k=0:1:N-1; %频域采样

WN=exp(-j*2*pi/N); nk=n'*k; WNnk=WN.^nk; Xk=xn*WNnk; subplot(2,1,1) stem(n,xn); subplot(2,1,2) stem(k,abs(Xk)); 2.矩形序列x(n)=R (n),求N分别取8，32时的DFT，最后绘出结果图形。 5 参考程序如下： function[Xk]=dft(xn,N) n=[0:1:N-1]; %n的行向量 k=[0:1:N-1]; %k的行向量 WN=exp(-j*2*pi/N); %旋转因子 nk=n'*k; %产生一个含nk值的N乘N维矩阵WNnk=WN.^nk; %DFT矩阵 Xk=xn*WNnk; %DFT系数的行向量调用上面函数解题。 N=8;x=[ones(1,5),zeros(1,N-5)]; n=0:N-1; X=dft(x,N); %N＝8点离散傅立叶变换 magX=abs(X);phaX=angle(X)*180/pi; k=(0:length(magX)'-1)*N/length(magX); subplot(2,2,1);stem(n,x);ylabel('x(n)'); subplot(2,2,2);stem(k,magX);axis([0,10,0,5]);ylabel('|X(k)|'); N=32;x=[ones(1,5),zeros(1,N-5)]; n=0:N-1; X=dft(x,N); %N＝32点离散傅立叶变换 magX=abs(X);phaX=angle(X)*180/pi; k=(0:length(magX)'-1)*N/length(magX);

DSP实验三和实验四

实验三离散时间系统及响应 [实验目的] 1、熟悉离散时间系统的时域和频域分析方法。 2、掌握稳定系统性的判别。 3、掌握利用MATLAB求解零极状态、差分方程和频率响应的方法。 [实验仪器] 计算机、Matlab6.5（或更高版本）软件 [实验参考书] ：自编实验指导书与本实验有关的MA TLAB函数： 1. conv.m用来实现两个离散序列的线性卷积。其调用格式是：y=conv(x,h) 2．filter.m求离散系统的输出y(n) 。若系统的h(n) 已知，可用conv.m文件可求出y(n) ；若系统的H(z) 已知，可用filter可求出y(n)，调用格式是: y=filter(b, a, x)；%其中x, y, a 和b都是向量。 3．impz.m在H(z) 已知情况下, 求系统的单位抽样响应h(n)。调用格式是： h = impz(b, a, N) 或 [h,t]=impz(b,a,N) N是所需的的长度。前者绘图时n从1开始，而后者从0开始。 4．freqz.m在H(z) 已知情况下, 求系统的频率响应。基本的调用格式是： [H,w]=freqz(b,a,N,'whole',Fs) N是频率轴的分点数，建议N为2的整次幂；w是返回频率轴座标向量，绘图用；Fs是抽样频率，若Fs＝1，频率轴给出归一化频率；’whole’指定计算的频率范围是从0～FS,缺省时是从0～FS/2. 5.zplane.m文件可用来显示离散系统的极－零图。其调用格式是： zplane(z,p), 或zplane(b,a), 前者是在已知系统零点的列向量z和极点的列向量p的情况下画出极－零图，后者是在仅已知H(Z)的A(z)、B(z) 的情况下画出极－零图。 [实验内容] 1、序列的卷积运算。下例是一个指数函数的卷积，结果如图所示。 clear all; fs=20; n =0:49; x=exp(-n/fs); subplot(211),stem(x) y=conv(x,x); subplot(212);stem(y);grid; 05101520253035404550 仿照上例，现在假设存在两个序列，()[,,,,,,] h n=54321,利用计 x n=9121041215()[,,,,]

DSP运行实验报告

DSP运行实验报告一、实验目的熟悉CCS软件仿真下，DSP程序的下载和运行；熟悉借助单片机的DSP程序下载和运行；熟悉借助仿真器的DSP程序下载和运行；熟悉与DSP程序下载运行相关的CCS编程环境。二、实验原理 CCS软件仿真下，借用计算机的资源仿真DSP的内部结构，可以模拟DSP程序的下载和运行。如果要让程序在实验板的DSP中运行、调试和仿真，可以用仿真器进行DSP程序下载和运行。初学者也可以不用仿真器来使用这款实验板，只是不能进行程序调试和仿真。在本实验板的作用中，单片机既是串口下载程序的载体，又是充当DSP 的片外存储器（相对于FLASH），用于固化程序。三、实验设备、仪器及材料安装有WINDOWS XP操作系统和CCS3.3的计算机。四、实验步骤（按照实际操作过程） 1、CCS软件仿真下，DSP程序的下载和运行。第一步：安装CCS，如果不使用仿真器，CCS 的运行环境要设置成一个模拟仿真器（软仿真）。

第二步：运行CCS，进入CCS 开发环境。第三步：打开一个工程。将实验目录下的EXP01目录拷到D:\shiyan下（目录路径不能有中文），用[Project]\[Open]菜单打开工程，在“Project Open”对话框中选 EXP01\CPUtimer\CpuTimer.pjt，选“打开”，第四步：编译工程。在[Project]菜单中选“Rebuild All”，生成CpuTimer.out文件。第五步：装载程序。用[File]\[Load Program]菜单装载第四步生成CpuTimer.out文件，在当前工程目录中的Debug 文件夹中找到CpuTimer.out文件，选中，鼠标左键单击“打开”。

DSP实验二要点

实验三 IIR 滤波器设计一、实验目的： 1.认真复习滤波器幅度平方函数的特性，模拟低通滤波器的巴特沃思逼近、切比雪夫型逼近方法；复习从模拟低通到模拟高通、带通、带阻的频率变换法；从模拟滤波器到数字滤波器的脉冲响应不变法、双线性变换法的基本概念、基本理论和基本方法。 2掌握巴特沃思、切比雪夫模拟低通滤波器的设计方法；利用模拟域频率变换设计模拟高通、带通、带阻滤波器的方法.。 3.掌握利用脉冲响应不变法、双线性变换法设计数字滤波器的基本方法；能熟练设计巴特沃思、切比雪夫低通、带通、高通、带阻数字滤波器。 4.熟悉利用MATLAB 直接进行各类数字滤波器的设计方法。二、实验内容 a. 设计模拟低通滤波器，通带截止频率为10KHz，阻带截止频率为16KHz，通带最大衰减1dB，阻带最小衰减20dB。 (1) 分别用巴特沃思、切比雪夫I、切比雪夫II 型、椭圆型滤波器分别进行设计，并绘制所设计滤波器的幅频和相频特性图。 (2) 在通带截止频率不变的情况下，分别用n=3,4,5,6 阶贝塞尔滤波器设计所需的低通滤波器，并绘制其相应的幅频响应和相频响应图。 %%%%%%%%%----巴特沃思-----%%%%%%% clc;clear all; omegap=10000*2*pi;omegas=16*10^3*2*pi; Rp=1;As=20; [N,omegac]=buttord(omegap,omegas,Rp,As,'s');%低通的节次 [b,a]=butter(N,omegac,'s'); [H,w]=freqs(b,a); %设计滤波器的幅频和相频特性图 subplot(211) plot(w/2*pi/1000,20*log10(abs(H)))

DSP实验指导书(DSP处理器原理与应用)

实验一：DSP软件集成开发环境CCS 一. 实验目的 1.了解DSP 软件开发集成环境Code Composer Studio (CCS) 的安装和配置过程； 2.熟悉并掌握CCS 的操作环境和基本功能，了解DSP软件开发的一般流程： (1) 学习创建工程和管理工程的方法; (2) 了解工程文件和映像文件的基本内容; (3) 了解基本的编译和调试功能; (4) 学习使用观察窗口。二. 实验设备计算机、Code Composer Studio 2.2 或以上版本三. 背景知识 1.DSP 应用系统开发一般需要借助以下调试工具： ―软件集成开发环境(CCS)：完成系统的软件开发，进行软件和硬件仿真调试； ―硬件开发及测试平台：实现系统的硬件仿真与调试，包括仿真器和评估模块。 https://www.wendangku.net/doc/b52559322.html,S 集成了适用于所有TI 器件的源码编辑、编译链接、代码性能评估、调试仿真等功能，为设计人员提供了一个完整的嵌入式系统软件开发环境。 3.系统软件部分可以通过CCS 建立的工程文件进行管理，工程文件一般包含以下几种文件： ―源程序文件：C 语言或汇编语言文件(*.ASM 或*.C) ―头文件(*.H) ―命令文件(*.CMD) ―库文件(*.LIB, *.OBJ)

四. 实验步骤 1．创建工程文件双击 CCStudio 图标，启动 CCS ，CCS 的初始界面包括工程管理区和工作区两部分。在菜单中选择“Project —>New…”，弹出“Project Create ”窗口：在“Project ”编辑框内填入项目名称，例如“hello1”；有必要的话，可以更改项目文件夹位置 (Location)；“Project ”下拉框用以确定输出文件类型；目标硬件类型在“Target ”下拉框中指定。点击“完成”后，CCS 自动在指定目录下生成“hello1.pjt ”工程文件。工程文件中保存了对应工程的设置信息和内部文件的引用情况。展开工程管理窗口“Projects ”下的“hello1.pjt ”，可以看到所有项目均为空。任务 1：在 notepad 中打开“hello1.pjt ”，了解并记录工程文件提供的设置信息若标题栏显示xxxx Emulator, 则需要先运行 Setup CCS, 选择合适的软件仿真环境

河海大学水力学2003年至2007年真题及答案

河海大学 2003 年攻读硕士学位研究生入学考试试题考试科目名称：水力学一、是非题（每题 2分，共 20 分） 1．紊流过渡粗糙区的沿程水头损失系数与雷诺数。（） 2．断面单位能量 E s 沿流程总是减小的。（） 3．计算消力池池深和池长时，设计流量一般选择使池深和池长最大的流量。（） 4．渗流的达西定律既适用于层流也适合于紊流。（） 5．在恒定均匀流中，沿程水头损失与速度的平方成正比。（） 6．两个液流在重力作用下的动力相似条件是它们的弗劳德数 F r 相等。（） 7．因为各并联支管的长度、直径及粗糙度可能不同，因此在各支管的水头损失也不同。（） 8．均匀流可以是恒定流，也可以是非恒定流。（） 9．流体质点做圆周运动不一定是有旋运动。（） 10．相对压强可以大于、等于或小于零。（）二、选择题（每题 2分，共 20 分） 1．当水流流态从层流转变为紊流后，过水断面上水流动能校正系数和动量校正系数将 2．溢流坝下泄水流收缩断面水深的跃后水深为，此时下游水深为，则形成远离水跃的条件是（）（1） h t h co2 （ 2） h t h co2 （3） h t h co2 （ 4）无法确定 3．渗流模型流速与实际渗流流速相比较（）（1）前者大于后者（ 2）二者相等（3）前者小于后者（ 4）无法比较 4．明渠的临界水深决定于（）（1）流量和底坡（ 2）断面形状、尺寸和底坡 1）均增大 3）增大、减小 2）均减小 4）减小、增大

3）流量和糙率（4）流量和断面形状、尺寸 5．均匀流断面和渐变流断面上的动水压强 1）均按静水压强规律分布 2）前者按静水压强规律分布，后者不静水压强规律分布 3）前者不静水压强规律分布，后者静水压强规律分布 4）均不按静水压强规律分布 6．下列哪个条件满足时，明渠水流是缓流 1) F r 12) F r 1 3) F r 14）无法确定 7．毕托管可以用来测 1）脉动流速（2）时均流速（3）脉动流速（4）脉动压强 8．满足dE s/ ds 0条件的流动是 1）非均匀渐变缓流非均匀渐变急流 3）均匀流4）临界流 9．有压管流，管轴线水平，管径逐渐增大且通过的流量不变，其总水头线沿流向应（1）逐渐升高（2）逐渐降低（3）与管轴线平行（4）无法确定 10．雷诺数的物理意义是 1）重力与惯性力之比2）重力与粘滞力之比 3）粘滞力与惯性力之比4）惯性力与粘滞力之比三、作图题（每题 5 分，共10 分） C 1．绘出下列两段不同糙率棱柱形渠道的水面曲线。（渠道每段都充分长） S2 i>i c ( > 2)

河海大学814水力学考研高分学长经验分享

河海大学水力学814高分学长考研经验水力学题量虽然大，其实只要理解透彻并不难的，静下心来仔仔细细的学明白教材，毕竟教材是根本。在熟悉书中知识结构的脉络的基础上理解知识点。做一些练习题巩固巩固，水力学考的是基础和细心，基础要扎实。说下我的复习过程：我是从5月份开始看书，那时候边看书边做课后习题，大概花了我两个月的时间，整理一些概念性的知识点，然后暑期参加了南京胜研考研的二期强化集训营，集训营上是水力学考研高分学长（专业课第一名，总分第二名）给我们把整个教材结合考试题内容深入讲了一遍，说实话讲的很透彻，总的来说感觉自己的知识理解上和熟练程度上比自己看的那些时间得到了很大提高，有很多讲的考试重难点和做题方法经验什么的确实自己悟不出来。完了我从开学的时候开始自己独立结合集训营老师讲的考试重点看第二遍，这一遍就是重点整理重点的知识点，大题常考的点，做一些相关的练习题，适当的结合真题，看书的时候碰到真题里的题目重点标注，重点整理，重点练习，我当时是用16开的白纸整理的，这一遍可能会花费的时间长些，做笔记是一个学习的过程，这遍下来我感觉自己当时的知识体系应该算建立起来了，我当时用了1个半月的时间，第三遍我就开始系统的做些练习和系统的做真题了，在做题的过程中进一步巩固知识点，做到拿起题目就知道它重点考的什么，反复的练习，做题的过程中可能会有些遗忘，多做多看做到温故而知新，然后就是第四遍，完全把课本看成一个整体，如果你的时间充裕，你可以拿着真题像考试那样做按时间做，这样效果更好，然后就可以上考场了。我总分三百八十多分，专业课挺好的，英语政治底子薄拉分了，要不能考的更高，总分算中等水平吧，也算考上了河海的好专业，希望你们也能参照我的学习方法，好好复习。

DSP实验

一：实验项目名称：TMS320F2812的GPIO实验二：实验目的：学会使用ccs5.3开发环境以及GPIO的基本操作，通过实验学习使用 TMS320F2812的通用输入输出管脚直接控制外围设备的方法。实验环境：pc一台，CCS5.3软件三：实验原理：四：实验步骤： 1：打开ccs5.3软件。 2依次单击File->New->CCS Project建立工程 3 在弹出的New Project对话框中单击Code Composer Studio，选择CCS Project，然后点击 Next 4 在弹出的对话框中，输入工程名字Led.c，在Family里选择C2000，在Variant的第二个复选框里选择TMS320F2812，在里选择empty project点击 Finish完成 5 将inc、source和F2812_nonBIOS_flash.cmd、rts2800_ml.lib拷贝到工程目录下； 6 在project面板里右键，New->Target Configuration File,在connect里面选择，然后在选择框中选择 TMS320F2812； 7 在main.c中编辑代码 8 连接好仿真器与开发板，按F11，下载并在线仿真五：记录与处理： #include"DSP281x_Device.h" #include"System.h" void main(void) { InitSysCtrl(); EALLOW; GpioMuxRegs.GPBDIR.all |=0x00FF; EDIS; GpioDataRegs.GPBDAT.all = 0x0000; DelayUs(60000); GpioDataRegs.GPBDAT.all = 0x00f0; while(1) {

河海大学水力学2019年至2019年考研真题试卷共26页word资料

河海大学2003年攻读硕士学位研究生入学考试试题考试科目名称：水力学一、是非题（每题2分，共20分） 1．紊流过渡粗糙区的沿程水头损失系数与雷诺数。（） 2．断面单位能量 E沿流程总是减小的。 s （） 3．计算消力池池深和池长时，设计流量一般选择使池深和池长最大的流量。（） 4．渗流的达西定律既适用于层流也适合于紊流。（） 5．在恒定均匀流中，沿程水头损失与速度的平方成正比。（） 6．两个液流在重力作用下的动力相似条件是它们的弗劳德数 F相等。 r （） 7．因为各并联支管的长度、直径及粗糙度可能不同，因此在各支管的水头损失也不同。（） 8．均匀流可以是恒定流，也可以是非恒定流。（） 9．流体质点做圆周运动不一定是有旋运动。（） 10．相对压强可以大于、等于或小于零。

（）二、选择题（每题2分，共20分） 1．当水流流态从层流转变为紊流后，过水断面上水流动能校正系数α和动量校正系数β将（1）均增大（2）均减小（3）α增大、β减小（4）α减小、β增大 2．溢流坝下泄水流收缩断面水深的跃后水深为，此时下游水深为，则形成远离水跃的条件是（）（1）2t co h h = （2）2t co h h > （3）2t co h h < （4）无法确定 3．渗流模型流速与实际渗流流速相比较（）（1）前者大于后者（2）二者相等（3）前者小于后者（4）无法比较 4．明渠的临界水深决定于（）（1）流量和底坡（2）断面形状、尺寸和底坡（3）流量和糙率（4）流量和断面形状、尺寸 5．均匀流断面和渐变流断面上的动水压强（）（1）均按静水压强规律分布

2008河海大学水力学试题及答案

河海大学水力学试卷水工港行农水一\选择题部分（1）在水力学中，单位质量力是指（） a、单位面积液体受到的质量力； b、单位体积液体受到的质量力； c、单位质量液体受到的质量力； d、单位重量液体受到的质量力。答案：c （2）在平衡液体中，质量力与等压面（） a、重合； b、平行 c、相交； d、正交。答案：d （3）液体中某点的绝对压强为100kN/m2，则该点的相对压强为 a、1 kN/m2 b、2 kN/m2 c、5 kN/m2 d、10 kN/m2 答案：b （4）水力学中的一维流动是指（） a、恒定流动； b、均匀流动； c、层流运动； d、运动要素只与一个坐标有关的流动。答案：d （5）有压管道的管径d与管流水力半径的比值d /R=（） a、8； b、4； c、2； d、1。答案：b （6）已知液体流动的沿程水力摩擦系数与边壁相对粗糙度和雷诺数Re都有关，即可以判断该液体流动属于 a、层流区； b、紊流光滑区； c、紊流过渡粗糙区； d、紊流粗糙区答案：c （7）突然完全关闭管道末端的阀门，产生直接水击。已知水击波速c=1000m/s，水击压强水头H = 250m，则管道中原来的流速v0为 a、1.54m b 、2.0m c 、2.45m d、3.22m 答案：c （8）在明渠中不可以发生的流动是（） a、恒定均匀流； b、恒定非均匀流； c、非恒定均匀流； d、非恒定非均匀流。答案：c （9）在缓坡明渠中不可以发生的流动是（）。 a、均匀缓流； b、均匀急流； c、非均匀缓流； d、非均匀急流。答案：b （10）底宽b=1.5m的矩形明渠，通过的流量Q =1.5m3/s，已知渠中某处水深h = 0.4m，则该处水流的流态为

DSP实验报告 (20)

汕头大学工学院三级项目报告课程名称：数字信号处理课程设计题目：FIR数字滤波器指导教师：姜永权系别：电子工程系专业：电子信息工程学号：08141017 姓名：邓桂杨合作者完成时间：2010 年11 月20 日至11 月26 日成绩：评阅人：姜永权

一项目意义与目标 FIR数字滤波器是一类广泛应用的数字滤波器。如自适应横向数字滤波器就采用FIR 数字滤波器，在GSM手机中就用到。在实际DSP系统中能够正确运用FIR数字滤波器，理解线性相位的含义和应满足的条件。通过恰当举例，通过亲自编程，给出设计结果，绘制频率响应曲线，掌握FIR数字滤波器设计算法。二项目内容 1.根据频率选择性数字滤波器的性能指标，计算线性相位FIR数字滤波器的理想单位冲激响应。选择Kaiser窗函数的阶数及形状控制参数，给出频率选择性线性相位FIR数字滤波器的设计结果，画出FIR数字滤波器的频率响应曲线。 2.根据任意频率响应（均衡器）性能要求，计算线性相位FIR数字均衡器的单位冲激响应。选择Kaiser窗函数的阶数及形状控制参数，给出均衡器的设计结果，画出均衡器的频率响应曲线。FIR数字滤波器、均衡器的性能分析。三项目报告正文 1.FIR滤波器的设计方法目前，FIR滤波器的设计方法主要是建立在对理想滤波器频率特征做某种近似的逼近。窗函数设计法是FIR滤波器设计的基本方法，它设计过程简单，性能也可满足常用的选频滤波器的技术要求。窗函数设计的基本思路是直接从理想滤波器的频率特性入手，通过积分求出对应的单位采样响应的表达式，最后通过加窗，得到满足要求的FIR滤波器的采样响应。窗函数在很大程度上决定了FIR滤波器的性能指标。

河海大学814水力学

814水力学 Gungi Peng 2016.11.22 一、绪论 1.牛顿内摩擦定律τ=μdu dy （du dy 为角变形率或流速梯度， μ为动力粘度系数=ρv ）（注意各物理量单位） 2.牛顿内摩擦定律的适用条件为牛顿液体层流运动 3.F =τA 。（F 为相邻层之间的内摩擦力） 4.圆盘F = τ2πrdr ；M = τ2πrrdr 5.宾汉塑性流体切应力与切应变率也成线性关系，但不过原点，牛顿流体过原点。（p7图） 6.表面力：压力、切力。质量力：重力、惯性力。 7.液体在流动(连续不断变形)过程中，其内部会出现抵抗，液体在流动状态下抵抗剪切的性质称为液体的粘性（又称粘滞性）二、水静力学 1.静水压强：垂直指向作用面，一点压强各方向数值相等。 2.平衡微分方程：dP =ρ（xdx +ydy +zdz ），当dP=0，即等压面，质量力与等压面垂直。 3.液体处于相对平衡：1.等压面水平，只有重力作用；2.等压面为斜面，重力和惯性力（产生水平加速度）作用；3.等压面为旋转抛物面，重力和惯性力（离心力）作用 4.标准大气压为101.4KN/m 2，绝对压强等于标准压强加相对压强 5.真空高度与真空压强的换算，P =ρgh ，真空压强=∣负压∣ 6.静止液体基本方程：z + P ρg =C ，静止液体中任一点相等。 7.解析法：平面上压力计算：P=AP c ，即受压平面形心压强。与放置方向无关，只与形心位置有关。yd =yc +Ic ycA 。矩形：Ic =1 12bh 3，圆：Ic =1 4πr4

8.压力图法：P=Ωb，压强分布图面积。压强分布：三角形e=1/3h，梯形：e=a 32 1+ 2 1+ 2 9.曲面：Px=PxcAx；Pz=ρgV，p=2+Pz2，θ=arctan Pz Px 10.虚压力体：水面和压力体在两侧，实压力体在同侧三、液体一元恒定总流基本原理 1.恒定流与非恒定流：液体流动时空间各点处的所有运动要素都不随时间而变化的流动称为恒定流，反之为非恒定流（与t是否有关，有无当地加速度） 2.均匀流与非均匀流：流场中位于同一流线上各质点的流速矢量不沿程变化（流线不相交）为均匀流，反之为非均匀流。非均匀流可以分为渐变流和急变流（有无迁移加速度） 3.元流：微小流管中的液流 4.总流：由无数元流集合而成的整股水流称为总流 5.一元流与三元流：液体的运动要素仅是一个坐标变量的函数，这种运动称为一元流，运动要素是三个坐标变量的函数，称三元流 6.流线：某一瞬间在流场中绘出的曲线，在此曲线上所有液体质点的速度矢量都和该曲线相切。 7.迹线：同一质点在一个时段的运动轨迹线 8.连续方程：V1A1=V2A2 9.能量方程：z1+P1 ρg +αv12 2g + P=z2+P2 ρg +αv22 2g + w；hw=hf（沿程）+hj（局部） 10.动量方程：P1?P2+R=ρQβ2v2?β1v1；R为固体边壁对水流的作用力 11.动能/动量校正系数，越不均匀越大，因此层流远大于紊流 12.同一管流断面上，动能校正系数α比动量校正系数β大 13.过水断面上的压强大小与基准面选择无关 14.均匀流断面上（非“流场内”）的压强分布规律与静水压强规律相同

河海大学水力学2019年至2019年考研真题试卷word资料19页

河海大学2019年攻读硕士学位研究生入学考试试题考试科目名称：水力学一、是非题（每题2分，共20分） 1．紊流过渡粗糙区的沿程水头损失系数λ与雷诺数。（） 2．断面单位能量s E 沿流程总是减小的。（） 3．计算消力池池深和池长时，设计流量一般选择使池深和池长最大的流量。（） 4．渗流的达西定律既适用于层流也适合于紊流。（） 5．在恒定均匀流中，沿程水头损失与速度的平方成正比。（） 6．两个液流在重力作用下的动力相似条件是它们的弗劳德数r F 相等。（） 7．因为各并联支管的长度、直径及粗糙度可能不同，因此在各支管的水头损失也不同。（） 8．均匀流可以是恒定流，也可以是非恒定流。（） 9．流体质点做圆周运动不一定是有旋运动。（） 10．相对压强可以大于、等于或小于零。（）二、选择题（每题2分，共20分） 1．当水流流态从层流转变为紊流后，过水断面上水流动能校正系数α和动量校正系数β将（1）均增大（2）均减小（3）α增大、β减小（4）α减小、β增大 2．溢流坝下泄水流收缩断面水深的跃后水深为，此时下游水深为，则形成远离水跃的条件是（）（1）2t co h h = （2）2t co h h > （3）2t co h h < （4）无法确定 3．渗流模型流速与实际渗流流速相比较（）（1）前者大于后者（2）二者相等（3）前者小于后者（4）无法比较 4．明渠的临界水深决定于（）（1）流量和底坡（2）断面形状、尺寸和底坡（3）流量和糙率（4）流量和断面形状、尺寸

5．均匀流断面和渐变流断面上的动水压强（）（1）均按静水压强规律分布（2）前者按静水压强规律分布，后者不静水压强规律分布（3）前者不静水压强规律分布，后者静水压强规律分布（4）均不按静水压强规律分布 6．下列哪个条件满足时，明渠水流是缓流（）（1）1r F > （2）1r F = （3）1r F < （4）无法确定 7．毕托管可以用来测（）（1）脉动流速（2）时均流速（3）脉动流速（4）脉动压强 8．满足/0s dE ds =条件的流动是（）（1）非均匀渐变缓流（2）非均匀渐变急流（3）均匀流（4）临界流 9．有压管流，管轴线水平，管径逐渐增大且通过的流量不变，其总水头线沿流向应（）（1）逐渐升高（2）逐渐降低（3）与管轴线平行（4）无法确定 10．雷诺数的物理意义是（）（1）重力与惯性力之比（2）重力与粘滞力之比（3）粘滞力与惯性力之比（4）惯性力与粘滞力之比三、作图题（每题5分，共10分） 1．绘出下列两段不同糙率棱柱形渠道的水面曲线。（渠道每段都充分长） 2．绘出下列二向曲面的水平压强分布图和压力体图。四、计算题（本题6小题，共100分） 1．粘度测量仪由两个圆筒组成，两筒的间隙充满某种液体。外筒与转轴连接，其半径为2r ，旋转角速度为ω。内筒悬挂于一扭丝上，其半径为1r ，扭丝所受的力矩M 可以测出。外筒与内筒底面间隙为a ，内筒高为H ，如图所示。试推出该液体的动力粘滞系数μ的计算式。 2．有一锥形容器如图所示，在容器的底部A 处接一U 形水银测压计，当容器没有水时，在A 以下注水后测压计压差读书为19cm ，当容器充满水后，测压计的压差读书又为多少？ 3．一泄洪管道直径4d m =，管长2l km =，上游水箱水面管道中心的水头25H m =，沿

DSP实验报告3

DSP 实验报告3 姓名学号时间 2015.6.9 地点实验大楼A631 实验题目卷积运算一、实验目的 1、学习如何用DSP 的C 语言程序实现卷积与相关的计算； 2、学习用CCS 调试程序的详细过程； 3、学习如何设置断点。二、实验内容及步骤 1、创建一个新项目工程将“maxminmath 文件夹”复制到“D:\CCStudio_v3.3\MyProjects ”下。 2、打开工程调试程序，并说明程序的功能； 3、编程实现以下功能：已知()()n R n x 4=，()()n R n h 4=，求()()()n h n x n y *= a 、采用子程序调用的方式实现，并且注释程序； b 、把x ，h 和y 添加到Watch 窗口中作为观察对象； c 、观察三个数组从初始化到卷积运算结束整个过程中的变化（可单击变量名前的“+”号把数组展开），并将变化过程记录下来。三、实验程序给出卷积运算的子程序： void main() { int x[7]={0,0,0,1,1,1,1},h[7]={1,1,1,1,0,0,0},y[7]={0}; int i,j; for(i=0;i<7;i++) for(j=0;j<7;j++) {if(i+j<7) y[i]=y[i]+x[i+j]*h[j]; y[i]=y[i]%2;} } 评阅

四、实验结果 y值的变化情况如下：五、实验思考题 1、试用汇编语言编写实现卷积运算的程序； .global start ;定义全局变量 .mmregs .data; 开始数据段 .bss m,1; 开设全局变量，每个变量占用一个字 .bss n,1 .bss z,1 .text; 开始文本段 start: ST #1000h,AR1 ST #3000h,AR2 ST #5000h,AR3 ST #0，@m RPT 999 RPT m+ LOOP: MVPD 6000h MVDD *AR1+,*AR4 MVDD *AR2+,*AR5 MACP *AR4,*AR5,src B LOOP ST A,*AR3+ xh: b xh;空循环 .end

DSP matlab实验三

实验三：IIR 数字滤波器的设计和实现学院：专业：班级：姓名：学号：班内序号：时间：【实验背景】数字信号处理的另一个重要应用是数字滤波器。数字滤波器是一个运算过程，将输入数列按既定的要求转换成输出数列。在数字信号的处理时只需利用数字相加、乘以常数和延时等运算，就可以完全达到传输特性的要求。数字滤波器分为无限冲激响应（IIR）和有限冲激响应（FIR）两大类。本实验用IIR 数字滤波器产生双音多频DTMF 拨号信号的产生，演示IIR 数字滤波器在通信系统中的应用。 DTMF 信号是将拨号盘上的0~F 共16 个数字，用音频范围的8 个频率来表示的一种编码方式。8 个频率分为高频群和低频群两组，分别作为列频和行频。每个字符的信号由来自列频和行频的两个频率的正弦信号叠加而成。频率组合方式如下图所示。【实验要求】根据ITU Q.23 建议，DTMF信号的技术指标是：传送/接收率为每秒10个号码，或每个号码100ms。每个号码传送过程中，信号存在时间至少45ms，且不多于55ms，100ms的其余时间是静音。在每个频率点上允许有不超过±1.5%的频率误差。任何超过给定频率±3.5%的信号，均被认为是无效的，拒绝承认接收。另

外，在最坏的检测条件下，信噪比不得低于15dB。DTMF 信号的编码: 把您的联系电话号码DTMF 编码生成为一个.wav 文件。其中关键是不同频率的正弦波的产生，要求采用滤波法生成所要求的DTMF信号。【实验内容】 1.实验代码： phone=input('请输入电话号码:','s'); >> sum=length(phone); fs=8000;total_x=[]; fH=[1336,1209,1336,1477,1209,1336,1477,1209,1336,1477]; fL=[941,697,697,697,770,770,770,852,852,852]; for i=1:sum num=phone(i)-48+1;n=1:400; x=sin(2*pi*n*fL(num)/fs)+sin(2*pi*n*fH(num)/fs); x=[x,zeros(1,400)]; total_x=[total_x,x];end audiowrite('num_sound.wav',total_x,8000) 警告: 数据在写入文件期间被裁剪。 > In audiowrite>clipInputData (line 396) In audiowrite (line 176) >> sound(total_x); >> plot(total_x);axis([0,10000,-2,2]);title('时域波形') >> k=[18,20,22,25,32,35,38];N=210; >> tm=[49,50,51;52,53,54;55,56,57;0,48,0]; >> for i=1:sum j=800*(i-1); X=goertzel(total_x(j+1:j+N),k+1); value=abs(X); figure(2) subplot(2,6,i); stem(k,value,'.','r'); title('FFT x(n)'); xlabel('k');ylabel('|X(k)|');limit=20; for i1=5:7 if value(i1)>limit break; end end for j1=1:4 if value(j1)>limit break;end end buffer(i)=tm(j1,i1-4);end j =0 j =800 j =1600 j =2400

河海大学研究生考试-水力学-2003真题答案

河海大学2003年攻读硕士学位研究生入学考试试题解析考试科目名称：水力学一、是非题（每题2分，共20分） 1． × 【参考答案】对于紊流： ①光滑区 λ仅仅是e R 的函数，λ～14 e R -，λ～14 v - ，f h ～14 v ； ②过渡区 λ是e R 和 d ?的函数，即(,)e R d λλ?=； ③粗糙区（又称阻力平方区） λ仅仅时d ?的函数，λ与e R 、v 均无关，f h ～2 v 。 2． × 【参考答案】由22221[1()]1()1c c dEs Fr dh k dEs i k ds k dh k i ds Fr ? =-???=-?-? =?-? 当, 1c c k h h k >< 即,c dEs h h ds ?>↑；同理,c dEs h h ds ?<↓ 亦即断面单位能量沿流程可能增大也可能减小。【点评】这种推导的方法要能熟悉运用。 3． × 【参考答案】根据通过泄水建筑物的几个不同的流量计算相应的co h 、2co h ，并确定t h ，使 2()co t h h -为最大值所对应的流量作为计算消力池深度d 的设计流量，而消力池长度则取最大流量作为设计流量。 4． × 【参考答案】达西定律： v k J =，式中：w h J l = —渗流水力坡度；w h —渗流水头损失； l —渗流长度；k —渗流系数。反映土壤渗流能力大小，其有速度的量纲。适用条件：均匀层流，砂土，具有线性规律。 5． √

【参考答案】2 2f l v h d g λ=??，故在恒定均匀流中，沿程水头损失与速度的平方成正比。 6． √ 【参考答案】动力相似?牛顿相似定律。两个液流在重力作用下的动力相似条件是它们的弗劳德数r F 相等。 7． × 【参考答案】几条管道在同一点分叉，然后又在另一点汇合的管道称为并联管道，其特点为： ①各条管道在分叉点和汇合点之间的水头损失相等； ②管道系统中的总流量等于各并联管道中的流量之和。 8． √ 【参考答案】这几组概念具体分类如下：?????????????? ?????? 均匀流恒定流急变流非均匀流渐变流均匀流非恒定流（仅限于管道）非均匀流【点评】?? ?层流紊流、 ???恒定流非恒定流、 ??? 均匀流非均匀流以及???急变流渐变流，对这四组概念要区别清楚，其中这四组概念之间并无必然的关系。 9． √ 【参考答案】流体质点做圆周运动不一定是有旋运动，流场中各点的旋转角速度0ω≠的液流运动称为有旋流动。 10． √ 【参考答案】相对压强可以大于、等于或小于零。二、选择题（每题2分，共20分） 1．（2）【参考答案】对于能量方程：22 1112221222w p v p v z z h g g g g ααρρ++=+++ α β ???与流速分布有关，分布越不均匀，α越大；层流的α、β比紊流的大；同是层（紊）流，