武漢理工大學(xué)數(shù)字信號處理第一次實驗報告完整答案

第一部分：例題調(diào)試（一）單位抽樣序列代碼如下：%先建立函數(shù)impseq(n1,n2,n0)function[x,n]=impseq(n1,n2,n0)n=[n1:n2];x=[(n-n0)==0];%編寫主程序調(diào)用該函數(shù)[x,n]=impseq(-2,8,2);stem(n,x);title('DSP_example01單位抽樣序列——徐櫻笑');（二）單位階躍序列代碼如下：%先建立函數(shù)stepseq(n1,n2,n0)function[x,n]=stepseq(n1,n2,n0)n=[n1:n2];x=[(n-n0)>=0];%編寫主程序調(diào)用該函數(shù)[x,n]=stepseq(-2,8,2);stem(n,x)title('DSP_example01單位階躍序列——徐櫻笑');（三）實數(shù)指數(shù)序列x(n)=an(運算符“.^”)代碼如下：n=[0:10];x=0.9.^n;stem(n,x);title('DSP_example03實數(shù)指數(shù)序列——徐櫻笑');（四）復(fù)數(shù)指數(shù)序列代碼如下：n=[-10:10];alpha=-0.1+0.3*j;x=exp(alpha*n);real_x=real(x);image_x=imag(x);mag_x=abs(x);phase_x=angle(x);subplot(2,2,1);stem(n,real_x);subplot(2,2,2);stem(n,image_x);subplot(2,2,3);stem(n,mag_x);subplot(2,2,4);stem(n,phase_x);suptitle('DSP_example04復(fù)數(shù)指數(shù)序列——徐櫻笑')（五）正、余弦序列（）代碼如下：n=[0:42];x=2*cos(0.1*pi*n+pi/3);stem(n,x);title('DSP_example05正余弦序列——徐櫻笑');

第二部分：課后習(xí)題（一）利用zeros函數(shù)生成單位抽樣序列：代碼如下：n=-10:10;y=[zeros(1,10),1,zeros(1,10)];stem(n,y)axis([-10,10,0,2]);title('DSP_test01單位抽樣序列——徐櫻笑');（二）利用zeros函數(shù)和ones函數(shù)生成單位階躍序列代碼如下：n=-5:10;y=[zeros(1,5),1,ones(1,10)];stem(n,y)axis([-5,10,0,2]);title('DSP_test02單位階躍序列——徐櫻笑');（三）已知x1(n)=u(n+1)(-3<n<5);x2(n)=u(n-3)(-4<n<7)。求：x(n)=x1(n)+x2(n)代碼如下：%建立移位函數(shù)sigshiftfunction[y,n]=sigshift(x,m,n0)n=m+n0;y=x;%建立疊加函數(shù)sigaddfunction[y,n]=sigadd(x1,n1,x2,n2)n=min(min(n1),min(n2)):max(max(n1),max(n2));y1=zeros(1,length(n));y2=y1;y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;y=y1+y2;%編寫主程序完成信號運算n1=-3:5;n2=-4:7;y1=[zeros(1,3),1,ones(1,5)];y2=[zeros(1,4),1,ones(1,7)];[y11,n11]=sigshift(y1,n1,-1);[y21,n21]=sigshift(y2,n2,3);[y3,n3]=sigadd(y11,n11,y21,n21);stem(n3,y3)axis([-5,15,0,4]);title('DSP_test03——徐櫻笑');（四）已知x1(n)=3e-0.25n(-2<n<8)，x2(n)=u(n+1)(-3<n<6)，求：x(n)=x1(n)*x2(n)將函數(shù)conv稍加擴展為函數(shù)conv-m，它可以對任意的序列求卷積。代碼如下：function[y,ny]=conv-m(x,nx,h,nh,p)%信號處理的改進(jìn)卷積程序nyb=nx(1)+nh(1);nyc=nx(length(x))+nh(length(h));ny=[nyb:p:nyc];y=conv(x,h);%編寫主程序完成信號運算p=1;alpha=-0.25;n1=[-2:p:8];n2=[-3:p:6];x1=3*exp(alpha*n1);x2=[zeros(1,2),1,ones(1,7)];[y,ny]=conv_m(x1,n1,x2,n2,p);subplot(3,1,1);stem(n1,x1)subplot(3,1,2);stem(n2,x2)subplot(3,1,3);stem(ny,y)suptitle('DSP_test04——徐櫻笑');（五）已知一個IIR數(shù)字低通濾波器的系統(tǒng)函數(shù)公式為輸入一個矩形信號序列x=square(n/5)(-2<n<10)，求該系統(tǒng)的響應(yīng)。function[y,ny]=conv-m(x,nx,h,nh,p)%信號處理的改進(jìn)卷積程序nyb=nx(1)+nh(1);nyc=nx(length(x))+nh(length(h));ny=[nyb:p:nyc];y=conv(x,h);%編寫主程序完成信號運算b=[0.1321,0.3963,