《实验一系统响应及系统稳定性实验报告.pdf》由会员分享，可在线阅读，更多相关《实验一系统响应及系统稳定性实验报告.pdf（3页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、一、实验目的 1掌握求系统响应的方法 2掌握时域离散系统的时域特性 3分析、观察及检验系统的稳定性 二、实验原理与方法 在时域中,描写系统特性的方法是差分方程和单位脉冲响应;已知输入信号,可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应,本实验仅在时域求解;在计算机上适合用递推法求差分方程的解,最简单的方法是采用MATLAB语言的工具箱函数filter函数;也可以用MATLAB语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积,求出系统的响应;系统的稳定性是指对任意有界的输入信号,系统都能得到有界的系统响应;或者系统的单位脉冲响应满足绝对可和的条件;系统
2、的稳定性由其差分方程的系数决定;实际中检查系统是否稳定,不可能检查系统对所有有界的输入信号,输出是否都是有界输出,或者检查系统的单位脉冲响应满足绝对可和的条件;可行的方法是在系统的输入端加入单位阶跃序列,如果系统的输出趋近一个常数包括零,就可以断定系统是稳定的;系统的稳态输出是指当n时,系统的输出;如果系统稳定,信号加入系统后,系统输出的开始一段称为暂态效应,随n的加大,幅度趋于稳定,达到稳态输出;注意在以下实验中均假设系统的初始状态为零;二、实验内容及步骤 1 编制程序,包括产生输入信号、单位脉冲响应序列的子程序,用 filter 函数或conv 函数求解系统输出响应的主程序;程序中要有绘制
3、信号波形的功能;程序代码 xn=ones1,32;hn=0.2 0.2 0.2 0.2 0.2;yn=convhn,xn;n=0:lengthyn-1;subplot2,2,1;stemn,yn,.titleayn波形;xlabeln;ylabelyn 输出波形 2 给定一个低通滤波器的差分方程为 输入信号)()(81nRnx 分别求出系统对)()(81nRnx和)()(2nunx的响应序列,并画出其波形;求出系统的单位冲响应,画出其波形;%内容 1：调用 filter 解差分方程,由系统对 un 的响应判断稳定性%=A=1,-0.9;B=0.05,0.05;%系统差分方程系数向量 B 和 A
4、 x1n=1 1 1 1 1 1 1 1 zeros1,50;%产生信号 x1n=R8n x2n=ones1,128;%产生信号 x2n=un hn=impzB,A,58;%求系统单位脉冲响应 hn subplot2,2,1;y=hn;tstemhn,y;%调用函数 tstem 绘图 titlea 系统单位脉冲响应 hn;box on y1n=filterB,A,x1n;%求系统对 x1n 的响应 y1n subplot2,2,2;y=y1n;tstemy1n,y;titleb 系统对 R8n 的响应 y1n;box on y2n=filterB,A,x2n;%求系统对 x2n 的响应 y2n
5、 subplot2,2,4;y=y2n;tstemy2n,y;titlec 系统对 un 的响应 y2n;box on 3 给定系统的单位脉冲响应为 用线性卷积法分别求系统 h1n 和 h2n 对)()(81nRnx的输出响应,并画出波形;%内容 3：调用 conv 函数计算卷积%=x1n=1 1 1 1 1 1 1 1;%产生信号 x1n=R8n h1n=ones1,10 zeros1,10;h2n=1 2.5 2.5 1 zeros1,10;y21n=convh1n,x1n;y22n=convh2n,x1n;figure2 subplot2,2,1;y=h1n;tstemh1n,y;%调用
6、函数 tstem 绘图 titled 系统单位脉冲响应 h1n;box on subplot2,2,2;y=y21n;tstemy21n,y;titlee h1n 与 R8n 的卷积 y21n;box on subplot2,2,3;y=h2n;tstemh2n,y;%调用函数 tstem 绘图 titlef 系统单位脉冲响应 h2n;box on subplot2,2,4;y=y22n;tstemy22n,y;titleg h2n 与 R8n 的卷积 y22n;box on 4 给定一谐振器的差分方程为 yn=1.8237yn-1-0.9801yn-2+b0 xn-b0 xn-2 令 b0=
7、1/100.49,谐振器的谐振频率为 0.4 rad;实验方法检查系统是否稳定;输入信号为 un 时,画出系统输出波形 y31n;给定输入信号为 xn=sin0.014n+sin0.4n,求出系统的输出响应 y32n,并画出其波形;%内容 4：谐振器分析%=un=ones1,256;%产生信号 un n=0:255;xsin=sin0.014n+sin0.4n;%产生正弦信号 A=1,-1.8237,0.9801;B=1/100.49,0,-1/100.49;%系统差分方程系数向量 B 和 A y31n=filterB,A,un;%谐振器对 un 的响应 y31n y32n=filterB,A
8、,xsin;%谐振器对 un 的响应 y31n figure3 subplot2,1,1;y=y31n;tstemy31n,y;titleh 谐振器对 un 的响应 y31n;box on subplot2,1,2;y=y32n;tstemy32n,y;titlei 谐振器对正弦信号的响应 y32n;box on 四、实验结果分析 由各实验结果的截图可看出,每个图都直观地反映了我们想要求得的单位脉冲响应、给定信号作用后的输出响应,都符合预期结果;五、思考题 1 如果输入信号为无限长序列,系统的单位脉冲响应是有限长序列,可否用线性卷积法求系统的响应 如何求 答：可以;把输入信号进行分段,分别进行卷积,最后将各段卷积结果相加即可;2 如果信号经过低通滤波器,把信号的高频分量滤掉,时域信号会有何变化 用前面第一个实验结果进行分析说明;答：时域信号的剧烈变化将被平滑,由实验内容 1 的内容可见,经过系统的低通滤波使输入信号和输出的阶跃变化变得缓慢上升与下降;六、实验心得及体会 通过本次实验我重新温习了MATLAB 这个软件的基本使用方法,运行环境;通过这款软件使我们的学习更加方便;实验中,我学会了filter和conv函数的基本用法,前者可计算知道输入信号的前提下求解输出响应的序列,后者则可通过输入信号和系统的单位脉冲响应的线性卷积,求出系统的响应;