• ベストアンサー

1次遅れ要素のベクトル軌跡: 1/2は何処から?

まずは添付画像をご覧ください。 1次遅れ要素のベクトル軌跡の計算の過程で { p - (1/2) }^2 + q^2 = (1/2)^2 と突然1/2が出てきます。 何処からその1/2が出てきたのかが分かりません。 まるで先にベクトル軌跡の図を知っていて、中心が 1/2 + j0 であるのを知っていたから1/2にしたかのようです。 因みに、問題の前半には、 1次遅れ要素の伝達関数は G(s) = 1/(1 + Ts) であるから、s=jωを代入する。 …とあります。 では、お願いします。

質問者が選んだベストアンサー

  • ベストアンサー
  • info33
  • ベストアンサー率50% (260/513)
回答No.1

>{ p - (1/2) }^2 + q^2 = (1/2)^2 >と突然1/2が出てきます。 >何処からその1/2が出てきたのかが分かりません。 すぐ上で導いたベクトル軌跡の式 p^2 + q^2 = p を変形して円の方程式に直せばいいだけです。 p^2 -p + q^2 = 0 p^2 -p +(1/2)^2 + q^2 = (1/2)^2 (p - 1/2)^2 + q^2 = (1/2)^2 (半径=1/2, 中心=(1/2, 0) の円 , q≦0))

futureworld
質問者

お礼

なるほど、そういう計算でしたか。 納得です。 ありがとうございました!

