0 引言

内弹道学是研究弹丸在膛内运动规律及其伴随射击现象的一门学科［1］。射击实验中，所经历的时间只有几毫秒到十几毫秒，由于膛内过程具有高温、高压、高速和瞬态的特点，给内弹道实验研究带来一定的困难。20世纪50年代和60年代，国内外的技术人员完成了微波干涉仪的研制，从而可以用实验方法间接测出弹丸在膛内的行程-时间曲线，为内弹道学的膛内运动规律研究提供一种必不可少的技术手段。随着科学技术的发展，微波、毫米波干涉仪［2-3］，激光干涉仪［4］等内弹道运动规律测试系统逐步完善与成熟，应用也越来越广［5-6］。在所有瞬时运动参数测量过程中，瞬时频率的估计精度决定了各参数的测量精度。微波干涉测量系统中瞬时频率估计的方法有多种，如极值法、短时傅立叶变换（STFT）、维格纳—威利变换（WVD）、S变换等。

1028 阿替普酶静脉溶栓治疗急性心源性缺血性脑卒中的预后影响因素分析 沈 芳，沈红健，邢鹏飞，姜 乐，黄石仁，张永巍，吴 涛，邓本强

极值法在时间上具有较高的分辨力，但是极易受到噪声的干扰使极值点的位置产生偏差，针对这一问题，研究者也提出了一些改进方法［7-8］，通过降低搜索的极值点位置产生的偏差来确保频率估计的精度，在一定程度上降低噪声，但当信号完全淹没在噪声中时，此方法频率的估计将无法保证；STFT是在时域信号上加窗函数进行傅立叶变换，是一种局部频谱分析方法，无法同时得到高的时域和频域测量精度，可以通过采用自适应改变窗口大小或者同其他方法联合弥补其不足［9-10］；WVD是信号时间自相关函数的傅立叶变换，具有良好的时间、频率分辨率和时频聚集性，但存在交叉项干扰，研究者也针对这个问题进行了抑制交叉项的研究 ［11-12］；文献［13］采用了改进的S变换来提高速度的测量精度。上述的方法各有其优缺点，使用者大多取其优点避其缺点。

而不断优化产品结构，着力打造多档次、多品牌、多规格的产品体系，是茅台一直以来努力的方向。高、中、低端产品的全覆盖，使之成为了与高端洋酒竞争的锐利武器。对于年轻人、外国人等消费群体，茅台则开发生肖纪念酒、珍品茅台酒等一系列符合消费者个性化消费潮流的新产品，在重大历史事件期间，更是推出了“奥运纪念酒”“开国盛世茅台酒”“上海世博会纪念酒”等，使茅台的产品更为丰富。

在时频联合分析中，时频聚集性反映了能量集中在瞬时频率周围的情况。针对信号模型构造匹配的变换核，通过选择恰当的参数，可得到具有较好的时频聚集性且没有交叉项干扰的时频分布。弹丸膛内运动的回波信号是非线性的调频信号，可采用多项式相位信号模型对其数学描述，本文采用多项式调频小波（PCT）［14-15］构造匹配变换核来对弹丸膛内运动回波信号进行变换，通过对叠加噪声的模拟信号和实际实验信号的分析，得到较高的时频聚集性，在一定频率分辨率的基础上，保障了瞬时频率估计精度。本文还针对实验信号在同等条件下对比了PCT和WVD分析的效果，结果说明PCT时频分布的优越性。

利用内弹道方程组［1］，假设微波干涉仪中心频率94 GHz，则模拟的弹丸膛内运动回波信号随时间变化的曲线如图1所示，是典型的非线性调频信号。

1 多项式调频小波变换原理

1.1 基本原理

参数化时频方法需根据信号的模型，构造匹配的变换核。如果变换核和信号模型较吻合，则可取得较好的时频聚集性，时频能量聚集性与变换核的参数选择有关。弹丸在膛内运动是一个变加速的过程，运动回波信号的频率随时间做非线性变化，是典型的非线性调频信号。频率随时间变化可采用多项式拟合逼近，因此，可构造多项式调频小波变换核来获取较好时频聚集性。

多项式调频解析信号模型为

 

其中，n为多项式的次数，为多项式相位系数，其瞬时频率为

 

