私はPython初心者なのです。複数のファイルを読み込み(今回はdata_0001とdata_x_7.6_0001というものを0001~0300まで)を行いたいと思うのですが、複数ファイルといっても数字の先頭に0という数字が残ってしまうためこのプログラムでは作動しません。
さらに、渦度を求めたく、このdata_0001とdata_x_7.6_0001というファイルには、x,y,z,u,v,w,pという順番でデータが入っており欲しいデータはx,y,u,vの4つで渦度の式は((vの変位)/(xの変位))ー((uの変位)/(yの変位))です。
申し訳ないのですが、どのようなプログラムが最適か教えてくださると幸いです。
import numpy as np
list_a= [None] * 300
for i in range(300):
with open('file/data_{}.txt'.format(i + 1), mode='r', encoding='utf-8') as f:
list_a[i] = f.read()
list_b= [None] * 300
for i in range(300):
with open('file/data_x_7.6_{}.txt'.format(i + 1),mode='r', encoding='utf-8') as f:
list_b[i] = f.read()
array_a=np.array(list_a)
array_b=np.array(list_b)
array_c=array_b-array_a
array_d=np.array(array_c,usecols=(1,2,3))
array_e=np.array(array_c,usecols=(0,2,3))
array_f=np.array(array_c,usecols=(0,1,3))
array_g=np.array(array_c,usecols=(0,1,2))
array_h=array_g/array_d
array_i=array_f/array_e
array_j=array_h-array_i
result = "{:.4f}\n".format(array_j)
print(array_j)
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。