1

私は二次微分方程式を持っています

 d^2(y(t))/dt^2+ B1*d(y(t))/dt+ c1*y(t)=A1*x(t)

t は秒単位で、0 より大きい値です。A1、B1、C1 は次の定数です。

A1= 3.8469x10^6
B1= 325.6907
C1= 3.8469x10^6

このシステムは線形で、時不変で、カジュアルです。このシステムは H1 と呼ばれます。Matlab を使用して、このシステムのインパルス応答関数 h1(t) と単位ステップ応答関数 g1(t) を計算してプロットしたいと考えています。

4

2 に答える 2

0

MATLAB には dsolve と呼ばれる便利な関数があり、文字列を引数として取ります (ezplot に似ています)。おそらく、定数を手動で入力する必要があります。しかし、それはDEを解決します。

通常、これらの出力、特にこのような出力は厄介なものになるため、少なくとも読みやすいようにするには、 pretty 関数を使用することをお勧めします。

入力してみてください：

解決を助ける

matlab に入って、何が得られるか見てみましょう :)

于 2012-12-03T03:12:13.983 に答える