机械零部件在工作状态下受各外力作用，内部会产生位错、微裂纹、滑移带等微损伤，常规超声方法由于其局限性，不能对上述微损伤进行有效检测. 非线性超声检测技术通过分析超声信号的频域特征可量化上述材料内部微损伤程度[1-3]，并通过时间反转技术[4]对损伤位置进行定位，其中超声非线性系数β是表征材料中微损伤程度的重要参数[5].

使用β表征微缺陷的大小时，β由超声信号的基波与二次谐波幅值决定[6]. 由于超声信号二次谐波幅值较小，容易受到各种软硬件干扰，故如何提高β的准确性是目前的研究热点. Thiele等[7]使用空气耦合传感器来降低检测中对被测件表面粗糙度的要求；高鹏等[8]设计和开发了非线性超声相控阵无损检测系统用于提高检测结果的空间分辨力；张荣繁等[9]研究了非线性声波衰减系数对非线性超声检测结果的影响.

除了硬件系统的影响外，对超声信号进行处理时，考虑到被测件的尺寸以及衰减等因素，发射的超声激励信号是一系列的脉冲串，而且由于采样率的限制，采集的可用于进行超声非线性系数计算的数据长度非常有限，导致超声信号频谱分析结果中，其频率分辨率较低. 由于超声在传播过程中出现了波形畸变，使得截断后的时域波形周期数与数据点数不协调，同时由于较低的频率分辨率，最终超声信号频谱分析结果中特征频率成分处可能不存在谱线[10]或引入较大的计算非线性，加窗虽能在一定程度上抑制信号边缘部分的非线性因素，但不能完全去除，同时加窗也会使得信号整体幅值谱幅值降低，后续需要进行一些幅值修正处理.

三是应坚持以消除突出瓶颈为先导。按照先易后难、先系统内后系统外、先干线后支线、先省内后区域的原则，将天然气主干管网、支线管网、LNG接收站、储气库和天然气生产设施等全部连接到一起，形成省内、省际、区域和全国性的天然气供气输送和配售系统，打造天然气基础设施“一张网”。

本文在通过对离散傅里叶变换(discrete Fourier transform，DFT)进行理论分析以及考虑超声信号不同截断长度、波形传播过程中的畸变等因素对超声非线性系数计算的影响的基础上，提出基于激励信号特征数据长度和端点数据特征的超声非线性系数计算方法，结果表明该方法能有效降低β的计算误差.

1 截断及较低频率分辨率对β的影响

1.1 超声非线性系数

为了量化固体中非线性超声传播规律，表征被测物体中的微缺陷，Breazeale等[11]从连续介质出发建立了固体介质一维纵波非线性波动方程，利用逐级近似微扰法，略去高次项，得到超声非线性系数β的表达式 [11]

(1)

式中：k为波数；x为波传播的距离；A1、A2分别为基波与二次谐波的幅值. 在实验中，k和x为固定的常数，β值由A1和A2决定，所以常使用来代替β值. 为了得到准确的β值，需要精确测量A1和A2的值. 一般使用透射法对被测件进行非线性超声检测，然后对接收到的信号进行幅值谱分析，得到A1和A2的值.

1.2 存在的问题

图1所示的非线性超声检测系统中，对施加227 MPa 拉伸应力的镁合金试件采用透射法，根据试件尺寸(60×φ13 mm)，确定发射信号频率为5 MHz，周期数为51，信号峰峰值为300 mV，类型为单频正弦脉冲串激励信号. 信号发生器产生的激励信号经功率放大器放大1 000倍达到300 V，使用高能低通滤波器滤除高阶谐波分量后，驱动发射传感器，向试件输入超声波. 超声波被接收传感器接收，并被示波器采集. 示波器采集到的原始信号如图2所示.

图1 超声非线性系数测量系统及试件 Fig.1 Ultrasonic nonlinear parameter measurement system and specimen

图2 原始超声信号 Fig.2 Original ultrasonic signal

