离散系统的时域及变换域分析.doc_第1页
离散系统的时域及变换域分析.doc_第2页
离散系统的时域及变换域分析.doc_第3页
离散系统的时域及变换域分析.doc_第4页
离散系统的时域及变换域分析.doc_第5页
已阅读5页,还剩3页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

实验1 离散系统的时域及变换域分析一、实验目的:1.加深对离散系统的差分方程、单位抽样响应和卷积分析方法的理解。2.加深对离散系统的频率响应分析和零、极点分布的概念理解。二、实验原理:1.时域离散系统Discrete-time systmex(n)y(n)其输入、输出关系可用以下差分方程描述:输入信号分解为冲激信号,系统单位抽样序列h(n),则系统响应为如下的卷积计算式:当时,h(n)是有限长度的(n:0,M),称系统为FIR系统;反之,称系统为IIR系统。 在MATLAB中,可以用函数y=filter(b,a,x)实现差分方程的仿真,也可以用函数 y=conv(x,h)计算卷积。2.变换域离散系统的时域方程为其变换域分析方法如下:系统函数为 分解因式 ,其中 和 称为零、极点。 在MATLAB中,可以用函数z,p,K=tf2zp(num,den)求得有理分式形式的系统函数的零、极点,用函数zplane(z,p)绘出零、极点分布图;也可以用函数zplane(num,den)直接绘出有理分式形式的系统函数的零、极点分布图。使用h=freqz(num,den,w)函数可求系统的频率响应,w是频率的计算点,如w=0:pi/255:pi, h是复数,abs(h)为幅度响应,angle(h)为相位响应。另外,在MATLAB中,可以用函数 r,p,k=residuez(num,den)完成部分分式展开计算;可以用函数sos=zp2sos(z,p,K)完成将高阶系统分解为2阶系统的串联。三 、实验内容1.时域(1.)编制程序求解下列系统的单位抽样响应,并绘出其图形。解 用MATLAB计算程序如下:N=15;n=0:N-1;b=1,-1;a=1,0.75,0.125;x=n=0;y=filter(b,a,x);subplot(3,2,1);stem(n,y,.);axis(0,N,-1,2);ylabel(y(n); (2.)给定因果稳定线性时不变系统的差分方程 对下列输入序列,求输出序列。解:MATLAB计算程序如下:N=80;n=0:N-1;b=1;a=1,-1,0.9;x=(n0&(n30);y=filter(b,a,x);subplot(3,2,3);stem(n,y,.);axis(0,N,-1,2);ylabel(y(n);2.变换域例1 求下列直接型系统函数的零、极点,并将它转换成二阶节形式 解 用MATLAB计算程序如下: num=1 -0.1 -0.3 -0.3 -0.2;den=1 0.1 0.2 0.2 0.5;z,p,k=tf2zp(num,den);disp(零点);disp(z);disp(极点);disp(p);disp(增益系数);disp(k);sos=zp2sos(z,p,k);disp(二阶节);disp(real(sos);zplane(num,den)输入到“num”和“den”的分别为分子和分母多项式的系数。计算求得零、极点增益系数和二阶节的系数:零点 0.9615 -0.5730 -0.1443 + 0.5850i -0.1443 - 0.5850i 极点 0.5276 + 0.6997i 0.5276 - 0.6997i -0.5776 + 0.5635i -0.5776 - 0.5635i 增益系数 1 二阶节 1.0000 -0.3885 -0.5509 1.0000 1.1552 0.6511 1.0000 0.2885 0.3630 1.0000 -1.0552 0.7679 系统函数的二阶节形式为: 极点图如右图。 例2 差分方程 所对应的系统的频率响应。 解:差分方程所对应的系统函数为用MATLAB计算的程序如下:k=256;num=0.8 -0.44 0.36 0.02;den=1 0.7 -0.45 -0.6;w=0:pi/k:pi;h=freqz(num,den,w);subplot(2,2,1);plot(w/pi,real(h);gridtitle(实部)xlabel(omega/pi);ylabel(幅度)subplot(2,2,2);plot(w/pi,imag(h);gridtitle(虚部)xlabel(omega/pi);ylabel(Amplitude)subplot(2,2,3);plot(w/pi,abs(h);gridtitle(幅度谱)xlabel(omega/pi);ylabel(幅值)subplot(2,2,4);plot(w/pi,angle(h);gridtitle(相位谱)xlabel(omega/pi);ylabel(弧度) 练习1.求系统 的零、极点和幅度频率响应和相位响应。要求:绘出零、极点分布图,幅度频率响应和相位响应曲线。 解:用MATLAB计算的程序如下:num=0.0528 0.0797 0.1295 0.1295 0.797 0.0528;den=1 -1.8107 2.4947 -1.8801 0.9537 -0.2336;z,p,k=tf2zp(num,den);disp(零点);disp(z);disp(极点);disp(p);零点 -1.5870 + 1.4470i -1.5870 - 1.4470i 0.8657 + 1.57795i 0.8657 - 1.5779i -0.0669 极点 0.2788 + 0.8973i 0.2788 - 0.8973i 0.3811 + 0.6274i 0.3811 - 0.6274i 0.4910 k=256;num=0.05280.07970.12950.12950.7970.0528;den=1-1.81072.4947-1.88010.9537-0.2336;w=0:pi/k:pi;h=freqz(num,den,w);subplot(2,2,1);plot(w/pi,real(h);gridtitle(幅度谱)xlabel(omega/pi);ylabel(幅值)subplot(2,2,4);plot(w/pi,angle(h);gridtitle(相位谱)xlabel(omega/pi);ylabel(弧度)四、实验结果分析1、系统函数的零、极点分别关于实轴和原点对称分布2、对于稳定的因果系统,H(z)的全部极点应落在单位圆内,所以描述的系统是稳定的因果系统3、通过Matlab,可以直观的看出系统函数的幅度和相位谱

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论