根据参数化时频方法［16］的基本思想，构造如下多项式变换核kp（t）

Step5：判断变换核参数是否最佳，如果最佳，结束程序；如果不是最佳，将变换核参数设置为初始变换核参数，返回Step3。

 

其中，为多项式变换核参数，因此，多项式调频小波变换为：

当前，大数据技术刚刚兴起，矿业生产逐渐对大数据技术概念引起重视，未来，随着现代化感知设备的不断开发，新的矿业软件引入，智能时代的开启，对数据的分析将会被智能设备自动检测、分析、调整等，这一系列流程操作将变得更加容易，大数据的应用价值也将会在工业生产中得到充分发掘。

安徽省淮南市大通区生活垃圾填埋场启用于1984年，于2009年6月停止使用，总占地面积110 000 m2，堆体顶部最大标高为50.5 m，填埋最大深度18.1 m，垃圾堆体约8.0×105m3，其天然密度为1.3 g/cm3。由于原有垃圾填埋场未做防渗，渗沥液导排设施导排不畅，而且填埋场一直未进行雨污分流，垃圾含水量很高，渗沥液向场地四周扩散，污染周边农田及地下水。

 

其中

 

根据参数化时频方法的基本原理，将多项式调频小波变换应用到弹丸膛内运动回波的瞬时频率估计。首先在模拟回波信号的基础上叠加不同信噪比的噪声干扰信号，分别利用WVD、PCT两种方法进行时频分布，可以看出PCT方法具有较强的抗噪声干扰能力。在频率窗口一致的条件下，即频率分辨率，频率估计精度一样的情况下，PCT能进行较好的瞬时频率估计，WVD由于自交叉干扰的问题使极大值提取的脊线有两条，必须先消除交叉干扰影响再进行瞬时频率估计。最后利用WVD、PCT对真实回波信号进行时频分析，完成瞬时频率的估计，发现在杂乱的炮口信号中，在特殊时段PCT仍能很好地完成瞬时频率估计，这一特殊时刻信号对精确预测炮口速度具有一定的贡献。但是PCT方法运算量大，算法需要进一步改进，提高效率。

 

是t0时刻附近时段内的多项式频率平移算子，t0起平移作用，gσ（t）是时宽为σ的高斯窗函数。

 

1.2 算法步骤

从图4中可以看出，WVD在分析单分量非线性频变信号也存在自交叉干扰，在利用极值提取脊线部分极值点在自交叉干扰时频特性曲线上，直接利用极值提取的脊线进行频率估计产生很大的误差。另外，在信噪比为-10 db情况下，信号原本的时频特征被淹没在噪声中，无法利用极值法提取脊线。

Step2：对脊线1用最小二乘法进行多项式拟合，拟合得到的参数作为变换核的初始参数；

Step3：根据初始变换核参数对信号进行PCT变换，从变换后的时频分布中提取脊线2；

Step4：对脊线2用最小二乘法进行多项式拟合，拟合得到的参数替换变换核参数；

丁香酚(E110640，99%，上海阿拉丁试剂有限公司)，酪蛋白酸钠(C8654，西格玛奥德里奇(上海)贸易有限公司)，海藻酸钠(180947，西格玛奥德里奇(上海)贸易有限公司)，磷酸氢二钠(5818118，上海麦克林生化科技有限公司)，吐温20、磷酸氢二钾、无水乙醇、氢氧化钠、柠檬酸、氯化钠等试剂均为分析纯。

前一观点把电机视为电路元件，对于定、转子上各有一套绕组的隐极电机，两绕组轴线之间的交角θm是随时间变化的。对磁链通过电流和电感表示，定、转子绕组之间的互感为Lsrcosθme，(θme=pθm是电角度)，由磁共能对θm求偏导数，可得电磁转矩

2 基于PCT的模拟回波信号瞬时频率估计

2.1 弹丸膛内运动的模拟回波信号

多年的天葬院生活，摒弃欲望，与世无争，置身于世态炎凉之外，恶魔的种子早已被包上了厚厚的壳，压在了心底的深处。然而，黑袍人的一句话，就让这颗种子生出了根。

  

图1 弹丸膛内运动的模拟回波信号

