14271039商都连续系统响应的MATLAB求解

《14271039商都连续系统响应的MATLAB求解》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、信号与系统》课程研究性学习报告姓名]j指导教师薛健时间2015.12学号14271039信号与系统的时域分析专题研讨【目的】(1)加深对信号与系统时域分析基本原理和方法的理解。(2)学会利用MATLAB进行信号与系统的分析。(3)培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。【研讨题目】题目1连续系统响应的MATLAB求解一连续LTI系统满足的微分方程为心)+3心)+2曲)=2")+班/)(1)已知x(r)=e_3/w(r)，试求该系统的零状态响应y』)。解：yzs(t)=-0.5e-t+3e~2t-2.5e~3t;(2)用lsim求出该系统的

2、零状态响应的数值解。利用(1)所求得的结果，比较不同的抽样间隔对数值解精度的影响。解：clear;sys=tf([21],[132]);%调用tf函数获得连续时间LTI模型te=0.01;t=0:te:5;%确定抽样间隔tel=0.2;tl=0:tel:5;te2=0.4;t2=0:te2:5;te3=0.5;t3=0:te3:5;y=lsim(sys,exp(-3*t).*(t>=0),t);%调用lsim函数获得系统的零状态响应数值解yl=lsim(sys,exp(-3*tl).*(tl>=0),tl);%改变抽样间隔y2=lsim(sys,exp(-3*

3、t2)・*(t2>=0),t2);y3=lsim(sys,exp(-3*t3).*(t3>=0),t3);yr=-0.5*exp(T*t)+3*exp(-2*t)-2.5*exp(-3*t);%系统零状态响应真实解plot(t,y);%数值解t-y图像axis([O5-0.10.6]);xlabel('Time(sec)');holdon;plot(t,yr);%真实解t-yr图像plot(tl,yl);不同抽样间隔的图像plot(t2,y2);plot(t3,y3);legend('y(t)','yr(t)','te=0.2','te=0.4','te=0.

4、5');holdoff;%由图像可知，抽样间隔增大,结果精度降低(3)用命令[x,Fsam,bits]=wavread(1Yourn1):将硬盘上的语音文件Yourn.wav读入计算机。用命令sound(x,Fsam)；播放该语音信号；解：clear;[x,Fsam,bits]=wavreadYourn');sound(x,Fsam);(1)用命令loadmodel01将磁盘文件modelOl.mat读入计算机后，MATLAB的workspace屮将会新增变量den和numoden表示微分方程左边的系数，变量num表示微分方程右边的系数。写出磁盘文件model

5、Ol.mat定义的微分方程；解：clear;loadmodelOl;%引入变量den;%显示数组dennum;%显示数组num(1)计算(3)中的信号通过(4)中系统的响应，播放系统输出的语音信号。与处理前的信号比较，信号有何不同？能用已学知识解释所得结果吗？解：clear;[x,Fsam,bits]=wavread(JYourn');%读取声音数据[数据，采样频率，采样深度]N=length(x);%计算抽取后数据点数t=(O:N-1)/Fsam;%计算数据时长loadmodelOl;%引入变量sys=tf(num,den);%建立连续时间LTI模型y=ls

6、im(sys,x,t);%获得音频信号通过系统的响应ysound(y,Fsam);subplot(211)；%绘制原信号图形plot(t,x);subplot(212);%绘制处理后信号图形plot(t,y);