从原始信号时域图形可看出，中间幅值较高部分的波形幅值明显产生了畸变，后边幅值较低部分为回波信号，利用Matlab对采集到的时域波形进行分析，实际采样率为500 MHz，采样点数为10 000点. 在进行频谱分析时，截取幅值较高部分的51个完整周期(图3所示波形实际数据点数为5 118点). 对图3所示信号进行快速傅里叶变换(fast Fourier transform，FFT)并求幅值谱，得到图4所示幅值谱图.

图3 一般截断信号 Fig.3 Truncation signal in general

图4 一般截断信号幅值谱 Fig.4 Amplitude spectrum of truncation signal in general

2) 通过补零，使得时域信号数据点数为采样率与基波比值的整数倍，可以保证基波与二次谐波的谱线不会偏离，同时，可以使信号边缘部分引入的计算非线性较小.

上述5 MHz峰值点右侧频谱谱线的频率值为4.89×106，10 MHz峰值点右侧频谱谱线的频率值为1.01×107，由此看出，按照完整周期数进行波形截取时，由于示波器采样率较高，而实际有效数据段长度较短，造成信号幅值谱实际频率分辨率较低，5 MHz与10 MHz处不存在谱线，信号幅值谱分析中已经不包含基波与二次谐波的准确幅值.

为了计算β值，一般使用最靠近5 MHz与10 MHz处的峰值点幅值代替基波与二次谐波幅值，由此得到的β值与实际β值之间存在一定的误差. 根据以上两峰值点幅值，计算得到的β值为0.005 05.

2 理论分析与解决方法

2.1 DFT分析

DFT的公式定义为

(2)

式中：X(j)为时域信号幅值数组；Y(k)为时域信号的DFT计算结果；Wn=e(-2πi )/n，n为时域信号幅值数组的长度，i=1,2,…,n；j，k=1, 2, 3，…，n.

通过对DFT算法以及信号时域波形特征的分析，提出以下超声信号截断方法：

fx(i)=iFs/n

(3)

式中：fx(i)为幅值谱第i点的频率值；Fs为采样率；n为输入的时域信号幅值数组长度.

当时域信号采样率Fs和基波频率f1为某一固定值时，从式(3)可以看出时域信号幅值数组长度n的大小将直接决定超声时域信号幅值谱的频率分辨率，不合适的时域信号幅值数组长度n将会导致没有幅值谱谱线会落在基波频率f1与二次谐波频率f2处.

将图3中的数据段实际点数5 118代入式(3)，当i分别等于51、52、102、103时，对应的fx(i)分别等于4.98×106、5.08×106、9.97×106、1.01×107 Hz，与图4中得到的结果一致，因此，在图4幅值谱图中，不存在频率成分刚好等于5 MHz与10 MHz的幅值谱谱线.

2.2 超声信号截断方法

信号段幅值谱的频率轴位置计算公式为

1) 按照激励信号周期数，在接收信号中幅值最高部分截取相同周期数的完整波形，要求首尾幅值一致且大小接近零.

2) 用采样率与基波频率的比值乘以激励信号周期数得到理想信号长度，将1)中实际截取的信号长度与理想信号长度进行比较，若实际截取信号的长度小于理想信号长度，则对实际信号两端进行对称补零以达到理想信号长度；若实际截取信号的长度大于理想信号长度，则对激励信号周期数加一再乘以采样率与基波频率的比值，从而得到次优的理想信号长度，然后对实际截取信号两端进行对称补零达到次优理想信号的长度.

3 实验研究

2) 当实际被处理数组长度n为Fs/f1的整数倍，但是长度超出激励脉冲串信号中正弦波形幅值数组长度时，超声非线性系数的计算结果呈增大趋势.

3.1 仿真信号

由式(1)可知，超声非线性系数的计算主要与基波和二次谐波有关，因此，在建立理想超声信号时可只考虑基波和二次谐波，将其他无关的信号分量忽略掉. 根据实际实验参数，将仿真信号的基波频率f1设置为5 MHz，幅值为1，周期n1为10，二次谐波频率f2设置为10 MHz，幅值为0.1，周期n2为20，采样率Fs设置为500 MHz. 利用正弦函数分别生成基波和二次谐波幅值数组(相位均为0)，然后将二者幅值叠加生成叠加信号.

