资源描述：
《数字信号处理(matlab版)上机实验操作》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验一离散时间信号与系统一、实验目的：1、熟悉常见离散时间信号的产生方法；2、熟悉离散时间系统的单位脉冲响应和单位阶跃响应的求解方法；3、熟悉离散时间信号经过离散时间系统的响应的求解方法。二、实验内容：已知离散时间系统差分方程为y(n)-0.5y(n-1)+0.06y(n-2)=x(n)+x(n-1)，求1、该系统的单位脉冲响应并绘图；2、该系统的单位阶跃响应并绘图；3、已知x(n)=可自己指定用filter函数经过系统的响应并绘图；4、用conv_m函数求系统响应并绘图。三、实验平台：MATLAB集成系统四、
2、设计流程：此处写个人自己的设计流程五、程序清单：此处写程序内容六、调试和测试结果：此处写程序的执行结果和实验过程中的调试经过、出现的错误和对应的解决方法七、教师评语与成绩评定此处由老师填写上机操作：实验一离散时间信号与系统实验内容：1．脉冲响应>>b=[1,1];a=[1,-0.5,0.06];n=[-10:25];>>impz(b,a,n);>>title('ImpulseResponse');xlabel('n');ylabel('h(n)')172.单位阶跃响应>>x=stepseq(0,-10,25);
3、s=filter(b,a,x);Warning:Functioncallstepseqinvokesinexactmatchd:MATLAB7workSTEPSEQ.M.>>stem(n,s)>>title('StepResponse');xlabel('n');ylabel('s(n)')3.>>a=[1,-0.5,0.06];b=[1,1];>>n=-20:120;>>x1=exp(-0.05*n).*sin(0.1*pi*n+pi/3);>>s1=filter(b,a,x1);>>stem(n,s1
4、);;xlabel('n');ylabel('s1(n)');174．>>a=[1,-0.5,0.06];b=[1,1];>>n=-20:120;>>h=impz(b,a,n);>>x1=exp(-0.05*n).*sin(0.1*pi*n+pi/3);>>[y,m]=conv_m(x1,n,h,n);Warning:Functioncallconv_minvokesinexactmatchd:MATLAB7workCONV_M.M.>>stem(m,y);title('系统响应');xlabel('m'
5、);ylabel('y(m)');17实验二离散信号与系统的连续频域分析一、实验目的：1、掌握离散时间信号的DTFT的MATLAB实现；2、掌握离散时间系统的DTFT分析；3、掌握系统函数和频率相应之间的关系。二、实验内容：1、自定义一个长度为8点的信号，信号幅度值也由自己任意指定，对该信号作DTFT，分别画出幅度谱和相位谱2、已知离散时间系统差分方程为y(n)-0.5y(n-1)+0.06y(n-2)=x(n)+x(n-1)，求出并画出其频率响应3、求该系统系统函数，并画极零点图，并通过freqz函数求频率响
6、应。三、实验平台：MATLAB集成系统四、设计流程：此处写个人自己的设计流程五、程序清单：此处写程序内容六、调试和测试结果：此处写程序的执行结果和实验过程中的调试经过、出现的错误和对应的解决方法七、教师评语与成绩评定此处由老师填写上机操作1>>n=0:7;>>x=(0.9*exp(j*pi/3)).^n;>>w=-pi:pi/200:pi;>>X=dtft(x,n,w);>>magX=abs(X);angX=angle(X);>>subplot(2,1,1);plot(w/pi,magX);xlabel('w/
7、pi');ylabel('幅度
8、X
9、');>>subplot(2,1,2);plot(w/pi,angX);xlabel('w/pi');ylabel('相位(rad/π)');172>>a=[1,-0.5,0.06];b=[1,1,0];>>m=0:length(b)-1;l=0:length(a)-1;>>w=0:pi/500:pi;>>num=b*exp(-j*m'*w);>>den=a*exp(-j*l'*w);>>H=num./den;>>magH=abs(H);angH=angle(H);>>H1=
10、freqz(b,a,w);>>magH1=abs(H1);angH1=angle(H1);>>subplot(2,2,2);plot(w/pi,angH/pi);grid;xlabel('w(frequencyinpiunits)');ylabel('相位(rad/π)');>>subplot(2,2,1);plot(w/pi,magH);grid;xlabel('w(frequenc
