ついに夏本番!さぁ、家族でキャンプに行くぞ! >>

パデ近似はテーラー展開を分母と分子にもつようなものですが、
・これの具体的な計算の仕方を教えて下さい。
・テーラー展開に対する利点を教えて下さい。

手元にある本には、「有限の級数展開を無限級数に置き換えることが出来る」、と書かれてあります。ネットで検索してみると「有理関数には基礎演算しか使われないので,数値的評価が非常に簡単である.」
と書かれていますが、一体パデ近似の御利益は何なのでしょうか?
どういうときにこれを使えば良いのでしょうか?

このQ&Aに関連する最新のQ&A

A 回答 (6件)

・これの具体的な計算の仕方を教えて下さい。



たとえば
f(x)=log(x) の(x-1)^6項までTaylor展開で関数近似をする
=x-1-(x-1)^2/2+(x-1)^3/3-(x-1)^4/4+(x-1)^5/5-(x-1)^6/6
+...
これを使ってwxMaximaを使ってパデ近似すると
(%i1) pade(f,6,4);
を実行するとパデ近似関数の5通り計算してくれる。その中の分子と分母の次数差が少ない近似関数を選んでやると
f(x)=(11*x^3+27*x^2-27*x-11)/(3*x^3+27*x^2+27*x+3)
が得られる。

・テーラー展開に対する利点を教えて下さい。
同じ次数のテーラー展開の項数を使用しても、テーラー展開より広い範囲で近似精度が良い。(log(x)のグラフ、x=1の周りのTaylor展開6次打切りとそれに対するPade近似の比較図参照)

> 一体パデ近似の御利益は何なのでしょうか?
パデ近似は、ある変数の広い変化範囲で近似度が一定で誤差が少ない。一方。テーラー展開では展開の中心xから離れるほど近似誤差が大きくなり
、パデ近似と同じ精度の数値計算をしようとするとテーラー展開の項数を相当沢山とらないと同じ程度の誤差にできない。つまり展開の中心から離れるとテーラー展開は急速に近似度が落ち収束が悪くなってしまうがパデ近似は近似次数を上げなくても近似度が落ちないという長所を持っている。

より少ない、項数の有理多項式で、高精度の計算が広い変域に渡って可能である。
グラフを参照して貰えば明らかですね。
「パデ近似の利点について教えて下さい。」の回答画像1
    • good
    • 0
この回答へのお礼

回答ありがとうございます。

大体のことは分かったのですが、やはり気になることがあります。

・wxMaximaを使って計算して下さったようなのですが、一体どういうアルゴリズムというか計算で展開しているのでしょうか?テイラー展開みたいに微分を繰り返していくのでしょうか?

・このグラフを見る限りパデ近似はもとのグラフと完全に重なっており、テイラー展開は用なしに思えてくるのですが、正直パデ近似を使っているものは、いままで2回ほどしか見かけたことがありません。
テイラー展開に対してこれほど有用に思えるのに、それほど利用されないのはなぜなのでしょうか?

・それとできればパデ近似に関して詳述してある書籍などがありましたら教えて下さい。一応学部レベルの物理と数学は勉強したつもりですが、これは習ったことがありません。


よろしくお願いいたします。

お礼日時:2009/05/24 13:28

> テイラー展開は用なしに思えてくるのですが、


> 正直パデ近似を使っているものは、いままで2回ほどしか見かけたことがありません。
> テイラー展開に対してこれほど有用に思えるのに、それほど利用されないのはなぜなのでしょうか?

恐らく、面倒くさいからです。
所詮は、少ない項数で近似精度を上げよう… というミミッちい話ですから、
多項式近似でも目的の精度が得られるなら、式が簡単なほうが好まれる訳です。

その際、テーラー近似を使うよりも、
一般の多項式近似の中から曲線の当てはめを行ったほうが、区間全体での
近似精度は上がります。テーラー近似は、テーラー展開に由来しているために、
どうしても、中心近傍での精度だけが偏重されてしまうのです。
    • good
    • 0

#1、#4です。


>体どういうアルゴリズムというか計算で展開しているのでしょうか?
以下にPade近似の求め方の具体的な手順が載っていますので参考にして下さい。

