Python计算防蓝光眼镜加权阻隔率
生活随笔
收集整理的這篇文章主要介紹了
Python计算防蓝光眼镜加权阻隔率
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
400nm以下是紫外線波段,
根據透光率可以看出是基本滿足uv400標準的。
?
?
?
?
###################################################################################################
?
############################################################################################################
透光率擬合.py
import numpy as np import matplotlib.pyplot as plt #定義x、y散點坐標(下面這個是波長) x = [380,390,400,410,420,430,440,450,460,470,480,490,500] x = np.array(x) print('波長范圍 is :\n',x) num = [0.1,0.1,0.1,2.9,8.1,10.3,11.2,13.1,15.7,19.4,26.3,33.4,46.8] y = np.array(num) print('y is :\n',y) #用3次多項式擬合 f1 = np.polyfit(x, y, 5) print('f1 is :\n',f1) p1 = np.poly1d(f1) print('p1 is :\n',p1)#也可使用yvals=np.polyval(f1, x) yvals = p1(x)#擬合y值 print('yvals is :\n',yvals) #繪圖 plot1 = plt.plot(x, y, 's',label='original values') plot2 = plt.plot(x, yvals, 'r',label='polyfit values') plt.xlabel('x') plt.ylabel('y') plt.legend(loc=4) #指定legend的位置右下角 plt.title('polyfitting') plt.show()?
加權阻隔率計算.py
import numpy as np from scipy.interpolate import InterpolatedUnivariateSplineλ = np.array([380,385,390,395,400,405,410,415,420,425,430,435,440,445,450,455,460,465,470,475,480,485,490,495,500])#波長λ W_λ=np.array([2 ,4 ,10 ,22 , 47,112,269,564,660,722,771,849,911,930,946,933,864,776,706,639,532,479,266,194,122])#WB(λ) # 下面是透射率(一部分是插值后得到的數據) τ=np.array([0.1*0.01,#3800.1*0.01,#3850.1*0.01,#3900.1*0.01,#3950.1*0.01,#4001.5*0.01,#4052.9*0.01,#4105.2*0.01,#4158.1*0.01,#4209.5*0.01,#42510.3*0.01,#43010.6*0.01,#43511.2*0.01,#44011.9*0.01,#44513.1*0.01,#45013.9*0.01,#45515.7*0.01,#46016.7*0.01,#46519.4*0.01,#47021.0*0.01,#47526.3*0.01,#48029.2*0.01,#48533.4*0.01,#49039.3*0.01,#49546.8*0.01#500])f = InterpolatedUnivariateSpline(λ, W_λ, k=3) #分母曲線擬合 a=f.integral(380, 500)#積分運算分母 print("分母:",a)f2=InterpolatedUnivariateSpline(λ,np.multiply(τ,W_λ),k=3)#分子曲線擬合 b=f2.integral(380, 500)#積分運算分子 print("分子:",b)print("加權透射率:",b/a) print("加權阻隔率:",1-b/a)?
上述是對外宣稱阻隔率70%的藍光眼鏡。
該品牌(為了通過csdn審核,已經取消名字)阻隔藍光99%的藍光眼鏡也只有30元不到的售價。
總結
以上是生活随笔為你收集整理的Python计算防蓝光眼镜加权阻隔率的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: ps性能首选项在哪里
- 下一篇: win10中Tiworker.exe占用