系统建模与仿真作业.docVIP

四川理工学院 系 统 建 模 与 仿 真 课程报告 姓 名：师兄 学 号：11041050118 专 业：自动化 班 级：自动化113 任课教师：李兆飞 四川理工学院自动化与电子信息学院 二O一三年十一月 目录 1.1 题目 1 1.2 解题 1 1.3 本题总结 4 2.1 题目 4 2.2 解题 4 2.3 本题总结 9 3.1 题目 9 3.2 解题 10 3.3本题小结 14 4.1 题目 14 4.2 解答 14 5.1 题目 14 5.2 解答 15 参考文献 15 1.1 题目 1、＝1+/6+/7+/10， ，采用面积法和最小二乘法把系统G等效成特性，求系统的k T 等参数。写出等效过程及相关程序。 最后把等效后的两个模型及原系统对阶跃输入的响应曲线绘制在一个图上进行比较，并分析优劣问题。 1.2 解题 ＝11 ＝3 ＝8 ＝1 ＝1 ＝1+/6+/7+/10＝3.78 clc? ruxuenian=input(请输入入学年：);? banji=input(请输入班级：);? xuehao=input(请输入学号：);? rxn=num2str(ruxuenian);? xh=num2str(xuehao);? a=ruxuenian;? b=;? c=str2num(ah(length(ah)));? d=str2num(ah(length(ah)-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??];?a=(@(a,s)a(1)./(a(2).*s+1).*exp(-a(3).*s),a0,s,gs) hold爋n?plot(s,a(1)./(a(2).*s+1).*exp(-a(3).*s),m);?sys1=([],[-1/x?1/x?1/x?1/x],2*x)牋? figure(2)?subplot(1,2,1)?step(sys1);?[,dent]=(a(3),5);?syst=(,dent);?num1=a(1);?den1=[a(2)?];?=(num1,den1);?sys2=series(,syst);?subplot(1,2,2)爏tep(sys2);?请输入入学年：2011 请输入班级：3 请输入学号：11041050118 ??a?=?? ?????0?????0?????0???? Zero/pole/gain:??? ?673.8857? -------------- ?(s+0.002968)^4 1.3 本题总结 最小二乘法（又称最小平方法）是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达。简单地说，最小二乘的思想就是要使得观测点和估计点的距离的平方和达到最小.这里的“二乘”指的是用平方来度量观测点与估计点的远近（在古汉语中“平方”称为“二乘”），“最小”指的是参数的估计值要保证各个观测点与估计点的距离的平方和达到最小. 2.1 题目 编写微分方程dy/dx＝xy, 当x＝0时y＝1+8/10+11/100＝1.91, x属于0～3之间，编写积分程序，包括欧拉数值积分程序，预报校正数字积分程序、4阶龙格库塔积分程序，它们的积分步长分别取0.01，0.1, 0.5，绘制积分结果曲线,比较在同一步长下不同算法的误差和同一算法在不同步长下的误差，得出结论。 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=