Loading [MathJax]/jax/output/CommonHTML/jax.js

2017年5月15日月曜日

学習環境

解析入門〈1〉(松坂 和夫(著)、岩波書店)の第5章(各種の初等関数)、5.4(三角関数(続き)、逆三角関数)、問題3.を取り組んでみる。

n

  1. eax(Acosbx+Bsinbx)aeax(Acosbx+Bsinbx)+eax(Absinbx+Bbcosbx)=eax(aAcosbx+aBsinbxbAsinbx+bBcosbx)=eax((aA+bB)cosbx+(aBbA)xsinbx)(aA+bB)cosbx+(aBbA)xsinbx=cosbxaA+bB=1aBbA=0B=1aAbaa2AbbA=0aa2Ab2A=0A=aa2+b2B=1a2a2+b2b=ba2+b2

コード(Emacs)

Python 3

#!/usr/bin/env python3
# -*- coding: utf-8 -*-

from sympy import Symbol, symbols, Derivative, sin, cos, solve, exp, pprint

a, b = symbols('a b', nonzero=True, real=True)
A, B, x = symbols('A B x', real=True)

expr1 = Derivative(A * exp(a * x) * cos(b * x) +
                   B * exp(a * x) * sin(b * x), x)

expr2 = exp(a * x) * cos(b * x)

pprint(expr1)
pprint(expr2)
pprint(expr1 - expr2)
pprint(solve((expr1 - expr2), A, B, dict=True))
expr11 = expr1.doit()
pprint(expr11)
pprint(expr11 - expr2)
pprint(solve((expr11 - expr2), A, B, dict=True))

入出力結果(Terminal, IPython)

$ ./sample3.py
∂ ⎛   a⋅x               a⋅x         ⎞
──⎝A⋅ℯ   ⋅cos(b⋅x) + B⋅ℯ   ⋅sin(b⋅x)⎠
∂x                                   
 a⋅x         
ℯ   ⋅cos(b⋅x)
   a⋅x            ∂ ⎛   a⋅x               a⋅x         ⎞
- ℯ   ⋅cos(b⋅x) + ──⎝A⋅ℯ   ⋅cos(b⋅x) + B⋅ℯ   ⋅sin(b⋅x)⎠
                  ∂x                                   
[]
     a⋅x                 a⋅x                 a⋅x                 a⋅x         
A⋅a⋅ℯ   ⋅cos(b⋅x) - A⋅b⋅ℯ   ⋅sin(b⋅x) + B⋅a⋅ℯ   ⋅sin(b⋅x) + B⋅b⋅ℯ   ⋅cos(b⋅x)
     a⋅x                 a⋅x                 a⋅x                 a⋅x           
A⋅a⋅ℯ   ⋅cos(b⋅x) - A⋅b⋅ℯ   ⋅sin(b⋅x) + B⋅a⋅ℯ   ⋅sin(b⋅x) + B⋅b⋅ℯ   ⋅cos(b⋅x) -

  a⋅x         
 ℯ   ⋅cos(b⋅x)
⎡⎧   -B⋅a⋅sin(b⋅x) - B⋅b⋅cos(b⋅x) + cos(b⋅x)⎫⎤
⎢⎨A: ───────────────────────────────────────⎬⎥
⎣⎩           a⋅cos(b⋅x) - b⋅sin(b⋅x)        ⎭⎦
$

0 コメント:

コメントを投稿