sympy.solve()関数の解が出てきません

Question

y = sympy.Symbol('y')

ans = sympy.solve(sympy.sin(y) - sympy.sin(y - mouseh) - hm)

という計算をしたいのですが、値が出ません。
mousehを予め、y=0のときのhmの値となるように設定していて、hmを変化させたときのyを計算したいです。
mouseh, hmを定数にしても(mouseh=0.4, hm=0.609375)変わりませんでした。

y = sympy.Symbol('y')
z = sympy.Symbol('z')

exp1 = sympy.sin(y) - sympy.sin(z) - hm #   縮むほうで計算
exp2 = y - mouseh -z
ans = sympy.solve([exp1,exp2])

としても変わりませんでした。

追記、

ansy = sympy.solve(((1 - sympy.sin(x)) / (sympy.sin(x) + sympy.cos(x))) - ratio)

という式は計算できていて、(ratioは0.6~0.7程度です）式の書き方がおかしいのかと思ったんですが、、、

Comments

[tiitoi]

mouseh、hm はなんでしょうか？
コピペして動くコードを質問欄に記載してください。

Answer 1

Sympy の問題ではなく、mouseh, hm の値によっては解が存在しないので、solve() で解が出てこないのは当然と言えます。

import matplotlib.pyplot as plt
import numpy as np

def f(x, mouseh=0.4, hm=0.609375):
    return np.sin(x) - np.sin(x - mouseh) - hm

xs = np.linspace(-3, 3, 100)
ys = f(xs)

fig, ax = plt.subplots(figsize=(8, 8))
ax.set_xlim(-3, 3)
ax.set_ylim(-3, 3)
ax.grid(True)
ax.plot(xs, ys)
plt.show()

Comments

[YukaU]

解がある前提で作っていたのと検算方法がわからなかったためとても助かりました。
解がない場合は別の処理をするように変更します。ありがとうございました。