RK4 Windkessel误差的方法

3视图(30天)

显示旧的评论

oumi k 2022年10月21日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1832608-rk4-method-for-windkessel-error

评论道: oumi k2022年10月22日

答:接受詹姆斯Tursa

我有这些错误

试图访问P (1.16667);索引必须是一个正整数或逻辑。

@误差(t, P, Q) (- P (t) / (R * C)) + (Q (t) / C)

错误RK4_WINDKESSEL(24)行

k2 = f (t(我)+ 0.5 * k1 * h P (i) + 0.5 * k1,问(我));

清除所有

关闭所有

clc

%的参数

R = 1.5;%的阻力

C = 4;%的合规

h = 1;%步长

t = 1: h: 100;%的时间向量

P (1) = 1;%初始条件

问(1)= 1;%初始条件

f = @ (t, P, Q) (- P (t) / (R * C)) + (Q (t) / C);

% RK4循环

i = 1:装天花板(100 /小时)

如果(我> = 1 | |我< = 5)

问(i) = i + 1;%坡道

结束

Q (+ 5) = Q(我);%坡道

t (i + 1) = t(我)+ h;

k1 = f (t (i)、P (i)、问(我));

k2 = f (t(我)+ 0.5 * k1 * h P (i) + 0.5 * k1,问(我));

k3 = f (t(我)+ 0.5 * k2 * h P (i) + 0.5 * k2);

k4 = f (t(我)+ k3 * h P (i) + 0.5 * k3);

P (i + 1) = P (i) + (1 / h) * 6 * (k1 + 2 * k2 + 2 * k3 + k4);

结束

情节(t, P)

包含(“时间”)

ylabel(压力)

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

接受的答案

詹姆斯Tursa 2022年10月21日

1
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/1832608-rk4-method-for-windkessel-error answer_1080923

编辑:詹姆斯Tursa 2022年10月21日

多个错误。当你处理调用你的函数f,它假定P和Q值在时间t。他们不是通过函数或函数句柄。所以就用P和Q,非P (t)和Q (t):

f = @ (t, P, Q) (- P / (R * C)) + (Q / C);

然后在RK4代码中,k的不是t计算的一部分。所以看起来应该像这样呼吁t部分:

                              k1 = f (t (i)等);
                             
                              k2 = f (t(我)+ 0.5 * h等);
                             
                              k3 = f (t(我)+ 0.5 * h等);
                             
                              k4 = f (t(我)+ h等);