lzth.net
当前位置:首页 >> mAtlAB中求微分方程组通解 >>

mAtlAB中求微分方程组通解

此题有误,两个方程的左侧应是dx/dt,dy/dt [x y]=dsolve('Dx=4*x-2*y', 'Dy=2*x-y','t') x = 2*C2*exp(3*t)+1/2*C1 y = C1+C2*exp(3*t) [x y]=dsolve('Dy=4*x-2*y', 'Dx=2*x-y','t') x = 1/4*C1+1/2*C1*t+1/2*C2 y = C1*t+C2

syms x(t) y(t) z(t) [x,y,z]=dsolve(diff(x)==4*x+3*y+2*z,diff(y)==6*x-3*y+2*z,diff(z)==6*x+5*y+4*z) 按上面试试

可这样: >> syms Vx Vy Vz b m y1(t) y2(t) y3(t) y4(t) y5(t) y6(t) y1=dsolve(diff(y1)==Vx) y2=dsolve(diff(y2)==-b/m*(Vx^2+Vy^2+Vz^2)*Vx) y3=dsolve(diff(y3)==Vy) y4=dsolve(diff(y4)==-b/m*(Vx^2+Vy^2+Vz^2)*Vy) y5=dsolve(diff(y5)==V...

数值求解法:先建立个函数文件里面写上方程组,如: function f=weifen(t,x) a=3.5;d=3;b=0;f=0; c=-0.7;e=-0.5; f=[(a+b*x(1)+c*x(2)),(d+e*x(1)+f*x(2))]'; 这里的f表示方程组dx/dt=a+bx+cy;dy/dt=d+ex+fy 然后用ode45函数求数值解即可代入初...

在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dxdt=zhidao_rk4_5(t,x) global a b c d r x1=x(1); x2=x(2); x3=x(3); x4=x(4); dxdt=[ a*(x2-x1)+x4; d*x1-x...

好像没什么特别好的办法: >> syms t q y; q=4*t^3;>> syms Dy>> eq=Dy-(q-2*t);>> y=dsolve(char(subs(eq)),'y(0)=5','t')y = t^4-t^2+5 也就是要先移项把微分方程写成等号右端为零的形式(需额外定义Dy),然后subs代入,再求解。

1、对于比较简单的二阶微分方程组,可以用dsolve()函数求得其特解,例如: syms y(t) z(t)%定义变量 Dy=diff(y);Dz=diff(z);%对y、z求一阶导数 s=dsolve(Dy==3*y+2*z-(2*t^2+1)*exp(2*t),Dz==4*y+z+(t^2+2*t-4)*exp(2*t),y(0)==1,z(0)==1) %求...

建立myfun.m文件,把以下代码复制进去 function dx=myfun(t,x)dx=zeros(2,1);dx(1)=4*x(1)+x(2).^2;dx(2)=8*x(1)+2*x(2).^2+5*x(2);运行以下代码 t=[0 0.1]; x0=[1 1]; [tt,xx]=ode45('myfun',t,x0); plot(tt,xx(:,1),'-',tt,xx(:,2),'-.')%画出...

说明 求解常微分方程一般使用ode系列函数(其中最常用的是ode45),基本做法就是把微分方程用函数表示出来,然后调用ode函数求解即可。对于本题而言,使用匿名函数描述微分方程更为方便一些。 参考代码 k0=8; k1=1; K=1; k2=5; n=2; % 常数定义d...

分段函数的一个简单方法就是分了几段就写几个函数: 定义一个主调函数、两个子函数 最后一个主调程序调用 function dx=fun(t,x) dx=zeros(2,1); f=; dc(t)=diff('f','t');%f就是你的c(t) dx(1)=x(2)*dc(t); dx(2)=x(1)+6*t^2;%0

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