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

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

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

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

Q&A

1回答

1382閲覧

微分方程式を解き,Excelに出力したいです.

nonnon615

総合スコア19

Python

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

0グッド

1クリップ

投稿2021/07/27 14:13

import sympy as sp import math import cv2 import numpy as np import matplotlib.pyplot as plt import csv import os import glob import pandas as pd in_excel_file = 'matome' file = in_excel_file pd.options.display.precision = 10 options = {} options['strings_to_formulas'] = False options['strings_to_urls'] = False with pd.ExcelWriter('yre_z_DH=0.4_.xlsx', mode='w') as writer: results = [] for i in range(35): print(i) def loadyre(file,i): yre = pd.read_excel(file+'.xlsx', skiprows=4, usecols=[i+3], nrows=5, header=None) #解析対象となるデータの読み込み return yre y0 = loadyre(file,i) print(y0) D = 0.016 H = 0.04 d = 0.015 h = 0.039 U = i*0.25 + 0.25 print(U) π = math.pi m = 0.0156 L = 0.3 fc = 9.0 ρ = 1.18 M = (m/(ρ*H*D*L))*(D/H) Vr = U/(fc*H) A = H*L sp.init_printing() sp.var('z') y = sp.Function('y')(z) #yはzの関数(Function) E = 4.28*10*10*10*10*10*10 I = (1/12)*(D*H*H*H-d*h*h*h) dCL = 0.3 #dCL/dα FL = -((ρ*U*A)/2*M)*dCL M = (FL*y)*(y/2) eq = sp.Eq( sp.diff(y, z ,2)+(M/(E*I))*sp.diff(y, z, 1), 0 ) yrez = sp.dsolve(eq, ics={y.subs(z,0):y0, y.diff(z,1).subs(z,0):0}) results.append(yrez.subs(z, 190)) df = pd.DataFrame(results).T df.to_excel(writer)

上記を実行して,運動方程式を解きたいのですが,どうもうまく解けていないように感じます.

Excelからデータを取り出すことはできています.

以下にエラーコードを示します.
よろしくお願いいたします。

0 3 0 0.0000000058 0.25 Traceback (most recent call last): File "C:\Users\片持ち.py", line 53, in <module> yrez = sp.dsolve(eq, ics={y.subs(z,0):y0, y.diff(z,1).subs(z,0):0}) File "C:\Users\Python\Python37\site-packages\sympy\solvers\ode\ode.py", line 631, in dsolve x0=x0, n=n, **kwargs) File "C:\Users\Python\Python37\site-packages\sympy\solvers\deutils.py", line 210, in _desolve n=terms, x0=x0, prep=prep) File "C:\Users\Python\Python37\site-packages\sympy\solvers\ode\ode.py", line 963, in classify_ode ics = sympify(ics) File "C:\Users\Python\Python37\site-packages\sympy\core\sympify.py", line 361, in sympify return conv(a) File "C:\Users\Python\Python37\site-packages\sympy\core\containers.py", line 310, in <lambda> converter[dict] = lambda d: Dict(*d.items()) File "C:\Python\Python37\site-packages\sympy\core\containers.py", line 234, in __new__ elements = frozenset(items) File "C:\Users\Python\Python37\site-packages\sympy\core\containers.py", line 106, in __hash__ return hash(self.args) TypeError: unhashable type: 'DataFrame' During handling of the above exception, another exception occurred: Traceback (most recent call last): File "C:\Users\片持ち.py", line 60, in <module> df.to_excel(writer) File "C:\Users\Python\Python37\site-packages\pandas\io\excel\_base.py", line 1020, in __exit__ self.close() File "C:\Users\Python\Python37\site-packages\pandas\io\excel\_base.py", line 1024, in close content = self.save() File "C:\Users\Python\Python37\site-packages\pandas\io\excel\_openpyxl.py", line 80, in save self.book.save(self.handles.handle) File "C:\Users\Python\Python37\site-packages\openpyxl\workbook\workbook.py", line 407, in save save_workbook(self, filename) File "C:\Users\Python\Python37\site-packages\openpyxl\writer\excel.py", line 293, in save_workbook writer.save() File "C:\Users\Python\Python37\site-packages\openpyxl\writer\excel.py", line 275, in save self.write_data() File "C:\Users\Python\Python37\site-packages\openpyxl\writer\excel.py", line 89, in write_data archive.writestr(ARC_WORKBOOK, writer.write()) File "C:\Users\Python\Python37\site-packages\openpyxl\workbook\_writer.py", line 148, in write self.write_views() File "C:\Users\Python\Python37\site-packages\openpyxl\workbook\_writer.py", line 135, in write_views active = get_active_sheet(self.wb) File "C:\Users\Python\Python37\site-packages\openpyxl\workbook\_writer.py", line 33, in get_active_sheet raise IndexError("At least one sheet must be visible") IndexError: At least one sheet must be visible

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

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

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

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

ppaul

2021/07/27 15:30

微分方程式を数式処理で解きたいのですか、それとも数値処理で解きたいのですか?
nonnon615

2021/07/27 15:49

どちらでも構いません。 微分方程式を数式処理したことがあったので、同様に組んでみたのですが上手くいきませんでした。 何らかの手段で数値解析が出来ればと考えています。
guest

回答1

0

微分方程式を数式処理で解きたいのであればsympyをマスターしましょう。
そのためにはWelcome to SymPy’s documentation!を読みましょう。
それを読めば、sp.dsolveの引数にDataFrameを渡してはいけない理由もわかるはずです。

微分方程式を数式処理で解きたいのであれば、運動方程式なら常微分方程式なのでscipyのodeintをマスターしましょう。
そのためにはscipy.integrate.odeintを読みましょう。

投稿2021/07/27 16:13

ppaul

総合スコア24670

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

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

アカウントをお持ちの方は

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

ただいまの回答率
85.35%

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

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

質問する

関連した質問