武汉科技大学《信号与系统》实验四实验报告参考模板

1/7专业班级电信(DB)1102班姓名XX学号 3505X验时间2013年6月18日指导教师陈华丽成绩1.掌握离散信号谱分析的方法：序列的傅里叶变换、离散傅里叶级数、离散傅里叶变换、快速傅里理解这些变换之间的关系；叶级数、离散傅里叶变换、快速傅里叶变换的Matlab实现；4.学习用FFT对连续信号和离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便1.对连续信号xa(t)=Ae-atsin(业0t)u(t)1 (A=444.128，a=502"，业0=502")进行理0.6想采样，可得采样序列a0.4图1给出了xa(t)的幅频特性曲线，由此图可以确00100200300400500定对xa(t)采用的采样频率。分别取采样频率为f/Hz(4)观察上述三种情况下，x(n)的幅度曲线是否一致？为什么？3.(1)编制信号产生子程序，产生以下典型信号供谱分析用。"x(n)=cosn24"x(n)=sinn384(2)对信号x(n)，x(n)，x(n)进行两次谱分析，FFT的变换区间N分别取8和16，观察两23 (3)连续信号x(n)的采样频率f=64Hz，N=16,32,64。观察三次变换的结果是否一致？为什4s武汉科技大学《信号与系统》实验四实验报告functionx=ff(T)%定义函数文件n=0:50;A=444.128;a=50*sqrt(2)*pi;w0=50*sqrt(2)*pi;x=A*exp(-a*n*T).*sin(w0*n*T);k=-250:250;w=pi/125*k;X=x*(exp(-j*pi/125)).^(n'*k);实验记录及个人小结(包实验源程注结果分析与讨论Y=abs(y);plot(w/pi,Y)ndy=ff(1/1000);title('f=1000Hz')200HZ：y=ff(1/300);title('f=300Hz')7000600050004000300020000002000800004002000-1.5-1.5-0.5-0.5f=1000Hz0f=300Hz00.50.5111.5.522300HZ：y=ff(1/200);title('f=200Hz')00512武汉科技大学《信号与系统》实验四实验报告functiony=x(n)%定义函数文件y=cos(0.48*pi*n)+cos(0.52*pi*n);end(1)n=0:10;f1=x(n);f2=(fft(f1));stem(n,f2)title('x(n)(n=0:10)')axis([0,10,-2,8])8765432100123x(n)(n=0:10)456789(2)f3=zeros(1,20);forn=0:10a=2*n+1;f3(a)=n;endf4=x(f3);f5=fft(f4)stem(f5)8642005x(n)(n=0:20)2025结(包括：(3)n=0:100;f1=x(n);f2=(fft(f1));stem(n,f2)002030405060708090x(n)的幅度曲线是否一致(1)N=8;n=0:N-1;x1=(n+1).*(u(n-4)-u(n))+(8-n).*(u(n-8)-u(n-4));x2=cos(pi*n/4);x3=sin(pi*n/8);%x4=cos(8*pi.*t)+cos(16*pi.*t)+cos(20*pi.*t);y1=fft(x1,N);subplot(2,2,1);stem(n,abs(y1));y2=fft(x2,N);subplot(2,2,2);武汉科技大学《信号与系统》实验四实验报告y2=fft(x2,N);subplot(2,2,2);stem(n,abs(y2));y3=fft(x3,N);subplot(2,2,3);stem(n,abs(y3));50064200224466884321002468(2)行周期延拓后的傅里叶级数变换的主值系列,里叶级数也不相同,所以变换后的结果也不相同.(3)n=0:15;f=64;x4=cos(8*pi*n/f)+cos(16*pi*n/f)+cos(20*pi*n/f);x4k16=fft(x4,16);结(包括：subplot(3,1,1)stem(n,abs(x4k16))title('x4FFT变换N=16')axis([0,14,0,15])n=0:31;f=64;x4=cos(8*pi*n/f)+cos(16*pi*n/f)+cos(20*pi*n/f);xsubplot(3,1,2)5stem(n,abs(x4k5title('x4FFT变换N=32')axis([0,31,0,20])0204681020551015510152025000武汉科技大学《信号与系统》实验四实验报告n=0:63;f=64;x4=cos(8*pi*n/f)+cos(16*pi*n/f)+cos(20*pi*n/f);x4FFT变换N=x4FFT变换N=16subplot(3,1,3)stem(n,abs(x4ktitle('x4FFT变换N=64')axis([0,63,0,40])502406824068x4FFT变换N=3255列也不相同,所以傅里叶级数也不相同,所以变换后的结果也不相同.0051505x4FFT变