信号的频谱分析实验报告.docVIP

实验四 信号的频谱分析（1天） 一．实验内容 1. 利用FFT分析连续周期，非周期信号的频谱，如周期，非周期方波，正弦信号等。理解CFS，CTFT与DFT（FFT）的关系； 2. 利用FFT分析离散周期，非周期信号的频谱，如周期，非周期方波，正弦信号等。理解DFS，DTFT与DFT（FFT）的关系，并讨论连续信号与离散信号频谱分析方法的异同。 二．实验步骤： 1.正余弦的fft程序： f0=50; t=0:0.001:0.02; w0=2*pi*f0; y1=sin(w0*t); figure(3),plot(t,y1); y2=cos(w0*t); figure(4),plot(t,y2); fy1=fft(y1); fy2=fft(y2); figure(1),subplot(211)，stem(abs(fy1)); title(正弦fft幅度) subplot(212)，stem(angle(fy1)); title(正弦fft相位’) figure(2),subplot(211)，stem(abs(fy2)); title(余弦fft幅度’) subplot(212)，stem(angle(fy2)); title(余弦fft相位) 正弦fft（幅度和相位）： 余弦fft（幅度和相位）： 2.半波对称正负方波： t=0:0.01:5; w0=2*pi*1; y=square(w0*t,50); figure(1),plot(t,y); fy=fft(y); figure(2),subplot(211)，stem(abs(fy)); subplot(212)，stem(angle(fy)); 3.方波： t=-5:0.02:5; y=((t=-1)-(t=1)); fy=fft(y); figure(2),subplot(211)，stem(abs(fy)); subplot(212)，stem(angle(fy)); 4．信号的频谱： t=0:0.1:5; y=exp(-t).*(t=0); figure(1),plot(t,y); figure(2),subplot(211)，stem(abs(fft(y))); subplot(212)，stem(angle(fft(y))); 5．做FFT变换验证: t=0:0.1:5; y=exp(-t).*(t=0); fy=fft(y)/10; f=0:0.1:5; w=2*pi*f; m=w.*w; fy1=1./sqrt(1+m); figure(1),plot(fy1,r); hold on;stem(abs(fy)); axis([0 50 0 1]); 6．Fft计算线性卷积验证实验三: 1.方波的脉冲响应验证： t=0:0.1:2; x=(t=0); m=-[(t=1)-(t=0)]; fy1=fft(x,41); fy2=fft(m,41); fy=fy1.*fy2; %时域卷积频域相乘 Z=ifft(fy); %做逆傅里叶变换 z=conv(x,m); figure(1),subplot(121)，stem(Z); axis([0 45 0 12]); title(ifft反卷积结果) subplot(122)，stem(z,r); axis([0 45 0 12]); title(conv函数卷积结果) 2.系统为exp（-2*t）u（t）的方波响应和阶跃响应验证： t=0:0.1:10; x=(t=0); m=-[(t=1)-(t=0)]; y=exp(-2*t).*(t=0); k=conv(x,y);h=conv(m,y); M=fft(m,200);Y=fft(y,200); fy1=M.*Y; Fy1=ifft(fy1); figure(1),subplot(121)，stem(h); title(方波响应conv实现) subplot(122)，stem(Fy1); axis([0 300 0 5]); title(ifft实现) X=fft(x,200); Y=fft(y,200); fy2=X.*Y; Fy2=ifft(fy2); figure(2),subplot(121)，stem(k); title(阶跃响应conv实现) subplot(122)，stem(Fy2); axis([0 300 0 6]); title(ifft实现) 三．实验结果： 1.实验曲线与理论曲线之间的逼近：实际曲线的横坐标需乘以截断长度（L）的倒数才能实现从K轴——f轴，纵坐标需除以N（L/Ts）。才能与理论曲对应。 2.经验证利用FFT也能实现线性卷积，即利用时域卷积，频域相乘的原理，将原信号做FFT，在频域做乘，经过逆傅里叶变换即可得出两