《系统响应及系统稳定性_论文-毕业文章.pdf》由会员分享，可在线阅读，更多相关《系统响应及系统稳定性_论文-毕业文章.pdf（11页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、.word.实验一及课堂作业 实验一：系统响应及系统稳定性 一、实验原理与方法 1、在时域求系统响应的方法有两种：第一种是通过解差分方程求得系统输出；第二种是系统的单位脉冲响应，通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出。2、检验系统的稳定性，其方法是在输入端参加单位阶跃序列，观察输出波形，如果波形稳定在一个常数值包括零上，系统稳定，否那么不稳定。3、系统的频域特性包括传输函数/特性系统单位脉冲响应的傅里叶变换幅频、相频、系统函数/特性系统单位脉冲响应的 Z 变换、零极点分布等。分析系统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果，因为零、极点分布对系统的频域特性有
2、影响，通过控制系统函数的零、极点分布就可以设计出不同特性需求的系统。二、实验容 1、编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用 filter 函数或 conv 函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。2、给定一个低通滤波器的差分方程为)1(9.0)1(05.0)(05.0)(nynxnxny 输入信号)()(),()(281nunxnRnx 1分别求出)(8)(1nRnx和)()(2nunx)的系统响应，并画出其波形。.word.2求出系统的单位脉冲响应，画出其波形。程序见附录 1.1、实验结果见图 1.1。3、给定系统的单位脉冲响应为 )()(101nRnh
3、)3()2(5.2)1(5.2)()(2nnnnnh 用线性卷积法求)(8)(1nRnx)分别对系统)(1 nh和)(2 nh)的输出响应并画出波形。程序见附录 1.2、实验结果见图 1.2。4、给定一谐振器的差分方程为 )2()()2(9801.0)1(8237.1)(00nxbnxbnynyny 令49.10010b，谐振器的谐振频率为rad4.0。1用实验方法检查系统是否稳定，输入信号为)(nu时，画出系统输出波形。2给定输入信号为 )4.0sin()014.0sin()(nnnx 求出系统的输出响应，并画出其波形。程序见附录 1.3、实验结果见图 1.3。三、实验结果和分析、讨论及结论
4、 1、实验结果：解差分方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果因为零极点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和
5、的输出响应并画出波形程.word.051015202530354000.020.040.060.080.1(a)系统单位脉冲响应h(n)0510152025303540455000.20.40.60.8(b)系统对R8(n)的响应y1(n)010203040506000.20.40.60.81(c)系统对u(n)的响应y2(n)图 1.1 依据差分方程求取系统脉冲响应和输出响应 实验分析、讨论及结论：（a）中 25 个点数和程序所写一致。Filter 函数实现线性常系数差分方程的递推求解，调用格式如下：Y=filter(B,A,x)*计算系统对输入信号 x 的零状态响应输出信号向量 Y，B、A
6、 是差分方程的系数向量。即 B=a1,a2am A=b1,b2bn 2、实验结果：解差分方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果因为零极点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见
7、附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应并画出波形程.word.0510152000.20.40.60.81(d)系 统 单 位 脉 冲 响 应 h1(n)05101520253002468(e)h1(n)与 R8(n)的 卷 积 y21(n)0246810121400.511.522.5(f)系 统 单 位 脉 冲 响 应 h2(n)051015202501234567(g)h2(n)与 R8(n)的 卷 积 y22(n)图 1.2 线性卷积求取二个不同系统输出响应 实验分析、讨论及结论：df单位脉冲响应点数与程序要求一致；eg卷积点数满足M+N-1 的要
8、求，图形也满足要求。Conv 函数用于计算两个有限长序列的卷积；C=convA,B计算两个有限长序列向量 A 和 B 的卷积 3、实验结果：解差分方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果因为零极点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的
9、系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应并画出波形程.word.050100150200250300-0.0500.05(h)谐 振 器 对 u(n)的 响 应 y31(n)050100150200250300-1-0.500.51(i)谐 振 器 对 正 弦 信 号 的 响 应 y32(n)图 1.3 依据差分方程求取系统脉冲响应和输出响应 实验分析、讨论及结论：在系统的输入端参加单位阶跃序列，如果系统的输出趋近一个常数包括零，就可以断定系统是稳定的，h中的输出显然趋近于零，所以系统是稳定的。在i中
10、，谐振器具有对某个频率进展谐振的性质，本实验中的谐振器的谐振频率是 0.4 rad,因此稳定波形为sin(0.4n)。四、思考题 1、如果输入信号为无限长序列，系统的单位脉冲响应是有限长序列，可否用线性卷积法求系统的响应?如何求？答：如果输入信号为无限长序列，系统的单位脉冲响应是有限长序列，可用解差分方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果因为零极
11、点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应并画出波形程.word.分段线性卷积法求系统的响应。具体步骤是对输入信号序列分段；求单位脉冲响应 hn与各段的卷积；将各段卷积结果相加。具体实现方法有第三章P91 介绍的重叠相加法和重叠保存法。2、如果信号经过低通滤波器，把信号的高频分量滤掉，时域信号会有何变化，用前面第一个实验结果进展分析说明。答：如果信号
12、经过低通滤波器，那么信号的高频分量将被过滤掉，时域信号的剧烈变化变的平滑，在有阶跃处附近产生过渡带。由实验容 1 的结果图可见，经过系统低通滤波使得输入信号)()()(81nRnxn、和)()(2nunx的阶跃变化变的缓慢上升与下降。五、总结与心得体会 实验总结即在实验原理中说明的两点：1、在时域求系统响应的方法有两种：第一种是通过解差分方程求得系统输出；第二种是系统的单位脉冲响应，通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出。2、检验系统的稳定性，其方法是在输入端参加单位阶跃序列，观察输出波形，如果波形稳定在一个常数值包括零上，系统稳定，否那么不稳定。实验的心得体会见下：在此次试验
13、中，通过课堂所留的三个例子，温习了关于 MATLAB 软件的操作及应用，根本使用方法和它的运行环境。又进一步地通过实验加深了对MATLAB 软件的了解，体会到了 MATLAB 具有完备的图形处理功能，实现计算结果和编程的可视化等功能。通过做实验的过程以及实验分析的结果，了解并学会了 filter 函数和 conv 函数的根本用法，前者可计算知道输入信号的前提下求解解差分方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系
14、统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果因为零极点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应并画出波形程.word.输出响应的序列，后者那么可以通过输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。当然，在实验过程中，也遇到了一些问题，比方课堂作业的第三个问题里面，由于程序里面缺少“；，导致少了一个结果图，通过检查并修改程序
15、，解决了问题。总得来说，实验还是比拟圆满的。通过这次的实验。极提升了自己对于程序编辑的熟练度，增加了对于书本里面知识点的应用，更深一层的加深了对 MATLAB 软件的使用。这对自己以后的实验积累了丰富的经历。六、附件：MATLAB 原程序清单 1.1调用 conv 函数计算卷积 close all;clear all A=1,-0.9;B=0.05,0.05;%系统差分方程系数向量B和A x1n=1 1 1 1 1 1 1 1 zeros(1,40);%产生信号x1(n)=R8(n)，用zeros用来加点的个数 x2n=ones(1,60);%产生信号x2(n)=u(n)hn=impz(B,A
16、,40);%求系统单位脉冲响应h(n)subplot(3,1,1);stem(hn);%调用函数stem绘图 title(a)系统单位脉冲响应h(n);y1n=filter(B,A,x1n);%求系统对x1(n)的响应y1(n)subplot(3,1,2);stem(y1n);title(b)系统对R8(n)的响应y1(n);y2n=filter(B,A,x2n);%求系统对x2(n)的响应y2(n)subplot(3,1,3);stem(y2n);title(c)系统对u(n)的响应y2(n);1.2利用线性卷积求取二个不同系统输出响应 close all;clear all x1n=one
17、s(1,8);%产生信号x1(n)=R8(n)h1n=ones(1,10)zeros(1,10);h2n=1 2.5 2.5 1 zeros(1,10);y21n=conv(h1n,x1n);y22n=conv(h2n,x1n);subplot(2,2,1);stem(h1n);title(d)系统单位脉冲响应h1(n);解差分方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系统的频域特性是为了知晓系统对不同频率的输入
18、信号所产生的响应结果因为零极点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应并画出波形程.word.subplot(2,2,3);stem(y21n);title(e)h1(n)与R8(n)的卷积y21(n);subplot(2,2,2);stem(h2n);title(f)系统单位脉冲响应h2(n);subplot(2,2,4);stem(y22n);t
19、itle(g)h2(n)与R8(n)的卷积y22(n);1.3 调用 filter 解差分方程求取系统单位脉冲响应和输出响应 close all;clear all un=ones(1,300);%产生信号u(n)n=0:299;xsin=sin(0.014*n)+sin(0.4*n);%产生正弦信号 A=1,-1.8237,0.9801;B=1/100.49,0,-1/100.49;%系统差分方程系数向量B和A y1n=filter(B,A,un);%谐振器对u(n)的响应y31(n)y2n=filter(B,A,xsin);%谐振器对u(n)的响应y31(n)subplot(2,1,1);
20、stem(y1n);title(h)谐振器对u(n)的响应y31(n);subplot(2,1,2);stem(y2n);title(i)谐振器对正弦信号的响应y32(n);课堂作业 1、绘出11)()(NNmmnnx的频谱。解：实验程序：实验结果：1.1 n=200;stept=2*pi/n;w=stept:stept:2*pi;y=sin(2.5*w)./sin(0.5*w);plot(w,y,w,zeros(size(w);axis(stept 2*pi-2 6);ylabel(y=sin(2.5*pi)/sin(0.5*pi);xlabel(w=02*pi);grid on;1.2 n
21、=200;stept=2*pi/n;w=stept:stept:2*pi;y=sin(2.5*w)./sin(0.5*w);plot(w,abs(y),w,zeros(size(w);axis(stept 2*pi 0 6);ylabel(y=sin(2.5*pi)/sin(0.5*pi);xlabel(w=02*pi);grid on;解差分方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系统的频域特性是为了知晓系统
22、对不同频率的输入信号所产生的响应结果因为零极点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应并画出波形程.word.123456-2-10123456y=sin(2.5*pi)/sin(0.5*pi)w=02*pi1234560123456y=sin(2.5*pi)/sin(0.5*pi)w=02*pi 图 1.1 图 1.2 2、输入5,4,3,2,1)
23、(nx，单位脉冲响应2,4,6,3,2,6)(nh，求输出序列，解：实验程序：N=5;M=6;L=N+M-1;x=1,2,3,4,5;h=6,2,3,6,4,2;y=conv(x,h);nx=0:N-1;nh=0:M-1;ny=0:L-1;subplot(2 3 1);stem(nx,x,.k);xlabel(n);ylabel(x(n);grid on;subplot(2 3 2);stem(nh,h,.k);xlabel(n);ylabel(h(n);grid on;subplot(2 3 3);stem(ny,y,.k);xlabel(n);ylabel(y(n);grid on;解差分
24、方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果因为零极点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应并画出波形程
25、.word.实验结果：024012345nx(n)050246nh(n)0510020406080ny(n)3、分析频谱a，b，c 保证幅频特性的最大值为 1。11111)(pzzazH，11211)(pzzbzH，)1)(1()1)(1()(11113zrezrezzczHjj 4,5.0,8.0rp 解：实验程序：图 2.1 图 2.2 图 2.3 p=0.8;r=0.85;alpha=pi/4;N=25;b1=1,1;a1=1-p;a=(1-p)/2;b1=b1*a;b2=1,-1;a2=1-p;b=(1+p)/2;b2=b2*b;b3=1 0-1;a3=1-2*r*cos(alpha)
26、r*r;c1=exp(j*2*alpha);c=abs(1-c1)/(1-r)*abs(1-r*c1);b3=b3/c;h1=impz(b1,a1,N);subplot(331);stem(h1,.k);subplot(332);zplane(b1,a1)H1,P=freqz(b1,a1,256,whole,1);subplot(333);plot(P,abs(H1);grid on;h2=impz(b2,a2,N);subplot(334);stem(h2,.k)hold on;plot(zeros(size(h2);subplot(335);zplane(b2,a2)H2,P=freqz(
27、b2,a2,256,whole,1);subplot(336);plot(P,abs(H2);grid on;h3=impz(b3,a3,N);subplot(337);stem(h3,.k)hold on;plot(zeros(size(h3);subplot(338);zplane(b3,a3)H3,P=freqz(b3,a3,256,whole,1);subplot(339);plot(P,abs(H3);grid on;解差分方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波
28、形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果因为零极点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应并画出波形程.word.实验结果：051015202500.050.10.150.2-1-0.5 00.5 1-1-0.500.51Real PartImaginary Pa
29、rt00.20.40.60.8100.20.40.60.810510152025-0.500.51-1-0.5 00.5 1-1-0.500.51Real PartImaginary Part00.20.40.60.8100.20.40.60.810510152025-0.2-0.100.10.2-1-0.5 00.5 1-1-0.500.51Real PartImaginary Part00.20.40.60.8100.511.5 图 3.1 图 3.2 图 3.3 图 3.7 图 3.8 图 3.9 图 3.4 图 3.5 图 3.6 解差分方程求得系统输出第二种是系统的单位脉冲响应通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出检验系统的稳定性其方法是在输入端参加单位阶跃序列观察输出波形如果波形稳定在一个常数值包括零上系统稳定脉冲响应的变换零极点分布等分析系统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果因为零极点分布对系统的频域特性有影响通过控制系统函数的零极点分布就可以设计出不同特性需求的系统二实验容编制程的功能给定一个低通滤波器的差分方程为输入信号分别求出和的系统响应并画出其波形求出系统的单位脉冲响应画出其波形程序见附录实验结果见图给定系统的单位脉冲响应为用线性卷积法求分别对系统和的输出响应并画出波形程