clzm.net
当前位置:首页 >> 用mAtlAB解微分方程 >>

用mAtlAB解微分方程

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

微分方程的话可以用dsolve函数,格式可以总结为dsolve('equation','condition',自变量)的形式。自变量默认为t。 这种方法解微分方程还好。如果是微分方程组的话建议用ode45函数(该方法也用于微分方程) 首先建立M文件写方程组的表达形式 functi...

楼主是要解析解吧 dsolve('Dy=A/V-B/V*sqrt(2*g*y)','y(0)=h0','t') ans = (2^(1/2)*A + 2^(1/2)*A*lambertw(0, (2^(1/2)*B*g^(1/2)*exp(-(A + B^2*g*(t/V - (A*(log(-(exp(-(A + 2^(1/2)*B*(g*h0)^(1/2))/A)*(2*B*(g*h0)^(1/2) + 2^(1/2)*A))/(...

1、对于求数值解的微分方程,你可以用ode45()函数求解。如求下列微分方程 func。m %自定义微分方程的函数 function z = func(t,y) z =[y(2);(1-y(1)^2)*y(2)-y(1)]; main。m %主程序 clear all;close all;clc y0 = [0.25;0]; h = 0.1; a = 0; b ...

matlab对微分方程的解析解并不特别人性化。可以的话列出你的微分方程。让我试试。

在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()函数,求得M(x)d的表达式 >> syms M(x) >> syms a q B11 c1 c2 >> D2M=diff(M,2); >> M=dsolve(a^2+D2M==M+q*x^2/(2*B11)-c1*x-c2) M =c1*x + C3*exp(-x/a) + C4*exp(x/a) + (- q*a^2 + B11*c2)/B11 - (q*x^2)/(2*B11)

可以用desolve方程: r=dsolve('eqn1','eqn2',...,'cond1','cond2',...,'var'). eqn1表示第1个微分方程,cond1表示第1个初始条件,var表示微分方程中的自变量,默认为t。

这是一个延迟微分方程; MATLAB可以解这类延迟微分方程,但是是数值解法;所以需要之到一个初始条件 x(0)的值; 你能给出x(0)的值我可以帮你解 首先编写关于延迟函数的M文件; function dx=yanchi(t,x,z) tau=z;%定义延迟时间 dx=x*(1-tau);%延...

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