syms xy=dsolve('Dy=x*y','y(0)=1','x')Y=subs(y,x,0:0.1:4);plot(0:0.1:4,Y)
>> syms r k positive>> dsolve('Dy = r*x*(1-x/k) - x^2/(1+x^2)')ans =C2 - (t*(k*x。
K3=f(t+h/2,X+h/2*K2);K4=f(t+h,X+h*K3);Y=X+h/6*(K1+2*K2+2*K3+K4) ;end以上是4阶龙格库塔法的代码:自己写函数,存为f. mfuncti。
如果是常微分方程,可以用dsolve函数。该函数可以解单变量常微分方程或者多变量常微分方程组,所以5个变量也不在话下。 调用格式如下: [y1,...,yN] = dsolve(eq。
好像没什么特别好的办法: >> syms t q y; q=4*t^3; >> syms dy >> eq=dy-(q-2*t); >> y=dsolv。
知道初始值就更好了.clc;clear dx=inline('(0.5-x/300)/(0.5-x/600)*x*(1-x/t)')[t,x]=ode45(dx,[1/600,1/30。
偏微分函数可用diff(z,x)。z表示函数,x是自变量 偏微分函数可用diff(z,x)。z表示函数,x是自变量
用Matlab编程计算二阶常微分方程组的数值解的方法有很多种,最常用有ode45函数。ode45函数的使用格式:【t,y】=ode45(odefun,tspan,y0)>> odefun=。
syms A B C x y;>> dsolve('Dy=A-B*y-C*x*y','x')ans =-1/2*A*erf(1/2*(-2*C)^(1/2)。
改动: function S=history(t) S=ones(3,1); 改为 S=ones(4,1); OVER!PS.多翻翻help!
猜猜你还想问: | ||
---|---|---|
MATLAB求解微分方程 | matlab代码大全 | matlab怎么保存代码 |
matlab基本代码 | python求微分方程 | matlab phase函数 |
matlab运行 | matlab求解矩阵方程 | matlab心形函数代码 |
回顶部 |