资源描述:
《数字信号处理实验答案.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、实验一熟悉Matlab环境一、实验目的1.熟悉MATLAB的主要操作命令。2.学会简单的矩阵输入和数据读写。3.掌握简单的绘图命令。4.用MATLAB编程并学会创建函数。5.观察离散系统的频率响应。二、实验内容认真阅读本章附录,在MATLAB环境下重新做一遍附录中的例子,体会各条命令的含义。在熟悉了MATLAB基本命令的基础上,完成以下实验。上机实验内容:(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=
2、[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);stem(n,c);xlabel('n');xlim([05]);ylabel('C');subplot(4,2,4);stem(n,d);xlabel('n');xlim([05]);ylabel('D');subplot(4
3、,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实现下列序列: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为
4、周期的函数x16(n)=x(n+16),绘出四个周期。e)将c)中的x(n)扩展为以10为周期的函数x10(n)=x(n+10),绘出四个周期。clearall;N=0:15;%a)x(n)=0.8n0≤n≤15xa=0.8.^N;-37-figure;subplot(2,1,1);stem(N,xa);xlabel('n');xlim([016]);ylabel('xa');%b)x(n)=e(0.2+3j)n0≤n≤15xb=exp((0.2+3*j)*N);subplot(2,1,2);stem(N,xb);xlabel('n');xlim([016]);ylabel('xb'
5、);figure;%c)x(n)=3cos(0.125πn+0.2π)+2sin(0.25πn+0.1π)0≤n≤15xc=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');%d)将c)中的x(n)扩展为以16为周期的函数x16(n)=x(n+16),绘出四个周期。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(
6、m)+0.2*pi)+2*sin(0.25*pi*n(m)+0.1*pi);endendsubplot(3,1,2);stem(n,x16);xlabel('n');ylabel('x16');%e)将c)中的x(n)扩展为以10为周期的函数x10(n)=x(n+10),绘出四个周期。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');-37-(3)x(n)=[1,-1,3,5],产
7、生并绘出下列序列的样本: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);subplot(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