系统建模及仿真课程报告.docVIP

系 统 建 模 与 仿 真 课程报告 姓 名： 学 号： 专 业：自动化 班 级： 任课教师： 二O一三年十一月 1.1题目 、＝1+/6+/7+/10， ，采用面积法和最小二乘法把系统G等效成特性，求系统的k T 等参数。写出等效过程及相关程序。 最后把等效后的两个模型及原系统对阶跃输入的响应曲线绘制在一个图上进行比较，并分析优劣问题。 1.2程序 clc ruxuenian=input(请输入入学年：); banji=input(请输入班级：); xuehao=input(请输入学号：); rxn=num2str(ruxuenian); xh=num2str(xuehao); a=ruxuenian; b=banji; c=str2num(xh(length(xh))); d=str2num(xh(length(xh)-1)); x=1+a/6+d/7+c/10; s=[1:0.1:10]; gs=2.*x./((x.*s+1).^4); figure(1) plot(s,gs); a0=[0 0 0]; a=lsqcurvefit(@(a,s)a(1)./(a(2).*s+1).*exp(-a(3).*s),a0,s,gs) hold on plot(s,a(1)./(a(2).*s+1).*exp(-a(3).*s),m); sys1=zpk([],[-1/x -1/x -1/x -1/x],2*x) figure(2) subplot(1,2,1) step(sys1); [numt,dent]=pade(a(3),5); syst=tf(numt,dent); num1=a(1); den1=[a(2) 1]; sysT=tf(num1,den1); sys2=series(sysT,syst); subplot(1,2,2) step(sys2); 请输入入学年：2011 请输入班级：5 请输入学号：11021010216 Optimization terminated: first-order optimality less than OPTIONS.TolFun, and no negative/zero curvature detected in trust region model. a = 0 0 0 Zero/pole/gain: 673.819 -------------- (s+0.002968)^4 2.1 题目 编写微分方程dy/dx＝xy, 当x＝0时y＝1+8/10+11/100=1.9, x属于0～3之间，编写积分程序，包括欧拉数值积分程序，预报校正数字积分程序、4阶龙格库塔积分程序，它们的积分步长分别取0.01，0.1, 0.5，绘制积分结果曲线,比较在同一步长下不同算法的误差和同一算法在不同步长下的误差，得出结论 说明（绿色线为欧拉法曲线，红色为预报校正法曲线，蓝色为4阶龙格库塔法曲线） 2.2 程序 欧拉数字积分程序 function [x,y]=naler(dyfun,xspan,y0,h) x=xspan(1):h:xspan(2); y(1)=y0; for n=1:length(x)-1; y(n+1)=y(n)+h*feval(dyfun,x(n),y(n)); end x=x;y=y; 预报校正数字积分程序 function[x,y]=zseuler(dyfun,xspan,y0,h) x=xspan(1):h:xspan(2); y(1)=y0; for n=1:length(x)-1 k1=feval(dyfun,x(n),y(n)); y(n+1)=y(n)+h*k1; k2=feval(dyfun,x(n+1),y(n+1)); y(n+1)=y(n)+h*(k1+k2)/2; end x=x;y=y; 4阶龙格库塔积分程序 function[x,y]=do(dyfun,xspan,y0,h) x=xspan(1):h:xspan(2); y(1)=y0; for n=1:length(x)-1 k1=feval(dyfun,x(n),y(n)); k2=feval(dyfun,x(n)+h/2,y(n)+h/2*k1); k3=feval(dyfun,x(n)+h/2,y(n)+h/2*k2); k4=feval(dyfun,x(n)+h,y(n)+h*k3); y(n+1)=y(n)+h*(k1+2*k2+2*k3+k4)/6; end x=x;y=y 步长为0.01, 不