根据多普勒效应的基本原理，多普勒频率与运动径向速度成正比，为满足弹丸膛内运动速度的误差要求，频率估计的精度是首要保障。用一个关于时间t的N阶多项式去逼近弹丸膛内径向运动多普勒频率曲线，分别采用4、5、6、8阶多项式对模拟内弹道方程得到的多普勒频率关于时间的曲线进行拟合逼近，得出随着阶数的增加，拟合误差越来越小，计算量越来越大，当阶数为6阶能基本满足拟合精度要求，因此，在这里采用6阶多项式拟合逼近内弹道运动多普勒频率关于时间的曲线。如图2所示。

  

图2 不同阶数拟合逼近的频率变化曲线

2.2 模拟回波信号的瞬时频率估计

在实际获取信号的过程中，信号会受到外界干扰和系统噪声的影响。因此，本文在采用PCT对模拟信号进行瞬时频率估计之前，在模拟回波信号上不同信噪比的噪声干扰信号来模拟更加真实的模拟实际弹丸运动回波信号。

图3为模拟理想回波信号分别叠加不同信号比噪声的PCT时频分布图，图中噪声采用均匀分布的高斯噪声。

从图 3（a）～（d）可以看出，采用 PCT 时频分析方法，在时频分布图上瞬时频率的时频聚集性很好，从图3（e）可以看出，PCT分析理想模拟回波信号、不同信噪比的回波信号都利用极大值就能够很好地提取脊线，瞬时频率估计结果基本一致，这说明PCT方法具有较强的抗均匀分布高斯噪声的能力。下页图4为模拟理想回波信号分别叠加不同信号比噪声的WVD时频分布图，图中噪声仍采用均匀分布的高斯噪声。

太湖新城地下空间工程位于太湖新城核心启动区地下，项目总建筑面积30万m2，占地面积6.8万m2，地下共3层，结合地铁4号线支线溪霞路站沿中轴大道南北向布置。地下1层为商业及公共配套设施，设置大量下沉式广场，进行地面地下交通转换。地下2层和3层为公共停车空间，共20万m2，可容纳停车2 800辆。该项目计划于2013年12月开工建设。

Step1：对信号进行STFT变换，从变换后的时频分布图中，获取脊线1；

  

图3 模拟信号加不同噪声的PCT时频分布

总的来说，PCT和WVD方法在一定噪声下都能取得很好的时频聚集性，但是WVD分析非线性调频信号中存在自交叉干扰项，必须消除自交叉干扰，否则利用极值法直接提取WVD时频分布的脊线进行频率估计会产生很大的误差。对于PCT方法，在SNR=-10 dB噪声干扰情况下所得的误差均值为0.008 0，误差方差为1.075 9，最大误差为2.928 6。总的来说，采用PCT来进行分析弹丸膛内运动模拟信号的瞬时频率估计是有效的。

(3) 当逆作法施工采用盆式开挖法挖土时，最大位移出现在地下连续墙墙顶处，此时应在确保将墙顶位移控制在合理范围的前提下，基于首层土的开挖深度，确定盆式挖土预留土体的经济宽度。针对逆作法盆式开挖不同坡肩宽度下地下连续墙的位移变化情况，采用有限元法进行了研究。以往的工程经验表明，当软土区首层土盆式开挖深度取6.5 m、坡肩宽度取6.0 m时，即可实现将地下连续墙的变形控制在合理范围内的目的。

  

图4 模拟信号加不同噪声的WVD时频分布

3 基于PCT的实际弹丸信号瞬时频率估计

图5（a）为某次试验获取的实际回波信号，图5（b）为（a）的局部信号，可以明显看到弹丸有效回波信号上叠加了缓变信号，需要采用一定的预处理方法消

  

图5 实际获取的回波信号

从图中可以看出PCT时频分布时频聚集性优于WVD。对PCT和WVD时频分布结果利用极大值进行脊线提取，瞬时频率估计结果如图7所示。

  

图6 实际信号经过预处理得到的时频分布

真实回波信号图5（a）经过预处理利用PCT、WVD 方法进行时频分析，得到如图 6（a）、（b）所示的时频分布图。

  

图7 PCT时频分布的瞬时频率估计

