微分方程式の補完された時変パラメータを確認する方法
2 views (last 30 days)
Show older comments
この質問で扱われている時変パラメータがinterp1によってどのように補完されたかを確認する方法はないでしょうか。
function dydt1 = myODE1(t1,y1)
ft = interp1(tv, ftv, t1,'spline'); % 時刻 t におけるパラメータの値を補間
dydt1 = - ft *y1; % 時刻 t における導関数の値を計算
end
このコードにおけるftを確認したいです。
確認というのはグラフでプロットするなど元の時変パラメータと比較したいということです。
0 Comments
Answers (1)
Naoya
on 10 Jun 2021
一例となりますが、以下のコマンドでどのように補間されるか確認できます。
% interp1 で利用する 代表点 を定義
tv = linspace(0,pi/2,10); % 時間ベクトル(tv)と対応するパラメータ値(ftv)を定義
ftv = sin(tv);
% interp1 で利用する 内挿する点群を定義
t1 = linspace(0,pi/2,1000);
% interp1 の実行 / 結果表示
ft = interp1(tv, ftv, t1,'spline');
plot(tv, ftv, 'o',t1, ft)
0 Comments
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!