• ベストアンサー

変数変換

下記の変数変換はあっているでしょうか?また、論文から引用した式なのですが、どうやって変数変換を行っているか、途中式が分かりません。分かる方がいましたら、教えてください! 1/y*(d^2x/dφ^2)+y*(d^2x/dψ^2)+d/dφ*(1/y)*(dx/dφ)+(dy/dψ)*(dx/dψ)=0 ξ=tanhaφを用いてφ=-∞~+∞をξ=-1~+1に変換すると上式から、次式が得られる。 -2a^2*ξ(1-ξ^2)*(dx/dξ)+(a(1-ξ^2))^2*(d^2x/dξ^2)+y*(d^2x/dψ^2)-(1/y)*(a(1-ξ^2))^2*(dy/dξ)*(x/dξ)+y*(dy/dψ)*(dx/dψ)=0

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

  • ベストアンサー
  • Akira_Oji
  • ベストアンサー率57% (45/78)
回答No.2

変数変換後の式にはエラーがあると思います。 ご存知のようにξ=tanh(aφ)を用いて(微分のChain ruleとで) (d/dφ)=(dξ/dφ).(d/dξ) と dξ/dφ=sech^2(aφ)a=a{1-tanh^2(aφ)}=a(1-ξ^2) を使って (d/dφ)=(dξ/dφ).(d/dξ)=a(1-ξ^2).(d/dξ) (d/dφ)^2=(d/dφ){a(1-ξ^2).(d/dξ)} =a(1-ξ^2).(d/dξ){a(1-ξ^2).(d/dξ)} =a^2(1-ξ^2).(d/dξ){(1-ξ^2).(d/dξ)} あるいは微分を実行して =a^2(1-ξ^2).{(-2ξ).(d/dξ)+(1-ξ^2).(d/dξ)^2} など d/dφ*(1/y)=(dξ/dφ).(d/dξ)(1/y)=a(1-ξ^2).(-1/y^2)(dy/dξ) など ---------- ξ=tanhaφを用いての変換では、ψの部分は変換されていませんので、第1式にある第2項と第4項は、変換後もそのまま変わらずに第3項と第5項にあるべきですが、 y*(d^2x/dψ^2) と (dy/dψ)*(dx/dψ) は変換後には y*(d^2x/dψ^2) と y*(dy/dψ)*(dx/dψ) になっています。変換前の第4項にyが余分に掛かって変換後の第5項になっています。 φに関係していた部分は、変換後の第1項、第2項、第4項ですが、それぞれ(1/y)の因子が足りません。もし、全体にyを掛けたとすれば、変換後の第3項にYがひとつ足りません。変換後の対応する順に並べれば以下のようになると思います。 -2(1/y)*a^2*ξ(1-ξ^2)*(dx/dξ) +(1/y)*(a(1-ξ^2))^2*(d^2x/dξ^2) +y*(d^2x/dψ^2) -(1/y)^2*(a(1-ξ^2))^2*(dy/dξ)*(x/dξ) +(dy/dψ)*(dx/dψ)=0

abc1031man
質問者

お礼

丁寧な解説ありがとうございます。 もう一度、見直してみます。

その他の回答 (1)

  • rabbit_cat
  • ベストアンサー率40% (829/2062)
回答No.1

普通に、 dx/dφ = dx/dξ・dξ/dφ d^2x/dφ^2 = d^2x/dξ^2・(dξ/dφ)^2 + dx/dξ・d^2ξ/dφ^2 といった関係を使えばいいのでは。

