资源描述:
《matlab实验文档》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、.实验一1.用Matlab语言求下列系统的状态方程、传递函数、零极点增益、和部分分式形式的模型参数,并分别写出其相应的数学模型表达式:(1)G(s)=(2)=y=[0202]X2.已知元件的实验数据如下,拟合这一数据,并尝试给出其特性方程。X0.01001.01002.01003.01004.0100Y2.54377.88849.624211.607111.9727X5.01006.01007.01008.01009.0100y13.218914.267914.613415.404515.0805num=[1144848]den=[1207010048][ABCD]=tf2
2、ss(num,den)[ZPK]=tf2zp(num,den)[RPH]=residue(num,den)A=[2.25-5-1.25-0.52.25-4.25-1.25-0.250.25-0.5-1.25-11.25-1.75-0.25-0.75]B=[4220]'C=[0202]D=[0];[numden]=ss2tf(A,B,C,D)[ZPK]=ss2zp(A,B,C,D)[RPH]=residue(num,den)X=0.01:1:9.01;Y=[2.54377.88849.624211.607111.972713.218914.267914.613415.4045
3、15.0805]p=polyfit(X,Y,5)Xi=0.01:1:9.01Yi=polyval(p,Xi)...plot(X,Y,Xi,Yi)gridon实验二1.分别用欧拉法、二阶龙格-库塔法、四阶龙格-库塔法求下面系统的输出响应y(t)在0≤t≤1上,h=0.1时的数值解。要求保留4位小数,并将三种方法的结果与真解进行比较。2.若为如何编程计算?h=0.1;disp('º¯ÊýµÄÊýÖµ½âΪ');disp('y=');y=1;fort=0:h:1m=y;disp(y);y=m-2*m*h;endh=0.1;disp('º¯ÊýµÄÊýÖµ½âΪ');disp(
4、'y=');y=1;fort=0:h:1disp(y);k1=-2*y;k2=-2*(y+k1*h);y=y+(k1+k2)*h/2;endh=0.1;disp('ËĽ×Áú¸ñ-¿âËþ·½·¨Çó½âº¯ÊýÊýÖµ½âΪ');disp('y=');y=1;fort=0:h:1disp(y);k1=-2*y;k2=-2*(y+k1*h/2);k3=-2*(y+k2*h/2);...k4=-2*(y+k3*h);y=y+(k1+2*k2+2*k3+k4)*h/6;endh=0.1;disp('º¯ÊýµÄÊýÖµ½âΪ');disp('y=');y=1;fort=0
5、:h:1y=exp(-2*t)endh=0.1;disp('º¯ÊýµÄÊýÖµ½âΪ');disp('y=');y=1;fort=0:h:1m=y;disp(y);y=m+m*m*h;endh=0.1;disp('º¯ÊýµÄÊýÖµ½âΪ');disp('y=');y=1;fort=0:h:1disp(y);k1=y*y;k2=(y+k1*h)^2;y=y+(k1+k2)*h/2;endh=0.1;disp('ËĽ×Áú¸ñ-¿âËþ·½·¨Çó½âº¯ÊýÊýÖµ½âΪ');disp('y=');y=1;fort=0:h:1disp(y);k1=y*y;...
6、k2=(y+k1*h/2)^2;k3=(y+k2*h/2)^2;k4=(y+k3*h)^2;y=y+(k1+2*k2+2*k3+k4)*h/6;end实验三1.分别使用解微分方程方法、控制工具箱、Simulink求解具有如下闭环传递函数的系统的阶跃响应。num=[10];>>den=[18364010];>>[ABCD]=tf2ss(num,den)A=-8-36-40-10100001000010B=1000C=00010D=0得到状态方程编写m文件求解微分方程组...functiondx=wffc(t,x)u=1;%阶跃响应,输入为1%dx=[-8*x(1)-36*x(
7、2)-40*x(3)-10*x(4)+u;x(1);x(2);x(3)];保存文件wffc.m%注意:保存文件的名字与函数名一致!%在命令行键入>>[t,x]=ode45('wffc',[0,8],[0;0;0;0]);>>y=10*x(:,4);>>plot(t,y);>>grid(2)控制工具箱:在matlab命令行中键入>>num=[10];>>den=[18364010];>>sys=tf(num,den);>>step(sys);>>grid(3)simulink求解:在simulink模型窗口中建立如下