How to integrate linear system of vectorial equations?

Question

Anshuman Pal 2021-8-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1437654-how-to-integrate-linear-system-of-vectorial-equations

评论： Star Strider 2021-8-22

Hello,

I have the following system of vectorial equations that describe a curve in space r, its tangent vector t and its normal $n$, parametrised by the arc length s:

Suppose I have the function

given as a vector of values, as well as initial values of r, t and n. Then is there some simple way (or package) for numerically integrating this system of vectorial equations?

Thank you very much!

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Wan Ji 2021-8-20

编辑：Wan Ji 2021-8-20

Use ode45 or other ode tools

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2021-8-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1437654-how-to-integrate-linear-system-of-vectorial-equations#answer_771499

在 MATLAB Online 中打开

syms kappa n(s) r(s) s t(s) r0 t0 n0 Y

Eqn = [diff(r) == t; diff(t) == -r + kappa*n; diff(n) == -kappa*t];

rtn = dsolve(Eqn, r(0)==r0, t(0)==t0, n(0)==n0)

rtn = struct with fields:

t: [1×1 sym] r: [1×1 sym] n: [1×1 sym]

t = simplify(rtn.t, 500)

t =

r = simplify(rtn.r, 500)

r =

n = simplify(rtn.n, 500)

n =

Alternatively:

[VF,Subs] = odeToVectorField(Eqn)

VF =

Subs =

rtnfcn = matlabFunction(VF, 'Vars',{s,Y,kappa})

rtnfcn = function_handle with value:

@(s,Y,kappa)[kappa.*Y(3)-Y(2);Y(1);-kappa.*Y(1)]

Use ‘rtnfcn’ with the approopriate numeric ODE solver (for example 0de45, ode15s) depending on the magnitude of κ.

For example:

sspan = linspace(0,10);                                     % Vector Of 's' Values
initconds = rand(3,1);                                      % Initial Conditions
[s,rtn] = ode15s(@(s,rtn,kappa), sspan, initconds);         % Integrate
figure
plot(s, rtn)
grid

.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Anshuman Pal 2021-8-22

Thank you! A couple more comments:

1) Is there a typo in `ode15s(@(s,rtn,kappa), sspan, initconds)`? Shouldn't it be `ode15s(@rtnfcn, sspan, initconds)`?

2) Can I use `rtnfcn` with a boundary-value solver like `bvp4c`?

Star Strider 2021-8-22

在 MATLAB Online 中打开

As always, my pleasure!

‘Is there a typo in `ode15s(@(s,rtn,kappa), sspan, initconds)`? ’

There is. It should be:

[s,rtn] = ode15s(@(s,rtn) rtnfcn(s,rtn,kappa), sspan, initconds);       % Integrate

‘Can I use `rtnfcn` with a boundary-value solver like `bvp4c`?’

Probably. One addition would be to create ‘kappa’ as an anonymous function, for example with ‘s’ as the independent variable and ‘kapa’ as the dependent variable:

kapamtx = [s(:) kapa(:)];
kappa = @(s) interp1(kapamtx(:,1), kapamtx(:,2), s);

It would be necessary either to keep ‘s’ within the limits of ‘kapamtx(:,1)’ in order to avoid either extrapolating or returning NaN values for ‘s’ outside that range.

.

请先登录，再进行评论。