関連するQ&A

  • ベクトル軌跡について

    周波数伝達関数 G(jω)=(1+j2ω)/(1+j0.125ω) のベクトル軌跡の書き方を教えて下さい。 G(jω)=K/(1+jωT)であれば、実部と虚部の関係からωを消去して、 (Re-K/2)^2+Im^2=(K/2)^2 と円の方程式を求めることができるのは分かります。 が、最初の問題を同様に考えても、ωが消去できずに困っています。 解説には、ゲインと偏角に適当なωの値を代入し図示すると書いてあり、答えは円でした。 出来るなら円の方程式を導き出し、円の中心・半径からベクトル軌跡を求めたいのです。 よろしくお願いします。

  • 円の方程式になるベクトル軌跡

    円の方程式になるベクトル軌跡 G(s)=K/(1+TS),G(jω)=K/(1+jωT) この時、ω→可変にするとG(jω)のベクトル軌跡は円の方程式になるということを証明したいのですが、 どういう手順で行えばいいのでしょうか。 出来れば詳しく、そうじゃなくてもヒントだけでも頂きたいです。

  • ベクトル軌跡(制御工学)について

    ベクトル軌跡(制御工学)について 制御工学を参考書を用いて勉強しているのですが、ベクトル軌跡でよくわからないことがあります。 伝達関数がG(jω)=1/(1+jω(L/R))のとき ∠G(jω)=∠-arctan(ωL/R) となっているのですが、∠G(jω)=∠-arctan(ωL/R)の-がなぜつくのかわかりません。 公式としてそうなっているのでしょうか。 授業ではG(jω)=k/(a+jb)のときの ∠G(jω)=arctan(b/a)と習ったのですが、これは間違いなのでしょうか。

  • ベクトル軌跡について

    大学で制御工学の勉強をしています。 ベクトル軌跡を考える際には、伝達関数G(s)を利用して周波数伝達関数G(jw)を求め実部と虚部に分ける。 ここまでは理解でき、実際にG(s)=1/(RCs+1)のようなsが一次の式の場合は ベクトル軌跡を図示できたのですがG(s)=4/(s^2+2s+4)のように次数が上がった場合には、 どのように考えればいいのでしょうか? 急ぎの質問で申し訳ないのですが、宜しくお願いします。

  • 二次遅れ要素の特性について

    制御工学の二次遅れ要素( 伝達関数G(s)=K/(s^2+2ζωs+ω^2) )に、時定数というものは存在しますか? 参考書にも載っておらず、インターネットで検索してもはっきりしたことが書かれていませんでした よろしくお願いします

  • 制御系の安定性について教えて下さい

     1.伝達関数L(s)のナイキスト線図を描いたとき、ベクトル軌跡が-1+j0を通ると持続振動で、この点より左側を通ると発散してしまいます。この-1+j0の点での周波数をω0とするとL(jω0)=-1と表すことが出来ます。  次に、L(s)の閉ループ伝達関数 G(s)=L(s)/(1+L(s))を考えて、ここにL(jω0)=-1を代入します。 この時、分母が0になるので、G(s)=∞となりこれは発散を意味していると思います。しかし周波数ω0では持続振動になるはずなのにどうして閉ループ伝達関数は∞になるのでしょうか? 2.また、別の伝達関数 P(s)がある周波数ω1の時にP(jω1)=-10となったとします(ナイキスト線図によると不安定)。これをG(s)=P(s)/(1+P(s))に代入するとG(s)=-10/(1+(-10))=1.111・・・となり閉ループ伝達関数は収束してしまいます。しかし一巡伝達関数を考えると位相が180°遅れて、ゲインが10ならば発散すると思うのですが、考え方の何処が間違っているのでしょうか?  説明不足な点は補足致しますので、どなかた教えて下さい。よろしくお願いします。

  • 制御工学

    開ループ伝達関数が L(s)=K/(s-1)だとします。 これのベクトル軌跡を描くときs=jωとして ゲインと位相を求めます。 ゲインは、k/√(ω^2+1) 位相は ∠k-∠(jω-1)ですが ωが0のときの位相はなぜ-180度なのでしょうか? -∠(jω-1)をみて、ベクトル軌跡を描くとき、 まず負の実軸に1進み、虚部は0になるので-180度になり、それれに-がついてるので180度となるのではないのでしょうか? おそらく自分が勘違いしてるのですがなぜかわかりません。 よろしくお願いします。

  • 伝達関数のゲイン図の描き方

    制御工学を勉強中です。 伝達関数G(s) = 1 / ( s * (1 + Ts) )のゲイン図を描きたいのですが、描き方がよくわかりません。 G(s) = G1(s) * G2(s) = {1/s} * {1/(1+Ts)} と分割し、s = jω を代入して Gain1(jω) = -20log10(ω) Gain2(jω) = -10log10(1 + ω^2T^2) となり、Gain(jω) = Gain1(jω) + Gain2(jω) = -20log10(ω)-10log10(1 + ω^2T^2) が得られましたが、これはどうやって図示してやればよいのでしょうか。 Gain2だけに関しては、ω^2T^2の値について (i) << 1 (ii) = 1 (iii) >> 1 の3通りを具体的に考えてプロットしましたが、Gain1はωTを持たないのでうまくGain2と合わせることができず困っています。 (i)(iii)はいいとして、(ii)をω=1/Tと考えてGain1を図示して2つの図を足し合わせればよいのでしょうか。 よろしくお願いします。

  • 周波数応答の計算の一部

    K / {ω(j-Tω)} = -KT / {1+(Tω)^2} - jK / [ ω{1+(Tω)^2} ] の計算過程が分かりません。 実部と虚部に分けたのは分かります。 これは G(s)=K/{s(1+Ts)} のベクトル軌跡を作成する問題です。 ω→+0の場合のための計算です。 必要であれば、補足します。

  • 2次遅れ要素の分解

    伝達関数の2次遅れ要素の分解で困っています。 2次遅れ要素は伝達関数として1/(T1s+1)(T2s+1)と記述されますが、これは1/(T1s+1)*1/(T2s+1) と書きなおせるので、入力Y(s)に対し、G1(s) {1/(T1s+1)}とG2(s){1/(T2s+1)}の直列結合になるかと思います。 1次遅れ要素は微分方程式T*dx/dt+x=yと書け、微分を(x(n)-x(n-1))/dt :dt=サンプル時間、で近似するとx(n)=(dt*y(n)+T*x(n-1))/(T+dt)と変換できるので、入力に対する応答がエクセル上で簡便にできます。 この応答を入力として2段目の伝達関数に入れてやると、すなわち u(n)=(dt*x(n)+T*u(n-1))/(T2+dt),ただしx(n)は1段目の出力 として計算すると2次遅れ要素の応答が計算されるのではないかと考えます。 実際、ステップ応答では上記2次遅れ要素の逆ラプラス変換の解である h(t)=1-T1/(T1-T2)*Exp(-t/T1)-T2/(T2-T1)*Exp(-t/T2) によく合致します。 ところが、これを実システムに置き換えるとずいぶんと実際の値と異なってしまいます。 実システムとは、ヒーターを用いたPIDコントローラーです。ヒーターのステップ入力(電力)に対しヒーター本体の温度(tmp1)と被加熱物質(tmp2)の応答を見てみると、tmp1はヒーターへのステップ入力に対し1次遅れ応答を、tmp2は2次遅れ応答を示し、Excelのソルバーを用いるとtmp1レスポンスの比例定数K1, 時定数T1が、またtmp2レスポンスのK2, T1, T2を求めることができます。 tmp2のパラメーターをソルバーを用いて求める際に、T1はtmp1レスポンスの値に固定してやると、(K2/K1)/(T2s+1)が被加熱物質単体の伝達関数として求められます。(と思っていますが、この辺りが間違い?) 被加熱物質の温度コントロールをしたいのですが、直接被加熱物質の温度をモニターすることが難しいので、何とかヒーターの温度から被加熱物質の温度を計算して、PIDの入力に使えないか、というのがやりたいことです。 すなわちヒーターの温度と上記で求めた被加熱物質単体の伝達関数から被加熱物質の温度を計算したいわけです。 実測値と計算値のグラフを添付いたします。 ここでの計算は次の式で行っています。 BlockTemp = (dt * (K2/K1) * HeaterTemp + T2 * BlockTemp_old ) / (T2 + dt) BlockTemp_old = BlockTemp dt:サンプリング時間 K1,K2:ヒーターとブロックのステップ応答から計算した比例定数 T2:上記方法でエクセルのソルバーを用いて導出した時定数 どこが間違っているのか、またどうしたら正しい値を得ることができるのか、ご指摘いただけると助かります。 よろしくお願いいたします。