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

2025-04-13 22:33:40
推荐回答(1个)
回答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——初值

运行结果