资源描述:
《数字信号处理实验(吴镇扬版)matlab程序》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、(1)数组的加、减、乘、除和乘方运算。输入A=[1234],B=[3456],求C=A+B,D=A-B,E=A.*B,F=A./B,G=A.^B并用stem语句画出A、B、C、D、E、F、G。clearall;a=[1234];b=[3456];c=a+b;d=a-b;e=a.*b;f=a./b;g=a.^b;n=1:4;subplot(4,2,1);stem(n,a);xlabel('n');xlim([05]);ylabel('A');subplot(4,2,2);stem(n,b);xlabel('n');xlim([05]);ylabel('B');subplot(4,2,3);
2、stem(n,c);xlabel('n');xlim([05]);ylabel('C');subplot(4,2,4);stem(n,d);xlabel('n');xlim([05]);ylabel('D');subplot(4,2,5);stem(n,e);xlabel('n');xlim([05]);ylabel('E');subplot(4,2,6);stem(n,f);xlabel('n');xlim([05]);ylabel('F');subplot(4,2,7);stem(n,g);xlabel('n');xlim([05]);ylabel('G');(2)用MATLAB实现
3、下列序列:a)x(n)=0.8n0≤n≤15b)x(n)=e(0.2+3j)n0≤n≤15c)x(n)=3cos(0.125πn+0.2π)+2sin(0.25πn+0.1π)0≤n≤15d)将c)中的x(n)扩展为以16为周期的函数x16(n)=x(n+16),绘出四个周期。e)将c)中的x(n)扩展为以10为周期的函数x10(n)=x(n+10),绘出四个周期。clearall;N=0:15;xa=0.8.^N;figure;subplot(2,1,1);stem(N,xa);xlabel('n');xlim([016]);ylabel('xa');xb=exp((0.2+3*j)*
4、N);subplot(2,1,2);stem(N,xb);xlabel('n');xlim([016]);ylabel('xb');figure;xc=3*cos(0.125*pi*N+0.2*pi)+2*sin(0.25*pi*N+0.1*pi);subplot(3,1,1);stem(N,xc);xlabel('n');xlim([016]);ylabel('xc');k=0:3;m=0;fori=1:4forj=1:16m=m+1;n(m)=N(j)+16*k(i);x16(m)=3*cos(0.125*pi*n(m)+0.2*pi)+2*sin(0.25*pi*n(m)+0.1*
5、pi);endendsubplot(3,1,2);stem(n,x16);xlabel('n');ylabel('x16');forj=1:10x10(j)=x16(j);endfori=1:3form=1:10x10(i*10+m)=x10(m);endendn=1:40;subplot(3,1,3);stem(n,x10);xlabel('n');ylabel('x10');(3)x(n)=[1,-1,3,5],产生并绘出下列序列的样本:a)x1(n)=2x(n+2)-x(n-1)-2x(n)b)clearalln=1:4;T=4;x=[1-135];x(5:8)=x(1:4);su
6、bplot(2,1,1);stem(1:8,x);grid;fori=1:4ifi-1<0x1(i)=2*x(i+2)-x(i-1)-2*x(i);elsex1(i)=2*x(i+2)-x(i-1+T)-2*x(i);endendx1(5:8)=x1(1:4);subplot(2,1,2);stem(1:8,x1);grid;(4)绘出下列时间函数的图形,对x轴、y轴以及图形上方均须加上适当的标注:a)x(t)=sin(2πt)0≤t≤10sb)x(t)=cos(100πt)sin(πt)0≤t≤4sta=0:0.05:10;xa=sin(2*pi*ta);subplot(2,1,1);
7、plot(ta,xa);xlabel('t');ylabel('幅度');tb=0:0.01:4;xb=cos(100*pi*tb).*sin(pi*tb);subplot(2,1,2);plot(tb,xb);xlabel('t');ylabel('幅度');(5)编写函数stepshift(n0,n1,n2)实现u(n-n0),n1