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

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

ただいまの
回答率

90.50%

  • Python

    8002questions

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

Scipyのodeintでエラーが出てしまいます。

受付中

回答 0

投稿

  • 評価
  • クリップ 0
  • VIEW 108

 前提・実現したいこと

Scipyのodeintを使って二次元ハミルトン系を解析したい

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

RuntimeWarning: invalid value encountered in double_scalars
ODEintWarning: Excess accuracy requested (tolerances too small). Run with full_output = 1 to get quantitative information.
  warnings.warn(warning_msg, ODEintWarning)

 該当のソースコード

m = 8.0
D = 10.0
z = 2.3
l = 1.95
E = 0.10

def model(f,t):

    p1 =f[0]
    q1 =f[1]
    q2 =f[2]

    dp1dt = -8*(2*(q1**3)-q1)*np.exp(-z*l*q2)
    dq1dt = (1/m)*p1
    dq2dt = ((2/m)*(E-(1/2*m)*(p1**2)-4*((q1**4)-(q1**2))*np.exp(-z*l*q2)+D*(1-np.exp(-l*q2))**2))**(1/2)
    dfdt = [dp1dt,dq1dt,dq2dt]

    return dfdt

#time
n1 = 101
t = np.linspace(0,2,n1)
p1 = np.empty_like(t)
q1 = np.empty_like(t)
q2 = np.empty_like(t)

#initial condition
n2 = 10
ini = np.linspace(-1.0,1.0,n2)

for h in range(0,n2):
    p1_0 = ini[h]

    for i in range(0,n2):
        q1_0 = ini[i]

        for j in range(0,n2):
            q2_0 = ini[j]

            f0 =[p1_0,q1_0,q2_0]
            p1[0] = f0[0]
            q1[0] = f0[1]
            q2[0] = f0[2]

            for k in range(1,n1):
                tspan = [t[k-1],t[k]]
                f = odeint(model,f0,tspan)
                p1[k] = f[1][0]
                q1[k] = f[1][1]
                q2[k] = f[1][2]
                f0 = f[1]

 試したこと

Python 初心者のためどうしたらいいのか全く分かりません

 補足情報(FW/ツールのバージョンなど)

AnacondaによりPython3系を使っています。

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

質問への追記・修正の依頼

  • RyoKaji

    2018/07/25 23:52

    エラーではなくWarningです。出力は意図したものにならないかもしれませんが、動作はします。現状どのような問題が生じているかご提示ください。

    キャンセル

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

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

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

関連した質問

  • 解決済

    ソート結果が正しくない

    前提・実現したいこと csvファイルdata1.csvおよびそれと同じフォーマットのcsvファイルdata2.csv, csvファイルdata3.csvのそれぞれに対して,温度(

  • 解決済

    Pythonにて

    例えばhttp://sample.com/a/b/c/d/e/1?ima=0000&cd=member  にブログの1ページ目があるとします。  2ページ目は  http:/

  • 受付中

    ターミナルで実行するのに時間がかかりすぎる

    ターミナルで実行するのに時間がかかりすぎます。 画像圧縮のアルゴリズムを書いています。 N × N ピクセルのグレースケール画像があり各ピクセルの画素値は 0 から 255

  • 解決済

    配列の中身を文字列として取得させたい

    配列の中身を文字列として取得させたいです。 image配列に [2, 5, 5, 2] という要素が入っています。 これを、 2 5 5 2 というように2行2列で

  • 解決済

    cythonの型定義

    cythonファイルに型を定義したいです もともとのpyxファイルと 自分で型を定義したときのコードをのせます gamma1.pyx from math import exp

  • 解決済

    nanの入れ替えをしたい

    python3で、配列内に計算上発生したnanを0で入れ替えるプログラムを書きたいと思っているのですが、入れ替えができません。 for i in range(0,10)

  • 解決済

    Pythonでバンドパスフィルタ

    Pythonで信号をバンドパスフィルタにかけたいと考えています。 fs = 250 # サンプリング周波数 nyq = fs / 2 # ナイキスト周波数

  • 解決済

    Pythonのtkinterのエラー

    2017-12-19 15:48:07.790 python[51777:1128073] *** Terminating app due to uncaught exceptio

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

  • Python

    8002questions

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