http://www.finetune.co.jp/~lyuka/technote/fract/ …
http://amath.colorado.edu/courses/7400/2008fall/ …
http://math.fullerton.edu/mathews/n2003/PadeAppr …
↑のsolutionのところをクリックすると詳細があります。
    • good
    • 0

#1です。


A#1の補足質問の回答
>・wxMaximaを使って計算して下さったようなのですが、一体どういうアルゴリズムというか計算で展開しているのでしょうか?テイラー展開みたいに微分を繰り返していくのでしょうか?

微分は繰り返さず未定係数法で係数の連立方程式を解くだけでいいですね。

f_n(x)をf(x)のTaylor展開をn次項までで打ち切った多項式とします。
パデ近似関数の計算原理はPade近似の定義に従って、求めたf_n(x)の係数とパデ近似の有理多項式g_n(x)=p_m(x)/q_r(x),(m+r=n)の係数の関係の連立方程式を解くだけです。つまり、q_r(x)を両辺に掛けた多項式方程式
f_n(x)q_r(x)=p_m(x)
を恒等式と見做して、未定係数法つまり、xの同じ係数間に成り立つ関係式から全ての係数が整数である有理多項式であるp_m(x),q_r(x)を求めてパデ近似式g_n(x)を決定してやれば良いですね。

>・このグラフを見る限りパデ近似はもとのグラフと完全に重なっており、
拡大すれば完全には重なりません。

>テイラー展開は用なしに思えてくるのですが、
>正直パデ近似を使っているものは、いままで2回ほどしか見かけたことがありません。
>テイラー展開に対してこれほど有用に思えるのに、それほど利用されないのはなぜなのでしょうか?

x=aの周りにTaylor展開近似(n次の項までの打ち切ったもの)は|x-a|が大きくなると急速に近似精度が落ちます。近似精度を同じに保つには|x-a|が増加するにつれ、打ち切り項数nを大きくしないといけません。逆に|x-a|が小さいxの範囲(近傍)では打ち切り項数n≦3でも十分実用になります。
つまり,f_3(x)=f(a)+f'(a)(x-a)+f''(a)(x-a)^2/2+f'''(a)(x-a)^3/6
で十分役立ち、手計算でもこの近似式は簡単に出せます。
実際の手計算での関数近似は3項位が限界ですね。
この近似計算は|x-a|<<R(Rは収束半径)という条件が必要です。
その範囲のxにおける近似計算ならTaylor展開近似が有利かと思います。

Pade近似はx=aの近傍でもより広い範囲を対象した近似式のため、その範囲での近似精度を確保するため、分子と分母の多項式の次数を多少多めに捕りますので連立方程式を解く手間がかかります。3次の微係数を求めるより3変数の連立方程式を解く方が手間がかかるかと思いませんか?
計算機で一旦パデ近似関数を求めておいても、有理多項式の分数形式の有理関数のの値を手計算で計算するのも大変ですね。Taylor級数近似式から近似値を計算した方が簡単かと思いますね。

>それとできればパデ近似に関して詳述してある書籍などがありましたら教えて下さい。
書籍でPade近似を勉強していません(多分大学で授業に使う数学や物理の教科書には載っていない)ので紹介できません。数式処理ソフトの関連マニュアルや高度な専門書やそれを取り扱う大学院の授業などでないと詳細には扱っていないでしょう。

内外の大学院の講義資料や数式ソフトの説明や専門家の論文資料などがネットに公開されていますので、#2さんも触れて見えますが「パデ近似」,「pade approximation」、「Pade 近似」などのキーワードでGoogle検索すれば、沢山のサイトが出てきますので、参考になるかと思います。
    • good
    • 0

御理解戴けなかったようですが、No.2 では、


「テーラー展開」は、一点近傍での無限級数による表示、
「テーラー近似」は、それを有限項で打ち切った近似、
「多項式近似」は、テーラー近似とは限らない一般の多項式による近似
の意味で使い分けています。
そのつもりで読み直して戴いたほうが、内容が伝わると思います。

ローラン展開も、テーラー展開も、近似ではなく、
収束円内ではもとの関数と一致するものですから、
誤差は0であり、「精度」という考えが意味を持ちません。

パデ近似とテーラー近似(を含む多項式近似)の比較については、
No.2 に述べました。
    • good
    • 0

以前、他の方の質問に沿って、


Pade 近似について調べてみたことのある者です。
http://oshiete1.goo.ne.jp/qa3885115.html
上記質問内に、参考サイトへのリンクがあります。

まず、
> パデ近似はテーラー展開を分母と分子にもつようなものですが、
というのが、少し違います。
パデ近似は、有理式による近似ですから、分子分母に多項式を持ちますが、
分子分母の多項式は、級数ではなく、「有限次の」多項式です。
この点は、パデ近似の利点と関連して、重要です。

パデ近似は、分子と分母の次数の和が一定以下という条件の下で最良の
有理式近似を与えます。条件内の有理式には、分母が定数式のものも
含まれますから、パデ近似の精度は、(有限次)多項式近似より悪いことは
ありえません。

テーラー近似は、数ある多項式近似の中でも、ある区間での最良ではなく、
一点(展開の中心)近傍での近似精度をねらったものですから、
二乗積分誤差など区間での誤差評価をすれば、当然、
テーラー近似の誤差 > 多項式近似の最良誤差 > パデ近似の最良誤差
となります。

では、どういう場合に、有理式近似のほうが真に良いのかというと、
目的の関数が特異点を持つ場合が、顕著に良いです。

テーラー級数の収束半径が、展開中心から一番近い特異点までの距離であること
を見ても解るように、べき級数展開は、特異点の近傍で収束が遅くなります。
それは、ベキ級数を有限次で打ち切ったものである多項式近似の誤差が大きくなる
ということ。
有理式近似であれば、近似式自体が極を持ちますから、目的関数の極まで
再現することができます。

繰り返しになりますが、有理式近似は、多項式近似をも含むものですから、
有限次の近似を考える限り、多項式近似に劣ることはありません。
多項式近似が最適な場合には、有理式近似を求めた結果が、たまたま
多項式になるだけです。

これに対して、テーラー展開の利点は、式形の単純さにあります。
そのために、理論的な取り扱いがし易い。
次数→∞ の極限を考える際には、有限次近似での精度など関係ありませんからね。
    • good
    • 0
この回答へのお礼

ありがとうございます。

結局はテーラー展開の利点は理論的な取り扱いがし易い、ということだけだということですね。

それと気になったのですが、特異点の持つ場合にテーラー級数が不利だという話ですが、ではパデ近似はローラン展開と比べるとその精度はどうなのでしょうか?

テーラー近似の誤差 > 多項式近似の最良誤差 > パデ近似の最良誤差

の関係いくとローラン展開はテーラー近似に入るのでしょうか?

お礼日時:2009/05/25 00:42

このQ&Aに関連する人気のQ&A

お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!

このQ&Aを見た人はこんなQ&Aも見ています

このQ&Aを見た人が検索しているワード

このQ&Aと関連する良く見られている質問

Qパデ(Pade)近似

皆さんよろしくお願いいたします。

表題の近似計算とはどのようなものでしょうか。
ネットを見ても理解できるものが有りません。
手持ちの教科書(微積分、解析学)には乗っていませんでした。
例題として小職手持ちの制御工学の教科書に、e^(Lx)のパデ近似は次式で表わされるとあります。
どのように導いたのでしょうか。
1次/1次の場合e^(Lx)≒(1-Lx/2)/(1+Lx/2)
2次/2次の場合e^(Lx)≒{1-(Lx)/2+(L^2)*(x^2)/12}/{(1+(Lx)/2+(L^2)*(x^2)/12}
3次/3次の場合e^(Lx)≒{1-(Lx)/2+(L^2)*(x^2)/12-(L^3)*(x^3)/120}/{(1+(Lx)/2+(L^2)*(x^2)/12+(L^3)*(x^3)/120}
当方工学出身ですが、初めてこの近似の存在を知りました。
数学的に、どのような定義でどう導いたかご教示いただければ幸いです。

Aベストアンサー

(1) 概ね賛成です。
r^(k)(0) - f^(k)(0) = 0 は、第 k 次項における係数を一致させること
であり、それを k = 0, 1, 2, …, n+m について行う。』 だと思います。
それが Pade近似の定義そのものです。

「 零点を一致させること 」 = 「 係数を一致させること 」
の方は、意味がよく分かりません。 私の理解としては、
「 近似すること 」 = 「 係数を一致させること 」
といった感じなのですが、コレで伝わるでしょうか?

x = a の近傍において f(x) を r(x) で近似するために、
両者のテーラー展開の係数を 低次(0~n+m 次)の項について一致させる。
それは、x = a を f(x) - r(x) の n+m 位の零点とするのと同じことである。』
とでも言えばよいのかな。

係数を一致させるのは f(x) と r(x) であり、
零点が一致するのは f(x) - r(x) と (Σa_i x^i) (Σq_i x^i) - Σp_i x^i です。
x = a が f(x) や r(x) の零点になる訳ではありません。

(2) そうかも知れません。
有理式の零点と、その分子の零点は、(有理式が既約である限り)同じですが。

(1) 概ね賛成です。
r^(k)(0) - f^(k)(0) = 0 は、第 k 次項における係数を一致させること
であり、それを k = 0, 1, 2, …, n+m について行う。』 だと思います。
それが Pade近似の定義そのものです。

「 零点を一致させること 」 = 「 係数を一致させること 」
の方は、意味がよく分かりません。 私の理解としては、
「 近似すること 」 = 「 係数を一致させること 」
といった感じなのですが、コレで伝わるでしょうか?

x = a の近傍において f(x) を r(x) で近似するために、
両者のテーラー...続きを読む

Q制御工学における無駄時間要素をパデ近似(3次/3次)したときのランプ応答について

皆さんよろしくお願いいたします。
標題の時間応答を求めようとしています。
無駄時間要素を無駄時間Lとすると伝達関数はG(s)=e^(-Ls)と表わせます。
この伝達関数をパデ近似(3次/3次)すると次式で表わせます。
e^(-Ls)={1-(Ls)/2+(L^2*s^2)/10-(L^3*s^3)/120}/{1+(Ls)/2+(L^2*s^2)/10+(L^3*s^3)/120}
この分母と分子の3次方程式を解くと1つの実根aと2つの虚根σ±jωが得られます。
すると上式は次式で表わせます。
G(s)={(s-a)(s-σ-jω)(s-σ+jω)}/{(s+a)(s+σ-jω)(s+σ+jω)}
この場合のランプ応答r(t)={G(s)/s^2}を求めようとしております。
このr(t)の計算結果を計算で求めたご経験のある方がいらっしゃいましたら、
計算方法と結果をご教示いただきたくお願いいたします。
小生の不足している経験でできることして、部分分数に展開したり
留数定理を用いて計算しようとしておりますがうまくいきません。

皆さんよろしくお願いいたします。
標題の時間応答を求めようとしています。
無駄時間要素を無駄時間Lとすると伝達関数はG(s)=e^(-Ls)と表わせます。
この伝達関数をパデ近似(3次/3次)すると次式で表わせます。
e^(-Ls)={1-(Ls)/2+(L^2*s^2)/10-(L^3*s^3)/120}/{1+(Ls)/2+(L^2*s^2)/10+(L^3*s^3)/120}
この分母と分子の3次方程式を解くと1つの実根aと2つの虚根σ±jωが得られます。
すると上式は次式で表わせます。
G(s)={(s-a)(s-σ-jω)(s-σ+jω)}/{(s+a)(s+σ-jω)(s+σ+jω)}
この場合のランプ応答r(t)={G(s)/s...続きを読む

Aベストアンサー

>多項式の正負をどのように判定されているのか.....

単純な目視に過ぎません。
有理式のままではゴチャゴチャして焦点が定まらないので、分子多項式だけ書き出してみましょう。

L=1のとき、原式を因数分解。
  -(s^3/120) + (s^2/10) - (s/2) + 1 = -(s-a)(s^2-bs+c)/120
この両辺に 120 を乗算すれば、
  -s^3 + 12s^2 - 60s + 120 = -(s-a)(s^2-bs+c)  (ac=120)

一方、mathstudy さんの因数分解の結果は、
  s^3 - 12s^2 + 60s - 120 = (s-a)(s-σ-jω)(s-σ+jω) = (s-a)(s^2-bs+c)
これは原式を正負逆転したものですね。
 

Q一巡伝達関数と開ループ伝達関数

一巡伝達関数と開ループ伝達関数は何が違うのでしょうか?
本によって定義がまちまちで、あまり正しい定義がないのかなと思ってしまいますが、ちゃんとした定義が存在するのでしょうか?
インターネットでは一巡伝達関数と開ループ伝達関数は同一視していますが、私の学校の教科書では開ループ伝達関数はフィードバック系を取り除いたときのもの(すなわちC(S)P(S))、一巡伝達関数は閉ループ系を一巡したときのもの(すなわちC(S)P(S)H(S))となっています。

ご存じの方がいたらご教授よろしくお願いします。

Aベストアンサー

教科書の定義が正しいです。

一巡伝達関数は、ループをどこかで切り開いた時に、ループ全体一周する伝達関数で、ループの安定性(位相余裕など)なんかを調べるときに使います。

開ループ伝達関数は、ループをどこかで切り開いた時に、入力と出力の比です。

つまり、ループを切り開いて考えるのは同じですが、一巡伝達関数がループを一周(フィードバックの要素も考える)のに対して、開ループ伝達関数は入力と出力の比です(したがってフィードバックの要素は考えない)。

フィードバックの要素がない場合には、2つは同じになります。

Q一次遅れ系の制御における時定数Tの求め方

計量士の資格を勉強していると自動制御の問題が出てきました。

単位ステップ応答は1-exp(t/T)である。
一次遅れ系の時定数Tの求め方として2つの方法がある。
一つは、1-exp(-t/T)が63.2%になったとき。
もう一つは、過渡応答曲線の原点での接線が定常値に交わるまでの原点からの時間を求める。
とあります。
ここで質問なのですが、この過渡応答曲線とは、1-exp(-t/T)の曲線のことでしょうか?

また単位インパルス応答はexp(-t/T)です。
これが36.8%になったとき時定数Tを求められることは知っているのですが、
同様に、過渡応答曲線の原点での接線が定常値に交わるまでの
原点からの時間を求めると時定数Tを求められるのでしょうか?

Aベストアンサー

>この過渡応答曲線とは、1-exp(-t/T)の曲線のことでしょうか?
そうです。
t=Tとおくと、このときの振幅v(T)=1-exp(-T/T)=1-exp(-1)≒0.6321
と定常値の振幅1に対して0.6321は63.21%にあたります。

>原点からの時間を求めると時定数Tを求められるのでしょうか?
求められます。

過渡応答曲線v(t)=exp(-t/T)に対して、t=0における接線は
u(t)=1-(t/T)ですので、u(t)=0(定常値)になる時間は
1-(t/T)=0からt=T(時定数)が求められます。
このときの振幅はv(T)=exp(-T/T)=exp(-1)≒0.3679
これは定常値(0)までの振幅1に対して36.79%にあたります。

Q固有値が複素数のときの固有ベクトルの求め方

固有値が複素数のときの固有ベクトルの求め方

( -7 -5 )
( 13 9 )

の2x2行列で固有値を求めると 1±2i になると思いますが

Av = λv の形で固有ベクトルを求めようとすると

( -8 + 2i ) x - 5 y = 0
13 x + ( 8 + 2i ) y = 0

の形になり、その先を求めることが出来ません。
何度も計算したので最後の2つの式は間違いは無いと思うのですが、
固有値が複素数の時は、Av = λv の方法で計算することは出来ないということでしょうか?
またどのように計算できるのでしょうか?
お知恵をお貸しいただければ幸いです。

Aベストアンサー

固有値は1±iになるかと…

そこから先の計算は普通に実数の時と同じ方法で計算できます.

Q制御工学における不安定零点の定義

皆さんよろしくお願いいたします。

制御工学において、ボードの定理を勉強中です。
ボードの定理の前提となっている最小位相推移系の1つの条件として
不安定零点を持たない系(もう1つは安定極であること)があります。
不安定零点とは何でしょうか?

教科書では、伝達関数
    K(s-z1)(s-z2)・・・(s-zm)
G(s)=----------     ただし(m≦n)
    (s-p1)(s-p2)・・・(s-pn)
において、不安定零点とはRe(zm)<0と定義していますが、なぜでしょう。
不安定という形容詞が付いているくらいなので、その名の
由来及び定義または、それらを掲載しているURLなどをご存知の方、
いらっしゃいましたらご教示いただきたく、お願いいたします。

Aベストアンサー

不安定零点とは連続系では、「右半平面に存在する零点」です。
つまり、実部が正の零点です。こいつがあると、位相が最小位相にならない。
Gp1(s)=( s+1)/(s^2 + s + 1) と
Gp2(s)=(-s+1)/(s^2 + s + 1) は
ゲイン特性は同じだけど位相特性は違う。
(Gp2の位相はGp1よりも180°も遅れる)

で、このときのGp2の零点を不安定零点と呼ぶ。
この不安定零点の実部は正になっています。

実部が正の極は不安定極と呼ばれますので、
零点の場合も合わせて不安定零点と呼んで
いるのだと思います。
もっと言えば、
不安定零点を持つシステムの逆は
当然ながら不安定極を持ちます。
「逆システムを構成すると“不安定”な極となってしまう“零点”」
ということだと思われます。

なお、このシステムにステップ入力などを入れても
不安定な挙動は示しませんのでご注意ください。
(安定性を支配しているのは極なので)

制御屋さんから見ると、不安定零点があると
アンダーシュートしたり、逆システムを組んだときに
不安定になってしまうことの方が問題になっています。
とはいっても、いろんな対処法があります。

はっきり言って、リファ可能な「由来」については・・わかりません。
この不安定零点は制御をやっている人なら誰でも
知っていますが、由来というのは過去に聞いたこと
無いですね。
数学屋さんに聞いたほうがいいかも知れませんね。

(参考までに)
不安定零点についてある程度書かれた教科書としては、
Goodwin著のControl system design がいいと思われます。

不安定零点とは連続系では、「右半平面に存在する零点」です。
つまり、実部が正の零点です。こいつがあると、位相が最小位相にならない。
Gp1(s)=( s+1)/(s^2 + s + 1) と
Gp2(s)=(-s+1)/(s^2 + s + 1) は
ゲイン特性は同じだけど位相特性は違う。
(Gp2の位相はGp1よりも180°も遅れる)

で、このときのGp2の零点を不安定零点と呼ぶ。
この不安定零点の実部は正になっています。

実部が正の極は不安定極と呼ばれますので、
零点の場合も合わせて不安定零点と呼んで
いるのだと思います。
もっと...続きを読む

QPID制御で 極と零点を置く場所について、

PID制御で 極と零点を置く場所について質問です。

まず1つの零点をてきとうにおいて、最適な痕軌跡を得て、
そのあと極を原点におき、のこりの零点を原点のすぐそばに置くことで、
痕軌跡のかたちをあまり変化させずにすべての条件を満たせる痕軌跡を作れる。

というような感じで授業で習ったのですが、

さいごの

“極を原点におき、のこりの零点を原点のすぐそばに置く”

にかんして、なぜ両方とも原点に置いたらいけないのでしょうか?

その場合、S/S=1になるとおもいますが、その場合、
極と零点はないものと考えるのでしょうか?


よろしくおねがいします。

Aベストアンサー

極零相殺が起きるからですね。
不安定極、不安定零点が相殺(不安定な極零相殺といいます)したら、内部安定性を確保できなくなります。(No.1さんのおっしゃっていることです)

また、そのモードに対して不可制御あるいは不可観測に陥ってしまうので望ましくありません。
(現代制御の話も含みますが、状態方程式の固有値と、伝達関数の極が一致しなくなります。)

Q制御工学の本を購入して勉強をしようと思いますが、種類が多くて、

制御工学の本を購入して勉強をしようと思いますが、種類が多くて、
どれを購入しようか迷っています。
制御工学は、大学時代に履修していましたが、
勉強はほとんどしていなく、初心者同然です。
制御にも電気、機械などその他多くの種類がありますが、電気の制御を勉強したく思っています。
良い参考書、演習問題がありましたら、ご教授ください。

Aベストアンサー

>ちなみに、これはっと思う書物はありますか?
>あれば教えてください。

数学の本を読んで一人で演習問題をやっても面白くないので

Scilabで学ぶディジタル信号処理 CQ出版 ISBN 978-4789830928

がお薦め。
サブタイトルが、「DSPシミュレータで試しながら理解できる」 とある通り、音声で遊べます。
制御とは関係なさそうですが、データの離散処理を学習するにはよいかも。
とにかくデジタル信号処理の、周波数特性、位相、遅れ の3要素を会得することです。そうすればボーデ線図などから、どういうふうに回路やソフトを組めば良いか見えてくることでしょう。

Silabは無料でダウンロードできますから、この本だけで遊べます。

Q伝達関数と周波数伝達関数の違い(関係)

「伝達関数と周波数伝達関数の違い(関係)」ってなんでしょうか?

伝達関数G(s):入力のラプラス変換と出力のラプラス変換したものの比
周波数伝達関数G(jw):入力の複素振幅と出力の複素振幅の比

ぐらいしかいえなんですが、これらがどう関係していてどう違うのかが説明できません。後はせいぜい「sにjwを代入」したら出てくるとか・・・

周波数伝達関数については、「加速度入力に対する加速度応答の比」という説明もあったのですが・・・

Aベストアンサー

#2です。補足します。
Q:これは、これでいいんでしょうか?わたしとしては、G(jw)のほうが複素数な気がするのですが。

A:G(jw)は虚軸で定義されたといった方が適切かも
しれません。(周波数wは実数です)
これを複素平面に拡張したのがG(s)です。
すなわち引数を純虚数以外にも拡張したといったほう
が適切でした。

Q:複素変数に拡張して、ナイキスト線図?とか書くと特異点が出てくるってことでしょうか?

A:G(jw)自体が複素数ということと引数が複素数(a+jw)であるということは違います。
ナイキスト線図はwを0から無限に動かしたとき
G(jw)がどうなるかを複素平面に書いているだけです。
特異点は引数をjwから一般のsに拡張したときに全て現われます。例えば1次遅れ要素はG(jw)=1/(Tjw+1)なので虚軸に特異点はありませんが、G(s)=1/(Ts+1)は実軸にs=-1/Tと
いう特異点があります。(もちろんこのようにGの解析形が分かっているときはG(s),G(jw)のありがたみはおなじです)
減衰のない系の共振の場合は、ある周波数w0で発散し虚軸でも特異的に振る舞いますが、一般には減衰があるので、その場合は周波数wを動かしても特異的には振舞いません。そのときの特異点はr+jw0という虚軸以外のところに隠れています(rは減衰率みたいなもの)。この特異点を(全て)求めれば、伝達関数の形が決定されるということです。

#2です。補足します。
Q:これは、これでいいんでしょうか?わたしとしては、G(jw)のほうが複素数な気がするのですが。

A:G(jw)は虚軸で定義されたといった方が適切かも
しれません。(周波数wは実数です)
これを複素平面に拡張したのがG(s)です。
すなわち引数を純虚数以外にも拡張したといったほう
が適切でした。

Q:複素変数に拡張して、ナイキスト線図?とか書くと特異点が出てくるってことでしょうか?

A:G(jw)自体が複素数ということと引数が複素数(a+jw)であるということは違います。
ナ...続きを読む

Qe^-2xの積分

e^-2xの積分はどうしたらよいのでしょうか…。e^xやe^2xsinxなどはのってるのですがこれが見つかりません。お願いします。

Aベストアンサー

いささか、思い違いのようです。

e^-2x は、 t=-2x と置いて置換してもよいけれど、牛刀の感がします。

e^-2x を微分すると、(-2)*( e^-2x )となるので、

e^-2x の積分は、(-1/2)*( e^-2x )と判明します。


このQ&Aを見た人がよく見るQ&A

人気Q&Aランキング