物理実験データのPythonコード簡略化処理

  物理実験データのコード処理

  ****

  物理実験(ヤング率の測定)データ処理コード

  コードが間違っている可能性があります。書き、試してみたいだけです。慎重に使用してください。

  プログラムの途中で機能を繰り返しているコードがいくつかあるので、コメントアウトされています

  時間:3/27/17:30。バグが見つかり、コードが少し変更されました

  数学のインポートから*

  #ワイヤー伸びにおけるスケールデータの処理について

  list1 = []

  #list2 = []

  list3 = []

  #list4 = []

  m = 0

  #t = 0

  sum0 = 0#Differenceメソッドを使用してデータを取得し、差の2乗を合計します

  sum = 0#読み取り差法で得られたデータの合計

  一方1:

  data = eval(input( "ワイヤーの重量が小さいものから大きいものへと増加するときのスケール画像の読み取り値の平均値(ni)を入力し、負の値を入力して停止してください"))#サイクル入力

  データ<0の場合:

  ブレーク

  list1.append(1)

  list1 [m] = data#ストレージソースのテストデータを一覧表示します

  m + = 1

  #while 1:

  #data = eval(input( "小さいものから大きいものへ、読み取りなどの軽量化でワイヤのスケールを入力してください。停止するには負の値を入力してください"))#循環入力

  #if data <0:

  #ブレーク

  #list2.append(1)

  #list2 [t] = data#ストレージソースのテストデータを一覧表示します

  #sum + = data

  #t + = 1

  range(m-1)のyの場合:

  list3.append(1)

  list3 [y] = list1 [y + 1] -list1 [y]#変化量を見つけるための元の測定データ差分法

  sum + = list3 [y]

  #for y in range(t-1):

  #list4.append(1)

  #list4 [y] = list2 [y + 1] -list2 [y]

  #list3.extend(list4)

  dn_ = sum /(m-1)

  deltn = dn_

  print( "スケール読み取り値の変化の平均値は次のとおりです:{}"。format(dn_))

  list3 [:] = [y-dn_ for y in list3]#平均値が各リストから差し引かれます

  list3のiの場合:

  sum0 + = i ** 2

  on =(sum0 /(m-2))** 0.5#標準偏差を計算する

  on_ = on /((m-1)**(0.5))

  print( "標準偏差は次のとおりです:{}"。format(on))

  #線径dのデータ処理

  list = []

  sum1 = 0#合計スコア

  n = 0#総数

  平均= 0#平均スコア

  sum2 = 0#差分二乗和

  一方1:

  data = eval(input( "複数のワイヤの修正された直径diを入力してください、停止するには負の値を入力してください"))#サイクル入力

  データ<0の場合:

  ブレーク

  list.append(1)

  list [n] = data#ストレージソースのテストデータを一覧表示します

  sum1 + =データ

  n + = 1

  平均= sum1 / n#dの平均値を計算します

  deltaverage =平均

  print( "線径の平均値は次のとおりです:{}"。format(average))

  list [:] = [y-リスト内のyの平均]#各リストは平均から差し引かれます

  リストにある私のために:棗荘ピープルフロー病院にとってどちらが良いかhttp://mobile.3635555.cn/

  sum2 + = i ** 2

  oi =(sum2 /(n-1))** 0.5#標準偏差を計算する

  hi_ = hi /(n **(0.5))

  print( "標準偏差odは:{}"。format(oi))

  print( "odを部首で割った値n = {}"。format(oi_))

  #テストした金属線のヤング率を計算します

  d = eval(input( "入力ミラー距離(mm):"))* 0.001

  l = eval(input( "入力ワイヤ長(mm):"))* 0.001

  k = eval(input( "入力光学レバー定数(mm):"))* 0.001

  dn = eval(input( "目盛りの読みの変化の平均値(mm)を入力してください"))* 0.001#国際単位系に変換

  平均* = 0.001#国際単位系へ

  from = 0

  ol = 0

  ok = 0#実験では1セットのデータしか測定されなかったため、計算を簡単にするために、当面は0として扱います。

  om = 0.003

  f = 9.79163

  of = om *(9.79163)

  E _ =(8 * d * l * f)/(pi *(平均** 2)* k * dn)

  Ee =((ol / l)** 2 +(of / f)** 2 +(od / d)** 2 +(oi_ / deltaverage)** 2 +(ok / k)** 2+(on_ / deltn)** 2)**(0.5)

  oe_ = Ee * E_

  print( "E = {} +/- {}"。format(E_、oe_))

  print( "Ee = {} * 100%"。format(Ee))

  #コードの終わり

おすすめ

転載: blog.51cto.com/14503791/2676571