请教一下
这个怎么用Matlab编程求解这个非线性常微分方程
用龙格库塔法
1条回答 默认 最新
- 一大岐 2022-06-02 11:54关注
调用这个方法
function [T,X,dX] = ODE_RK4( Hfun,t,h,x0 ) T = t(1):h:t(2); % 计算 N = length(T); x0 = x0(:); x0 = x0'; % 初值变为行向量 m = length(x0); % 状态量维数 X = zeros(N,m); % 初始化状态量 dX = zeros(N,m); % 状态导数 X(1,:) = x0; for k = 2:N h = T(k) - T(k-1); K1 = Hfun( T(k-1) , X(k-1,:)' ); K2 = Hfun( T(k-1)+h/2 , X(k-1,:)'+h*K1/2 ); K3 = Hfun( T(k-1)+h/2 , X(k-1,:)'+h*K2/2 ); K4 = Hfun( T(k-1)+h , X(k-1,:)'+h*K3 ); X(k,:) = X(k-1,:)' + (h/6) * ( K1 + 2*K2 + 2*K3 + K4 ); dX(k-1,:) = (1/6) * ( K1 + 2*K2 + 2*K3 + K4 ); end dX(N,:) = Hfun( T(N),X(N,:) ); end
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 无用评论 打赏 举报
悬赏问题
- ¥15 2020长安杯与连接网探
- ¥15 关于#matlab#的问题:在模糊控制器中选出线路信息,在simulink中根据线路信息生成速度时间目标曲线(初速度为20m/s,15秒后减为0的速度时间图像)我想问线路信息是什么
- ¥15 banner广告展示设置多少时间不怎么会消耗用户价值
- ¥16 mybatis的代理对象无法通过@Autowired装填
- ¥15 可见光定位matlab仿真
- ¥15 arduino 四自由度机械臂
- ¥15 wordpress 产品图片 GIF 没法显示
- ¥15 求三国群英传pl国战时间的修改方法
- ¥15 matlab代码代写,需写出详细代码,代价私
- ¥15 ROS系统搭建请教(跨境电商用途)