質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.84%

  • Python 3.x

    4843questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

pythonで偏微分の結果を描画したい

受付中

回答 1

投稿

  • 評価
  • クリップ 0
  • VIEW 102
退会済みユーザー

退会済みユーザー

 前提・実現したいこと

x0=3.0 x1=4.0 のときに、x0について偏微分した場合、
その結果のmatplotlibでの表示方法が知りたいです。

現在のファイルの記述は以下の通りです。
以下のページで引数の二乗和を出力する記述を真似ました。
https://endoyuta.com/2017/01/09/python-matplotlib%E3%81%A73%E6%AC%A1%E5%85%83%E3%82%B0%E3%83%A9%E3%83%95%E3%82%92%E6%9B%B8%E3%81%8F/

import numpy as np
import matplotlib.pylab as plt
from mpl_toolkits.mplot3d import Axes3D

def function(x0, x1):
    return x0**2 + x1**2

x0 = np.arange(-3.0, 3.0, 0.25)
x1 = np.arange(-3.0, 3.0, 0.25)
X0, X1 = np.meshgrid(x0, x1)
Y = partial_diff(X0)
fig = plt.figure()
ax = Axes3D(fig)
ax.plot_wireframe(X0, X1, Y)
plt.show()

 試したこと

def numerical_diff(f, x):
    h = 1e-4
    return (f(x+h) - f(x-h)) / (2*h)

def function(x0):
    return x0*x0 + 4.0**2.0

print(numerical_diff(partial_diff, 3.0))

上記のコードで偏微分することはできていると思うのですが、
matplotlibにどのようにその値を適用させ、出力するのかがわかっていません。

数学に対する知識も鈍いので、初歩的なミスに陥っているのかもわかりませんが、
どうぞよろしくお願いいたします。

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

+1

上記のコードで偏微分することはできていると思うのですが

失礼ながらできているようには見えません。

(1) partial_diffの定義が提示されてない
(2) 2つ目のコード内のfunction(x0)の内容x0*x0 + 4.0**2.0が間違っている
(3) 微分計算に数値計算アプローチ・解析的アプローチの2つがあるがどちらを取ろうとしているか見えない。
両方を不自然に混ぜて解こうとしているかのようにも見えます。((1)が提示されてないので実際にどうかは不明です)

解析的に解くなら以下になるでしょうし

∂f/∂x1 = 2x1

数値計算的に解くなら以下になると思います。

∂f/∂x1 = (f(x0, x1 + h) - f(x0, x1))/h


さて元の質問に戻ります。

x0=3.0 x1=4.0 のときに、x0について偏微分した場合、その結果のmatplotlibでの表示方法が知りたいです。 

特定の位置(x0=3.0 x1=4.0)の偏微分の値はスカラでしかありません。それを3Dグラフとしてプロットしても興味のあるようなものにはならない気がします。だって3D空間といってもプロットする点が1つしかないのですから・・・実際の質問者さんの意図は「fをx0, x1が-3~3の間のx1についての偏微分値をプロットすること」ではないのでしょうか?もしそうなら以下のようなグラフになると思います。

イメージ説明


と、ここまでコメントしておいてなんですが「プログラミングの質問になっていない」気がします。

質問の際には「何を求めているか」「プログラミング上の課題が何か」を明確にすべきですが、「何を求めているか」は問題領域(本件では数学の微分法)について第三者にもわかるように明確に述べる必要があります。しかしご質問ではそれが少々欠けているようです。それゆえ「プログラミング上の課題が何か」を閲覧者が推測できない気がするのです。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 90.84%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

同じタグがついた質問を見る

  • Python 3.x

    4843questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。