主权项 |
一种导波信号包络线绘制方法,其步骤如下:步骤1从所采集并经过滤波的导波回波信号中提取I个数据x<sub>i</sub>,i=1,2,…I,作为待分析数据,所述的I个待分析数据为第一次待检回波信号到达前后由导波传感器获取的导波信号,该信号必须包括至少一次激励波和待检回波,引入循环计数变量P,用于记录循环次数,令P=1,之后执行步骤2;步骤2将步骤1所述的I个待分析数据x<sub>i</sub>,i=1,2,…I的每一个元素取绝对值,得到一新的数据序列|x<sub>i</sub>|,i=1,2,…I,并同时记录对应的采集时间序列为t<sub>i</sub>,i=1,2,…I,之后执行步骤3;步骤3寻找步骤2所得数据序列|x<sub>i</sub>|,i=1,2,…I中所有极大值点A<sub>j</sub>,j=1,2……J,并同时记录下各极大值所对应的时间,得对应的时间序列t<sub>j</sub>,j=1,2,…J,总分析区间记为t<sub>总</sub>=t<sub>J</sub>‑t<sub>1</sub>,之后执行步骤4;步骤4计算步骤3所得极大值序列A<sub>j</sub>,j=1,2……J中各点的导数,形成新的数据序列A’<sub>j</sub>,j=1,2……J,其对应的时间序列不变,仍为t<sub>j</sub>,j=1,2,…J,之后执行步骤5;步骤5对步骤3所得极大值序列A<sub>j</sub>,j=1,2……J进行插值运算,包括如下步骤:1)将总分析时间t<sub>总</sub>划分为J‑1个独立的插值区间,如[t<sub>j</sub>,t<sub>j+1</sub>]为第j个独立插值区间,之后执行步骤2);2)在第j个插值区间内,以原采样周期T为插值步长,将插值区间长度[t<sub>j</sub>,t<sub>j+1</sub>]除以采样周期T,得第j个区间的插值点数为M<sub>j</sub>=(t<sub>j+1</sub>‑t<sub>j</sub>)/T‑1,j=1,2......J‑1,每个插值点处对应的时间为t<sub>m</sub>=t<sub>j</sub>+m*T,m=1,2......M<sub>j</sub>,之后执行步骤3);3)在第j个插值区间内,计算每个插值点处的基函数α<sub>jm</sub>(t),β<sub>jm</sub>(t),j=1,2......J‑1,m=1,2......M<sub>j</sub>;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>α</mi><mrow><mi>j</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>t</mi><mi>m</mi></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><mrow><msub><mi>t</mi><mi>j</mi></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mn>2</mn><mfrac><mrow><msub><mi>t</mi><mi>m</mi></msub><mo>-</mo><msub><mi>t</mi><mi>j</mi></msub></mrow><mrow><msub><mi>t</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>t</mi><mi>j</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><msub><mi>t</mi><mi>m</mi></msub><mo>∈</mo><mrow><mo>(</mo><msub><mi>t</mi><mi>j</mi></msub><mo>,</mo><msub><mi>t</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000794710320000021.GIF" wi="1261" he="175" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>β</mi><mrow><mi>j</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>t</mi><mi>m</mi></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><mrow><msub><mi>t</mi><mi>j</mi></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow></mfrac><msup><mo>)</mo><mn>2</mn></msup><mrow><mo>(</mo><msub><mi>t</mi><mi>m</mi></msub><mo>-</mo><msub><mi>t</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>t</mi><mi>m</mi></msub><mo>∈</mo><mrow><mo>(</mo><msub><mi>t</mi><mi>j</mi></msub><mo>,</mo><msub><mi>t</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000794710320000022.GIF" wi="1076" he="171" /></maths>4)在第j个插值区间内,计算每个插值点的值<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>g</mi><mrow><mi>j</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mi>j</mi></msub></munderover><mo>[</mo><msub><mi>α</mi><mrow><mi>j</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>A</mi><mi>j</mi></msub><mo>+</mo><msub><mi>β</mi><mrow><mi>j</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msubsup><mi>A</mi><mi>j</mi><mo>′</mo></msubsup><mo>]</mo><mo>;</mo></mrow>]]></math><img file="FDA0000794710320000023.GIF" wi="907" he="189" /></maths>5)将所有的插值点与原极大值点按时间顺序组成新的数据序列y<sub>i</sub>,i=1,2,…I,之后执行步骤6;步骤6判断P是否等于2,若P不等于2,令P=P+1,对步骤5所得数据序列y<sub>i</sub>,i=1,2,…I,重复步骤2至步骤5,得到数据序列z<sub>i</sub>,i=1,2,…I,之后执行步骤7;若P等于2,执行步骤7;步骤7平滑连接步骤6所得数据序列z<sub>i</sub>,i=1,2,…I,得到一平滑曲线,即为待检导波信号的包络曲线。 |