从图7（a）可以看出，由于信号自交叉干扰，WVD时频分布图用极大值法抽取的脊线除了信号本身还有干扰项的脊线，无法直接利用极大值的结果进行瞬时频率估计，必须消除干扰项的脊线。PCT时频分布结果用极大值提取脊线，在出炮口处由于气体泄漏等因素造成信号质量下降，部分极大值点偏离了信号本身的时频特征，剔除粗大点后实现的频率估计结果如图7（b），PCT在实际信号中也能较好地完成瞬时频率估计。

4 结论

由参数组定义的多项式频率旋转算子，

参考文献：

由于采用RS-485的硬件协议，上位机为主机，而单片机作为下位机也就是从机，从机不会主动向主机发送信息，在接收到主机的信息后作出响应，因此从机在空闲状态时处于接收状态。

［1］张小兵.枪炮内弹道学［M］.北京：北京理工大学出版社，2014.

［2］陈征祥，施力琳.基于LabVIEW的微波干涉仪数据处理系统的开发［J］.科技资讯，2007，5（24）：6-7.

［3］于洋.电磁炮膛内弹丸测速系统设计［J］.四川兵工学报，2011，32（11）：43-46.

［4］王德田，彭其先，刘俊，等.激光干涉测速技术在内弹道弹丸速度测量中的应用研究［J］.高压物理学报，2011，25（2）：133-137.

［5］张玉成，张江波，严文荣，等.基于弹丸膛内速度微波测量的发射药燃烧规律［J］.火炸药学报，2010，33（4）：74-77.

［6］程诚，宋盛义，关永超，等.多普勒雷达测量电磁轨道炮内弹道速度［J］.强激光与粒子束，2014，26（4）：83-86.

［7］魏媛媛，姚金杰.基于小波去噪和改进极值法的动目标速度测量［J］.微计算机信息，2010，26（19）：217-218，200.

［8］刘丹，张丕状，马春燕.一种精确判定多普勒信号极值点位置的方法 ［J］. 激光与光电子学进展，2015，52（7）：53-57.

［9］刘丹，张丕状，马春燕.基于最小二乘法和STFT的弹丸回波信号估计分析 ［J］. 电子设计工程，2015，23（17）：131-134.

［10］姚金杰，韩焱.时频分析在变加速运动目标速度测量的应用［J］.火力与指挥控制，2011，36（5）：133-135，139.

［11］陈曦，许建中.Wigner高阶谱在弹丸膛内运动分析中的应用［J］.南京理工大学学报，2012，36（3）：442-447.

［12］杨健，陈曦，李学连.基于多项式Wigner-Ville分布的膛内弹丸回波信号处理［J］.四川兵工学报，2012，33（8）：112-114，123.

结合图8和表5数据可知,受绝热加速量热仪炉体材质影响,在样品反应剧烈阶段,炉体的追踪能力有限,无法做到百分百绝热。但是,补偿后的动力学参数更接近表1文献值,补偿前偏差过大,这反映了通过对反应过程的动态特性补偿能很大程度上使系统满足式(2),使得在动力学参数求解过程中误差减小。

［13］JIANG L，LI L，ZHAO G Q，et al.Instantaneous frequency estimation of nonlinear frequency-modulated signals under strong noise environment［J］.Circuits，Systems，and Signal Processing，2016，35：3734.

［14］LU W L，XIE J W，WANG H M，et al.Signal component extraction method based on polynomial chirp Fourier transform ［J］.Acta Physica Sinica ，2016，65（08）：9-18.

［15］XU L J，YANG Y X，YU S D.Analysis of moving source characteristics using polynomial chirplet transform［J］.Journal of the Acoustical Society of America ，2015，137（4）：EL320-6.

中国人习惯于将大多数蔬菜炒熟后食用，而在西方饮食中，蔬菜大多以生食为主。由于受西方生食蔬菜习俗时尚的影响，近年来，生食饮食与自制酵素（自制酶）等饮食概念与习惯，受到不少人追捧与效仿。不少喜欢这么吃的人都崇尚生吃蔬菜，认为这样可以最有效地保证维生素的摄入。与之持相反意见的人则认为从安全卫生的角度来说，应当将蔬菜烹熟了吃。

［16］杨扬.参数化时频分析理论、方法及其在工程信号分析中的应用［D］.上海：上海交通大学，2013.