由于非线性超声检测中，实际发射的是一系列脉冲串，为了尽可能模拟真实数据，将上述叠加后的幅值数组前后分别补100个零点. 由(Fs/f1)n1可知，最终生成的正弦波形幅值数组长度为1 000点，加上两侧一共补的200个零点，图5中最终生成的超声仿真信号幅值数组的长度为1 200点.

最终生成的超声信号仿真波形如图5所示.

图5 仿真信号 Fig.5 Simulation signal

针对图5中的超声仿真信号，图6给出了仿真信号幅值数组在1 080、1 040、1 000、960、920点5种长度时的幅值谱图. 在以上述点数对信号进行截断时，以1 080点为例，先用信号总长(1 200点-1 080点=120点)得到总共应该被去除的数据点数，然后在仿真时域信号两端各减去60个数据点，取中间剩余部分为被分析信号.

术中：①手术视野较小，没有清晰的显露手术视野：脊柱外科手术中出血较多，视野不清，以及椎板“开窗”术或半椎板切除术术野局限，导致粘连松解或椎间盘摘除时对硬脊膜造成损伤；②术中操作不当：手术医生没有足够的实践经验、术中操作过程不熟练或者术中动作粗暴，都有可能会对硬膜造成损伤进而引发脑脊液漏，特别是在棘突斜坡进行减压时钳咬椎板或黄韧带时角度及深浅度不当，造成硬脊膜被咬损或撕破，椎管内肿瘤摘除切开硬脊膜后硬脊膜缝合方式不对等。

图6 不同截断点数的仿真信号频谱 Fig.6 Spectrum of the simulation signal for different truncated points

对比图6中1 080、1 040、1 000、960、920点5组信号分别在5 MHz与10 MHz处的谱线可以发现，数据长度为1 000点时，5 MHz与10 MHz处频谱谱线与频率轴上5 MHz与10 MHz处的偏离最少. 其他4组信号的频谱分析结果在5 MHz与10 MHz处均发生不同程度偏离.

什么？经血发黑就是体内有毒？你是不是武侠小说看多了啊？那是不是每次换卫生巾，还要拿银针试一下啊？女人每个月排毒一次，那男人不是早就应该中毒身亡了？

为了进一步验证信号截断对超声非线性系数的影响，依次对上述仿真超声波信号两侧的数据进行删除，每次两侧各删除10个点(单次共计删除20个点)，考察信号幅值数组长度从1 200点降到800点的情况，对应的超声非线性系数有21组. 超声非线性系数随信号段长度的变化情况如图7所示，其中实际值是由仿真超声信号截断后做FFT(不加窗)，根据FFT结果，计算得到. 其中标准值是由基波与二次谐波的初始设置参数直接计算得出的，标准值为一固定值，用于衡量实际值的偏差情况.

人们对美好生活的向往，不仅要有充足的食物，更需要有丰富、安全的食物，更需要有优质的食物，更需要有自己喜欢的、舒适的环境，不能每人都穿着航空服，吃着所谓的“有机食品”，这样真的幸福吗？充满污染的土壤、水、空气、还有人为增加的那些化学品，还能生产出有机食品，谁信呀！

图7 不同截断点数下超声非线性系数的计算结果 Fig.7 Calculation results of ultrasonic nonlinear coefficients at different cutoff points

由图7可以看出，超声信号幅值数组的长度为800～1 000点，幅值数组的长度刚好为100的整数倍时，超声非线性系数的实际值与标准值的偏差较小，当幅值数组的长度由k×100(k=8,9,10)点逐渐变成(k×100+50)点且k=8,9时，超声非线性系数实际值会逐渐偏离标准值，且实际值与标准值的偏差在幅值数组的长度为(k×100+50)点且k=8,9时达到最大.

超声信号幅值数组的长度为1 000～1 200点时，幅值数组的长度等于k×100(k=8,9,10)点处的实际值最接近标准值. 但是考察点数在1 000、1 100、1 200点三处的实际超声非线性系数可以发现，实际值与标准值的偏差逐渐增加.

