Python计算防蓝光眼镜加权阻隔率

心已赠人 2023-01-18 13:23 25阅读 0赞

400nm以下是紫外线波段，

根据透光率可以看出是基本满足uv400标准的。

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70][]

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 1][]

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 2][]

＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 3][]

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 4][]![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 5][]

＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃

**透光率拟合.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,#380
    	0.1*0.01,#385
    	0.1*0.01,#390
    	0.1*0.01,#395
    	0.1*0.01,#400
    	1.5*0.01,#405
    	2.9*0.01,#410
    	5.2*0.01,#415
    	8.1*0.01,#420
    	9.5*0.01,#425
    	10.3*0.01,#430
    	10.6*0.01,#435
    	11.2*0.01,#440
    	11.9*0.01,#445
    	13.1*0.01,#450
    	13.9*0.01,#455
    	15.7*0.01,#460
    	16.7*0.01,#465
    	19.4*0.01,#470
    	21.0*0.01,#475
    	26.3*0.01,#480
    	29.2*0.01,#485
    	33.4*0.01,#490
    	39.3*0.01,#495
    	46.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元不到的售价。

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70]: /images/20221021/a65a3d5a37204c2ca87c22fbb4a6a301.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 1]: /images/20221021/d520e025856c4f67ba6bef1735def507.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 2]: /images/20221021/4e9aadf09af24b58a3f7f2aa64fccf21.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 3]: /images/20221021/814ed74958154e0d8430396132fb7e94.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 4]: https://img-blog.csdnimg.cn/2021050616381445.jpg?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk=,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2FwcGxleXVjaGk_size_16_color_FFFFFF_t_70 5]: /images/20221021/71376b28d1af4a9089c0225c0bfe5872.png