5实验报告汇总.docxVIP

实验五：连续信号的采样与重构已编的小程序段模拟低通滤波器设计程序代码：clf;Fp=3500;Fs=4500;Wp=2*pi*Fp;Ws=2*pi*Fs;[N,Wn]=buttord(Wp,Ws,0.5,30,s);[b,a]=butter(N,Wn,s);wa=0:(3*Ws)/511:3*Ws;h=freqs(b,a,wa);plot(wa/(2*pi),20*log10(abs(h)));gridxlabel(Frequency,Hz);ylabel(Gain,dB);title(Gain response);axis([0 3*Fs -60 5]);运行结果：时域过采样程序代码：n=0:50;x=sin(2*pi*0.12*n);y=zeros(1,3*length(x));y([1:3:length(y)])=x;subplot(2,1,1)stem(n,x);title(输入数列);subplot(2,1,2)stem(n,y(1:length(x)));title(输出数列);运行结果：时域欠采样程序代码：clf;n=0:49;m=0:50*3-1;x=sin(2*pi*0.042*m);y=x([1:3:length(x)]);subplot(2,1,1)stem(n,x(1:50));axis([0 50 -1.2 1.2]);title(输入数列);subplot(2,1,2)stem(n,y);axis([0 50 -1.2 1.2]);title(输出数列);运行结果：正弦信号采样程序代码：clf;t=0:0.0005:1;f=13;xa=cos(2*pi*f*t);subplot(2,1,1)plot(t,xa);grid;xlabel(时间，msec);ylabel(幅值);title(连续时间信号 x-{a}(t));axis([0 1 -1.2 1.2])subplot(2,1,2);T=0.1;n=0:T:1;xs=cos(2*pi*f*n);k=0:length(n)-1;stem(k,xs);grid;xlabel(时间，mesec);ylabel(幅值);title(连续时间信号 x[n]);axis([(length(n)-1)-1.2 1.2]);运行结果：程序设计实验（注：两种情况直接在一段程序中给出结果）clf;t=0:0.0000005:0.02;f1=150;f2=3000;xa1=3*sin(2*pi*f1*t);xa2=3*sin(2*pi*f2*t);fs=5120;T=1/fs;nn1=-1:T:1;nn2=-1:T:1;xs1=3*sin(2*pi*f1*nn1);xs2=3*sin(2*pi*f2*nn2);k1=0:length(nn1)-1;k2=0:length(nn2)-1;subplot(3,2,1);plot(t,xa1);gridxlabel(时间,msec);ylabel(幅值);title(连续时间信号 x_{a1}(t));axis([0 0.02 -3 3]);subplot(3,2,2);plot(t,xa2);gridxlabel(时间,msec);ylabel(幅值);title(连续时间信号 x_{a2}(t));axis([0 0.002 -3 3]);subplot(3,2,3);stem(k1,xs1);grid;xlabel(时间,msec);ylabel(幅值);title(离散时间信号 x1[n]);axis([0 100 -3 3]);subplot(3,2,4);stem(k2,xs2);grid;xlabel(时间,msec);ylabel(幅值);title(离散时间信号 x2[n]);axis([0 5 -3 3]);subplot(3,2,5);N1=length(xs1);fx1=fft(xs1);df1=fs/N1;n1=0:N1/2;f1=n1*df1;plot(f1,abs(fx1(n1+1))*2/N1);grid;运行结果：