根据仿真信号的处理结果来看，主要有以下结论：

我省是国家重要的大豆生产和供给基地，种植面积和产量均位居全国之首。全省各种组织形式生产大豆单产相差较大，黑龙江农垦系统种植大豆平均亩产175kg，农垦系统大豆亩产约比地方合作社高15kg，而农民散户大豆亩产又比合作社低约25kg，主要是整体标准化栽培水平和技术到位率的差异，具体表现为以下几个方面：

钢筋安装应该对受力钢筋的牌号、规格、搭接长度等相关指标满足施工设计的规范，而关于钢筋绑扎的细部构造，也应该顺应施工质量相关验收指标，满足技术要求。钢筋施工属于较为隐蔽的工程，在进行混凝土浇筑过程中，应该顾及振捣扰动等诸多不利因素的产生，以此造成钢筋骨架与混凝土内出现偏移。钢筋进行绑扎过程中，应该着重绑扎的质量问题，必须满足相关规范要求以及技术指标，并且绑扎必须要牢固，不会出现松动[1]。

分别采用仿真信号和实测信号来验证超声时域信号长度对超声非线性系数的影响以及上述时域波形截断规则.

3.2 实测信号

针对图2中的原始超声信号，在图3所示的截断信号基础上进行进一步分析，给出了超声信号幅值数组为5 140、5 120、5 100、5 080、5 060点5种长度时如图8与图9所示的幅值谱图形.

图8 超声信号5 MHz处幅值谱 Fig.8 Amplitude spectrum of the ultrasonic signal at 5 MHz

图9 超声信号10 MHz处幅值谱 Fig.9 Amplitude spectrum of the ultrasonic signal at 10 MHz

由图9、10可以发现，当超声信号幅值数组长度等于5 100点时，5 MHz与10 MHz频率成分的谱线基本落在5 MHz与10 MHz处，而当仿真信号幅值数组长度为5 140、5 120、5 080、5 060点时，幅值谱在5 MHz与10 MHz频率成分处均发生了偏离.

通过分级才能按级包装、定价、收购和销售。分级不仅可以贯彻优质优价的政策，而且可以推动果树栽培管理技术的发展和提高。通过分级、剔除病虫果和机械损伤果，既可使产品按大小分级便于包装标准化，又可减少在贮运中的损失，减轻一些危险病虫害的传播，并将这些残次产品及时销售或加工处理，以降低成本和减少浪费。总之，苹果的分级是苹果商品化生产中的一个重要环节，应引起高度的重视。

对采集到的信号依次进行截断，做频谱，求超声非线性系数，最终得到如图10所示的超声非线性系数随实验数据长度变化曲线. 其中超声非线性系数标准值(β=0.005 05)是由图4中频谱分析结果计算得出的，在此仅用于评价实际值的偏差情况.

1) 当被处理信号的幅值数组长度与激励脉冲串信号中正弦波形幅值数组长度接近，且实际被处理数组长度n基本为Fs/f1的整数倍时，超声非线性系数的计算结果比较准确.

图10 超声非线性系数随实验数据长度变化 Fig.10 Variation of ultrasonic nonlinear parameter with experimental data length

道路设计车速50km/h，双向6车道。红线宽40m，标准横断面：3m人行道（含树穴）+3m非机动车道+3m绿化带+22m机动车道+3m绿化带+3m非机动车道+3m人行道（含树穴）=40m（见图1），其中机动车道采用沥青混凝土路面，非机动车道采用水泥混凝土路面，人行道采用陶瓷步转。

由于回波等其他因素的影响，用于计算标准值的超声信号的实际长度为5 118点，与实验中Fs/f1≈100的整数倍依然有一定的偏差. 为了得到较准确的幅值谱结果，以截取的5 118点数据为基础，分为2种情况对其进行进一步处理.

1) 在后端减去18点得到5 100点数据，对5 100点数据进行分析得到其在5 MHz处的频谱谱线的坐标为(x:5.00×106 y:0.775)，10 MHz处频谱谱线坐标为(x:1.00×107 y:0.004)，以此为基础的超声非线性系数计算结果为0.007 20.

