資源簡介
西交大計算方法上機可用。RK-4
該算法用Runge-Kutta 法解初值問題。a,b—求解區間的左右端點,y 0 —初
值,步長選用h。依次由公式計算出K1,K2,K3,K4 的值并計算出出每步的y(n)。
同時每算依次就更新依次t(n)的取值,最終輸出t i 和y
代碼片段和文件信息
%%Runge?Kutta
clear;
h=1/128;???%步長?1/81/321/128
n=1/h;
%初值
y=-1;y_2=3;y_3=2;
y=[y;y_2;y_3];
i=1;
for?t=0:h:1-h
???t1=t;
???y1=y;
???f1=[0?1?0?0;0?0?1?0;-1?1?1?2]*[y1;t1]+[0?0?-3]‘;
???K1=h*f1;????%計算k1
???t2=t+h/2;
???y2=y+1/2*K1;???
???f2=[0?1?0?0;0?0?1?0;-1?1?1?2]*[y2;t2]+[0?0?-3]‘;
???K2=h*f2;????%計算k2
???t3=t+h/2;
???y3=y+1/2*K2;
???f3=[0?1?0?0;0?0?1?0;-1?1?1?2]*[y3;t3]+[0?0?-
評論
共有 條評論