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

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

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

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

解決済

測定データを積分関数でフィッティングしたいです。

saka_kei
saka_kei

総合スコア0

Python

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

2回答

0評価

1クリップ

4972閲覧

投稿2017/08/01 15:35

編集2022/01/12 10:58

###前提・実現したいこと
scipy.optimizeのlastsqを用いたフィッティングについて質問です。
フィッティングする関数に、積分関数を使いたいと思ってます。
残渣を求める行にintegrate.quad関数を用いて、その引数にパラメーターを指定してフィッティングをしようと思いましたが、下記のエラーが出てつまづいています。

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

--------------------------------------------------------------------------- ValueError Traceback \(most recent call last\) <ipython-input-32-02de193ee868> in <module>\(\) 25 26 param0 = \[1\.2, 1, 10\.5, -24\.3\] ---> 27 param = optimize\.leastsq\(fanction_gauss, param0, args=\(x, y\)\) ~/\.pyenv/versions/anaconda3-2\.4\.1/lib/python3\.5/site-packages/scipy/optimize/minpack\.py in leastsq\(func, x0, args, Dfun, full_output, col_deriv, ftol, xtol, gtol, maxfev, epsfcn, factor, diag\) 375 if not isinstance\(args, tuple\): 376 args = \(args,\) --> 377 shape, dtype = _check_func\('leastsq', 'func', func, x0, args, n\) 378 m = shape\[0\] 379 if n > m: ~/\.pyenv/versions/anaconda3-2\.4\.1/lib/python3\.5/site-packages/scipy/optimize/minpack\.py in _check_func\(checker, argname, thefunc, x0, args, numinputs, output_shape\) 24 def _check_func\(checker, argname, thefunc, x0, args, numinputs, 25 output_shape=None\): ---> 26 res = atleast_1d\(thefunc\(\*\(\(x0\[:numinputs\],\) \+ args\)\)\) 27 if \(output_shape is not None\) and \(shape\(res\) != output_shape\): 28 if \(output_shape\[0\] != 1\): <ipython-input-32-02de193ee868> in fanction_gauss\(param, x, y\) 21 22 def fanction_gauss\(param, x, y\): ---> 23 residual = y - integrate\.quad\(\(param\[0\]\*np\.exp\(-\(x-param\[2\]\)\*\*2/\(2\*param\[1\]\)\)\+param\[3\]\), x-0\.25, x\+0\.25\) 24 return residual 25 ~/\.pyenv/versions/anaconda3-2\.4\.1/lib/python3\.5/site-packages/scipy/integrate/quadpack\.py in quad\(func, a, b, args, full_output, epsabs, epsrel, limit, points, weight, wvar, wopts, maxp1, limlst\) 321 if \(weight is None\): 322 retval = _quad\(func, a, b, args, full_output, epsabs, epsrel, limit, --> 323 points\) 324 else: 325 retval = _quad_weight\(func, a, b, args, full_output, epsabs, epsrel, ~/\.pyenv/versions/anaconda3-2\.4\.1/lib/python3\.5/site-packages/scipy/integrate/quadpack\.py in _quad\(func, a, b, args, full_output, epsabs, epsrel, limit, points\) 370 def _quad\(func,a,b,args,full_output,epsabs,epsrel,limit,points\): 371 infbounds = 0 --> 372 if \(b != Inf and a != -Inf\): 373 pass # standard integration 374 elif \(b == Inf and a != -Inf\): ValueError: The truth value of an array with more than one element is ambiguous\. Use a\.any\(\) or a\.all\(\)

###該当のソースコード

python

import sys, csv import numpy as np from scipy import integrate sys\.path\.append\(""\) import matplotlib\.pyplot as plt cnt = 0 x = np\.array\(\[\]\) y = np\.array\(\[\]\) da_int = np\.array\(\[\]\) f = open\("beamsize_y_up_45mm\.csv", "r", encoding="utf_8_sig"\) try: reader = csv\.reader\(f\) for row in reader: x = np\.append\(x, float\(row\[0\]\)\) y = np\.append\(y, float\(row\[1\]\)\) cnt = cnt \+ 1 finally: f\.close\(\) def fanction_gauss\(param, x, y\): residual = y - integrate\.quad\(\(param\[0\]\*np\.exp\(-\(x-param\[2\]\)\*\*2/\(2\*param\[1\]\)\)\+param\[3\]\), x-0\.25, x\+0\.25\) return residual param0 = \[1\.2, 1, 10\.5, -24\.3\] param = optimize\.leastsq\(fanction_gauss, param0, args=\(x, y\)\)

###試したこと
フィッティング関数を単なるgaussianにするとうまくいきました。

python

def gaussian\(param, x, y\): residual = y - \(\(\(param\[0\]/\(np\.sqrt\(2\*np\.pi\)\*param\[1\]\)\)\*np\.exp\(-\(x-param\[2\]\)\*\*2/\(2\*param\[1\]\*\*2\)\)\)\+param\[3\]\) return residual

###補足情報(言語/FW/ツール等のバージョンなど)
Python 3.5.3 :: Anaconda 2.4.1 (x86_64)を使っています。
python初心者です。よろしくお願いいたします。

良い質問の評価を上げる

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

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

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

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

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

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

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

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

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

まだ回答がついていません

会員登録して回答してみよう

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

ただいまの回答率
87.20%

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

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

質問する

関連した質問

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

Python

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