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

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

ただいまの
回答率

88.59%

python 微分でxを定義する方法

解決済

回答 1

投稿 編集

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

退会済みユーザー

pythonにて微分の問題を解くプログラムを組んでいるのですがxの定義の仕方が分かりません
下のコードの一番最後です。

class Expression:
    def eval(self, x_value):
        '''変数 x の値が x_value だった時の式の値を返すメソッド'''
        raise NotImplementedError

    def diff(self):
        '''式を変数 x で微分した式に対応する Expression オブジェクト
        を返すメソッド'''
        raise NotImplementedError

class Number(Expression):
    '''変更の必要なし'''
    def __init__(self, number):
        self.number = number

    def eval(self, x_value):
        return self.number

    def diff(self):
        return Number(0)

class BinalyExpression:
    def __init__(self, a, b):
        self.a=a
        self.b=b
class Add(BinalyExpression):
    def eval(self, x_value):
        return self.a.eval() + self.b.eval()
    def diff(self):
        return self.a.diff() + self.b.diff()
    pass

class Sub(BinalyExpression):
    def eval(self, x_value):
        return self.a.eval() - self.b.eval()
    def diff(self):
        return self.a.diff() - self.b.diff()
    pass


class Mul(BinalyExpression):
    def eval(self, x_value):
        return self.a.eval() * self.b.eval()
    def diff(self):
        return self.a.diff() * self.b.diff()
    pass

class Div(BinalyExpression):
    def eval(self, x_value):
        return self.a.eval() / self.b.eval()
    def diff(self):
        return self.a.diff() / self.b.diff()
    pass


class X(Expression):

    pass


このテストファイルを実行できるようにしたいです
↓テストファイル↓

import unittest
from expr import Expression, Number, X, Add, Sub, Mul, Div

class TestExpr(unittest.TestCase):
    def test_x_is_expression(self):
        self.assertIsInstance(X(), Expression)

    def test_x_eval_10(self):
        self.assertEqual(X().eval(10), 10)

    def test_x_diff_is_a_number(self):
        self.assertIsInstance(X().diff(), Number)

    def test_x_diff(self):
        self.assertEqual(X().diff().eval(0), 1)

    def test_2x_3_eval(self):
        f = Add(Mul(Number(2), X()), Number(3))  # f(x) = 2x + 3
        self.assertEqual(f.eval(10), 23)

    def test_2x_3_diff(self):
        f = Add(Mul(Number(2), X()), Number(3))  # f(x) = 2x + 3
        self.assertEqual(f.diff().eval(100), 2)

    def test_xx_3x_2_eval(self):
        x = X()
        f = Add(Add(Mul(x, x), Mul(Number(3), x)), Number(2))
        # f(x) = x**2 + 3x + 2
        self.assertEqual(f.eval(10), 132)

    def test_xx_3x_2_diff(self):
        x = X()
        f = Add(Add(Mul(x, x), Mul(Number(3), x)), Number(2))
        # f(x) = x**2 + 3x + 2
        self.assertEqual(f.diff().eval(10), 23)

    def test_sub(self):
        x = X()
        f = Sub(Mul(x, x), Mul(Number(2), x))  # f(x) = x**2 - 2x
        self.assertEqual(f.eval(2), 0)

    def test_sub_diff(self):
        x = X()
        f = Sub(Mul(x, x), Mul(Number(2), x))  # f(x) = x**2 - 2x
        self.assertEqual(f.diff().eval(2), 2)

    def test_div(self):
        x = X()
        f = Div(x, Add(Mul(x, x), Number(1)))  # f(x) = x / (x**2 + 1)
        self.assertAlmostEqual(f.eval(2), 0.4)

    def test_div_diff(self):
        x = X()
        f = Div(x, Add(Mul(x, x), Number(1)))  # f(x) = x / (x**2 + 1)
        self.assertAlmostEqual(f.diff().eval(1), 0)
        self.assertAlmostEqual(f.diff().eval(2), -3/25)

if __name__ == '__main__':
    unittest.main(verbosity=2)
  • 気になる質問をクリップする

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

回答 1

checkベストアンサー

+1

根本的なところからコメントを試みます。

計算機における微分には概ね2つのアプローチがあります。

  • (1) 数値微分
    任意のf(x)に対してxがある値X0のとき、f(X0)とf(X0+delta)を求めその傾きの近似値
    (f(X0 + delta) - f(X0))/delta
    により数値計算的に求める手法です。例えばf(x)=x**2のとき、x=1における微分をdelta=1E3ぐらいにとって計算すると2.0009999999996975ぐらいの値が求まります。正確な値には遠いですがf(x)がどのような複雑な式であってもこういう具体値を簡単に求められるのが利点です。

  • (2) 代数的な微分
    f(x) = x**2
    のような代数式を代数式のまま微分する手法です。ご存じのとおり
    f(x)/dx = 2 * x
    になりますが、これと同様のことを計算で求めるのですね。得られた結果 2 * x という式そのものに任意のxを与えると具体的な微分値がわかります。例えばx=1を代入すれば結果は2となり、(1)の数値微分と比べて「代数的な微分式が求まりさえすれば」より正確な値を求めることができます。微分値の具体的な値を求めるのではなく代数的な微分結果を求めること自体が目的のこともあるでしょう。

さて、ご質問の課題(?)は(2)の方にあたります。Expressionのdiffの定義をよくごらんください。

'''式を変数 x で微分した式に対応する Expression オブジェクト を返すメソッド'''

質問者さんはこのメソッドの仕様を表す文章が「代数式を微分してそれを代数式として返すのがExpressionの機能(目的)」であることを充分把握できていなかったようです。

Xの定義の前に四則演算(Add, Sub, Mul, Div)のdiffの定義をどうすべきかもう一度よく考えなおしてください。例えばAddですが

class Add(BinalyExpression):
    ...

    def diff(self):
        return self.a.diff() + self.b.diff()


これでは目的を果たしていないことがおわかりでしょうか?この定義ではdiffの結果がPythonの数値であると仮定してしまってます。そうではなく元の式を微分した代数式を表すようなExpressionを返さなければなりません。つまり本来diffメソッドは次のようになっているべきです。

class Add(BinalyExpression):
    ...

    def diff(self):
        return Add(self.a.diff(), self.b.diff())


また上記は

f(x) = f1(x) + f2(x)

のときの微分の公式

f'(x) = f1'(x) + f2'(x)

に基づいているという点に注意してください。質問者さんのMul, Divのdiffメソッドの定義は微分の定義そのものからかけ離れています。例えば乗算の微分は

f(x) = f1(x) * f2(x)

のときの公式

f'(x) = f1(x) * f2'(x) + f1'(x) * f2(x)

に従ってません。

以上の点に注意し代数的な微分公式を適用して四則演算の全てを考えなおせば、自ずとXをどう定義すべきかがおわかりになると思います。

投稿

編集

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

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

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

関連した質問

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