资源描述:
《5实验报告汇总》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实验五:连续信号的采样与重构一、已编的小程序段(1)模拟低通滤波器设计程序代码:elf;Fp=3500;Fs=4500;Wp=2*pi*Fp;Ws=2*pi*Fs;[NaWn]=buttord(Wp,0•5了30)1s1);[b,a]=butter(N,Wn,'s');wa=0:(3*Ws)/511:3*Ws;h=freqs(bza,wa);plot(wa/(2*pi),20*logl0(abs(h)));gridxlabel(1FrequencyfHz1);ylabe丄(1Gain,dB1);tit
2、le(1Gainresponse1);axis([03*Fs-605]);运行结果:(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('输入数列1);subplot(2,1,2)stem(n,y(1:length(x)));title('输岀数列T;偷入数列输出数列(1)时域欠采样程序代码:elf;n=0:49;m=0:50*3
3、-1;x=sin(2*pi*0.042*m);y=x([1:3:length(x)]);subplot(2,1,1)stem(nzx(1:50));axis([050-1.21.2]);title(1输入数列1);subplot(2flf2)stem(n^y);axis([050-1.21.2]);titleC输出数列T;输入数列(1)正弦信号采样程序代码:elf;t=0:0.0005:1;f=13;xa=cos(2*pi*f*t);subplot(2W)plot(t,xa);grid;xlabel(
4、1时间,msec');ylabel「幅值f);title(f连续时间信号x-{a}(t)1);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(n)-1;stem(kzxs);grid;xlabel(1时间,mesec!);ylabel(*幅值1);title(f连续时间信号x[n]1);axis([(length(n)-1)-1.21.2]);运行结果:连续时间信号)C・3(t)连续时间信号x[n]05
5、005•12674时间,mesec■!1(••••ff•••••i••••1••:<1>■•■•1‘<••>:18910二、程序设计实验(注:两种情况直接在一段程序中给出结果)elf;t=0:0.0000005:0.02;fl=150;f2=3000;xal=3*sin;xa2=3*sin(2*pi*f2壮);fs=5120;T=l/fs;nnl=-l:T:1;nn2=-l:T:1;xsl=3*sin(2*pi*fl*nnl);xs2=3*sin(2*pi*f2*nn2);kl=O
6、:length(nnl)-1;k2=0:length(nn2)-1;subplot(3,2,1);plot(t,xal);gridxlabel('时间,msec');ylabel('幅值');title('连续时间信号x_{al}(t)*);axis([00.02-33]);subplot(3,2,2);plot(t,xa2);gridxlabe丄(*时间,msec*);ylabel(*幅值');title('连续时间信号x_{a2}(t)');axis([00.002-33]);subplot(3,
7、2,3);stem(klzxsl);grid;xlabel('时间,msec*);ylabel(M畐值');title('离散时间信号xl[n]');axis([0100-33]);subplot(3,2Z4);stem(k2,xs2);grid;xlabel('时间,msec');ylabel('幅值');title('离散时间信号x2[n]');axis([05-33]);subplot(3,2,5);Nl=length(xsl);fxl=fft(xsl);dfl=fs/Nl;nl=0:Nl/2;
8、fl=nl*dfl;plot(f1zabs(fxl(nl+1))*2/Nl);grid;运行结果: