• ベストアンサー

非線形分布定数システム→線形時変システム

非線形分布定数システム dx  d --=--[A(t,z)x(t,z)]+B(t,z)u(t,z) dt dz を,いじくりやすいように線形化したいのですが,z(位置)についてのみ差分化し, dx’ --=A’(t)x’(t)+B’(t)u’(t) dt として線形時変システムとして扱うことは,数学的に良いのでしょうか? 差分化に伴い,xやAなどはx’やA’などに記述し直されるとします. 難しい数学でいじくる非線形分布定数システムと,線形システムと, 別個に解説する本は何冊か持っているのですが,具体的な線形化の妥当性については 解説がなく,上記の手続きが不安になりました. 宜しくご教授下さいませ.

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

  • ベストアンサー
  • mmky
  • ベストアンサー率28% (681/2420)
回答No.1

mmkyです。難しすぎてちょっとだめだけど、参考になるかもと思い書き込みます。参考にならなかったらごめん。 dx/dt=d{[A(t,z)x(t,z)+B(t,z)u(t,z)}/dz =dI/dz とすれば、(dx'/dt)=(dI'/dz) は可能性がある? (dz/dz)=1、(dt/dt)=1 とすると、 (dx/dt)(dz/dz)=dI/dz =(dx/dz)(1/dt)dz=dI/dz=I' x'/dt=I'/dz {x'/dt}={I'/dz} おのおのを新しい関数と考えれば、 {F}={G} d{F}=d{G} d{F}={d{F}/dx'}dx'+{d{F}/d(dt)}d(dt) =(1/dt)dx'+(-x'/(dt)^2)d(dt) d{G}=(1/dz)dI'+(-I'/(dz)^2)d(dz) d{F}-d{G}=0 (1/dt)dx'+(-x'/(dt)^2)d(dt)-(1/dz)dI'-(-I'/(dz)^2)d(dz)=0 {(1/dt)dx'-(1/dz)dI'}+{(I'/(dz)^2)d(dz)-(x'/(dt)^2)d(dt)}=0 ここで、{(I'/(dz)^2)d(dz)-(x'/(dt)^2)d(dt)}→0 (なぜならd(dz)=(dz/dz)dz=dz, d(dt)=(dt/dt)dt=dt, I'/dz=x'/dt ゆえ) だから、 (1/dt)dx'-(1/dz)dI'=0 (dx'/dt)=(dI'/dz) 参考にならないか?

First_Noel
質問者

お礼

雑用多々で返事が遅くなってしまいました.スミマセン.. ご意見をもとに進めていきたいと思います. お教え頂きましてありがとうございました.<(__)>

関連するQ&A

  • 線形化 状態方程式

    非線形システムdx/dt=f(x,u), f(x,u)=(2・x1・x2+u ,x1+x2^2)T x=(x1,x2)T テーラー展開一次近似により、平衡点(x',u'), x'=(x1',x2')T に関して局所線形化したX=x-x',U=u-u'に関する状態方程式 dX/dt=AX+BUを求めよ。 上記の問題で平衡点をx0=(-1,1)T u0=2としてその近傍でテーラー展開しようとしたのですが、AとBの計算の仕方がわかりません。 お手数ですが、教えていただけないでしょうか? 似たような問題があるサイトがあれば教えてください

  • 非線形対象システムの線形化について

    非線形対象システムの線形化をおこなってください.ただし入力u(t) 出力y(t) dy(t)/dt + 10[y(t)]^2 + 20 = 30u(t) という問題を解く事が出来ません.分かる方 解答お願いします.

  • (信号処理で) 線形微分方程式で係数が定数のシステム について

    (信号処理で)  線形微分方程式で係数が定数のシステム について 1.それが 線形であること 2.そして、時不変であること これらはどのようにして証明というか、確認できるのでしょうか? たとえば入力がx(t)で出力がy(t)とすると 一階の線形微分方程式で係数が定数の場合の例: dy(t)/dt + ay(t) = bx(t) このシステムを考えたとき、どのようにして線形 そして時不変であることを確認できるのでしょうか? 上の式の形のまま確認する方法はあるのでしょうか? それとも、 先に解いて、y(t) = ~~~ の形にしてから確認するのでしょうか? (とくに、時不変になるという確認方法が知りたいです。) よろしくおねがいします。

  • 不定積分の計算で出た定数は捨てて良いのでしょうか

     46歳の会社員です。思うところがあって、1 年前から数学を独学で勉強しています。  非常にレベルが低い質問をしているのかもしれませんが、周りに聞ける人がいないのでここに質問をすることにしました。  不定積分の計算で出てきた定数は積分定数と扱って捨ててよいのでしょうか ?  例えば、 ∫(x + 1)^2 dx ((x + 1)の 2乗を積分) を ∫(x^2 + 2 * x + 1) dx に変形すると、 x^3 / 3 + x^2 + x になりますが、 x + 1 = t とおいて ∫t^2 dt に変形すると、 x^3 / 3 + x^2 + x + 1 / 3 となり、定数 1 / 3 が出てきます。  また、 ∫{2 / (2 * x + 2)} dx を ∫{1 / (x + 1)} dx に変形すると、 log|x + 1| になりますが、 2 * x + 2 = t とおいて ∫(2 / t) * (1 / 2) dt に変形すると、 log|2 * x + 2| になります。  これを log|2 * x + 2| = log|(x + 1) * 2| = log|x + 1| + log|2| と変形すると、定数 log|2| が出てきます。  これらの定数は積分定数として扱って捨ててよいのでしょうか ?

  • 線形システムの問題で・・・

     先日、学校で線形システムのこんな問題が出されたのですが、よくわかりません。 どなたか教えていただけませんか。よろしくおねがいします。    [x1(k+1)] [a 0][x1(k)] [1]  [x2(k+1)]=[0 b][x2(k)]+[1]u(k) という2次線形システムついて、状態推移行列と インパルス応答x(k)を求めよ(a,bは定数で0<a<1,0<b<1)という問題なんです。 x1とか書いてるのは、一乗、二乗という意味ではなく、ただのx(k)を区別する添え字です。 ちなみに、縦に並んでいるかぎ[は、2次の行列なのでひとつとしてみてください。

  • 統計学、同時分布、密度関数の和について

    統計学、同時分布、密度関数の和について 統計の問題で解答を読んで、よく分からない箇所があります。ぜひご教授ください。問題、解答は以下の通りです。 [問題]X,Yは、ともに正の値のみをとる連続型の確率変数でそれらの同時分布はf(x,y)である。このときU=X+Yの密度関数を求めよ。 [解答]P(U<=u) = ∬[x+y<=u]f(x,y)dxdy であるが、積分変数をxとyからxとt=x+yに変換すると P(U<=u) = ∫[0,u]dt∫[0,t]f(x,t-x)dx      …(1) これをuで微分することにより f(u) = ∫[0,u]f(x,u-x)dxが得られる。 …(2) ∫[,]で[]は積分の区間を表しています。 質問内容は2つです。 (1)のところで変数変換をすると、 ∫[0,u]∫[0,t]f(x,t-x)dxdtとなります。なぜ∫[0,u]dtが∫[0,t]f(x,t-x)dxと別々にできるのかが分かりません。 また(2)のところでも なぜ∫[0,u]dt∫[0,t]f(x,t-x)dxを微分して∫[0,u]f(x,u-x)dxが得られるのかも良く理解できていません。 どちらか一方でも分かる方はご教授お願いいたします。

  • 微分方程式

    dx/dt=a^2-x^2 (aは実数の定数) (1)この微分方程式は1階の線形同次・線形非同次・非線形のどれにあてはまるか。 (2)この微分方程式の一般解を変数分離法で求めよ。 考えたことは(1)は非線形だと思いますが、合っていますか? (2)はdx/(x^2-a^2)=-dtと変形し、両辺積分します。  すると、1/(2a)log(|x-a|/|x+a|) = -t + C このあとx=が分からないです。 教えてください。お願いします

  • 分布関数

    参考書の解説一部抜粋 F(x^2+y^2+z^2)=F(x^2)F(y^2)F(z^2)のとき、この式の関係は exp(a+b+c)=(e^a)(e^b)(e^c) のように満たされる。したがった、分布関数は F(x^2)=Aexp(-kx^2) の関数形となる。Aとkは定数 定数Aは規格化条件 ∫[-∞,∞]f(x)dx=1 となる。 F(x^2)=Aexp(-kx^2)を規格化すると A∫[-∞,∞]exp(-kx^2)dx=A(π/k)^(1/2)=1 A=(k/π)^(1/2) 質問 (1)なぜ、F(x^2)=Aexp(-kx^2)の関数形となるのですか? そもそも、なぜF(x^2+y^2+z^2)=F(x^2)F(y^2)F(z^2)のとき、この式の関係はexp(a+b+c)=(e^a)(e^b)(e^c)のように満たされるのですか? (2)なぜ、A∫[-∞,∞]exp(-kx^2)dx=A(π/k)^(1/2)=1となるのですか? πはどこからでてきたのですか?

  • 2変量の確率分布について

     統計学の勉強をしています。一様分布における2変量の確率変数についてわからなくなったので質問させてください。  一様分布の確率密度関数はfx(x)=1/(b-a)ですが、b=1,a=0とするとfx(x)=1となりますよね。  このことを踏まえて2変量t=x+y(yもxと同様の一様分布の確率でxとyは独立)を定義して、その確率密度関数はf(t)=∫fx(x)fy(t-x)dxで与えられますよね(ここで間違っていたならすみません…)  そこでこの関数にfx(x)=1,fy(y)=1を代入して∫範囲を0から1(dxで積分ですから)に設定して積分をするとf(t)=1となってしまいました。  このままtにおいてtの期待値を求めると∫(0,2)tf(t)dt=2となりました。(積分範囲はdtについてですから0から2までとしました)  しかし、よく考えてみると0から2までの範囲の一様分布でその期待値となるのは普通1じゃないかと思います。  計算が間違っているのか、そもそも考え方が違うのか、わかる方がいらっしゃったら、ご教授していただけませんでしょうか?よろしくお願いします。

  • 非斉次な2階の線形微分方程式

    fを実数として、非斉次な2階の線形微分方程式を積分定数を用いて解け。(xはtの関数) (d^2x/dt^2)+9x=fe^(2it) x=Ae^(2it)とすると、 A=f/5 ここからどうすればよいかわかりません。 詳しい解説お願いします。