matlab解偏微分方程组如何编程 matlab中用obe45求解微分方程组?

[更新]
·
·
分类:行业
3195 阅读

matlab解偏微分方程组如何编程

matlab中用obe45求解微分方程组?

matlab中用obe45求解微分方程组?

ode45是最常用的求解非刚性微分方程(组)的指令,它采用变步长四、五阶Runge-Kutta法,求解精度比ode23、ode15要高。但ode45的计算量比较大,而ode23计算量小,且误差大
调用格式:[t,y]ode45(odefun,tspan,y0)
odefun 用以表示f(t,y)的函数句柄或inline函数,t是标量,y是标量或向量;
tspan 如果是二维向量[t0,tf],表示自变量初值t0和终值tf;如果是高维向量[t0,t1,t2,。。。,tn],则表示输出结点列向量;
y0 表示初值向量;
t 表示结点列向量(t0,t1,t2,。。。,tn)^T;
y 表示数值解矩阵,每一列对应y的一个分量。
实例说明:
例1 yy-2t/y,y(0)1,0t4
odefuninline(y-2*t/y,t,y);
[t,y]ode45(odefun,[0,4],1)
plot(t,y,o-)
例2 x-x^3-y,x(0)1;yx-y^3,y(0)0.5,0t30
首先,建立自定义函数文件myfun.m
function fmyfun(t,x)
f(1)-x(1)^3-x(2);
f(2)x(1)-x(2)^3;
ff(:)
然后,在命令窗口中输入
[t,y]ode45(myfun,[ 0 30],[1:50]])
plot(t,x(:,1) ,t,x(:,2),o-)
figure
plot(x(:,1) ,x(:,2),:)

matlab中偏微分表达式如何表示?

偏微分函数可用diff(z,x)。z表示函数,x是自变量

matlab编程求解微分方程?

示例:
1. 求解微分方程 y 2xy xe-x2
syms x y ydsolve(Dy 2*x*yx*exp(-x^2),x)
2.. 求微分方程 xy y - e x 0 在初始条件 y (1) 2e 下的特解并画出解函数的图形.
syms x y ydsolve(x*Dy y-exp(1)0,y(1)2*exp(1),x)ezplot(y)