数字信号处理上机.doc

《数字信号处理上机.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、2.38(4).x(n)=,-10≤n≤10在给定的区间上产生信号，使用stem()函数画图，分别画出幅度、相位、实部和虚部.n=-10:0.5:10;m=-0.2+j*0.4;x=exp(m*n);subplot(2,2,1);stem(n,real(x));gridon;title('Realpart');xlabel('n');ylabel('Amplitude');subplot(2,2,2);stem(n,imag(x));gridon;title('Imaginarypart');xla

2、bel('n');ylabel('Amplitude');subplot(2,2,3);stem(n,abs(x));gridon;title('Absolutepart');xlabel('n');ylabel('y');subplot(2,2,4);stem(n,(180/pi)*angle(x));gridon;title('Anglepart');xlabel('n');ylabel('y');2.40.已知x(n)={1,2,3,4,5},0≤n≤4h(n)={1,-2,1,3},0≤n≤

3、3nx=[0:4];x=[1,2,3,4,5];nh=[0:3];h=[1,-2,1,3];[y,ny]=conv_m(x,nx,h,nh)%ny为y在横坐标轴上的位置取值%扩展x(n)和y(n)的长度n=length(ny);x1=zeros(1,n);h1=zeros(1,n);x1(find((ny>=min(nx))&(ny<=max(nx))==1))=x;h1(find((ny>=min(nh))&(ny<=max(nh))==1))=h;subplot(3,1,1);stem(ny,x

4、1);xlabel('n');ylabel('x(n)');axis([07-36]);subplot(3,1,2);stem(ny,h1);xlabel('n');ylabel('h(n)');subplot(3,1,3);stem(ny,y);xlabel('n');ylabel('y(n)');function[y,ny]=conv_m(x,nx,h,nh)ny0=nx(1)+nh(1);nyf=nx(end)+nh(end);y=conv(x,h);ny=[ny0:nyf];end2.42.

5、已知某线性移不变系统由下列差分方程描述：y(n)-y(n-1)+0.9y(n-2)=x(n);(1)计算并画出在-20≤n≤100内的冲激响应h(n)；(2)计算并画出在-20≤n≤100内的单位阶跃响应；(3)判断系统是否稳定。本题使用impseq()、stepseq()和filter()函数求解。n=[-20:100];b=[1];a=[1,-1,0.9];x=impseq(0,-20,100);h=filter(b,a,x);subplot(2,1,1);stem(n,h);title('冲激

6、响应');xlabel('n');ylabel('h(n)');x=stepseq(0,-20,100);s=filter(b,a,x);subplot(2,1,2);stem(n,s);title('阶跃响应');xlabel('n');ylabel('s(n)');z=roots(a);magz=abs(z);function[x,n]=impseq(np,ns,nf)if((np

7、(np>nf)

8、(ns>nf))error('参数必须满足ns<=np<=nf')endn=ns:nf;x

9、=[(n-np)==0];function[x,n]=stepseq(np,ns,nf)if((np

10、(np>nf)

11、(ns>nf))error('参数必须满足ns<=np<=nf')endn=ns:nf;x=[(n-np)>=0];2.45(3)x(n)=,0≤n≤7,计算下列序列的傅里叶变换(DTFT)X()，并画出其幅度和相位函数，使用plot()函数画图。figure()n=0:7;x1=exp(-j*0.3*pi*n);w=[-600:1:600]*4*pi/600;x=x1*ex

12、p(-j*(n'*w));subplot(2,1,1);plot(w,abs(x));gridon;axis([-15,15,-1,10]),xlabel('w');title('幅度');subplot(2,1,2);plot(w,angle(x));gridon;axis([-15,15,-5,5]),xlabel('w');title('相位');第二次上机3.22x(n)={1,1,1,1},0≤n≤3,计算(DTFT)X(exp(jw)),并画出幅度相位。x