系统建模与仿真作业课程报告.docVIP

四川理工学院 系 统 建 模 与 仿 真 课程报告 姓 名： 学 号：11021010416 专 业：自动化 班 级：自动化112 任课教师： 四川理工学院自动化与电子信息学院 二O一三年十一月 目录 第1题 3 1.1 题目 3 1.2 解题 3 1.3 本题总结 5 第2题 5 2.1 题目 5 2.2 解题 5 2.3 本题总结 5 第3题 5 3.1 题目 5 3.2 解题 5 3.3 本题总结 5 第4题 5 4.1 题目 5 4.2 解题 5 第5题 5 5.1 题目 5 5.2 解题 5 参考文献 5 第1题 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 请输入班级：2 请输入学号：11021010416 a = 0 0 0 Zero/pole/gain: 673.819 -------------- (s+0.002968)^4 1.3 本题总结 第二题 2.1 题目 编写微分方程dy/dx＝xy, 当x＝0时y＝1+6/10+11/100=1.71, x属于0～3之间，编写积分程序，包括欧拉数值积分程序，预报校正数字积分程序、4阶龙格库塔积分程序，它们的积分步长分别取0.01，0.1, 0.5，绘制积分结果曲线,比较在同一步长下不同算法的误差和同一算法在不同步长下的误差，得出结论 说明（绿色线为欧拉法曲线，红色为预报校正法曲线，蓝色为4阶龙格库塔法曲线） 2.2 程序 解：(1)4阶龙格库塔 Funval.m function fv = Funval(f,varvec,varval) var = findsym(f); varc = findsym(varvec); s1 = length(var); s2 = length(varc); m =floor((s1-1)/3+1); varv = zeros(1,m); if s1 ~= s2 for i=0: ((s1-1)/3) k = findstr(varc,var(3*i+1)); index = (k-1)/3; varv(i+1) = varval(index+1); end fv = subs(f,var,varv); else fv = subs(f,varvec,varval); end LGKT4.m function y=LGKT4(f,h,a,b,y0,varvec) format long; N=(b-a)/h; y=zeros(N+1,1); y(1)=y0; x=a:h:b; var=findsym(f); for i=2:N+1 K1=Funval(f,varvec,[x(i-1) y(i-1)]); K2=Funval(f,varvec,[x(i-1)+h/2 y(i-1)+K1*h/2]); K3=Funval(f,varvec,[x(i-1)+h