主权项 |
一种表征腐蚀速度的电化学噪声处理方法,其特征在于,按照下述步骤进行:(1)采用三电极体系传感器采集电化学噪声信号;得到采样条件下的电压值和电流值,分别做两者的标准差,再将两者相除即可得到腐蚀电阻;(2)利用小波变换的多分辨分析对电压、电流信号按预设层数进行分解;(3)将分解的各层电压、电流信号进行重构;(4)将各层对应的频率换算成实际信号的频率;(5)计算不同频率下的电压、电流标准差,得到小波噪声电阻;采用Mallat算法和Sym4小波函数进行分解和重构,用式(1)、(2)对所述的电压、电流信号按预设层数进行多分辨分析分解:<maths num="0001"><math><![CDATA[<mrow><msub><msup><mi>ca</mi><mi>m</mi></msup><mi>n</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mo>-</mo><mi>∞</mi></mrow><mi>∞</mi></munderover><mover><msub><mi>h</mi><mrow><mi>k</mi><mo>-</mo><mn>2</mn><mi>n</mi></mrow></msub><mo>‾</mo></mover><msup><msub><mi>ca</mi><mi>k</mi></msub><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001024233740000011.GIF" wi="1670" he="161" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><msup><mi>cd</mi><mi>m</mi></msup><mi>n</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mo>-</mo><mi>∞</mi></mrow><mi>∞</mi></munderover><mover><msub><mi>g</mi><mrow><mi>k</mi><mo>-</mo><mn>2</mn><mi>n</mi></mrow></msub><mo>‾</mo></mover><msup><msub><mi>ca</mi><mi>k</mi></msub><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001024233740000012.GIF" wi="1632" he="156" /></maths>分解表示成:<img file="FDA0001024233740000013.GIF" wi="1554" he="300" />其中ca<sup>j</sup>和cd<sup>j</sup>分别是2<sup>j</sup>分辨率下的离散逼近信号和离散细节信号,其中ca<sup>j</sup>针对频率不超过2<sup>‑j</sup>部分,cd<sup>j</sup>针对频率在2<sup>‑j</sup>和2<sup>‑j+1</sup>中间的部分;m=1,2…n,n为多分辨分析的预设层数;k为电压、电流信号的数据点数;<img file="FDA0001024233740000014.GIF" wi="154" he="70" />为低通滤波器系数,根据选取的小波波形确定;<img file="FDA0001024233740000015.GIF" wi="150" he="60" />为高通滤波器系数,根据选取的小波波形确定;用式(4)对生成的新生成的信号进行重构:<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>c</mi><mi>n</mi><mi>m</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><msqrt><mn>2</mn></msqrt></mfrac><mrow><mo>(</mo><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mo>-</mo><mi>∞</mi></mrow><mi>∞</mi></msubsup><mover><msub><mi>h</mi><mrow><mi>k</mi><mo>-</mo><mn>2</mn><mi>n</mi></mrow></msub><mo>‾</mo></mover><msubsup><mi>ca</mi><mi>k</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mo>-</mo><mi>∞</mi></mrow><mi>∞</mi></msubsup><mover><msub><mi>g</mi><mrow><mi>k</mi><mo>-</mo><mn>2</mn><mi>n</mi></mrow></msub><mo>‾</mo></mover><msubsup><mi>ca</mi><mi>k</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001024233740000021.GIF" wi="1687" he="100" /></maths>其中,<img file="FDA0001024233740000022.GIF" wi="70" he="60" />为重构后的电压、电流信号;小波噪声电阻由式(5)、(6)所得:Rcd(j)=S(Ecd(j))/S(Icd(j)) (5)Rca=S(Eca)/S(Ica) (6)其中j是频率因子,S表示信号的标准偏差。被分解的信号频率低于j的增长,当信号频率是f时,则j尺度下的信号频率f<sub>j</sub>=f*2<sup>‑j</sup>。 |