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

四川理工学院 系 统 建 模 与 仿 真 课程报告 学 生： 学 号： 专 业：自 动 化 班 级：自动化112班 任课教师： 老师 四川 二O一三年十一月 目录 第一章………………………………………………………3 1.1 题目……………………………………………………3 1.2 解答……………………………………………………3 第二章………………………………………………………3 2.1 题目………………………………………………………3 2.2 程序………………………………………………………3 2.3 小结……………………………………………………8 第三章………………………………………………………11 3.1 题目………………………………………………………8 3.2 解答………………………………………………………8 第四章………………………………………………………16 4.1 题目……………………………………………………16 4.2 解答……………………………………………………16 第五章………………………………………………………16 5.1 题目……………………………………………………17 5.2 解答……………………………………………………17 第一章 1.1题目 、＝1+/6+/7+/10， ，采用面积法和最小二乘法把系统G等效成特性，求系统的k T 等参数。写出等效过程及相关程序。 最后把等效后的两个模型及原系统对阶跃输入的响应曲线绘制在一个图上进行比较，并分析优劣问题。（x=3.27） 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 请输入学号：11021010212 a = 0 0 0 Zero/pole/gain: 673.8857 -------------- (s+0.002968)^4 第二章 2.1 题目 编写微分方程dy/dx＝xy, 当x＝0时y＝1+8/10+10/100=1.9, 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 LG