2) 在两端补83个零点得到5 200点数据，对5 200点数据进行分析，得到其在5 MHz处的频谱谱线的坐标为(x:5.00×106 y:0.761)，10 MHz处频谱谱线坐标为(x:1.00×107 y:0.003)，以此为基础的超声非线性系数计算结果为0.005 50.

商业人像摄影自从90年代初盛行全国以来，在摄影风格上大致经历了三个阶段。最初是国营照相馆的中规中矩的固定模式，人物造型比较古朴，不注重化妆造型，虽然讲究用光技法，但照片没有时代感，缺乏新意。而后，随着台湾影楼的大举进入，其平光重妆，讲究人物包装与造型多变的摄影风格成为全国影楼效仿的样板。90年代末，“前卫”、“时尚”摄影又成为潮流。迎合时尚是名利双收的一条捷径。创意摄影运用特殊手段以夸张的色彩、反常规构图、两底叠放等表现形式和手法，构筑作品的情调，刻画被摄者的性格和心态。

将以上2种情况得到的超声非线性系数与图10中标准值(β=0.005 05)进行比较，得到各自的相对偏差分别为42.57%和8.91%，由此可以判断通过补零得到5 200点数据时，最终得到的超声非线性系数与标准值的偏差更小.

3.3 实验分析

从仿真信号不同截断程度下的超声非线性系数来看，当对信号的截断不满足2.2节的要求时，计算得到的超声非线性系数最大约为标准值的200%；当满足要求时，实际计算得到的超声非线性系数与标准值之间的误差最多不超过标准值的5%，由此可以看出，2.2节提出的截断方法能有效降低仿真信号中超声非线性系数的计算误差.

实测信号中，当对信号的截断不满足2.2节的要求时，实际超声非线性系数与标准值的偏差最大为标准值的260%；通过2.2节的方法对时域信号进行处理后，实际超声非线性系数与标准值的偏差最多不超过标准值的10%.

4 结论

1) 结合仿真信号和实验信号的分析结果，在对超声信号进行截断处理时，根据激励信号周期数，在接收信号中幅值最高部分截取相同周期数的完整波形，可保证超声信号基波与二次谐波幅值的准确性.

由图10可以发现，实际值与标准值在信号幅值数组长度为k×100且k=50,51,52,53,54处的偏差值与信号幅值数组长度为k×100+50且k=50,51,52,53处的偏差值相比，后者的偏差明显比前者大. 图10中实验信号实际超声非线性系数的变化规律与仿真信号实际超声非线性系数的变化规律基本一致.

从图4中可以发现，基波频率5 MHz附近的峰值点坐标为x: 4.98×106,y: 0.718 (均保留3位有效数字)；二次谐波频率10 MHz附近的峰值点坐标为x:9.97×106,y: 2.59×10-3.

参考文献：

[1] 张剑锋, 轩福贞, 项延训. 材料损伤的非线性超声评价研究进展[J]. 科学通报, 2016, 61(14): 1536-1550.

纵观当前中小学教师数学阅卷，智学网阅卷早已是大势所趋。智学网作为我们日常各类考试的网络辅助平台，主要依托于计算机网络技术及其图像处理技术[1]。其中，选择题部分采用OMR技术，可智能识别学生涂卡信息，该部分由电脑自动判定得分。主观题部分，则采用图像切割技术，将考生所有主观题按照题序逐题完整剪切出来，同时将相应数据以特殊的命名方式上传至系统。待教师开始阅卷，系统会将学生答题数据进行还原，将剪切好的答题图像分发到教师阅卷终端，由教师进行人工判分。教师判分完毕后，由智学网系统自动统计每位学生的总分，并给出不同级别的学业报告与分析[2]。

ZHANG J F, XUAN F Z, XIANG Y X. Evaluation of material damage using nonlinear ultrasonic wave[J]. Science China Press, 2016, 61(14): 1536-1550. (in Chinese)

[2] 吴斌, 李佳锐, 颜丙生, 等. LY12铝合金早期性能退化下超声非线性系数测量和金相观察[J]. 北京工业大学学报， 2012， 38(1)： 22-27.

