问ode45内插
EN

Stack Overflow用户

提问于 2014-10-29 18:47:51

回答 1查看 1.3K关注 0票数 0

我有一个具有离散值的向量，需要传递到我的ODE系统中，我想使用ode45命令。当我使用这个向量时，它需要在求解器内进行插值。有什么办法可以做到吗？

我有一个线性的耦合系统。

dxdt = f(t)*x + g(t)*y
dydt = g(t)*x + h(t)*y

我有三个向量f( t )，g(t)和h(t)作为t的函数，我需要把它们传递给求解者。

我可以用C或C++编写解算器的代码。有人建议，在Matlab中这样做会更快。有人能建议一个办法吗？

matlab

ode

回答 1

Stack Overflow用户

回答已采纳

发布于 2014-10-29 18:56:47

当然，您可以使用interp1。假设向量fvec、gvec和tvec分别包含f、g的值以及这些值的时间点，则可以将导数函数定义为：

dxydt = @(t,x) [interp1(tvec, fvec, t) * x(1) + interp1(tvec, gvec, t) * x(2)
                interp1(tvec, gvec, t) * x(1) + interp1(tvec, hvec, t) * x(2)];

并在ode45中使用

[T,Y] = ode45(dxydt, tspan, [x0; y0]);

票数 2

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/26638326

复制

相似问题

问ode45内插
EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问ode45内插EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问ode45内插
EN