《《信号与系统》实验四》由会员分享，可在线阅读，更多相关《《信号与系统》实验四（6页珍藏版）》请在金锄头文库上搜索。
1、 . 实验容1.对连续信号，进展理想采样，可得采样序列。图1给出了的幅频特性曲线，由此图可以确定对采用的采样频率。分别取采样频率为 1KHz、300Hz和200Hz，画出所得采样序列的幅频特性。并观察是否存在频谱混叠。 图1 连续信号2. 设1取时，求的FFT变换，并绘出其幅度曲线。2将(1)中的以补零方式加长到，求并绘出其幅度曲线。3取，求并绘出其幅度曲线。4观察上述三种情况下，的幅度曲线是否一致？为什么？3. 1编制信号产生子程序，产生以下典型信号供谱分析用。2对信号，进展两次谱分析，FFT的变换区间N分别取8和16，观察两次的结果是否一致？为什么？3连续信号的采样频率，。观察三次变换的结
2、果是否一致？为什么？实验记录及个人小结包括：实验源程序、注释、结果分析与讨论等题一:function x=caiyang(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);h=1,2.5,2.5,1;m=0:3;H=h*(exp(-j*pi/125).(m*k);y=H.*X;Y=abs(y);plot(w/pi,Y)endy=caiyang(1/1000);title(f
3、=1000Hz)y=caiyang(1/300);title(f=300Hz)y=caiyang(1/200);title(f=200Hz)通过观察图像,可以看出采样频率为300Hz和200Hz时存在频谱混叠.题二:function y=x(n)%绘制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)(2):f3=zeros(1,20);for n=0:10 a=2*n+1; f3(a)=n;endf4=x(f3
4、);f5=fft(f4)stem(f5)(3):n=0:100;f1=x(n);f2=(fft(f1);stem(n,f2)题三(1)X1(n)函数:function y=x(n)if n=0&n=4&n=7 y=8-n;stem(n,y)else y=0; stem(n,y) hold onendendX2(n)函数:function y=x(n)y=cos(pi*n/4)endX3(n)函数:function y=x(n)y=sin(pi*n/8)endX4(n)函数:function y=x(t) y=cos(8*pi*t)+cos(16*pi*t)+cos(20*pi*t);end题三
5、(2):x1=1 2 3 4 4 3 2 1;subplot(3,1,1)stem(x1)title(x1原图像)x1k8=fft(x1);subplot(3,1,2)stem(abs(x1k8)title(x1(n)FFT变换 N=8)x2=x1,zeros(1,8)x1k16=fft(x2);subplot(3,1,3)stem(abs(x1k16)title(x1(n)FFT变换 N=16)X2(n)de 8,16次FFT变换n=0:7;x4n=cos(pi*n/4);subplot(3,1,1)stem(n,abs(x4n)title(x4n图像)x4k8=fft(x4n);subpl
6、ot(3,1,2)stem(n,abs(x4k8)title(X3(n)FFT变换 N=8)m=0:15;x4k16=fft(x4n,16);subplot(3,1,3)stem(m,abs(x4k16)title(X3(n)FFT变换 N=16)X3(n)de 8,16次FFT变换n=0:7;x4n=sin(pi*n/8);subplot(3,1,1)stem(n,abs(x4n)title(x4n图像)x4k8=fft(x4n);subplot(3,1,2)stem(n,abs(x4k8)title(X3(n)FFT变换 N=8)m=0:15;x4k16=fft(x4n,16);subpl
7、ot(3,1,3)stem(m,abs(x4k16)title(X3(n)FFT变换 N=16)结论:FFT的变换区间分别取8,16时,同一函数的两次结果不一样,因为FFT变化可以看成是将序列进展周期延拓后的傅里叶级数变换的主值系列.N不同,进展的周期延拓序列也不一样,所以傅里叶级数也不一样,所以变换后的结果也不一样.题三(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);x4k32=fft(x4,32);subplot(3,1,2)stem(n,abs(x4k32)title(x4FFT变换 N=32)axis(0,31,0,20)n=0:63;f=64;x4=cos(8*pi*n/f)+cos(16*pi*n/f)+cos(20*pi*n/f);x4k64=fft(x4,64);subplot(3,1,3)stem(n,abs(x4k64)title(x4FFT变换 N=64)axis(0,63,0,40) /
