资源描述：
《实验五：连续信号的采样与重构》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
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('Gainresponse');axis([03*Fs-605]);运行结
2、果：(1)时域过采样程序代码：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('输出数列');运行结果：(1)时域欠采样程序代码：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([050-1
3、.21.2]);title('输入数列');subplot(2,1,2)stem(n,y);axis([050-1.21.2]);title('输出数列');运行结果：(1)正弦信号采样程序代码：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([01-1.21.2])subplot(2,1,2);T=0.1;n=0:T:1;xs=cos(2*pi*f*n);k=0:length(
4、n)-1;stem(k,xs);grid;xlabel('时间，mesec');ylabel('幅值');title('连续时间信号x[n]');axis([(length(n)-1)-1.21.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
5、);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([00.02-33]);subplot(3,2,2);plot(t,xa2);gridxlabel('时间,msec');ylabel('幅值');title('连续时间信号x_{a2}(t)');axis([00.002-33]);subplot(3,2,3);stem(k1,xs1);grid;xlabel('时间,mse
6、c');ylabel('幅值');title('离散时间信号x1[n]');axis([0100-33]);subplot(3,2,4);stem(k2,xs2);grid;xlabel('时间,msec');ylabel('幅值');title('离散时间信号x2[n]');axis([05-33]);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;运行结果：
