实验五连续系统分析

《实验五连续系统分析》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验五连续系统分析一、实验目的深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性屮的重要作用及意义，掌握根据系统阑数的零极点设计简单的滤波器的方法。掌握利用MATLAB分析连续系统的时域响应、频响特性和零极点的基本方法。二、实验原理MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。1.连续系统的时域响应连续时间LTI系统可用如下的线性常系数微分方程来描述：+U(”"(，)+…+(V，⑴+aoy(t)=bmx(,nt)+u(z

2、)+…+办，(z)+bQx(t)o已知输入信号40以及系统初始状态y((T),j/((r)，…,^"4(0"),就可以求出系统的响应。MATLAB提供了微分方程的数值il•算的函数，可以il•算上述n阶微分方程描述的连续系统的响应，包括系统的单位冲激响应、单位阶跃响应、零输入响应、零状态响应和完全响应。在调用MATLAB函数时，需要利用连续系统对应的系数函数。对微分方程进行Laplace变换即可得系统函数：X（5）67^5+Cln_ySn+•••+在MATLAB中可使用向量和叫量分别保存分母多项式和分子多项式的系

3、数：“丨，••Wo]6哦,，bm-，…，々I，办o]这些系数均按s的降幂直至沁排列。(1)连续系统的单位冲激响应A(/)的计算impulse(sys)计算并画出系统的冲激响应。参数：sys可由函数tf(b,a)获得。K屮：“=b=[bm，bm_，…,b'，b0]h=impulse(sys,t):计算井画出系统在向量t定义的区间上的冲激响应，向量h保存对应区间的系统冲激响应的输出值。(2)连续系统的单位阶跃响应以0的计算step(sys):计算并画出系统的阶跃响应。参数：sys可由函数tf(b，a)获得。其中:

4、a=[an,an.{,•••,«,,«01/?=H，…，AA1g=step(sys,t):计算并画出系统在向量t定义的区间上的阶跃响应，向量g保存对应区间的系统阶跃响应的输出值。(1)连续系统的零状态响应的计算lsim(sys,x,t)计算并M山系统的零状态响应。参数:sys可由函数tf(b,a)获得，X为输入信号，t为定义的时间向量。1.连续系统的系统函数零极点分析连续LTI系统的系统阑数//⑻可以表示为部分分式形式："⑴==k(5-4)(卜z2)•••&-、)D(S)(5-p,)(s-p2)".O-)设m

5、,KHCs)的极点pf全部为单极点，则:k.nH⑴=Z——Kt)=XWOi=s一Pii=]系统函数//⑴的极点p/决定丫冲激响应/价)的基本形式，而零点和极点共同确定丫冲激响应/2⑺的幅值欠。MATLAB屮提供了roots函数计算系统的零极点，提供了pzmap函数绘制连续系统的零极点分布图。2.连续系统的频率响应若连续因果LTI连续系统的系统函数//⑴的极点企部位于S左半平面，则系统的频率响应可由H⑴求出，即MATLAB中freqs函数可以分析连续系统的频响，格式如下：H=freqs(b，a，w):计算系统在指

6、定频率点肉量w上的频响H;w为频率点向量。[H,w]=freqs(b,a):自动选取200个频率点计算频率响应。三.实验内容1.己知描述连续系统的微分方程为#^+10>，(Z)=2x(Z),输入=w⑺，初始状态dty(0_)=l,计算该系统的响应，并与理论结果比较，列出系统响应分析的步骤。实验代码：a=fl101;b=[2];[ABCD]=tf2ss(b,a);sys=ss(A，B，C，D);t=0:0.001:5;xt=t>0;sta=[l];y=lsim(sys，xt，t，sta);subplot(3，l，l

7、);plot(t，y);xlabel('t’)；title(’系统完全响应y(t)’);subplot(3，l，2);plot(t,y；-b');holdonyt=4/5*exp(-10*t)+l/5;plot(t，yt，’:r’)；legendf数值计算7理论计算’);holdoffxlabel(’t’)；subplot(3，l，3);k=y’-yt;plot(t，k);k(l)titlef误差’)；实验结果:系统完全响应y(t)00.511.522.533.544.550结果分析：求输入x(r)分別为w(r)

8、，理论值y(t)=0.8*exp(-10t)*U(t)+0.2。程序运行出的结果与理论预期结果相差较大误差随吋间增大而变小，初始值相差最大，终值基本相同。1.已知连续时间系统的系统函数为=―———，s"I-2»ysintu(t),时，系统地输出j，0),并与理论结果比较。实验代码：a=[l,3,2,0];b=「4，ll;sys=tf(b,a);t=0:0.001:5;xl