- 积分
- 0
- 注册时间
- 2009-3-2
- 仿真币
-
- 最后登录
- 1970-1-1
|
求助:
想做X(1)与dx(1)/dt的平面图。
程序中,X(6)=dx(1)/dt
,采用了微分代数方程求解的方法。
但做不出来,从输出的图形来看X(1)和X(6)基本上是同相位,没有看到微分的效果。
不知道程序有啥错误。
请指教!
谢谢!
定义函数程序
function xdot = fun(t,x)r1=0.0003;r2=0.0003;r3=0.0003;g=1e-9;a1=41.3*1e-6;a2=41.3*1e-6;a3=41.3*1e-6;b1=0.246*1e-6;b2=0.246*1e-6;b3=0.246*1e-6;c1=0;c2=0;c3=0;C=0.3334;Em=8100;xdot=[-r1*(a1*x(1)+b1*x(1)^3+c1*x(1)^5)+x(4)+Em*sin(314*x(5));-r2*(a2*x(2)+b2*x(2)^3+c2*x(2)^5)+x(4)+Em*sin(314*x(5)-2.1);-r3*(a3*x(3)+b3*x(3)^3+c3*x(3)^5)+x(4)+Em*sin(314*x(5)+2.1);-(a1*x(1)+b1*x(1)^3+c1*x(1)^5+a2*x(2)+b2*x(2)^3+c2*x(2)^5+a3*x(3)+b3*x(3)^3+c3*x(3)^5)/(3*C)-g*x(4)/C;1;-r1*(a1*x(1)+b1*x(1)^3+c1*x(1)^5)+x(4)+Em*sin(314*x(5))-x(6)];
运行程序程序名K1
clear,clc,clf;M=[1,0,0,0,0,0;0,1,0,0,0,0;0,0,1,0,0,0;0,0,0,1,0,0;0,0,0,0,1,0;0,0,0,0,0,0];options=odeset;option.Mass=M;tfinal=0.3;x0=[2,-0.2,-1,0,0,0];[t,x]=ode45('fun',[0,tfinal],x0);figure(1)hold onplot(t,x(:,1),'-r')xlabel('t')ylabel('va')figure(2)plot(t,x(:,2),'-.b')xlabel('t')ylabel('vb')figure(3)plot(t,x(:,3),'k')xlabel('t')ylabel('vc')figure(4)%plot3(x(:,1),x(:,2),x(:,3))%plot(x(:,5),x(:,1))plot(x(:,1),x(:,6)) |
|