数字信号处理(MATLAB版)上机实验操作..doc
文本预览下载声明
实验一 离散时间信号与系统
一、实验目的：
二、实验内容：
三、实验平台： MATLAB集成系统
四、设计流程：五、程序清单：、结果：
七、教师评语与成绩评定
此处由老师填写
上机操作：实验一 离散时间信号与系统
实验内容 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));
实验二 离散信号与系统的连续频域分析
一、实验目的：
MATLAB实现；
2、掌握离散时间系统的DTFT分析；
3、掌握系统函数和频率相应之间的关系。
二、实验内容：
三、实验平台： 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）
一、实验目的：
MATLAB实现；
2、熟悉DTFT和DFT之间的关系。
3、了解信号不同变形的DFT与原信号DFT
显示全部