matlab上机实验微盘下载,数字信号处理(MATLAB版)上机实验操作.doc

数字信号处理(MATLAB版)上机实验操作.doc

(18页)

本资源提供全文预览，点击全文预览即可全文预览,如果喜欢文档就下载吧，查找使用更方便哦！

9.90 积分

实验一 离散时间信号与系统一、实验目的：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集成系统四、设计流程：此处写个人自己的设计流程五、程序清单：此处写程序内容六、调试和测试结果：此处写程序的执行结果和实验过程中的调试经过、出现的错误和对应的解决方法七、教师评语与成绩评定此处由老师填写上机操作：实验一 离散时间信号与系统实验内容：1．脉冲响应>> b =[1,1]; a = [1,-0.5,0.06];n = [-10:25];>> impz(b,a,n);>> title('Impulse Response'); xlabel('n'); ylabel('h(n)')2.单位阶跃响应>> x = stepseq(0,-10,25); s = filter(b,a,x);Warning: Function call stepseq invokes inexact match d:\MATLAB7\work\STEPSEQ.M.>> stem(n,s)>> title('Step Response'); 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);;xlabel('n');ylabel('s1(n)'); 4． >> 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: Function call conv_m invokes inexact match d:\MATLAB7\work\CONV_M.M.>> stem(m,y);title('系统响应');xlabel('m');ylabel('y(m)');实验二 离散信号与系统的连续频域分析一、实验目的：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函数求频率响应。三、实验平台： 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/pi');ylabel('幅度|X|');>> subplot(2,1,2);plot(w/pi,angX);xlabel('w/pi');ylabel('相位(rad/π)');2>> 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=freqz(b,a,w);>> magH1=abs(H1);angH1=angle(H1);>> subplot(2,2,2);plot(w/pi,angH/pi);grid; xlabel('w (frequency in pi units)');ylabel('相位(rad/π)');>> subplot(2,2,1);plot(w/pi,magH);grid;xlabel('w (frequency in pi units)');ylabel('幅度|H|');>> subplot(2,2,3);plot(w/pi,magH1);grid;xlabel('w (frequency in pi units)');ylabel('幅度|H1|');>> subplot(2,2,4);plot(w/pi,angH1/pi);grid; xlabel('w (frequency in pi units)');ylabel('相位(rad/π)');axis([0,1,-0.8,0]);>> figure(2);zplane(b,a);实验3 离散信号与系统的离散频域分析(DFT)一、实验目的：1、掌握离散时间系统的DFT的MATLAB实现；2、熟悉DTFT和DFT之间的关系。3、了解信号不同变形的DFT与原信号DFT之间的关系二、实验内容：自定义一个长为8点的信号；1、对信号分别做8点、16点、32点DFT，分别与DTFT合并作图并比较DFT与DTFT之间的关系。2、在信号后加零扩展至八点，每两个相邻样本之间插入一个零值，扩充为16点序列，作DFT，画出幅度谱和相位谱，并与原序列的DFT进行比较。3、将信号以8为周期扩展，得到长为16的两个周期，作DFT，画出幅度谱和相位谱，并与原序列的DFT进行比较。三、实验平台： MATLAB集成系统四、设计流程：此处写个人自己的设计流程五、程序清单：此处写程序内容六、调试和测试结果：此处写程序的执行结果和实验过程中的调试经过、出现的错误和对应的解决方法七、教师评语与成绩评定此处由老师填写上机操作1.>> n=0:7;>> x=(0.9*exp(j*pi。省略部分。设计指标：wp=0.3*pi;ws=0.4*pi;Rp=1;As=40;要点：设计指标要与结果图一致七、教师评语与成绩评定此处由老师填写上机操作模拟原型Butterworth滤波器：>> wp=0.4*pi;ws=0.6*pi;Rp=0.5;As=50;>> [b,a]=afd_butt(wp,ws,Rp,As);*** Butterworth Filter Order = 17 >> [db,mag,pha,w]=freqs_m(b,a,pi);>> [ha,x,t]=impulse(b,a,100);>> subplot(2,2,1);plot(w/pi,mag);xlabel('w');ylabel('|H|');grid on>> subplot(2,2,2);plot(w/pi,pha);xlabel('w');ylabel('相位'); grid on>> subplot(2,2,3);plot(w/pi,db);xlabel('w');ylabel('db'); grid on>> subplot(2,2,4);plot(t,ha);xlabel('t');ylabel('ha(t)'); axis([0,100,-0.2,0.4]);grid on双线性变换法求得的数字Butterworth低通滤波器：>> wp=0.4*pi;ws=0.6*pi;Rp=0.5;As=50;>> T=2;OmegaP=(2/T)*tan(wp/2);OmegaS=(2/T)*tan(ws/2);>> N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(OmegaP/OmegaS)));>> OmegaC=OmegaP/((10^(Rp/10)-1)^(1/(2*N)));wn=2/pi*atan(OmegaC*T/2);>> [b,a]=butter(N,wn);>> [db,mag,pha,grd,w] = freqz_m(b,a);>> n=0:100;>> h=impz(b,a,n);>> subplot(2,2,1);plot(w/pi,mag);xlabel('w');ylabel('|H|'); axis([0,1,0,1.15]);grid on>> subplot(2,2,2);plot(w/pi,pha);xlabel('w');ylabel('相位'); axis([0,1,-4,4]); grid on>> subplot(2,2,3);plot(w/pi,db);xlabel('w');ylabel('db'); axis([0,1,-140,20]); grid on>> subplot(2,2,4);stem(n,h);xlabel('n');ylabel('h(n)'); axis([0,100,-0.2,0.4]);grid on由前面的数字Butterworth低通滤波器映射得到的高通滤波器：接上面的程序>> wp1=ws;ws1=wp;>> [N1,wn1]=buttord(wp1/pi,ws1/pi,Rp,As);>> [b1,a1]=butter(N1,wn1,'high');>> [db1,mag1,pha1,grd1,w1] = freqz_m(b1,a1);>> h1=impz(b1,a1,n);>> subplot(2,2,1);plot(w1/pi,mag1);xlabel('w1');ylabel('|H1|');axis([0,1,0,1.15]);grid on>> subplot(2,2,2);plot(w1/pi,pha1);xlabel('w1');ylabel('相位'); axis([0,1,-4,4]); grid on>> subplot(2,2,3);plot(w1/pi,db1);xlabel('w1');ylabel('db1'); axis([0,1,-140,20]); grid on>> subplot(2,2,4);stem(n,h1);xlabel('n');ylabel('h1(n)'); grid on信号x(n)通过数字Butterworth低通滤波器：接上面程序>> w1=0.3*pi;w2=0.7*pi;>> x=cos(w1*n)+cos(w2*n);>> [y,ny]=conv_m(x,n,h',n);>> X=dtft(x,n,w);>> Y=dtft(y,ny,w);>> magX=abs(X);magY=abs(Y);figure(2);>> subplot(2,2,1);stem(n,x);xlabel('n');ylabel('x(n)');grid on>> subplot(2,2,2);plot(w/pi,magX);xlabel('w');ylabel('幅度|X|'); grid on>> subplot(2,2,3);stem(ny,y);xlabel('n');ylabel('y(n)'); grid on>> subplot(2,2,4);plot(w/pi,magY);xlabel('w');ylabel('幅度|Y|'); grid on(1)x=exp(-0.9*n);w=0:pi/500:2*pi;X=x*exp(-j*n'*w);magX=abs(X);angX=angle(X);x1=[x zeros(1,8)];x2=[x1 zeros(1,16)];Xk=dft(x,8);Xk1=dft(x1,16);Xk2=dft(x2,32);magXk=abs(Xk);angXk=angle(Xk);magXk1=abs(Xk1);angXk1=angle(Xk1);magXk2=abs(Xk2);angXk2=angle(Xk2);subplot(4,2,1);plot(w/pi,magX);subplot(4,2,2);plot(w/pi,angX);subplot(4,2,3);plot(4*w/pi,magX); hold ; stem(0:7,magXk);subplot(4,2,5);plot(8*w/pi,magX); hold ; stem(0:15,magXk1);subplot(4,2,7);plot(16*w/pi,magX); hold ; stem(0:31,magXk2);subplot(4,2,4);plot(4*w/pi,angX); hold ; stem(0:7,angXk);subplot(4,2,6);plot(8*w/pi,angX); hold ; stem(0:15,angXk1);subplot(4,2,8);plot(16*w/pi,angX); hold ; stem(0:31,angXk2); 关 键 词： 数字信号 matlab 操作 实验 处理 上机

 天天文库所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。

关于本文