% 廢話(huà)不多說(shuō)，直接復(fù)制就可以運(yùn)行，文章底部有參考例子，不明白就看看

function [ yy ] = R_K_4( f , y0 , x0 , xn , hh )

% f 是inline function的句柄

% y0 是初始值

% [x0 xn]是計(jì)算范圍

% hh 是步長(zhǎng)

%% ===輸入判斷===

if (4 == nargin) % == 當(dāng)沒(méi)有給步長(zhǎng)的時(shí)候使用默認(rèn)步長(zhǎng)==

h = 0.1 ;

elseif (5 == nargin)

h = hh ;

end

%% ===參數(shù)初始化(可以修改，可從網(wǎng)上找相關(guān)參數(shù)的值)，以下是通用參數(shù)===

c = [1/6 1/3 1/3 1/6] ;

r = [0.5 0.5 1] ;

u = [0.5 0 0.5 0 0 1] ;

%% ===

x = x0:h:xn ; % == 計(jì)算節(jié)點(diǎn)==

y = y0 ;

N = length(x) ; % == 節(jié)點(diǎn)數(shù)==

for k = 1:N-1 % == 以下就直接按書(shū)本公式寫(xiě)的，應(yīng)該容易理解== kk(1) = f( x(k) , y(k) ) ;

kk(2) = f( x(k)+r(1)*h , y(k) + u(1)*h*kk(1) ) ;

kk(3) = f( x(k)+r(2)*h , y(k) + u(2)*h*kk(1) + u(3)*h*kk(2) ) ;

kk(4) = f( x(k)+r(3)*h , y(k) + u(4)*h*kk(1) + u(5)*h*kk(2) + u(6)*h*kk(3) ) ;

y(k+1) = y(k) + h* ( c(1)*kk(1) + c(2)*kk(2) + c(3)*kk(3) + c(4)*kk(4) ) ;

end

總結(jié)

以上是生活随笔為你收集整理的四阶龙格库塔matlab计算例题,四阶龙格库塔法matlab实现的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。