• 締切済み
  • すぐに回答を!

2変数の2次導関数の離散化(差分)

関数f(x,y)の2次導関数である∂f(x,y)/∂x∂yを数値計算により, より高精度で計算したいと考えているのですが,∂f(x,y)/∂x∂x のように同じ変数に対する偏導関数の離散化方法しか見つけることが できませんでした. そこで,ご存じの方がいらっしゃいましたら,∂f(x,y)/∂x∂yのような2変数の導関数の高次計算法をご教示頂けないでしょうか. 宜しくお願い致します. またそれらについて詳細に記述された文献等もあれば,教えて頂ければ幸いです.

共感・応援の気持ちを伝えよう!

  • 回答数2
  • 閲覧数1372
  • ありがとう数1

みんなの回答

  • 回答No.2

No.1です。私の経験では、境界条件が高次差分に対応している場合が少なく、適当に処理すると後で大きな誤差となってしまうため、使わないことにしています。それよりもメッシュを細かく切る方が良いと思います。4次差分にした場合、境界条件は2階微分を与える必要があると思いますが、そちらの検討結果も勘案してください。

共感・感謝の気持ちを伝えよう!

質問者からのお礼

ご回答ありがとうございます.高次にすると境界におけるゴーストポイントが増えるので境界条件が複雑になり,高次差分による高精度化よりも格子間隔自体を小さくしたほうがよいということですね.参考にさせて頂きます.

質問者からの補足

たびたび申し訳ありません. 上記の内容についてですが,計算メモリで制限を受けている場合は手法の改善というよりは並列化というのが,現在の数値計算のセオリーになっているのでしょうか.

関連するQ&A

  • 回答No.1

2階偏微分方程式を差分法を用いて離散化する時は境界内をx方向にh、y方向にkのメッシュを掛けて、その交点(i,j)における値f(i,j)と隣接する点との関係に書きなおす必要があります。ご質問の,∂f(x,y)/∂x∂yは   ∂f(x,y)/∂x∂y=(f(i+1,j+1)-f(i+1,j-1)-f(i-1,j+1)+f(i-1,j-1))/(4hk) で与えられます。さらに隣までの関係を取り入れれば差分式としての精度は上がりますが、私の経験では面倒なだけで上式が扱いやすく、メッシュの切り方が適切ならば非線形の場合も収束もが早かったように思います。 境界条件も差分化して適切に処理すれば数値解としての正解に辿りつくはずです。

共感・感謝の気持ちを伝えよう!

質問者からの補足

ご回答ありがとうございます. 上記の方法の高精度化についてお尋ねしたいのですが, 上記の方法は2次精度の中心差分法をx方向偏導関数に対してy方向に用いることで算出される式だと思うのですが,精度が現実の現象に対して大きく問題となる場合には,同様に4次精度の中心差分法をx方向偏導関数に対してy方向に用いればよいのでしょうか.(あるいはその逆) もし宜しければご回答を宜しくお願い致します. ∂f(x,y)/∂x∂y ≒ (fx(i,j+1)-fx(i,j-1))/(2h) ≒ (f(i+1,j+1)-f(i+1,j-1)-f(i-1,j+1)+f(i-1,j-1))/(4hk)

関連するQ&A

  • 第5次導関数の問題です

    この解き方であっているか、わかる方よろしくお願いします。 関数f(x)=x^5+2mp第5次導関数f^(5)(x)を求めよ。 (f^(5)(x)の(5)の部分だけが指数です。) 答え 1次導関数:5x^4 2次導関数:20x^3 3次導関数:60x^2 4次導関数:120x 5次導関数:120 よって、f^(5)(x)=120

  • 高次導関数について

    log(1-x^2)の高次導関数を求めよという問題なんですが どうしてもわかりません 教えていただきたいです。

  • 2次の導関数についての質問

    導関数と2次関数について、下記の考え方であっていますか? もし間違えていたらご指摘おねがいします!

  • 2変数関数の2次導関数のことです。

    2回連続微分可能で、z=f(x,y),x=x(t),y=y(t)の関係があって、このときのzのtに関する2次導関数を求めるという問題なんですが、1次の導関数は dz/dt=(∂z/∂x)(dx/dt)+(∂z/∂y)(dy/dt) だと思うんですが、2次の場合は d^2z/dt^2=(d/dt)((∂z/∂x)(dx/dt)+(∂z/∂y)(dy/dt)) となって、それぞれの項を積の微分法で解けばいいのでしょうか?できたらその形も教えて下さい。お願いします。

  • 公式より導関数を求める

    lim h→0 f(a+h)-f(a)/h の公式より導関数を求めたいと思いますが 計算手順がわからないので、教えてください。宜しくお願いします。 普通に微分したほうが早いのですけど、式を定義にして解こうとすると分かりません。宜しくお願いします。 【問題】 y=1/ x^2 の導関数を求めよ。 

  • 2次関数の導関数についての質問

    下の2次関数の導関数についての設問の正誤について教えてください。もし間違っているなら解説していただけないでしょうか?

  • 導関数の問題です

    解析学を習い始めたばかりの初心者です。 以下の問題の答え(できれば途中の計算も)を教えていただけないでしょうか? 次の導関数を求めよ。 (1) f(x)=(2x+1)^3 (2) g(x)=1/(x^2+x+1) わかる方教えてください。m(_ _)m

  • 指数が含まれる第5次導関数の問題の解き方

    以下の問題を自分で解答してみましたが、 ぜんぜん自信がありません。 わかる方、いらっしゃいましたらご指導お願いします。 【問題】 関数f(x)=xe^(5x)の第5次(階)導関数f^(5)(x)を求めよ。 【解答】 1次導関数:5e^(4x) 2次導関数:5e^(3x) 3次導関数:5e^(2x) 4次導関数:5e^(1x) 5次導関数:5e^(0)=5・1 よって、f^(5)(x)=5 以上、よろしくお願いします。

  • 導関数の計算

    f(x)=x-1(マイナス一乗の意)の導関数を定義に従って求める問題なんですが、どうしてよいかわかりません(計算の仕方)。 どなたか教えて下さい。

  • 導関数を定義に従って求める問題

    数学の得意な方にとっては基本的な問題なのかもしれませんが、自分で問題を解きながら、以下どのようにしたらよいかわからなくなってしまいました。 問題は -5x^2+4xの導関数を定義に従って求めるというものです。 定義に従って、問題を解いていくと、今までの導関数の問題では分母のhが綺麗に消えて計算が出来たのですが、今回はhが残ってしまいます。 limit[(-5x^2-10xh)/h,h->0] + Limit[4x/h,h->0] までは出来ました。ここからやり方がわかりません。 いろいろ教科書や参考書を見てみたのですが、 どうも、自分が悩んでいる部分の箇所が見つからなくて…。 どうか、宜しくお願いいたします。

専門家に質問してみよう