原始信号的谱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
可见两者的振幅谱完全一致