信号与系统软件实验.ppt

jsu@uestc.edu.cn qguan@uestc.edu.cn 网络工程系 电子科大通信学院 软件实验 实验一 MATLAB介绍 表示信号、系统的MATLAB函数、工具箱 （实验指导书中实验五）P28 实验二 离散系统的冲激响应、卷积和 （实验指导书中实验六）P34 信号的频谱及频域分析 离散系统的转移函数，零、极点分布和模拟 （实验指导书中实验七）P38 表示信号、系统的MATLAB函数、工具箱 实验目的 加深对常用离散信号的理解； 熟悉表示信号的基本MATLAB函数； cos，sin，exp，imag，real，ones，pi，rand，randn，zeros，sinc，sum 实验原理 实验原理 Ex 已知t的取值范围，画出y=sin(t)的图型。 参考程序： t= -2*pi:0.05:2*pi; y=sin(t); plot(t,y) Ex 已知n的取值范围，画出y=sin(n)的图型。 参考程序： n= -2*pi:1:2*pi; y=sin(t); stem(n,y) Ex 产生单位脉冲序列 和单位阶跃序列u(n-3) 实验内容 离散系统的冲激响应、卷积和 实验目的 加深对离散系统冲激响应、卷积和分析方法的理解 利用MATLAB函数conv、filter计算卷积及系统输出 实验原理 实验原理 实验原理 matlab 中卷积运算的指令是c=conv(a,b),其中a和b是有限长的序列。 ex6:计算下面两个序列的卷积 a=[-2 0 -3 1 –1 3]; b=[1 2 0 -1]; 序列a的非0值从-2开始，序列b的非0值从-1开始 实验内容 用解析方法计算 实验内容 实验内容 信号的频谱及频域分析 实验目的： 1、加强对离散时间傅里叶变换（DTFT）的理解； 2、熟悉离散时间信号的频谱分析方法。 fft，ifft，abs，angle， fftshift，freqz fs=100; n=-1:1/fs:1; x=sin(2*pi*5*n)+2*sin(2*pi*20*n); figure,subplot(2,1,1),plot(n,x) title(时域波形) xlabel(时间 t) N=512; X=fft(x,N); f=(0:N/2-1)/N*fs; absX=abs(X); subplot(2,1,2),plot(f,absX(1:N/2)) title(幅度谱) xlabel(频率 Hz) [H，w]=freqz (b,a,n) 区间上n点等分的频率响应（没有n时，自动给出512点）；w为频率坐标，H为频率响应函数，可以利用plot函数画出频率响应曲线。 计算结果 离散系统的转移函数，零、极点分布和模拟 实验目的 加深对离散系统转移函数、零极点概念的理解 会根据系统转移函数求系统零极点分布 利用MATLAB函数tf2zp、zplane求系统零极点及绘制零极点图 实验原理 实验原理 实验内容 MATLAB软件实验 实验仿真系统的软件界面 表示信号、系统的MATLAB函数、工具箱 离散系统的冲击响应、卷积和 离散系统的转移函数，零、极点分布和模拟 附：评分标准 A：优秀 80～100分 符合实验报告格式要求。掌握系统零极点图的画法；在给定系统常系数差分方程或转移函数的的前提下，能够正确绘制系统零极点图；根据零极点分布，能够大致绘出系统的频率响应。页面整洁、语言通顺、图表清晰、结论正确。能够总结实验结论，得出心得体会。 B：中等 70～80分 符合实验报告格式要求。在给定系统常系数差分方程或转移函数的的前提下，能够正确绘制系统零极点图；根据零极点分布，基本上能够大致绘出系统的频率响应。页面整洁、语言通顺、图表清晰、结论正确。能够总结实验结论，得出心得体会。 C：合格 60～70分 符合实验报告格式要求。在给定系统常系数差分方程或转移函数的的前提下，能够正确绘制系统零极点图。页面整洁、语言通顺、图表清晰、结论正确。有实验结论及心得体会。 ? D：不合格 60分以下 实验报告格式不规范；程序错误较多，零极点图有误。 频域响应为 z域方程为 转移函数为 离 散 系 统 的 时域方程为 将转移函数分解因式： 得到系统的零、极点图。 在MATLAB中，可以用函数[z,p,K]=tf2zp（num,den）求得有理分式形式的系统转移函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统转移函数的零、极点分布图。 num=[0.008 -0.03