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...

在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...

dsolve()求解微分方程后,使用ezplot()作图,参考代码: >> y = dsolve('Dy=y*cos(x)/(1+2*y^2)','y(0)=1','x');>> ezplot(y,[-10,10]) % [-10,10]为x区间

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

好像没什么特别好的办法: >> 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代入,再求解。

>> syms x(t) y(t) r d a b r='1';d='0.5';a='0.1';b='0.02'; [x,y]=dsolve(diff(x)==(r-a*y)*x,diff(y)==-(d-b*x)*y,x(0)==25,y(0)==2) Warning: Explicit solution could not be found. > In dsolve at 194 x = [ empty sym ] y = [] 说明无解。

在Matlab下输入:edit zhidao_feiying.m,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_feiying(t,x) a=1; b=2; c=3; d=4; %%比方说 %f=a*x+y-z; %g=b*sin(x*y)-c*cos...

这说明此微分方程没有显示解,应用数值方法求解,例如ode45等函数 就一个例子(选自matlab帮助文档) 解入下微分方程组 1,建立方程组函数 function dy = rigid(t,y)dy = zeros(3,1); % a column vectordy(1) = y(2) * y(3);dy(2) = -y(1) * y(3);d...

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