• ベストアンサー

5次以上の特性多項式の根(伝達関数の極)の求め方

解析的に解けない5次のある代数方程式でも,それを特性方程式に持つ伝達関数の極はMATLABで求めることができます.なぜでしょうか?特性方程式が5次以上の伝達関数の極はどうやって求めるのでしょうか?

  • 科学
  • 回答数1
  • ありがとう数1

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

  • ベストアンサー
  • f272
  • ベストアンサー率46% (8027/17156)
回答No.1

解析的に解かなくても数値的に解けばいいだけですよね。 https://ja.wikipedia.org/wiki/%E6%B1%82%E6%A0%B9%E3%82%A2%E3%83%AB%E3%82%B4%E3%83%AA%E3%82%BA%E3%83%A0 色々な手法があるので,勉強してください。

octopass
質問者

お礼

なるほど数値的に求めているのですね!ありがとうございました!

関連するQ&A

  • 特性方程式に指数関数が入る場合の極を求める方法

    MATLABで特性方程式に指数関数が入る場合の極を求めたいと考えています。 例えば、s^2+s+1+exp(-s)=0のような方程式です。 特性方程式が多項式のみで構成されていないので、 poleコマンド、eigコマンド、rootsコマンドなどは使えません。 ですから、特性方程式を数値計算で解かなければなりません。 ニュートン法で解く場合、初期値を指定して、その周りにある解を探索します。 しかし、特性方程式に指数関数が入りますと、解が無限個あり、 しかも、原点に行くほど極配置が密になるので、 どのように初期値を設定すれば良いか分かりません。 現在は、複素平面上においてメッシュ状に初期値を指定しているため、 探索に非常に時間がかかるうえ、極の取りこぼしがあるかもしれません。 また、多次元のシステムを考える場合、極がどのあたりにあるのか見当をつけること難しいため、 複素平面上をしらみつぶしに調べなければならず、得られた極配置の信頼性が乏しく感じます。 何か、良いアイデアはないでしょうか。 ぜひともご教示頂ければ幸いです。

  • 伝達特性の意味

    あるシステムの入出力特性を表す伝達特性について質問がございます. あるブラックボックス(伝達関数は未知)のシステムに入力信号が入り、出力信号が得られる場合. 伝達関数は、「出力信号のラプラス変換/入力信号のラプラス変換」で示されますよね. 求めた伝達関数のゲイン特性、位相特性、またインパルス応答もここから導き出せる. ここで質問なんですが、これら3つから、そのシステムについてどんな事がいえるのかが分からないんです. ゲイン特性は、LPFのようです. 位相特性は一意に定まらず、振動している. インパルス応答もよくわからない波形を示している. ゲイン特性、位相特性、インパルス応答を求めた後、どうその結果を解釈すればよいか教えてください.

  • 伝達関数の近似方法について

    周波数応答測定器のデータを利用して、伝達関数へ近似する方法を探しています。 一部のアナライザには、極と零を計測してくれる機能があるようですが、あいにく持ち合わせておりません。 MATLABを利用する方法など、手段がありましたら教えてください。 よろしくお願いします。

  • 伝達関数の極の求め方について

    次の伝達関数の極をすべて求めて安定判別をしてみたのですが、やり方と答えがあっているのかがわからないので教えてください。 G(s)=2/{(s+1)(s+2)} 極を求めるためには(s+1)(s+2)=0となればよいので s=-1、-2となる。 ここで、上で求めた値で正の値が1つでもあれば不安定であるため今回は安定となる。 G(z)=1/(1+2.5z^-1+z^-2) 極を求めるためには1+2.5z^-1+z^-2=0となればよいので 両辺にz^2をかける。 z^2+2.5z+1=0 両辺に2をかける。 2z^2+5z+2=0 (z+2)(2z+1)=0 z=-2、-0.5となる。 ここで、上で求めた値で単位円の外に値が1つでもあれば不安定であるため今回は不安定となる。 これであっていますか? 誰かわかる人がいたら教えてください。

  • 伝達関数の周波数特性

    ある伝達関数の入出力の位相差は、周波数伝達関数の虚数部を実数部で割ったものの逆正接関数、とある書物にありました。 正接関数、つまりtanですから、この逆関数の取り得る値の範囲は-π/2からπ/2までです。 つまり、伝達関数の周波数特性としての入出力の位相差の絶対値は、高々π/2ということになりますが、私の理解は正しいでしょうか? 高々π/2ではなく、これを超えることがありえますか?

  • PLL 伝達特性

    PLL(Phase Locked Loop)について、勉強しているのですが、その中で、「伝達特性」という言葉がでてきます。「伝達特性」って何か?わかりません。「伝達特性」、何kHzってでてきます。「伝達特性」、何kHzは何を表しているのか、ご教示の程、宜しくお願い致します。 以上、宜しくお願い致します。

  • 伝達特性がどうやったら求まるのかわかりません、助けてください。

    伝達特性がどうやったら求まるのかわかりません、助けてください。 伝達特性がどうやったら求まるのかわかりません、助けてください。 次のRCの一時のローパスフィルタについて設問に答えよ。 この回路の伝達特性(周波数特性)を示しなさい。 周波数範囲;10Hz~10MHz RCフィルタ回路で30kΩの抵抗と1000pFの直列回路です。

  • DSP:伝達関数H(z)を求める目的は?

    差分方程式だけ立てればScilab上で、 伝達関数H(z)を使わず、差分方程式だけで 信号処理プログラムが書けるがわかりました。 なのに何故にZ変換で伝達関数H(z)を求める 必要があるのでしょうか? 今書いているプログラムではScilab上で、 伝達関数の式は出てこず、差分方程式だけ 使っています。 利得の式や周波数特性の式を計算するのに 伝達関数H(z)が必要ということですか?

  • 伝達関数

    伝達関数G(s)=2/{s(s+1)}の 周波数伝達関数G(jω) = -2(ω+j)/[ω{ω^(2)+1}] というのは理解できるのですが、 ゲイン特性|G(jω)| = 2/ω√{ω^(2)+1} 位相特性∠G(jω) = arctan(1/ω) というのは、どのように導出したのでしょうか。

  • 極と零点から伝達関数をつくる際の質問です。

    極と零点から伝達関数をつくる際の質問です。 z平面上で 極が二つ z = 2と4 に在り 零点が一つ z = 3 に在るケースを考えた場合、 二つの伝達関数: 伝達関数H(z) = (係数)* (z-3 )/ {(z-2)*(z-4)} ---<1> 伝達関数H(z) = (係数)* (1-3/z) / {(1-2/z)*(1-4/z)} --- <2> がおもいつくのですが、 この2つの式<1>と<2>は同じものですか? 式<1>=式<2> とならないのですか? 式<1>を式<2>の形に変形しようとすると H(z) = (係数)* Z*(1-3/z) / { z^2*(1-2/z)*(1-4/z) } = (係数)* (1-3/z) / { z*(1-2/z)*(1-4/z) } となり,分母にzが残ってしまい 式<1>は <2>の式とはちがうものになってしまいます。 そもそも式1 と式2は根本的に違うものなのでしょうか? その場合、どちらで考えるのがいいのでしょうか? よろしくおねがいます。