请教一下
这个怎么用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 outlook无法配置成功
- ¥15 Pwm双极模式H桥驱动控制电机
- ¥30 这是哪个作者做的宝宝起名网站
- ¥60 版本过低apk如何修改可以兼容新的安卓系统
- ¥25 由IPR导致的DRIVER_POWER_STATE_FAILURE蓝屏
- ¥50 有数据,怎么建立模型求影响全要素生产率的因素
- ¥50 有数据,怎么用matlab求全要素生产率
- ¥15 TI的insta-spin例程
- ¥15 完成下列问题完成下列问题
- ¥15 C#算法问题, 不知道怎么处理这个数据的转换