资源描述：
《系统建模与仿真作业》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、四川理工学院系统建模与仿真课程报告姓名：孙昌刚学号：11021010416专业：自动化班级：自动化112任课教师：李兆飞四川理工学院自动化与电子信息学院二O一三年十一月20目录第1题31.1题目31.2解题31.3本题总结5第2题62.1题目62.2解题62.3本题总结12第3题123.1题目123.2解题123.3本题总结18第4题184.1题目184.2解题18第5题195.1题目195.2解题19参考文献2020第1题1.1题目＝1+/6+/7+/10，，采用面积法和最小二乘法把系统G等效成特性，求系统的kT等参数。写出等效过程及相
2、关程序。最后把等效后的两个模型及原系统对阶跃输入的响应曲线绘制在一个图上进行比较，并分析优劣问题。1.2程序(1)(2)20clcruxuenian=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=
3、[1:0.1:10];gs=2.*x./((x.*s+1).^4);figure(1)plot(s,gs);a0=[000];a=lsqcurvefit(@(a,s)a(1)./(a(2).*s+1).*exp(-a(3).*s),a0,s,gs)holdonplot(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);sys
4、t=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请输入学号：11021010416a=00020Zero/pole/gain:673.819--------------(s+0.002968)^41.3本题总结20采用系统模型的方法比较简单可靠，并且使用起来非常的方便。而等效的程序写起来比较麻烦，而且还需要记忆很多的指令，用起来不如系统建模的
5、方法简便、快捷。因此，采用系统建模的方法更好、更实用。第二题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.mfunctionfv=Funval(f,
6、varvec,varval)var=findsym(f);varc=findsym(varvec);s1=length(var);s2=length(varc);m=floor((s1-1)/3+1);varv=zeros(1,m);ifs1~=s2fori=0:((s1-1)/3)k=findstr(varc,var(3*i+1));index=(k-1)/3;varv(i+1)=varval(index+1);endfv=subs(f,var,varv);elsefv=subs(f,varvec,varval);endLGKT4.m2
7、0functiony=LGKT4(f,h,a,b,y0,varvec)formatlong;N=(b-a)/h;y=zeros(N+1,1);y(1)=y0;x=a:h:b;var=findsym(f);fori=2:N+1K1=Funval(f,varvec,[x(i-1)y(i-1)]);K2=Funval(f,varvec,[x(i-1)+h/2y(i-1)+K1*h/2]);K3=Funval(f,varvec,[x(i-1)+h/2y(i-1)+K2*h/2]);K4=Funval(f,varvec,[x(i-1)+hy(i-1
8、)+h*K3]);y(i)=y(i-1)+h*(K1+2*K2+2*K3+K4)/6;endformatshort;plot(x,y,'r:')legend('h=0.01','h=0.1',
