数字信号处理上机实验.doc

数字信号处理实验报告班级： 测控0701 姓名： 黄小宝 学号： 30703020191-1clearclose alla=1,2,3,4,b=3,4,5,6c=a+bd=a-be=a.*bf=a./bg=a.bsubplot(2,4,1):stem(a)subplot(2,4,2):stem(b)subplot(2,4,3):stem(c)subplot(2,4,4):stem(d)subplot(2,4,5):stem(e)subplot(2,4,6):stem(f)subplot(2,4,7):stem(g)1-2clearclose alln=0:15x1=0.8.nx2=abs(exp(0.2+3*i)*n)x3=angle(exp(0.2+3*i)*n)x4=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)subplot(2,2,1):stem(n,x1)subplot(2,2,2):stem(n,x2)subplot(2,2,3):stem(n,x3)subplot(2,2,4):stem(n,x4)1-3clearclose allt1=0:0.01:10x1=sin(2*pi*t1)t2=0:0.01:4x2=cos(100*pi*t2).*sin(pi*t2)subplot(2,1,1);stem(t1,x1)subplot(2,1,2);stem(t2,x2)2-1clearclose allb=0.05,0.05a=1,-0.9n=0:100x1=(n=0)=1)-(n-8=0)=1)x2=(n=0)x3=(n=0)=1)-(n-1=0)=1)y1=filter(b,a,x1)y2=filter(b,a,x2)y3=filter(b,a,x3)subplot(2,2,1):stem(y1)subplot(2,2,2):stem(y2)subplot(2,2,3):stem(y3)2-2clearclose alln=0:20h1=(n=0)=1)-(n-10=0)=1)h2=(n-0)=1)+2.5*(n-1)=1)+2.5*(n-2)=1)+(n-3)=1)x1=(n=0)=1)-(n-8=0)=1)y1=conv(x1,h1)y2=conv(x1,h2)subplot(3,2,1):stem(n,h1)subplot(3,2,2):stem(n,h2)subplot(3,2,3):stem(n,x1)n=0:(length(y1)-1)subplot(3,2,4):stem(n,y1)n=0:(length(y2)-1)subplot(3,2,5):stem(n,y2)2-3clearclose alln=0:1000b=1/100.49,0,1/100.49a=1,-1.8237,0.9801u=(n=0)=1)-(n-1=0)=1)q=(n=0)x=sin(0.014*n)+sin(0.4*n)y1=filter(b,a,q)y2=filter(b,a,u)y3=filter(b,a,x)subplot(2,2,1):stem(n,y1)subplot(2,2,2):stem(n,y2)subplot(2,2,3):stem(n,y3)3-1clear all; close all;figure(1);A=444.128;a=50*1.414*pi;w=50*1.414*pi;Fs1=1000;T1=1/Fs1;Fs2=300;T2=1/Fs2;Fs3=200;T3=1/Fs3;Tp=64/1000;N1=Tp*Fs1;N2=ceil(Tp*Fs2);N3=ceil(Tp*Fs3);n1=0:N1-1;n2=0:N2-1;n3=0:N3-1;x1=A*exp(-a.*n1*T1).*sin(w.*n1*T1).*(n1=0);x2=A*exp(-a.*n2*T2).*sin(w.*n2*T2).*(n2=0);x3=A*exp(-a.*n3*T3).*sin(w.*n3*T3).*(n3=0);X1=T1*fft(x1,N1);X2=T2*fft(x2,N2);X3=T3*fft(x3,N3);f1=n1/Tp; f2=n2/Tp;f3=n3/Tp; subplot(3,2,1),stem(n1,x1,.);axis(0,N1,1.2*min(x1),1.2*max(x1)subplot(3,2,2),plot(f1,abs(X1);axis(0,Fs1,0,1.2*max(abs(X1)subplot(3,2,3),stem(n2,x2,.);axis(0,N2,1.2*min(x2),1.2*max(x2)subplot(3,2,4),plot(f2,abs(X2);axis(0,Fs2,0,1.2*max(abs(X2)subplot(3,2,5),stem(n3,x3,.);axis(0,N3,1.2*min(x3),1.2*max(x3)subplot(3,2,6),plot(f3,abs(X3);axis(0,Fs3,0,1.2*max(abs(X3);3-2figure(2);M=27,N1=32;N2=16;n=0:M-1;xa=1:floor(M+1)/2);xb=ceil(M/2)-1:-1:1;x=xa,xb;Xk=fft(x,1024);X32=fft(x,N1);X16=X32(1:2:N1);x32=ifft(X32,N1);x16=ifft(X16,N2);subplot(3,2,2),stem(n,x,.);axis(0,26,0,1.2*max(x)subplot(3,2,1),plot(0:1023,abs(Xk);axis(0,512,0,1.2*max(abs(Xk);subplot(3,2,4),stem(0:N1-1,x32,.);axis(0,N1,0,1.2*max(x)subplot(3,2,3),stem(0:N1-1,abs(X32),.);axis(0,N1/2,0,1.2*max(abs(X32);subplot(3,2,6),stem(0:N2-1,x16,.);axis(0,N2,0,1.2*max(x16)subplot(3,2,5),stem(0:N2-1,abs(X16),.);axis(0,N2/2,0,1.2*max(abs(X16);4-1x1=1,1,1,1x2=1,2,3,4,4,3,2,1x3=4,3,2,1,1,2,3,4n1=8n2=16y1=fft(x1,n1);z1=abs(y1)y2=fft(x1,n2);z2=abs(y2)y3=fft(x2,n1);z3=abs(y3)y4=fft(x2,n2);z4=abs(y4)y5=fft(x3,n1);z5=abs(y5)y6=fft(x3,n2);z6=abs(y6)subplot(3,2,1),stem(0:length(z1)-1,z1)subplot(3,2,2),stem(0:length(z2)-1,z2)subplot(3,2,3),stem(0:length(z3)-1,z3)subplot(3,2,4),stem(0:length(z4)-1,z4)subplot(3,2,5),stem(0:length(z5)-1,z5)subplot(3,2,6),stem(0:length(z6)-1,z6)4-2n=0:32x1=cos(pi*n)/4);x2=cos(pi*n)/4)+cos(pi*n)/8);n1=8;n2=16;y1=fft(x1,n1);z1=abs(y1);y2=fft(x1,n2);z2=abs(y2);y3=fft(x2,n1);z3=abs(y3);y4=fft(x2,n2);z4=abs(y4);subplot(2,2,1),stem(0:length(z1)-1,z1);subplot(2,2,2),stem(0:length(z2)-1,z2);subplot(2,2,3),stem(0:length(z3)-1,z3);subplot(2,2,4),stem(0:length(z4)-1,z4);5-1function st=mstgN=1600 Fs=10000;T=1/Fs;Tp=N*T; t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;fc1=Fs/10; fm1=fc1/10;fc2=Fs/20;fm2=fc2/10; fc3=Fs/40;fm3=fc3/10;xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); st=xt1+xt2+xt3; fxt=fft(st,N);subplot(3,1,1)plot(t,st);grid;xlabel(t/s);ylabel(s(t);axis(0,Tp/8,min(st),max(st);title(a) s(t)的波形)subplot(3,1,2)stem(f,abs(fxt)/max(abs(fxt),.);grid;title(b) s(t)的频谱)axis(0,Fs/5,0,1.2);xlabel(f/Hz);ylabel(幅度)5-2clear all;close allFs=10000;T=1/Fs;st=mstg;fp=280;fs=450;wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; N,wp=ellipord(wp,ws,rp,rs); B,A=ellip(N,rp,rs,wp);H,W=freqz(B,A,1000);y1t=filter(B,A,st); figure(2);subplot(2,1,1);plot(W/pi,20*log10(abs(H)/max(abs(H);grid on;axis(0,1,min(20*log10(abs(H)/max(abs(H),5); yt=y_1(t);subplot(2,1,2);n=0:length(st)-1;t=n*T;plot(t,y1t);grid on; fpl=440;fpu=560;fsl=275;fsu=900;wp=2*fpl/Fs,2*fpu/Fs;ws=2*fsl/Fs,2*fsu/Fs;rp=0.1;rs=60; N,wp=ellipord(wp,ws,rp,rs); B,A=ellip(N,rp,rs,wp);H,W=freqz(B,A,1000);y2t=filter(B,A,st);figure(3);subplot(2,1,1);plot(W/pi,20*log10(abs(H)/max(abs(H);grid on;axis(0,1,min(20*log10(abs(H)/max(abs(H),5); yt=y_2(t);subplot(2,1,2);n=0:length(st)-1;t=n*T;plot(t,y2t);grid on; fp=890;fs=600;wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; N,wp=ellipord(wp,ws,rp,rs); B,A=ellip(N,rp,rs,wp,high);H,W=freqz(B,A,1000);y3t=filter(B,A,st);figure(4);subplot(2,1,1);plot(W/pi,20*log10(abs(H)/max(abs(H);grid on;axis(0,1,min(20*log10(abs(H)/max(abs(H),5); yt=y_3(t);subplot(2,1,2);n=0:length(st)-1;t=n*T;plot(t,y3t);grid on;6-1function xt=xtgN=1000,Fs=1000;T=1/Fs;Tp=N*Tt=0:T:(N-1)*Tfc=Fs/10;f0=fc/10mt=cos(2*pi*f0*t)ct=cos(2*pi*fc*t)xt=mt.*ctnt=2*rand(1,N)-1fp=150;fs=200;Rp=0.1;As=70fb=fp,fs;m=0,1dev=10(-As/20),(10(Rp/20)-1)/(10(Rp/20)+1)n,fo,mo,W=remezord(fb,m,dev,Fs)hn=remez(n,fo,mo,W)yt=filter(hn,1,10*nt)xt=xt+ytfst=fft(xt,N);k=0:N-1;f=k/Tpsubplot(2,1,1);plot(t,xt);grid;xlabel(t/s);ylabel(x(t)axis(0,Tp/5,min(xt),max(xt);title(a)信号加噪声波形)subplot(2,1,2);plot(f,abs(fst)/max(abs(fst);grid;title(b) 信号加噪声频谱)axis(0,Fs/2,0,1.2);xlabel(f/Hz);ylabel(幅度)6-2clear all;close all;N=1000;xt=xtg