1

実はまた同じ質問をしている

Dx=y

Dy=-k*yx^3+9.8*cos(t)

inits=('x(0)=0,y(0)=0')

これらは、私がプロットしたかった微分方程式です。

まず、微分方程式を解いてグラフをプロットしてみました。

Dsolve('Dx=y','Dy=-k*yx^3+9.8*cos(t)', inits)

このように、しかし、このシステムには明確な解決策はありませんでした。

今私は立ち往生しています:(

方程式を解かずに、このシステムをどのようにプロットできますか?

Chris Taylor は、y 対 t および x 対 t のグラフをプロットして、私の質問を解決しました。

しかし、y対xグラフをプロットする方法を知りたかったのです（これは正しいですか？縦軸y、横軸xのグラフを意味しています）

私を助けてください

4

1 に答える 1

3

以下のように Chris Taylor の回答のコードを使用します。

k = 1;
f = @(t,x) [x(2); -k * x(2) - x(1)^3 + 9.8 * cos(t)];
tspan = [0, 10];
xinit = [0, 0];
[t, y] = ode45(f, tspan, xinit);

y の 1 列目と 2 列目にそれぞれ格納された 0 から 10 までの t の関数として値 x と y を与えます。したがって、単に呼び出すだけで y(x) をプロットできます

plot(y(:,2),y(:,1));
于 2013-04-16T18:27:12.597 に答える