Raspberry Pi Picoで標準偏差を求める。

Question

Raspberry Pi Picoで電圧を64回取得して、その標準偏差を出力したいと思い、下記のmicropythonのプログラムを書きました。

コード
import machine
import utime

sensor_smoke = machine.ADC(0)
conversion_factor = 3.3 / (65535)

reading =  sensor_smoke.read_u16() * conversion_factor

def data():
  result = []
  for reading in range(64):
      reading =  sensor_smoke.read_u16() * conversion_factor
      print(reading)
      utime.sleep(1)
      result.append(reading)
  return result

list1=data()
print(list1)

number = len(list1)
print("len:"+str(number))
avg = sum(list1)/len(list1)
print("avg:" + str(avg))


def func(x):
    return pow(2,x-avg)

resultk=0

for i in list1:
    resultk = resultk + func(i)

resultkk = resultk /number
s = pow(0.5, resultkk)

print("standard_deviation:" + str(s))
       

上記のプログラムがあっているかどうか確かめるために、for reading in range(64)の64を3に代えると、
1.300664
1.292607
1.293413
[1.300664, 1.292607, 1.293413]
len:3
avg:1.295561
standard_deviation:0.499999
が出力されました。

windowsでnumpyを使って、下記のプログラムで確かめると、標準偏差は、0.0036231031205614826となり間違っているようですが、micropythonのプログラムのどこが間違っているのでしょうか？

コード
import numpy as np

x  = [1.300664, 1.292607, 1.293413]

standard_deviation =np.std(x)

print(standard_deviation)

Answer 1

pow()の使い方が間違っています。

return pow(2,x-avg)
s = pow(0.5, resultkk)

pow(x, y)はxをy乗するという意味です。
前者は2を(x-avg)乗するのではなく、(x-avg)を2乗するのが正しいかと思います。
後者は0.5をresultkk乗するのではなく、resultkkを0.5乗するのが正しいかと思います。
これらの引数を逆にすると、正しく0.0036231031205614826と表示されるはずです。

前者の修正例:

def func(x):
    return pow(x-avg, 2)

後者の修正例:

s = pow(resultkk, 0.5)