Question

matlab解二阶微分方程 方程：y''=680000(y-60)/(46.05*(1.5-x)) 边界条件：y(0)=70; y'(1.5)=0;

Answer 1

matlab解二阶微分方程 方程：y''=680000(y-60)/(46.05*(1.5-x))， 边界条件：y(0)=70; y'(1.5)=0。可以用数值方法来求解，如变步长的Runge——Kutta法。求解思路：

1、自定义微分方程函数

function f = ode_fun(x,y)

f=[y(2);680000*(y(1)-60)/(46.05*(1.5-x))];

end

2、调用ode45函数

ode45函数使用格式:

[X,Y] = ode45(@ode_fun,tspan,y0)

tspan——自变量初值

y0——初值

运行结果