常微分⽅程初值问题:多步预测-修正⽅法[MATLAB] #先上代码后补笔记#
关注安#可以直接复制粘贴调⽤的MATLAB函数代码!#
港琪月饼
1. 亚当斯(Adams)预测-修正算法山东巡抚
由亚当斯-巴什福特(Adams-Bashforth)显式预测公式和亚当斯-莫顿(Adams-Moulton)隐式修正公式组成的预测-修正(PECE)对。function [ YMat ] = Adams( func, tvec, y_init, order )
% Adams预测-修正算法,⽤于求解常微分初值问题
头像男动漫呆萌可爱% 输⼊四个参数:函数句柄func(接收列向量、返回列向量),积分时间列向量tvec,初值⾏向量y_init,阶数order;
% 输出⼀个参数:数值解,每⼀⾏对应积分时间列向量的⼀⾏,各列为变量⼀个分量。
switch order
ca '4'工作焦虑
row = size(tvec, 1); col = size(y_init, 2);
YMat = zeros(row, col);
YMat(1:4, :) = Runge_Kutta(func, tvec(1:4), y_init, '4');
for i=4:row - 1
stepsize = tvec(i + 1) - tvec(i);
ydiff0 = func(tvec(i), YMat(i, :).');
ydiff1 = func(tvec(i - 1), YMat(i - 1, :).');
ydiff2 = func(tvec(i - 2), YMat(i - 2, :).');
ydiff3 = func(tvec(i - 3), YMat(i - 3, :).');
y_predict = YMat(i, :).' + (55*ydiff0 - 59*ydiff1 + 37*ydiff2 - 9*ydiff3)*stepsize/24;
y_corrector = YMat(i, :).' + (9*func(tvec(i + 1), y_predict) + 19*ydiff0 - 5*ydiff1 + ydiff2)*stepsize/24;
YMat(i + 1, :) = y_corrector.';当兵的歌
英语形容词有哪些
end
end
end
国庆快乐