资源描述：
《计算机仿真实验报告7.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、山东工商学院计算机仿真及应用实验报告实验七MATLAB的基本应用(二)及Simulink仿真（验证性实验）学院：专业班级：实验时间：学号：姓名：一、实验目的1、掌握连续信号的仿真和傅里叶分析方法2、掌握连续系统的分析方法（时域分析法，拉氏变换法和傅里叶分析法）；3、掌握离散信号的仿真和分析运算方法4、掌握离散系统的分析方法（时域分析法）；5、掌握符号运算方法；6、掌握Simulink仿真工具；二、实验原理1、连续信号的仿真和分析法，参考教材第6.1节，重点：单位冲激信号的仿真方法；单位阶跃信号的仿真方法；复指数信号的仿真方法2、连续系统的分析方法，参考教材第6.1节，重点：例6.2，LTI系
2、统的零输入响应的求解方法；例6.3，LTI系统的冲激响应的求解方法例6.5，LTI系统的零状态响应的求解方法例6.6，系统中有重极点时的计算3、系统的频域分析方法，参考教材第6.2节，重点：例6.7，方波分解为多次正弦波之和例6.8：全波整流电压的频谱例6.10：调幅信号通过带通滤波器例6.12：用傅里叶变换计算滤波器的响应和输出4、离散信号的仿真和分析法，参考教材第6.3节，7.1节，重点：单位脉冲序列impseq，单位阶跃序列stepseq例7.1：序列的相加和相乘例7.2：序列的合成与截取例7.3：序列的移位和周期延拓运算三、实验内容（包括内容，程序，结果）以自我编程练习实验为主，熟悉
3、各种方法和设计，结合课堂讲授，实验练习程序代码。1、根据教材第6.1节的内容，练习连续信号和系统的时域分析和拉氏变换方法。q602clear,clca=input('输入分母系数向量a=[a1,a2,...]=');n=length(a)-1;Y0=input('输入初始条件向量Y0=[y0,Dy0,D2y0,...]=');p=roots(a);V=rot90(vander(p));c=VY0';dt=input('dt=');tf=input('tf=');t=0:dt:tf;y=zeros(1,length(t));fork=1:ny=y+c(k)*exp(p(k)*t);endpl
4、ot(t,y),gridholdon输入分母系数向量a=[a1,a2,...]=[3571]输入初始条件向量Y0=[y0,Dy0,D2y0,...]=[100]dt=0.2tf=8Warning:ImaginarypartsofcomplexXand/orYargumentsignored.>Inq602at9输入分母系数向量a=[a1,a2,...]=[3571]输入初始条件向量Y0=[y0,Dy0,D2y0,...]=[010]dt=0.2tf=8Warning:ImaginarypartsofcomplexXand/orYargumentsignored.>Inq602at9输入分母系
5、数向量a=[a1,a2,...]=[3571]输入初始条件向量Y0=[y0,Dy0,D2y0,...]=[001]dt=0.2tf=8Warning:ImaginarypartsofcomplexXand/orYargumentsignored.>Inq602at9q603clear,clca=input('多项式分母系数向量a=');b=input('多项式分子系数向量b=');[r,p]=residue(b,a),disp('解析式h(t)=r(i)*exp(p(i)*t)')disp('给出时间数组t=[0:dt:tf]')dt=input('dt=');tf=input('tf=')
6、;t=0:dt:tf;h=zeros(1,length(t));fori=1:length(a)-1h=h+r(i)*exp(p(i)*t);endplot(t,h),grid多项式分母系数向量a=poly([0-1+2i-1-2i-2-5])多项式分子系数向量b=[831]r=0.62000.1300-0.3900i0.1300+0.3900i-0.90000.0200p=-5.0000-1.0000+2.0000i-1.0000-2.0000i-2.00000解析式h(t)=r(i)*exp(p(i)*t)给出时间数组t=[0:dt:tf]dt=0.2tf=81、根据教材第6.2节练习傅
7、里叶分析方法。q607t=0:.01:2*pi;y=sin(t);plot(t,y),figure(gcf),pausey=sin(t)+sin(3*t)/3;plot(t,y),pausey=sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/9;plot(t,y)z=zeros(10,max(size(t)));x=zeros(size(t));fork=1:
