1、实验五：连续信号的采样与重构一、 已编的小程序段(1) 模拟低通滤波器设计程序代码：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);运行结果：(2) 时域过采样程序代码：
2、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(输出数列);运行结果：(3) 时域欠采样程序代码：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
3、,2)stem(n,y);axis(0 50 -1.2 1.2);title(输出数列);运行结果：(4) 正弦信号采样程序代码：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);ylab
4、el(幅值 );title(连续时间信号 xn);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
5、);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(离散时间信号 x1n);axis(0 100 -3 3);subplot(3,2,4);stem(k2,xs2);grid;xlabel(时间 ,msec);ylabel(幅值 );title(离散时间信号 x2n);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;运行结果：
