资源描述：
《数字信号处理.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、实验一熟悉MATLAB环境一、实验目的(1)熟悉MATLAB的主要操作命令。(2)学会简单的矩阵输入和数据读写。(3)掌握简单的绘图命令。(4)用MATLAB编程并学会创建函数。(5)观察离散系统的频率响应。二、实验内容上机实验内容:(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。A=[1234];stem(A)B=[3356];stem(B)A=[1234];B=[3356];C=A+Bstem(C)A=[1234];B=[34
2、56];D=A-Bstem(D)A=[1234];B=[3456];E=A.*Bstem(E)A=[1234];B=[3456];F=A./Bstem(F)A=[1234];B=[3456];G=A.^Bstem(G)(2)用MATLAB实现下列序列:n=[0:15];x1=0.8.^n;stem(x1)n=[0:15];x2=exp((0.2+3*j)*n);stem(x2)n=[0:1:15];x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);stem(x3)d)将(c)中的x(n)扩展为以16为周期的函数,绘出四
3、个周期。n=[0:1:63];x4=3*cos(0.125*pi*rem(n,16)+0.2*pi)+2*sin(0.25*pi*rem(n,16)+0.1*pi);stem(x4)e)将(c)中的x(n)扩展为以10为周期的函数，绘出四个周期。n=[0:1:39];x5=3*cos(0.125*pi*rem(n,10)+0.2*pi)+2*sin(0.25*pi*rem(n,10)+0.1*pi);stem(x5)(3)产生并绘出下列序列的样本:n=0:3;x=[1-135];x1=circshift(x,[0,-2]);x2=circshift(x,[0,1]);x3
4、=2*x1-x2-2*x;stem(x3);n=0:3;x=[1-135];x1=circshift(x,[0,1]);x2=circshift(x,[0,2]);x3=circshift(x,[0,3]);x4=circshift(x,[0,4]);x5=circshift(x,[0,5]);xn=1*x1+2*x2+3*x3+4*x4+5*x5;stem(xn)(4)绘出下列时间函数的图形，对x轴、y轴以及图形上方均须加上适当的标注t=0:0.001:10;x=sin(2*pi*t);plot(t,x)title('sin(2*pi*t)');xlabel('x');
5、ylabel('t');t=0:0.01:4;x1=cos(100*pi*t);x2=sin(pi*t);x=x1.*x2;plot(t,x)title('x=cos(100*pi*t)*sin(pi*t)');xlabel('x');ylabel('t');(5)编写函数实现，绘出该函数的图形，起点为n1，终点为n2。>>n1=input('请输入起点')请输入起点3n1=3>>n2=input('请输入终点')请输入终点9n2=9>>n0=input('请输入移动距离')请输入移动距离6n0=6>>y(1:n0)=0;y(n0:n2)=1;stem(y)(6)给定一因
6、果系统求出并绘制H(z)的幅频响应与相频响应。b=[1sqrt(2)];a=[1-0.670.9];w=0:pi/100:pi;h=freqz(b,a,w);subplot(211)plot(abs(h));gridonsubplot(212);plot(angle(h))gridon(底表格)(7)计算序列和序列的离散卷积，并作图表示卷积结果。A=[8-2-123]B=[23-1-1];C=conv(A,B)stem(C)(8)求以下差分方程所描述系统的单位脉冲响应N=51;num=[1-2];den=[10.1-0.62];x=[1zeros(1,N-1)];n=0:
7、1:N-1;G=filter(num,den,x);stem(n,G)AXIS([-153-2.51.2])三、思考题(1)比较实验内容第(2)题中d)和e)两小题的结果，试说明对于周期性信号，应当如何采样，d能保证周期扩展后与原信号保持一致?答：对于周期性信号，在进行采样时，其采样周期必须满足采样定理，即采样频率应该大于信号最高频率的两倍，这样才能避免迭混，以便采样后仍能准确的恢复原信号。(2)对于有限长序列，如何用MATLAB计算其DTFT?答：用函数freqz可以计算序列在给定的离散频率点上的DTFT，该变换序列是以形如
