专利名称：一种基于小波变换和支持向量机的涡结构识别方法
技术领域：
本发明涉及一种基于小波变换和支持向量机的涡结构识别方法，可用于分析大气中高 速飞行的导弹和飞机的光学头罩与来流之间形成的气动光学效应，为图像校正、控制湍流 提供基础。
背景技术：
高速成像制导导弹和飞机等在大气层内高速飞行时，光学头罩周围会形成高速复杂流 场，致使目标光波通过高速流场时光波相位发生改变，目标图像产生偏移、模糊、抖动和 能量衰减，这种现象称为气动光学效应。气动光学效应的研究重点之一是畸变波前同湍流 涡动力学特性之间的关系。按照湍流的涡旋学说，湍流的脉动与混合主要由大尺度的涡结 构造成，具有高度的各向异性，且其流动、演变具有重复性和可预测性；小尺度涡几乎是 各向同性的，而且不同流动的小尺度涡有许多共性。大尺度涡造成了大部分的密度脉动， 在光学传输中造成了大部分的相位脉动，是影响气动光学传输效应的主要因素。因此要按 照不同的涡结构尺度建立精确的光学传输模型，就必须进行涡结构的识别。涡结构识别的关键在于特征提取和分类器设计。传统的涡结构识别方法，在特征提取 方面通常是依据小波变换具有多分辨分析特性，直接进行小波分解，利用小波系数矩阵作 为特征量使得特征冗余多，计算量大在分类器设计方面传统方法通常直接根据经验确定 小波分解水平来区分大、小尺度涡结构，识别率低。 发明内容本发明的目的是克服传统的涡结构识别方法的不足，提出一种基于小波变换和支持 向量机的涡结构识别方法。本发明的技术方案是把机载光学窗口处湍流折射率场经小波变换后得到各频带的系 数矩阵，以其统计量作为涡结构的特征量进行识别。提取的特征有小波系数矩阵的均值、 方差、能量、信息熵、对角线惯性矩、绝对均匀性和聚类重要性。再利用支持向量机方法 建立分类器，根据结构风险最小原则使分类面不仅能将涡结构正确分幵，而且分类间隔最 大。
具体步骤如下(1)利用小波变换对描述涡结构的折射率进行小波分解，得到不同尺度上的小波系数矩阵。
(2) 对小波系数矩阵进行特征提取，小波系数矩阵的统计量作为涡结构的特征值。(3) 根据所提取特征值，利用支持向量机的方法建立最佳分类判别函数，其建立准则 采用利用风险最小化原则，精确识别大小涡结构。(4) 利用步骤(3)建立的最佳分类判别函数，重复步骤(1) 、 (2)对涡结构样本 进行涡结构识别。(5) 利用小波反变换进行大、小尺度涡重建。
本发明的基本原理是导弹和飞机周围的绕流流场中存在不同尺度的涡结构，而具有 "数字显微镜"之称的小波4H斤具有良好吋、频特性以及多尺度多分辨分析特性，是研究精 细结构的理想工具，可以很好地反映涡结构。传统的方法直接利用小波变换的系数作为涡 结构的特征量，这样造成特征冗余大，存储和计算资源的浪费。表示涡结构的矩阵是一个 随机变量，用统计量来描述更为科学，因此本发明将折射率场经小波变换后得到各频带的 系数矩阵的统计量作为涡结构特征量，并以此进行识别，提取的特征有小波系数矩阵的均 值、方差、能量、信息熵、对角线惯性矩、绝对均匀性和聚类重要性，这些涡结构的特征 就构成了涡结构模式的样本，简称涡结构样本，记为X,;涡结构的类别记为 ^e{-l,+l}，/-l，2J， /为样本个数。已知涡结构类型的样本包括大尺度涡结构样本和小尺度涡结构样本，大尺度涡结构样本由流场的折射率求得，流场的折射率是通过利用大涡模 拟的方法求解Navier-Stokes方程(简写N-S方程)，得到密度场数据再经过Gladstone-Dale 变换(简称GD变换)求得，小尺度涡结构的样本通过N-S方程的精确数值解求得，N-S 方程的精确数值解中既包含大尺度涡模式又包含小尺度涡模式，从N-S方程的精确数值解 中减去大尺度涡模式便得到小尺度涡模式的密度场描述，在经过GD变换便可得到小尺度 涡结构。但由于建立N-S方程和求解精确数值运算量巨大，需要多节点计算机同时并行计 算多天，因此本发明提出基于小波变换和支持向量机方法识别方法，大大提高计算效率。 利用N-S方程精确解出的涡结构作为已知类别的样本建立最佳分类判别函数，其原理是因 为涡结构样本属于二维线性不可分问题，通过核函数进行非线性变换，将低维不可分的情 况转化为高维可分空间，即通过函数映射^ —把对利用特征x的二维线性不可分 问题转化为高维线性可分问题，因此定义高维线性判别函数/ = >5^|>/*2: + 6]， w、6为 权值向量及偏置，对特征x识别问题转换为对特征z的识别问题。如图3所示，表示了涡结 构分类问题的基本原理，支持向量机的主要思想是利用核函数将输入向量x映射到一个高 维特征空间z，并在该空间内构造一个最优分类超平面来逼近分类函数尸=>/.2+6 = 0。因 为原；c特征空间为非线性空间，映射到z线性空间，线性空间的线性函数表达式为^ ^ + 6 = 0，最优分类超平面的构造最终归结为在原空间上求解一个凸二次规划问题。恩3 中O和口分别代表两类涡结构样本，尸与^之间的间隔为i^1， F/与尸2之间的距离为分类间隔j^i，若使分类间隔j^f最大，即要^IMI2最小。因为样本不完全线性可分，可求解4w/w+;^^,的最小值，其中/i为惩罚参数，；z越大表示对错误分类的惩罚越大;《为分类2 ,'=1误差。采用拉格朗日乘子法求解这个具有线性约束的二次规划问题，目标函数表示为111严—{￡(>1^,2)-^卜||2+/^《-^>^,(>/'2, + 6)-1]}，进一步转化得w-t义j^，由此构成的最佳分类判别函数为<formula>formula see original document page 5</formula>虽然将特征x映射高维特征空间，维数有所提高，但也只考虑在高维特征空间的点积运算，故不必明确知道映射函数，构造内积核函数为〖Ov力-z^z。就可得到非线性支持向量机的 判决函数<formula>formula see original document page 5</formula>,Nnsv表示支持向量的个数，从最佳分类判别函 数可以看出，z为样本序号，/表示所有样本的个数，即所有特征值A所对应的样本。而经过推导只有不为零的A才对分类有贡献，它所对应的样本就是支持向量xnsv.在计算分类函数时，不需计算/个样本，只需计算NNsv样本即可，也就是支持向量:cNsv决定了分类超平面，而其他的样本(非支持向量)并不影响分类。因此支持向量机分类判决函数是一种用支 持向量集合对所有样本数据的压缩表示，降低了维数，这是其他机器学习技术(例如神经网 络技术)所不具备的优点，最优分类面就是使分类面不但能将两类正确分开(训练错误率为0)，而且使分类间隔最大，即分类器具有最好的泛化能力。 本发明与现有技术相比的优点在于(1) 在特征提取方面，把折射率场数据经小波变换后系数矩阵的统计量作为涡结构 特征量，包括信息熵、绝对均匀性和聚类重要性等七个特征，比传统的直接利用小波系数 矩阵作为特征量，不仅剔出了大量的冗余信息，降低了计算量。(2) 在分类器设计方面，提出了基于支持向量机的分类方法，利用结构风险最小原 则使分类面不仅能将涡结构正确分开，训练错误率为0，而且使分类间隔最大，比传统的经 验风险最小原则分类方法更具有泛化能力。本发明能更加准确地表征和区分湍流涡结构， 为导弹、飞机等机载光学设备的光学效应精确建模奠定了基础。
图l为本发明的涡结构识别方法流程图。图2为塔式小波变换图，D为低频部分，H为高频，下标l， 2表示图像的第一级和第 二级分解。塔式小波变换是不断地对每一级的低频子带进行滤波和采样处理，保留每级中 的高频子带部分。图3为最优分类超平面。O和口分别代表两类涡结构样本，F为分类超平面，A和F2分 别为各类中距离分类超平面最近的涡结构样本且平行于分类超平面的平面，它们之间的距 离即记为分类间隔。
具体实施方式
如图1 2所示，本发明的具体方法如下1利用Daubechines小波变换对折射率进行小波分解，得到不同尺度上的小波系数矩阵。 原始数据为光学窗口附近的折射率场数据。小波变换方式采用塔式信号分解方式，如 图2所示，D为低频部分，H为高频，下标1， 2表示图像的第一级和第二级分解。塔式 小波变换是不断地对每一级的低频子带进行滤波和下采样处理，保留每级中的高频子带部 分。小波函数采用Daubechies小波，小波分解水平W根据具体问题进行选择，与所要求的 细节有关，iV越大分解水平越高，得到的细节越多。假定对于某个分解水平M(0《MSAO ， W为最高的分解水平，分解从水平W到W-M，由综合考虑本实施例分解的折射率数据时AT , M都取为3，即分解水平分别为3， 2， 1， 0。小波分解过程为根据Daubechies小波函数和尺度函数确定( ^， -^^"^0， -Z2SwS0， (6"J是无限序列的截断，这两个序列可通过现有的小波分析方法得到。A, A相当于滤波窗口的大小，A = A=2W-i。对已知的折射率场数据 [Cb,J,"W-l,…W-M;-Z3("《"V0 ;—A'"M《A")，其中A为分解的水平，Z3，A分别表示数据的起始和终止的行数；4，丄6分别表示数据的起始和终止的列数。本实施例 = =丄,)="w = 128 。
令 《)=《+" /2， 《)=(Z( +f ))/2 ;Zf =《+1)/2，《> 一I2+ZT))/2。小波分解的系数表示为<formula>formula see original document page 6</formula>
2对小波系数矩阵进行特征提取。
因为大尺度涡是会产生较大的折射率变化，是气动光学效应的主要原因，故只计算低 频系数^ . 的特征值，包括均值、方差、能量、信息熵、对角线惯性矩、绝对均匀性和聚 类重要性作为涡结构的特征。1) 均值2) 方差3) 能量<formula>formula see original document page 7</formula>(2)<formula>formula see original document page 7</formula>(3)<formula>formula see original document page 7</formula>4) 信息熵<formula>formula see original document page 7</formula>5) 对角线的惯性矩<formula>formula see original document page 7</formula>6) 绝对均匀性<formula>formula see original document page 7</formula>7) 聚类重要性<formula>formula see original document page 7</formula>(8)3根据所提取的特征和己知涡结构样本类别，利用支持向量机的方法建立最佳分类判别 函数。因为折射率表示的涡结构:c是非线性的，要是使其线性可分，就必须对;c作非线性变 换z--(x),使其成为线性可分，则在映射后z空间中求取权值向量w及偏置6,使其满足线性可分。涡结构识别中输入A映射成z,，输出模式表示为^，分别表示大尺度涡结构和小尺度涡结构k,少,U-l,2J，^"^,.eFl,+lK ;， J,分别表示第/个样本所对应的特征和is样本类别，/为样本总数，表示z,.属于^维实数空间，y,样本类别有两类分别表示大涡结构和小涡结构。。而非线性变换的映射函数并不需要知道显式形式，只需计算核函数。定义大涡和小涡分别表示为若 >^=1， Or.z,) + 6 2l 若 (t/.z,) + "-1以上二式合写为袭、)+ 6]-12 0 (9) 最优分界面方程为F=wr.z+6-0 (10)由此可见求解最优分界面F应在式(9)的约束下，即使两类涡结构样本满足xo/^+6)^u-i,2j，这样分类间隔为j^j，要使分类间隔最大，即要^Mi2最小，因为样本不完全线性可分，可求解^v/m;+A力《的最小值，采用拉格朗日乘子法求解这个具有线性约束的二次规划问题，目标函数表示为maxmiii雖6,;i)4llMf +《《-土程w、 } ,，""1 " 2 ,=1 (ii)s厶 A^;i^:0其中A为惩罚参数，Ae[O，l]，越大表示对错误分类的惩罚越大，通过|^ = 0求得；《为 分类误差；义,为拉格朗日乘子，y"为;ii的最大值。由此得到 3丄 '沩 tT可得w-S^V,; (12)至此可求得w,6,;i，式中A可能是(DA=o; (Do< A<>" ; A= a。只有后两者所对应的A对w有贡献，对应的样本成为支持向量，计为xww，也就是对最优分界面、判别函数有贡献，所对应的学习方法称之为支持向量机。在支持向量中满足③所对应的2,称为边界支持向量，实际上是错分的训练样本点;满足②所对应的4称为标准支持向量。
将式(12)带到式(11)，进一步处理得111严[必(义)=2^-22义;A 7,力(；.勺)],'='户' (13)W, t為乂 -0， 》0 由式(13)可以求得各4，最优分界面方程为6= 0 (14)由此构成的最佳分类判别函数为涡结构是一种线性不可分的情况，无法知道映射函数z--(;c)的显式形式，可以利用内 积核函数的非线性映射将输入向量映射到高维空间，支持向量机能在此高维空间中给出最 佳分类超平面。构造内积核函数为. ;c)=z『.z (16)它表示特征空间中的一个内积。最佳分类判别函数进一步表示为—〈/=w"〖5] A乃尤"，+6] )/-尽管通过非线性函数将样本数据映射到具有高维的特征空间，并在特征空间中构造最 优分类超平面，但在求解最优化问题和计算判别函数时并不需要显式计算该非线性函数，而只需计算核函数，从而避免特征空间维数灾难问题。核函数的选择必须满足Merce条件。 常用的核函数有线性核函数 /:(x,,j:) = v^ (18)多项式核函数〖Oc,，"-(v" + l)" (19)径向基核函数^(《，;0-exp(20)径向基核函数式(20)具有较好的分类性能，比多项式核函ir和线性核函数性能要好，因此本发明选用径向机核函数式(20)构造分类器。4利用步骤3建立的最佳分类判别函数，重复步骤l、 2对涡结构样本进行涡结构识别。 输入为未知涡结构样本的特征量;c，经过最佳分类判别函数/进行判别，输出为l的为大涡模式，输出为-l的为小涡模式。5利用小波反变换进行大、小尺度涡重建。在重构中需对低频带的小波系数矩阵进行特征量的提取，大尺度涡结构的重构也只是 利用低频分量进行重构。高频分量重构后得到小尺度涡结构。根据Daubiches小波预先存储的有限序列{/\ }和,这两个序列可以通过查找有关小波分析的书得到。c'i+1; ,m《+1; , 分别为重构的大尺度涡结构和小尺度涡结构，其他具体的符号定义同步骤l。l A+l，rt,m 一 2/，m-2I/C*;/,_/ , ，m P (21)本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
权利要求
1、一种基于小波变换和支持向量机的涡结构识别方法，其特征在于包括下列步骤(1)利用小波变换对描述涡结构的折射率进行小波分解，得到不同尺度上的小波系数矩阵。(2)对小波系数矩阵进行特征提取，小波系数矩阵的统计量作为涡结构的特征值。(3)根据所提取特征值，利用支持向量机的方法建立最佳分类判别函数，其建立准则采用利用风险最小化原则。(4)利用步骤(3)建立的最佳分类判别函数，重复步骤(1)、(2)对涡结构样本进行涡结构识别。(5)利用小波反变换进行大、小尺度涡重建。
2、 根据权利要求书l所述的一种基于小波变换和支持向量机的涡结构识别方法，其特 征在于所述步骤(1)中小波分解采用塔式分解。
3、 根据权利要求书l所述的一种基于小波变换和支持向量机的涡结构识别方法，其特 征在于所述步骤(2)特征提取的是小波系数矩阵的方差、能量、信息熵、对角线惯性矩、绝对均匀性和聚类重要性。
全文摘要
本发明涉及一种基于支持向量机(SVM，SupportVectorMachine)的涡结构识别方法，在特征提取方面，提取折射率场数据经小波变换后系数矩阵的统计量作为涡结构特征量，与传统的直接利用小波变换后系数矩阵相比，剔出了大量的冗余信息，降低了计算量。在分类器设计方面，提出了基于支持向量机的分类方法，利用结构风险最小原则使分类面不仅能将涡结构正确分开，而且使分类间隔最大，比传统的经验风险最小原则分类方法从原理上降低了误识率。本发明能更加准确地表征和区分湍流涡结构，为导弹、飞机等机载光学设备的光学效应精确建模奠定了基础。
文档编号G01M9/00GK101158617SQ20071017756
公开日2008年4月9日 申请日期2007年11月16日 优先权日2007年11月16日
发明者冯浩楠, 琳 吴, 房建成, 杨照华 申请人:北京航空航天大学