《数字信号处理》上机实习报告

原始信号的谱151050050100150处理后的信号210-1-20510窗函数:

源程序见matlab文件夹“chuang”

运行结果:

带阻滤波器的频谱10.50050100150处理后的信号的谱1050050100150

六、设计一个二维滤波处理程序(分别做低通、高通等处理)。 程序思路:

二维滤波在地震中就是在时间上滤一遍得到频率,在空间域内在滤一遍得到波数。

选择的谱也要满足对称共轭,经过二维滤波后,再反变换回去时得到的信号一定是没有虚数的,为实信号;

源程序见matlab文件夹“erweilvbo”: 运行结果: 原始信号:

在行方向做高通滤波后的信号图像:

在列方向做低通滤波后的信号图像:

七、验证时间域的循环褶积对应的是频率域的乘积;线性褶积则不然。

由前面的程序(第三大题)可以看出对于循环褶积,不管是在时间域还是在频率域的计算结果都是一样。

a=[4 5 6 0];

>> b=[2 3 4 5]; 运行结果为 8 22 43 58 49 30

但是对于线性褶积,由于当两个信号非零个数不同时不能进行相乘,所以必须得先补零,然后做谱的乘积,

a=[4 5 6 0]; >> b=[2 3 4 5]; >> c=fft(a)

d=fft(b) e=c.*d ifft(e) ans =

57 52 43 58 可见两者并不相等

八、请用通俗、易懂的语言说明数字信号处理中的一种性质、一条定理或一个算例(顺便利用Matlab对其进行实现)。

时移性质:好比一个人在北京,时移相当于他去了上海,谱的相位发生了改变,相当于他的位置发生了变化,但是谱的模并没有改变,相对的是这个人本身并没有发生任何改变,谱的重要信息没有随时移,人的位置虽然改变了,但是人的性质也没有改变。 如:a=[1 2 3 4 5];

b=[2 3 4 5 1];%b是a的向左时移一位 subplot(1,2,1); plot(abs(fft(a))); subplot(1,2,2); plot(abs(fft(b)));

运行结果:

1616141412121010886644212345212345

可见两者的振幅谱完全一致

联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4