関連するQ&A

  • 変数分離の規則について

    変数分離の規則について 例えば、3y^2(dy/dx)=y^3+xですが、 y^3=α 3y^2=dα/dy より dα/dy・dy/dx=α+x dα/dx=α+x とαとxの式に変換して構わないのでしょうか? ご指導願います。

  • 変数変換?

    まずx^2y''+xy'-y=0という問いがあるのですがy1=xで、y=y1zとおきyが解であるときのzの満たす方程式とはどういうことなのでしょうか?またx^2y"+4xy'+2y=1,y(1)=y'(1)=0という問いで変数変換を考え(d^2y/dt^2)をy"=(d^2y/dx^2),y=(dy/dx),(dy/dt)で表し、元の方程式に代入してy(t)の満たす方程式にしたのですがそこからどう解けばいいのかわかりません。また元の方程式を解く場合(変数変換をしない場合?)と上記の作業をした場合では解き方、解等違ってくるのでしょうか?どうかご教授お願いします。

  • 変数変換と微分

    よろしくお願いします。ある本で次の図と記述があります。 ※vxのxはvの右下に付く添え字です。 Fxx=μ(dvx/dy)+μ(dvy/dx) これを以下の式でxy座標からXY座標に変換 X=(x+y)/√2,Y=(-x+y)/√2 上記の両式をt,x,yで微分すると VX=(dX/dt)=(vx+vy)/√2 VY=(dY/dt)=(-vx+vy)/√2 (∂X/∂x)=1/√2,(∂X/∂y)=1/√2,(∂Y/∂x)=-1/√2,(∂Y/∂x)=1/√2 が得られます。 さらに整理して vx=(VX-VY)/√2,vy=(VX+VY)/√2 これらを使って,一番最初の式に出てくる微分(dvx/dy)と(dvy/dx)を変数X,Yによる微分に変数変換します。 (dvx/dy)=(∂X/∂y)(∂vx/∂X)+(∂Y/∂y)(∂vx/∂Y) (dvy/dx)=(∂X/∂x)(∂vy/∂X)+(∂Y/∂x)(∂vy/∂Y) ここで質問です。  上記の式の最後の2行 つまり, (dvx/dy)=(∂X/∂y)(∂vx/∂X)+(∂Y/∂y)(∂vx/∂Y) (dvy/dx)=(∂X/∂x)(∂vy/∂X)+(∂Y/∂x)(∂vy/∂Y)  この式が出てきた経緯が分かりません。(dvxをdyで微分するとなぜ右辺のようになるのかが分かりません。同様にdvyをdxで微分するとなぜ右辺のようになるのかが分かりません。)  どなたか解説をお願いできないでしょうか。

  • 2次微分の変数変換

    dy/dx=(dy/du)(du/dx) とかけて、dy/dxからdy/duの関係に変換することは積分でよくあります。 ですが、2次微分 d^2y/dx^2 をdy/duの関係に書き換えるとどうなりますか。 たとえば、sinx=uとしますと、dy/dx=(dy/du)cosxになりますが、 d^2y/dx^2はどうでしょう。 うまく説明できていないかもしれませんが、 どなた分かる方がいらっしゃいましたら、ご教示お願いします!

  • 媒介変数の関数で(d^2y)/(dx^2)をお願いします

    媒介変数の関数で(d^2y)/(dx^2)をお願いします x=1-4t+t^2 y=1-t^2 どう計算しても 2/(t-2)^3になりますが 答えの分子が1でした dy/dx=-t/(t-2) 途中式をお願いいたします

  • 重積分の変数変換について

      u = x + y, v = x - y という変換をしたとき   du = (∂u/∂x)dx + (∂u/∂y)dy = dx + dy.   dv = (∂v/∂x)dx - (∂v/∂y)dy = dx - dy.   ∴dudv = dx^2 - dy^2 = 0 となってしまい   dudv = 2dxdy という関係式が導けません。  上記を行列で表すと   ┌ ┐ ┌        ┐┌ ┐ ┌  ┐┌ ┐   |du|=|∂u/∂x ∂u/∂y||dx|=|1 1||dx|   |dv| |∂v/∂x ∂v/∂y||dy| |1 -1||dy|   └ ┘ └        ┘└ ┘ └  ┘└ ┘.   abs|1 1|= 2     |1 -1| となって、たしかにヤコビアンは 2 になるのですが。  上の考え方はどこがおかしいのでしょう?

  • 変数分離形微分方程式の導入における式の展開について

    変数分離形微分方程式の導入において、参考書に次のような式の展開がありました。 y=f(x)において、         Δy/Δx≒dy/dx         Δy≒dy/dx・Δx         dy=dy/dx・Δx(近似式)・・・(1) ここで、y=x とすると、         dx=dx/dx・Δx より         dx=Δx・・・(2) (2)を(1)に代入すると、         dy=dy/dx・dx 上の式で、y=x としていますが、なぜなのでしょうか? y=x は y=f(x) の関数の1つに過ぎないと思うのです。一般化しているように思えないのですが・・・ 宜しくお願い致します。

  • 完全形でない3変数関数の微分方程式の解法

    全微分方程式A(x,y,z)dx+B(x,y,z)dy+C(x,y,z)dz=0がある。この式をPとおく。ここで、ベクトル値関数f=[A,B,C]とおき、f・(rotf)=0となるならばPは積分可能でその一般解は下記の手順により求まる。 手順1:Pについてdz=0とすると、Adx+Bdy=0となる。この式をQとおく。これが(∂A/∂y)=(∂B/∂x)を満たすとき、また満たさないときは積分因子μをかけることによりこのQの一般解ξ(x,y,z)=E (Eは定数)が得られる。 手順2:Pの両辺にλをかけたものの一般解を求める。するとλAdx=(∂ξ/∂x)となる。これから、λの値を求める。 手順3:ξの全微分はdξ=(∂ξ/∂x)dx+(∂ξ/∂y)dy+(∂ξ/∂z)dzとなり、このうち(∂ξ/∂x)dx+(∂ξ/∂y)dyはλAdx+λBdyとなるが、最後の(∂ξ/∂z)dzだけはλRdzとなるかは不明である。 dξ=(∂ξ/∂x)dx+(∂ξ/∂y)dy+(∂ξ/∂z)dzと(∂ξ/∂x)dx+(∂ξ/∂y)dy=λAdx+λBdyより、λAdx+λBdy=dξ-(∂ξ/∂z)dzとなる。 するとPの両辺にλをかけた式は、λAdx+λBdy+λCdz=dξ+{λC-(∂ξ/∂z)}dz=0となる。 ここで、λC-(∂ξ/∂z)=ηとおくと、λAdx+λBdy+λCdz=dξ+ηdz=0となり、2変数の全微分方程式dξ+ηdz=0が得られる。この解が結局全微分方程式Pの一般解となる。 ここで質問です。 手順1でdz=0とした式Adx+Bdy=0 (∂A/∂y)=(∂B/∂x)、またはμAdx+μBdy=0 (∂μA/∂y)=(∂μB/∂x)を解くとこの一般解、ξ(x,y,z)=Eが得られ、この関数ξの全微分はdξ=(∂ξ/∂x)dx+(∂ξ/∂y)dy=Adx+Bdy=0、またはdξ=(∂ξ/∂x)dx+(∂ξ/∂y)dy=μAdx+μBdy=0になるのが分かります。 手順2,3でλAdx+λBdy+λCdz=0という式が出てきますが、これはλをかける事により完全形になっていると思われます。しかしなぜλAdx=(∂ξ/∂x)となるのかが分かりません。ξはAdx+Bdy=0の解として現れる関数なので、λAdx+λBdy+λCdz=0を満たす関数は別にあり、例えばこれをσとすると、この関数の全微分はdσ=(∂σ/∂x)dx+(∂σ/∂y)dy+(∂σ/∂z)dz=λAdx+λBdy+λCdz=0となり、λAdx=(∂σ/∂x)dxとなるのではないのでしょうか? それともこの関数σがξと一致すると仮定しているのでしょうか? それからもう1つ気になるのですが、手順3で「最後の(∂ξ/∂z)dzだけはλRdzとなるかは不明である。」とありますが、これもよく意味が分かりません。なぜ(∂ξ/∂z)dzだけλRdzとはなるか分からないのでしょうか? おそらく私が根本的に間違っていると思いますので、詳しい方教えてください。お願いします。

  • 45°傾いた座標系への変数変換

    ある本で添付の図が示してあり,xy座標系から45°傾いたXY座標系へ変数を変換するという式が掲載されています。 その式は,X=(x+y/√2)),Y=((-x+y)/√2)と書かれています。 そこで私は,A(1,2) B(2,2) C(2,1) D(1,1)という点を仮定し上記の式を使ってA',B',C',D'への変換を試みました。 A→A’(1,2)→(√2,1) B→B’(2,2)→(2√2,0) C→C’(2,1)→(3√2/2,-√2/2) D→D’(1,1)→(√2,0)  ところが,上記のA’,B’,C’、D’を座標の上にプロットしても45°傾いた正方形になりません。 そもそも式がおかしいような気がするのですが,私がどこか勘違いをしているのでしょうか。  御回答をいただけたら助かります。

  • dx を変数として扱える理由

    高校の数学では、微分を次のように習いました。 y=f(x) ...f(x)はxの関数 yをxで微分することを次のように書く。 dy/dx=df(x)/dx 例えば y=f(x)=x^2+3x+4 なら dy/dx=2x+3 高校の授業で数学の先生の漏らした言葉に、 dx は、ひとつの変数と扱って計算してよい。 とすると dy=(2x+3)dx と書ける。 ここで積分の魔法をかけると ∫dy=∫(2x+3)dx y+A=x^2+3x+B (A,Bは定数) なんと、これはA,Bを B-A=4とすれば 最初のf(x)と一致してます。 このようなめちゃくちゃな話をそのまま信じるのも あれなので、こんな計算が許される理由を教えてください。