bycj.net
当前位置:首页 >> 大神求解动力学微分方程,用mAtlAB中oDE45编程实现... >>

大神求解动力学微分方程,用mAtlAB中oDE45编程实现...

参考代码: % 参数定义m1 = 1; m2 = 2; m3 = 3;k1 = 4; k2 = 5; k3 = 6; k4 = 7;P0 = 8; w = 9;% 微分方程A = [k1+k2 -k1 0; -k2 k2+k3 -k3; 0 -k3 k3+k4];dx = @(t,x)[x(4:6); (-A*x(1:3)+[P0*sin(w*t);0;0])./[m1;m2;m3]];% 初始条件(x4~x6 为 x1'~x3')x0 =

这个方程要解出来还要初始条件吧 不妨设(x(0),y(0))=(a,b); 还有用ode45只能得到数值解,不能给出解析式表达 那么先建立一个函数文件 test_fun.m function dx=tsst_fun(x,t) dx(1)= c-b*x(1)-w*x(1)*x(2); dx(2)=w*x(1)*x(2)-(b+m)*x(2); 那么调用ode45 x0=(0,0); t0=0:0.01:10;;%根据你所想要的求得值设定t0,间隔是任意的,与求解所用的步长无关,[x,t]=ode45(@test_fun,t0,x0); 得到了 这样可以么?

<p>你把lgkt中的x改成x(1)试试看,我试过应该没问题的</p> <p>function xx=lgkt(t,x)</p> <p>xx=3*x(1)+cos(t)+sin(2*t);</p> <p>[t,x]=ode45(@lgkt,[0 10],[2]);</p> <p>plot(t,x,'-o');</p> <p></p>

使用ode45解微分方程组,可以用下列格式来求解.[t,x]=ode45(odefun,tspan,x0) 式中:odefun给复定微分方程组的自定义函数;tspan自变量t的范围,制即【t0,tf】;x0x的初值向量,即x1(0),x1'(0),x2(0),x2'(0),x3(0),x3'(0) 由于没有给出自zd变量t的范围和x的初值,所以不能得到其x1,x2,x3的值.

大神求解动力学微分方程? 我要提问 大神求解动力学微分方程,用matlab中ode45编程实现 匿名 分享到微博 提交回答 1 问: VB编程如何实现文字的检索 答: 建一个TextBox=Private Sub Form_Load()

大神求解动力学微分方程,用matlab中ode45编程实现. 回答 2 5 matlab中的ode45怎样设置定步长解微分方程? 回答 2 数学 相关知识 教育科学 答: 数学:甲数、乙数与丙数的和是1400,甲数是乙数的

ode45是最常用的求解非刚性微分方程(组)的指令,它采用变步长四、五阶Runge-Kutta法,求解精度比ode23、ode15要高.但ode45的计算量比较大,而ode23计算量小,且误差大 调用格式:[t,y]=ode45(odefun,tspan,y0) odefun 用以表示f(t,y

要画成二元一价微分方程组

function dy=myfunc(x,y)a=b=c=dy=zeros(2,1)dy(1)=y(2)dy(2)=(b*dy(2)^2-c*sin(y(1)))/a这里的dy(1)是一阶导数,所以你的初值需要给出一阶导数值

主函数和子函数均没有错误将主函数中第一行[T,Y]=ode45('Mx1',20,[ 1; 1; 1; 2;3]);中的积分范围改为500,即可!

网站首页 | 网站地图
All rights reserved Powered by www.bycj.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com