发明名称 采用遗传算法和非线性凸规划理论优化参数的滤波方法
摘要 本发明涉及移动通信领域,旨在提供一种采用遗传算法和非线性凸规划理论优化参数的滤波方法。相对于传统的数字无线中继放大设备应用技术,本发明是在将数字基带GSM传输信号送至数字GSM滤波器,进行选频滤波处理,输出滤波后GSM基带信号的过程中,根据遗传算法和非线性凸规划理论,预先优化计算好滤波器参数,使用配有以上GSM新型滤波器对GSM低倍率数字基带信号进行选频滤波处理。本发明不对参数对称性,线性相位特性进行严格要求,而是以一系列代价函数将其进行折中,放宽了参数优化时的受限程度,同时引入解决非线性规划的遗传算法,进一步提升了优化的效果,因此在设计结果上大大优于现有技术设计结果。
申请公布号 CN102664646A 申请公布日期 2012.09.12
申请号 CN201210146243.7 申请日期 2012.05.10
申请人 杭州畅鼎科技有限公司 发明人 蔡斌;梅其灵
分类号 H04B1/26(2006.01)I;H04B1/00(2006.01)I 主分类号 H04B1/26(2006.01)I
代理机构 杭州中成专利事务所有限公司 33212 代理人 金祺
主权项 1.一种采用遗传算法和非线性凸规划理论优化参数的滤波方法,其特征在于,包括以下步骤:(1)以模拟超外差混频结构为基础,对天线接收信号进行模拟下变频,得到模拟中频信号;由采样频率为f<sub>sam</sub>的时钟控制的模数转换器ADC对模拟中频信号进行数字化处理,得到数字中频信号;然后由数字DDC进行数字下变频和数字降采样处理,得到数字基带GSM传输信号;(2)将数字基带GSM传输信号送至数字GSM滤波器,进行选频滤波处理,输出滤波后GSM基带信号;(3)将滤波后GSM基带信号送至数字DUC进行上变频,经数模转换、模拟上变频,放大至射频,输送到发射天线;所述步骤(2)是通过采用遗传算法和非线性凸规划理论对GSM滤波器参数进行优化设计得以实现的,具体包括以下步骤:(A)根据不同系统的设计指标要求,对滤波器性能进行分配,将滤波器性能指标具体划分为幅度响应要求、相位响应要求和时延要求;)(B)对步骤(A)中的各个要求在数字频带[0,2π)内进行量化,得到理想滤波器频响函数H<sub>d</sub>(e<sup>jω</sup>),以其对应的幅度响应|H<sub>d</sub>(e<sup>jω</sup>)|、相位响应arg (H<sub>d</sub>(e<sup>jω</sup>))、群时延响应(相位变化率)τ(e<sup>jω</sup>)作为优化的目标;(C)设置向量变量h=[h(0),h(1)…h(N-1)]<sup>T</sup>表示实际可得滤波器的参数,N表示目标滤波器的阶数,T表示向量转置;则滤波器的频域响应表示为下式:<maths num="0001"><![CDATA[<math><mrow><mi>H</mi><mrow><mo>(</mo><msup><mi>e</mi><mi>j&omega;</mi></msup><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>=</mo><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>h</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>&times;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j&omega;n</mi></mrow></msup><mo>=</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow></mrow></math>]]></maths>其对应的幅度响应为|H(e<sup>jω</sup>)|,相位响应为arg (H(e<sup>jω</sup>)),群时延响应为<img file="FDA0000162633882.GIF" wi="300" he="138" />;(D)设理想低通滤波器通带为[0,ω1],阻带为[ω2,π],设置代价函数<maths num="0002"><![CDATA[<math><mrow><mi>&phi;</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mo>=</mo><mi>&Delta;</mi><mn>0</mn><mo>&times;</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mi>&omega;</mi><mn>1</mn></mrow></msubsup><msup><mrow><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mi>j&omega;</mi></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mi>d&omega;</mi><mo>+</mo><mi>&Delta;</mi><mn>1</mn><mo>&times;</mo><msubsup><mo>&Integral;</mo><mrow><mi>&omega;</mi><mn>2</mn></mrow><mi>&pi;</mi></msubsup><msup><mrow><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mi>j&omega;</mi></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mi>d&omega;</mi><mo>+</mo><mi>&Delta;</mi><mn>2</mn><mo>&times;</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mi>&omega;</mi><mn>1</mn></mrow></msubsup><msup><mrow><mo>|</mo><mfrac><mrow><mi>d</mi><mi>arg</mi><mrow><mo>(</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mi>d&omega;</mi></mfrac><mo>-</mo><mi>&tau;</mi><mrow><mo>(</mo><msup><mi>e</mi><mi>j&omega;</mi></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mi>d&omega;</mi></mrow></math>]]></maths>△0,△1和△2为比例因子,用于调节通带、阻带和群时延的优化程度;在寻找h时,使上述代价函数最小,以使设计滤波器的频率响应和理想滤波器的频域响应最大程度的接近;(E)根据数字计算机的计算特点,对上述连续代价函数进行离散化,具体为将频带[0,2π)均匀量化为K个样值点<img file="FDA0000162633884.GIF" wi="162" he="129" />,i=0,1,……,K-1;则步骤(D)中代价函数中的各项分别离散化为:<maths num="0003"><![CDATA[<math><mrow><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mi>&omega;</mi><mn>1</mn></mrow></msubsup><msup><mrow><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mi>j&omega;</mi></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mi>d&omega;</mi><mo>&RightArrow;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>k</mi><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>&times;</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi></mrow></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mo>&Integral;</mo><mrow><mi>&omega;</mi><mn>2</mn></mrow><mi>&pi;</mi></msubsup><msup><mrow><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mi>j&omega;</mi></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mi>d&omega;</mi><mo>&RightArrow;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi><mn>2</mn></mrow><mrow><mi>K</mi><mo>/</mo><mn>2</mn></mrow></munderover><msup><mrow><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>&times;</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi></mrow></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mi>&omega;</mi><mn>1</mn></mrow></msubsup><msup><mrow><mo>|</mo><mfrac><mrow><mi>d</mi><mi>arg</mi><mrow><mo>(</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mi>d&omega;</mi></mfrac><mo>-</mo><mi>&tau;</mi><mrow><mo>(</mo><msup><mi>e</mi><mi>j&omega;</mi></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mi>d&omega;</mi><mo>&RightArrow;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><mi>arg</mi><mrow><mo>(</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>arg</mi><mrow><mo>(</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>&tau;</mi><mrow><mo>(</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>&times;</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi></mrow></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>k</mi><mn>1</mn><mo>=</mo><mfrac><mrow><mi>&omega;</mi><mn>1</mn></mrow><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mo>&times;</mo><mi>K</mi></mrow></math>]]></maths>,<maths num="0007"><![CDATA[<math><mrow><mi>k</mi><mn>2</mn><mo>=</mo><mfrac><mrow><mi>&omega;</mi><mn>2</mn></mrow><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mo>&times;</mo><mi>K</mi></mrow></math>]]></maths>(F)根据步骤(E)中离散化的形式,应用遗传算法和非线性优化方法求解代价函数<maths num="0008"><![CDATA[<math><mfenced open="" close=""><mtable><mtr><mtd><mi>&phi;</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mo>=</mo><mi>&Delta;</mi><mn>0</mn><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>k</mi><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>&times;</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi></mrow></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>&Delta;</mi><mn>1</mn><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi><mn>2</mn></mrow><mrow><mi>K</mi><mo>/</mo><mn>2</mn></mrow></munderover><msup><mrow><mo>|</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>&times;</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi></mrow></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo></mtd></mtr><mtr><mtd><mi>&Delta;</mi><mn>2</mn><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><mi>arg</mi><mrow><mo>(</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>arg</mi><mrow><mo>(</mo><mi>F</mi><mrow><mo>(</mo><mi>h</mi><mo>,</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>&tau;</mi><mrow><mo>(</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>&times;</mo><mi>i</mi><mo>/</mo><mi>K</mi><mo>&times;</mo><mn>2</mn><mi>&pi;</mi></mrow></msup><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mtd></mtr></mtable></mfenced></math>]]></maths>的最优解,其结果即为设计所需的滤波器;(G)利用步骤(F)所得的滤波器,对步骤(1)中输出的数字基带GSM传输信号进行滤波处理,输出滤波后数字基带GSM传输信号。
地址 310011 浙江省杭州市拱墅区祥园路45号东5楼