WU B, LI J R, YAN B S, et al. Ultrasonic nonlinear coefficient measurement and optical observation of LY12 aluminum alloy at early degradation stages[J]. Journal of Beijing University of Technology, 2012, 38(1): 22-27. (in Chinese)

[3] YAN B S, TANG B P. Simulation and testing of a non-linear ultrasonic system for testing the debonding of the adhesive layers of magnesium alloys[C]∥2016 13th International Conference on Ubiquitous Robots and Ambient Intelligence. New York: Institute of Electrical and Electronics Engineers, Inc. , 2016: 26-30.

[4] 张望, 吴维亮, 肖黎, 等. 闭合型裂纹时间反转损伤识别[J]. 振动与冲击, 2015, 34(9): 95-101.

ZHANG W, WU W L, XIAO L, et al. Time reversal based damage detection for closed crack[J]. Journal of Vibration and Shock, 2015, 34(9): 95-101. (in Chinese)

[5] CANTRELL J H. Quantitative assessment of fatigue damage accumulation in wavy slip metals from acoustic harmonic generation[J]. Philos Mag, 2006, 86(11): 1539-1554.

[6] WU B, YAN B S, HE C F. Nonlinear ultrasonic characterizing online fatigue damage and microscopic observation in situ[J]. Transactions of Nonferrous Metals Society of China, 2011, 21(12): 2597-2604.

(8) 车载安全设备中的安全通信模块接收到RSSP-I安全通信包后，进行安全协议解析和验证，并将验证通过的安全数据提交至安全应用模块进行处理。

[7] THIELE S, KIM J Y, QU J M, et al. Air-coupled detection of nonlinear Rayleigh surface waves to assess material nonlinearity[J]. Ultrasonics, 2014, 54(4): 1470-1475.

而对B班同学，请同学们观察下面这个图形，说说它有几条边，有几个角，并且测量一下边的长度，角的度数，说说这些边和角有什么大小关系。

[8] 高鹏, 李法新. 非线性超声相控阵无损检测系统及实验研究[J]. 实验力学, 2014, 9(1): 1-11.

GAO P, LI F X. Nonlinear ultrasonic phased array nondestructive testing: system and performance study[J]. Jounal of Expermental Mechanics, 2014, 9(1): 1-11. (in Chinese)

研究发现，在人际知觉或是印象形成过程中，个体对社群性维度的认知加工要优先于对能动性维度的加工，因为判断他人的意图(是否会伤害自己)要比判断他人是否有能力完成该意图更重要(Wojciszke & Abele，2008)。Ybarra等人发现社群性词语能够较快的被个体识别，个体对社群性词语进行褒贬分类的时间也较短(Ybarra，Chan， & Park，2001)。Abele和Bruckmüller(2011)也发现了同样现象，当要求个体根据给出的行为描述进行社群性和能动性特质推断时，个体对社群性特质的推断明显较快(Abele & Bruckmüller，2011)。

[9] 张荣繁, 李雄兵, 张书增, 等. 材料非线性衰减系数的二次谐波测量方法研究[J]. 声学学报, 2017, 42(2): 223-229.

ZHANG R F, LI X B, ZHANG S Z, et al. Research on nonlinear acoustic attenuation coefficients of material based on second harmonic measurement[J]. Acta Acustica, 2017, 42(2): 223-229. (in Chinese)

[10] 温和, 滕召胜, 王永, 等. 改进加窗插值FFT动态谐波分析算法及应用[J]. 电工技术学报, 2012, 27(12): 270-277.

WEN H, TENG Z S, WANG Y, et al. Improved windowed interpolation FFT algorithm and application for power harmonic analysis[J]. Transactions of China Electrotechnical Society, 2012, 27(12): 270-277. (in Chinese)

[11] BREAZEALE M A, PHILIP J. Determination of third order elastic constants from ultrasonic harmonic generation measurements[C]∥Physical Acoustics: Principles and Methods. Orlando: Academic Press, Inc., 1984: 1-60.