- 1、本文档共16页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
演示文稿演讲PPT学习教学课件医学文件教学培训课件
实验19 离散系统状态方程的求解 一、实验目的 二、实验原理 三、涉及的MATLAB函数 四、实验内容与方法 五、实验要求 六、思考题 一、实验目的 (1) 了解离散系统状态方程求解方法; (2) 了解离散系统信号流图化简的方法; (3) 了解函数ode45的调用方法。二、实验原理 离散系统状态方程的一般形式为 x(k+1)=Ax(k)+Bf(k) 在此只对单输入的n阶离散系统的状态方程求解。一般采用递推迭代的方式求解,由初始条件x(0)和激励f(0)求出k=1时的x(1),然后依次迭代求得所要求的x(0),…,x(n)的编程时应注意,MATLAB中变量下标不允许为零,则初始点的下标只能取1,第n步的x的下标为n+1。 三、涉及的MATLAB函数 采用函数ode45可以求解微分方程。其调用格式如下: [t,y]=ode45(odefun,tspan,y0) 其中,odefun指状态方程的表达式,tspan指状态方程对应的起止时间[t0,tf],y0指状态变量的初始状态。 MATLAB程序:%离散系统状态求解%A=input(′系数矩阵A=′)%B=input(′系数矩阵B=′)%x0=input(′初始状态矩阵x0=′)%n=input(′要求计算的步长n=′) %f=input(′输入信号f=′) %要求长度为n的数组clear allA=[0.5 0; 0.25 0.25]; B=[1; 0]; x0=[-1; 0.5]; n=10; f=[0 0.5*ones(1,n-1)]; x(:,1)=x0; for i=1:n x(:,i+1)=A*x(:,i)+B*f(i); endsubplot(2,1,1); stem([0:n],x(1,:)); subplot(2,1,2); stem([0:n],x(2,:));离散系统状态方程的求解结果如图19.1所示。 图 19.1 离散系统状态方程的求解 (2) 离散系统状态求解。MATLAB程序:A=[0 1; -2 3]; B=[0; 1]; %方程输入C=[1 1; 2 -1]; D=zeros(2,1); x0=[1; -1]; %初始条件N=10; f=ones(1,N); sys=ss(A,B,C,D,[]); y=lsim(sys,f,[],x0); k=0:N-1; subplot(2,1,1); stem(k,y(:,1),′b′); subplot(2,1,2); stem(k,y(:,2),′b′);离散系统状态方程的求解结果如图19.2所示。 图 19.2 离散系统状态方程的求解 实验19 离散系统状态方程的求解
文档评论(0)