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

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

新規登録して質問してみよう
ただいま回答率
85.69%
NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

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

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

2回答

1480閲覧

scipyのcurve_fitを使って累乗近似をしたいが、エラーが出てしまう

tukituki_84

総合スコア26

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

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

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2023/05/02 07:01

実現したいこと

Excelの累乗近似機能のように、次のx, yの数値をPython(Numpy, Scipy)を使って「y = bx^a + c」の形に累乗近似したいと考えています。

x = (38602.61004, 41670.86138, 45333.78572, 49376.79694, 53348.51068) y = (158489.3192, 125892.5412, 100000, 79432.82347, 63095.73445)

発生している問題・エラーメッセージ

コードを作成し実行しましたが、次のようなエラーが出てしまいます。

Traceback (most recent call last): File "/Users/202305/0502/0502b.py", line 10, in <module> popt, cov = curve_fit(funk_fit, x, y) File "/Users/opt/anaconda3/lib/python3.9/site-packages/scipy/optimize/minpack.py", line 794, in curve_fit raise RuntimeError("Optimal parameters not found: " + errmsg) RuntimeError: Optimal parameters not found: Number of calls to function has reached maxfev = 800.

エラーの意味が分からず、途方に暮れている状況です。
訳すると「最適なパラメータが見つからない」と書いているので、そもそも大事なところで間違っているということでしょうか。。

該当のソースコード

Python

1import numpy as np 2from scipy.optimize import curve_fit 3 4def funk_fit(x, a, b, c): 5 return b*x**a + c 6 7x = np.array([38602.61004, 41670.86138, 45333.78572, 49376.79694, 53348.51068]) 8y = np.array([158489.3192, 125892.5412, 100000, 79432.82347, 63095.73445]) 9 10popt, cov = curve_fit(funk_fit, x, y) 11popt

どの部分にこのエラーの原因があるのでしょうか。どなたか助けていただけると幸いです。

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

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

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

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

  • 質問になっていない投稿
  • スパムや攻撃的な表現を用いた投稿

適切な質問に修正を依頼しましょう。

また依頼した内容が修正された場合は、修正依頼を取り消すようにしましょう。

melian

2023/05/02 07:36

a を整数に限定してみてはどうでしょうか。 return b*x**(int(a)) + c この場合、結果は以下になります。 [ 1.00000000e+00 -6.31642507e+00 3.93831191e+05]

回答2

1

ベストアンサー

ChatGPTに聞いてみたところ、関数評価回数maxfevを増やしてみるとよいとの回答でした。
以下、それに従ったコードです。

Python

1import numpy as np 2from scipy.optimize import curve_fit 3import matplotlib.pyplot as plt 4 5def funk_fit(x, a, b, c): 6 return b*(x**a) + c 7 8x = np.array([38602.61004, 41670.86138, 45333.78572, 49376.79694, 53348.51068]) 9y = np.array([158489.3192, 125892.5412, 100000, 79432.82347, 63095.73445]) 10 11popt, cov = curve_fit(funk_fit, x, y, maxfev=20000) 12print(popt) # [ 1.13142281e-03 -2.53388891e+08 2.56587018e+08]

投稿2023/05/02 07:36

can110

総合スコア37977

hotjyu65_👍を押しています

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

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

このような回答には修正を依頼しましょう。

また依頼した内容が修正された場合は、修正依頼を取り消すようにしましょう。

1

今回の件によらず、Pythonで起きるエラーの大半は、エラーメッセージをそのままコピペするといい感じの先行事例行きつきやすいです。


今回のキーワードはすでに察しておられると思いますが、RuntimeError: Optimal parameters not found: Number of calls to function has reached maxfev = 800.になります。これをキーワードに探すと、StackOverflowに行きつきます。

Excelの累乗近似機能のように、次のx, yの数値をPython(Numpy, Scipy)を使って「y = bx^a + c」の形に累乗近似したいと考えています。

エクセルやPythonを抜きにした大枠の話になりますが、最適化問題を扱う際のアプローチとして、一般的に「何回も試して一番いいもの(落としどころ)を答える、無ければエラーを返す」が一般的に取られます。今回の場合も恐らくこれで、StackOverflowの回答にもあるように、(試行回数を増やす、制約を増やす等の)答えに行きつきやすくなるような策をして様子を見るが一番の解決策の様です。具体的には、curve_fit(略, maxfev=8000) のように指定できるはずです。


※試行数を大きくし過ぎると、計算できない場合にずっと待たされるのでご注意ください。もし既定の試行回数が一瞬で済むようであれば、手動で決める試行回数をx10やx100にしても問題ないとは思います。

投稿2023/05/02 07:28

編集2023/05/02 07:30
fourteenlength

総合スコア3029

hotjyu65_👍を押しています

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

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

このような回答には修正を依頼しましょう。

また依頼した内容が修正された場合は、修正依頼を取り消すようにしましょう。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.69%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

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

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

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

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。