《数字信息处理》课程教学资源（PPT实验）DFT及互相关的初步应用

实验 DFT及互相关的初步应用 195758

19:57:58 3 实验一 DFT及互相关的初步应用

姆一、实验目的 1.熟悉频谱(幅度谱)的特性 2熟悉互相关的初步应用 195758

19:57:58 4 一、实验目的 1. 熟悉频谱（幅度谱）的特性 2. 熟悉互相关的初步应用

姆三、实验内容 1利用 FFT/IFFT实现简单的频域滤波 要求与思路: 1)构造一低频离散信号f 2)在此低频信号上加随机噪声得到混合信号g; 3)对信号g进行FFT变换得到其幅度谱和相位谱 4)对高于某个频率的幅度谱全部置为0; 5)利用FT对频谱做反变换得到滤波后的信号f 6)画图比较f和f。 -部分提示见下页 195758

19:57:58 5 三、实验内容 1.利用FFT/IFFT实现简单的频域滤波 要求与思路： 1）构造一低频离散信号f； 2）在此低频信号上加随机噪声得到混合信号g； 3）对信号g进行FFT变换得到其幅度谱和相位谱； 4）对高于某个频率的幅度谱全部置为0； 5）利用IFFT对频谱做反变换得到滤波后的信号f’; 6）画图比较 f 和 f’。 --------部分提示见下页…

nx=100 x1=[l:nx]; y=sin(x1*0.3);%自己构造有用的信号序列 noise=rand(l, nx)-05 y2=yl noise; figure(1); plot(xl, y1,r, xl,y2, b: ) XXX=fit(y2); figure(2); plot(abs(XX);tite('信号的幅度谱); %改造频谱,实现滤波 fori=10:nx-10%自己选择合适漶波范围 XXXO=0 en f1= ifft(XXx) figure(3) lotx12yl,r,xl,,b:);itc滤波的效果); 195758

19:57:58 6 nx = 100; x1 = [1:nx]; y1 = sin(x1*0.3); %自己构造有用的信号/序列 noise = rand(1,nx) - 0.5; y2 = y1 + noise; figure(1); plot(x1,y1,'r', x1,y2,'b:'); XXX = fft(y2); figure(2); plot(abs(XXX)); title('信号的幅度谱'); %改造频谱，实现滤波 for i = 10 : nx-10 %自己选择合适滤波范围 XXX(i) = 0; end f1 = ifft(XXX); figure(3); plot(x1,y1,'r',x1,f1,'b:'); title('滤波的效果');

5 . 5 lrl 5 70

19:57:58 7

信号的幅度谱 40 25 20 15 设置为0 0 0102034570g 100

19:57:58 8 设置为0

滤波的效果 0.8 0.6 0.4 0.2 0.2 0.4 0.6 0.8 01020④405670890100

19:57:58 9

姆2互相关的一个简单应用 ■目的:利用互相关运算寻找被噪声淹没 的有用信号的位置。 互相关:rn=∑xmym- m=-00 195758

19:57:58 10 2. 互相关的一个简单应用 ◼目的：利用互相关运算寻找被噪声淹没 的有用信号的位置。 [ ] [ ] [ ] xy m r l x m y m l  =− 互相关： = − 

yln ay出现在y2m的某个位置 y∠n -0.5 20 10 在y回上加随机噪声→y3可 100 140

19:57:58 11 y1[n] y2[n] 在y2[n]上加随机噪声➔ y3[n] y1[n]出现在y2[n]的某个位置➔

yln n 2 120 140 rn= xcorr(y3, y1); 5 -5 150 100 100

19:57:58 12 y1[n] y3[n] r[n] = xcorr(y3,y1);