版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、实验2 离散系统的时域分析一、实验目的1、熟悉并掌握离散系统的差分方程表示法;2、加深对冲激响应和卷积分析方法的理解。 二、实验原理在时域中,离散时间系统对输入信号或者延迟信号进行运算处理,生成具有所需特性的输出信号,具体框图如下:其输入、输出关系可用以下差分方程描述:输入信号分解为冲激信号,记系统单位冲激响应,则系统响应为如下的卷积计算式:当时,hn是有限长度的(),称系统为FIR系统;反之,称系统为IIR系统。三、实验内容1 / 71、用MATLAB求系统响应1) 卷积的实现线性移不变系统可由它的单位脉冲响应来表征。若已知了单位脉冲响应和系统激励就可通过卷积运算来求取系统响应,即程序:x=
2、input(Type in the input sequence=); %输入xh=input(Type in the impulse response sequence=); %输入hy=conv(x,h); % 对x,h进行卷积N=length(y)-1; %求出N的值n=0:1:N; %n从0开始,间隔为1的取值取到N为止disp(output sequence=); disp(y); %输出ystem(n,y); %画出n为横轴,y为纵轴的离散图xlabel(Time index n); ylable(Amplitude); % 规定x轴y轴的标签输入为: x=-2 0 1 -1 3
3、h=1 2 0 -1 图形:2) 单位脉冲响应的求取线性时不变因果系统可用MATLAB的函数filter来仿真 y=filter(b,a,x); 其中,x和y是长度相等的两个矢量。矢量x表示激励,矢量a,b表示系统函数形式滤波器的分子和分母系数,得到的响应为矢量y。例如计算以下系统的单位脉冲响应y(n)+0.7y(n-1)-0.45y(y-2)-0.6y(y-3)=0.8x(n)-0.44x(n-1)+0.36x(n-2)+0.02x(n-3) 程序: N=input(Desired impuse response length=);b=input(Type in the vector b=)
4、;a=input(Type in the vector a=);x=1 zeros(1,N-1);y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel(Time index n); ylable(Amplitude);输入:N=41b=0.8 -0.44 0.36 0.02a=1 0.7 -0.45 -0.6图形: 2、以下程序中分别使用conv和filter函数计算h和x的卷积y和y1,运行程序,并分析y和y1是否有差别,为什么要使用xn补零后的x1来产生y1;具体分析当hn有i个值,xn有j个值,使用filter完成卷积功能,需要如何补零?程序:clf;h
5、= 3 2 1 -2 1 0 -4 0 3; %impulse responsex = 1 -2 3 -4 3 2 1; %input sequencey = conv(h,x);n = 0:14;subplot(2,1,1);stem(n,y);xlabel('Time index n'); ylabel('Amplitude');title('Output Obtained by Convolution'); grid;x1 = x zeros(1,8);y1 = filter(h,1,x1);subplot(2,1,2);stem(n,y1
6、);xlabel('Time index n'); ylabel('Amplitude');title('Output Generated by Filtering'); grid; 图形: 因为在y=filter(b,a,x)中,利用给定矢量a和b对x中的数据进行滤波,结果放入y矢量中,y与x长度要相等,所以要使用xn补零后的x1来产生y1。若hn有i个值,xn有j个值,则x1 = x zeros(1,i-1)3、编制程序求解下列两个系统的单位冲激响应,分别用filter 和 impz实现,并绘出其图形。给出理论计算结果和程序计算结果并讨论。第
7、一题:filter实现:程序:N=input('Desired impuse response length=');b=input('Type in the vector b=');a=input('Type in the vector a=');x=1 zeros(1,N-1);y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel('Time index n');ylabel('Amplitude');图形: impz实现:程序:b=input('Type in the v
8、ector b=');a=input('Type in the vector a='); N=25y=impz(b,a,N);k=0:1:N-1;stem(k,y);xlabel('Time index n'); ylabel('Amplitude');图形: 第二题:filter实现:程序:N=input('Desired impuse response length=');b=input('Type in the vector b=');a=input('Type in the vector a
9、=');x=1 zeros(1,N-1);y=filter(b,a,x);k=0:1:N-1;stem(k,y);xlabel('Time index n');ylabel('Amplitude');输入:Type in the vector b=0.25 0.25 0.25 0.25Type in the vector a=1N =25图形: impz实现:程序: b=input('Type in the vector b=');a=input('Type in the vector a='); N=30y=impz(b,a,N);k=0:1:N-1;stem(k,y);xlabel('Time index n'); ylabel('A
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 《电场强度》参考课件
- 3.1农业区位因素及其变化课件2024-2025学年高中地理 人教版(2019)必修第二册
- 2024年上半年资产证券化发展报告-18正式版
- 12.2 第2课时“边角边”
- 商用车保险市场深度分析与投资战略研究报告模板
- 中石油某分公司环氧乙烷乙二醇装置环氧乙烷系统改造项目可行性研究报告
- 北师大版小学数学四年级下册 第六单元测试题(共8套)
- XX集团有限责任公司2024年上半年总结会讲话材料
- 2024年XX镇预警线索核实研判会会议记录
- 低压配电系统运行维护操作手册
- 综合实践课程标准
- 铝材有限公司贴膜机安全风险分级管控清单
- 世界地图(每个国家边界标出-可任意填充颜色)
- 16 开设幼儿保育专业的论证报告
- 藏语常用语及简单对话
- 毕业论文写作指导讲义
- 主轴精度要求及其检测方法(V18下载版)
- 六年级安全教育试题
- 妇产科学女性生殖系统解剖
- WS/T 512-2016医疗机构环境表面清洁与消毒管理规范
- GB/T 22153-2008国际货运代理通用交易条件
评论
0/150
提交评论