滤波器广泛地应用于无线通信和数字信号处 理技术领域， 其特性由滤波 器结构和滤波系数决定。 在无线通信系统和数字信号处理领域中， 滤波器应 用的场景在很多时候是时变的， 也就是说一组滤波器系数组往往不能满足所 有场景的需求。 于是， 人们研究出了自适应的滤波器， 即滤波器的系数组是 实时计算并调整的。 这样， 自适应滤波器可以灵活地应对各种情况。 然而， 自适应滤波器在实际应用中也有其不可忽视的 缺点， 自适应滤波器的计算量 往往比较大， 这对设备的省电性能是非常不利的， 同^还要占用额外的计算 资源和存储资源， 增加了硬件成本。

实际上， 较为常用的一种折中的方法是， 预存储凡组不同的滤波器系数 组， 根据不同的场景选择相应的一组系数组， 这样既降低了计算量又能满足 不同场景的需求。 应用这种方法所需要解决两个主要问题是： ί) 预存储哪凡 组系数； 预存储系数与各种场景之间的映射关系。

在实际选择过程中， 如果采用比较直观的遍历式检测的方法进行选 择， 其运算量是非常巨大的。 以 LTE系统中的频域信道估计为例， 在给定滤波器 结构的情况下， 最优的滤波器系数取决于当前信道的^延扩展� �信噪比， 也 即， 其应用场景由时延扩展和信噪比共同决定， 如果将所有可能的时延扩展 和信噪比分别进行量化： i) 按照时延扩展对应的 RMS Delay (均方根时延） 量化为 50:50: (表示以 50ns为起始点， 50ns为步长进行量化）， 则 1000 ns 的时延扩展范围量化后得到一共 21个值。 ii) 信噪比 SNR量化为 - 5:2„5： (表 示以- 5dB为起始点， 2.5为步长进行量化）， 贝 ύί Β的信噪比的范围量化后 得到一共 15个值。由任一时延扩展的量化参数和任一信� ��比的量化参数组成 的组合， 可以定义一个量化场景， 这样， 一共可以得到 21 X 15= 315个量化 场景。 这 315个量化场景对应着最优的 315组滤波系数组。 如果在实际应用 中预存储 8组滤波器系数组， 就需要从 3】5组中选择出 8组滤波器系数组来 覆盖所有的场景， 那么可能的组合共有 ₅ =2, 19811236238094e ^{+ 15} 种， 将超过 200 万亿种情况。 对于如此巨大的数字， 如果采用直观的遍历式检验的方式 进行选择， 巨量运算将耗费大量时间， 在实际应^中几乎是不可行的。 发明内容

本发明所要解决的技术问题是提供一种应^于� �变系统的预存储滤波器 系数组的选择方法和装置， 能够以较低的滤波性能损失来覆盖所有的应用 场 景， 并能够极大地降低选择预存储滤波器系数组时 的运算量， 提高选择效率。

为解决上述技术问题， 本发明提供方案如下：

一种时变系统中预存储的滤波器系数组的选择 方法， 包括：

步骤 A, 对所述 ^变系统中定义^变场景的每个参数进行量化， 获得由 量化参数组所定义的量化场景， 每个量化参数组包括有所述每个参数量化后 得到的量化参数；

步骤 B, 根据所获得的量化场景和预定的滤波器结构， 确定每一量化场 景对应的滤波器系数组， 以及采用具有该滤波器系数组的滤波器对所有 量化 场景分别进行滤波所引起的滤波性能损失；

步骤 C, 根据所述滤波性能损失， 选择出数量大于预存储的滤波器系数 组的第一组数的待选择的量化场景， 将该待选择的量化场景对应的滤波器系 数组作为待选择的滤波器系数组；

步骤 D， 通过遍历方式， 从待选择的滤波器系数组中选择出针对所有量 化场景具有最佳滤波性能的、 组数等于所述第一组数的滤波器系数组， 作为 所述预存储的滤波器系数组。

优选地， 上述选择方法中， 所述步骤 D之后还包括：

步骤 E， 针对每一量化场景， 通过遍历方式， 获得采用所述预存储的滤 波器系数组中的每一滤波器系数组对应的滤波 器， 对该量化场景分别进行滤 波所能够获得的最小滤波性能损失， 进而建立该最小性能损失对应的滤波器 系数组与该量化场景之间的映射关系。 优选地， 上述选择方法中， 所述步骤 A具体包括：

根据所述时变系统中定义时变场景的每个参数 的取值范围， 对每个参数 分别进行量化， 得到量化后的量化参数；

根据量化后的量化参数， 确定由量化参数组所定义的量化场景， 其中， 每个量化参数组唯一对应于一个量化场景。

优选地， 上述选择方法中， 所述步骤 B具体包括：

对于预定的滤波器结构， 根据每个量化场景对应的量化参数组， 确定每 个量化场景对应的滤波器系数组；

根据每个量化场景对应的滤波器系数组， 确定每个量化场景对应的滤波 器；

针对每一量化场景， 确定采用任一量化场景对应的滤波器对该每一 量化 场景进行滤波， 相对于采^该每一量化场景对应的滤波器对该� �一量化场景 进行滤波所引起的滤波性能损失。

优选地， 上述选择方法中，

所述确定采用任一量化场景对应的滤波器对该 每一量化场景进行滤波， 相对于采用该每一量化场景对应的滤波器对该 每一量化场景进行滤波所引起 的滤波性能损失， 具体是通过滤波仿真或根据滤波器系数生成过 程通过理论 计算得到的。

优选地， 上述选择方法中， 所述步骤 C具体包括：

生成一滤波性能损失的损失矩阵， 所述损失矩阵的元素 ^表示第 i个量 化场景对应的滤波器对第 个量化场景进行滤波，相对于第 /个量化场景对应 的滤波器对第 个量化场景进行滤波所引起的性能损失；

对选择基础进行一次以上的选择处理， 直至选择基础中符合预定的性能 损失容忍门限要求的元素的数目为零， 其中， 每次选择处理具体包括： 从损 失矩阵中剔除列号集合所涉及的所有列上的元 素， 以及剔除与所述列号集合 的值相等的所有行上的元素， 得到本次选择处理的选择基础； 在本次选择处 理的选择基础上， 选择出具有最多的、 符合预定的性能损失容忍门限要求的 元素的行， 将该行的行号增加到行号集合中， 同时将该行上符合所述性能损 失容忍门限要求的所有元素的列号增加到列号 集合中； 在所述行号集合中行号的数量大于预存储的滤 波器系数组的组数时， 确 定所述行号集合中所有行号所对应的量化场景 的集合， 进而根据所述量化场 景的集合中所有量化场景的量化参数组， 确定对应的滤波器系数组， 得到一 组待选择的滤波器系数组。

优选地， 上述选择方法中，

在所述行号集合中行号的数量小于或等于预存 储的滤波器系数组的组数

B寸， 减小所述预定的性能损失容忍门限后， 返回所述对选择基础进行一次以 上的选择处理， 直至选择基础中符合预定的性能损失容忍门限 要求的元素的 数目为零的步骤。

优选地， 上述选择方法中， 所述步骤 D具体包括：

根据预存储的滤波器系数组的组数， 丛所述待选择滤波器系数组选择出 预存储的滤波器系数组的所有可能的组合；

针对每一组合， 通过遍历方式， 确定该每一组合中的每一滤波器系数组 对应的滤波器对每一量化场景分别进行滤波所 能够获得的最小性能损失， 计 算该每一组合针对所有量化场景所能获得的最 小性能损失的和值；

选择出最小的所述和值所对应的组合， 作为预存储的滤波器系数组。 本发明还提供了一种时变系统中预存储的滤波 器系数组的选择装置， 包 场景量化器，用于对所述时变系统中定义时变 场景的每个参数进行量化， 获得由量化参数组所定义的量化场景， 每个量化参数组包括有所述每个参数 量化后得到的量化参数；

性能损失生成器， 用于根据所获得的量化场景和预定的滤波器结 构， 确 定每一量化场景对应的滤波器系数组， 以及采用具有该滤波器系数组的滤波 器对所有量化场景分别进行滤波所引起的滤波 性能损失；

滤波器系数组初始选择器， ^于根据所述滤波性能损失， 选择出数量大 于预存储的滤波器系数组的第一组数的待选择 的量化场景， 将该待选择的量 化场景对应的滤波器系数组作为待选择的滤波 器系数组；

滤波器系数组二次选择器， 用于通过遍历方式， 从待选择的滤波器系数 组中选择出针对所有量化场景具有最佳滤波性 能的、 组数等于所述第一组数 的滤波器系数组， 作为所述预存储的滤波器系数组。

优选地， 上述的选择装置中， 还包括：

映射表生成器， 用于针对每一量化场景， 通过遍历方式， 获得采用所述 预存储的滤波器系数组中的每一滤波器系数组 对应的滤波器， 对该量化场景 分别进行滤波所能够获得的最小滤波性能损失 ， 进而建立该最小性能损失对 应的滤波器系数组与该量化场景之间的映射关 系。

优选地， 上述的选择装置中，

所述场景量化器具体用于：

根据所述时变系统中定义时变场景的每个参数 的取值范围， 对每个参数 分别进行量化， 得到量化后的量化参数；

根据量化后的量化参数， 确定由量化参数组所定义的量化场景， 其中， 每个量化参数组唯一对应于一个量化场景。

优选地， 上述的选择装置中，

所述性能损失生成器具体用于：

对于预定的滤波器结构， 根据每个量化场景对应的量化参数组， 确定每 个量化场景对应的滤波器系数组；

根据每个量化场景对应的滤波器系数组， 确定每个量化场景对应的滤波 器；

针对每一量化场景， 确定采 ^任一量化场景对应的滤波器对该每一量化 场景进行滤波， 相对于采用该每一量化场景对应的滤波器对该 每一量化场景 进行滤波所引起的滤波性能损失。

优选地， 上述的选择装置中，

优选地， 所述性能损失生成器进一步用于通过滤波仿真 或根据滤波器系 数生成过程通过理论什算， 得到所述滤波性能损失。

优选地， 上述的选择装置中，

所述滤波器系数组初始选择器具体用于：

生成一滤波性能损失的损失矩阵， 所述损失矩阵的元素 ^表示第 i个量 化场景对应的滤波器对第 个量化场景进行滤波，相对于第_/个量化场景� ��应 的滤波器对第 _/个量化场景进行滤波所引起的性能损失； 对选择基础进行一次以上的选择处理， 直至选择基础中符合预定的性能 损失容忍门限要求的元素的数目为零， 其中， 每次选择处理具体包括： 丛损 失矩阵中剔除列号集合所涉及的所有列上的元 素， 以及剔除与所述列号集合 的值相等的所有行上的元素， 得到本次选择处理的选择基础； 在本次选择处 理的选择基础上， 选择出具有最多的、 符合预定的性能损失容忍门限要求的 元素的行， 将该行的行号增加到行号集合中， 同时将该行上符合所述性能损 失容忍门限要求的所有元素的列号增加到列号 集合中；

在所述行号集合中行号的数量大于预存储的滤 波器系数组的组数时， 确 定所述行号集合中所有行号所对应的量化场景 的集合， 进而根据所述量化场 景的集合中所有量化场景的量化参数组， 确定对应的滤波器系数组， 得到一 组待选择的滤波器系数组。

优选地， 上述的选择装置中，

所述滤波器系数组初始选择器进一步用于在所 述行号集合中行号的数量 小于或等于预存储的滤波器系数组的组数时， 减小所述预定的性能损失容忍 门限后， 返回所述对选择基础进行一次以上的选择处理 ， 直至选择基础中符 合预定的性能损失容忍门限要求的元素的数目 为零的步骤。

优选地， 上述的选择装置中，

所述滤波器系数组二次选择器具体用于：

根据预存储的滤波器系数组的组数， 从所述待选择滤波器系数组选择出 预存储的滤波器系数组的所有可能的组合；

针对每一组合， 通过遍历方式， 确定该每一组合中的每一滤波器系数组 对应的滤波器对每一量化场景分别进行滤波所 能够获得的最小性能损失， 计 算该每一组合针对所有量化场景所能获得的最 小性能损失的和值；

选择出最小的所述和值所对应的组合， 作为预存储的滤波器系数组。 从以上所述可以看出， 本发明提供的预存储滤波器系数组的选择方法 和 装置， 通过两次选择的方式， 首先选择出一组待选择的滤波器系数组， 进而 利用遍历方式， 从待选择的滤波器系数组中选择出预存储的滤 波器系数组， 由于待选择的滤波器系数组的数量要小于量化 场景的数量， 从而可能的预存 储的滤波器系数组的组合也就大大减少， 因此本发明实施例可以极大地降低 遍历运算量， 提高预存储的滤波器系数组的选择效率。 并且， 本发明实施例 在首次选择 B寸根据滤波性能损失是否满足容忍门限来选� �出待选择的滤波器 系数组， 再基于该待选择的滤波器系数组选择出预存储 的滤波器系数组， JA 而能够较低的性能损失来覆盖所有的量化场景 ， 并且能够给出预存储滤波器 系数组和所有量化场景之间的映射关系。 附图说明

图 1为本发明实施例所述的预存储滤波器系数组� �选择方法的流程图； 图 2为本发明实施例中引用的通用的滤波器结构� �示意图；

图 3 为本发明实施例所述预存储滤波器系数组的选 择装置的结构示意 图。 具体实施方式

本发明提供了一种预存储滤波器系数组的选择 方法， 通过两次筛选的方 式进行预存储滤波器系数组选择， 在完成相同任务的情况下， 极大地降低了 运算量， 并能够给出 M组滤波器系数和所有场景之间的映射关系， 以较低的 性能损失来覆盖所有的时变场景。

为不失一般性和便于后文描述， 首先进行以下定义或说明：

1、 假设最终保留的预存储的滤波器系数组的数目 为 M。

2、 假设时变场景由由 ^，… 共 P个参数定义， 最优的滤波器系数将 根据这 P个参数生成。 即这 P个参数既唯一地定义了一个特定的时变场景� � 也唯一地定义了对应于这一个时变场景的一组 最优滤波器系数。

3、如何根据尸个参数生成滤波器系数组可以� �用现有技术的方法， 为节 约篇幅， 本文中不再赘述。

请参照图 1， 本发明实施例所述的应用于时变系统的预存储 滤波器系数 组的选择方法， 具体包括以下步骤：

步骤 11， 对时变系统中定义时变场景的每个参数进行量 化， 获得由量化 参数组所定义的量化场景， 每个量化参数组包括有所述每个参数量化后得 到 的量化参数。 步骤 12, 根据所获得的量化场景和预定的滤波器结构， 确定每一量化场 景对应的滤波器系数组， 以及采用具有该滤波器系数组的滤波器对所有 量化 场景分别进行滤波所引起的滤波性能损失。

步骤】 3， 根据所述滤波性能损失， 选择出数量大于预存储的滤波器系数 组的第一组数的待选择的量化场景， 将该待选择的量化场景对应的滤波器系 数组作为待选择的滤波器系数组。

步骤 14， 通过遍历方式， 从待选择的滤波器系数组中选择出针对所有量 化场景具有最佳滤波性能的、 组数等于所述第一组数的滤波器系数组， 作为 所述预存储的滤波器系数组。

通过以上步骤， 本实施例采^两次筛选的方式， 最终确定预存储的滤波 器系数组， 其中第一次筛选是从所有量化场景对应的滤波 器系数组中， 筛选 出数量大于所述第一组数 （假设为 M) 的待选择的滤波器系数组； 然后， 通 过第二次筛选， 从待选择的滤波器系数组中， 筛选出数量为所述第一组数的 预存储的滤波器系数组。 可以看出， 第二次筛选， 是基于第一次筛选所得到 的待选择的滤波器系数组进行， 其可能的组合共有 个。 而现有技术从 所有量化场景对应的滤波器系数组中选择， 其可能的组合共有 = 个， Ν 表示所有量化场景的数量。 由于待选择的滤波器系数组的数量要小于量化 场 景的数量， 从而可能的预存储的滤波器系数组的组合也就 大大减少， 因此本 实施倒可以极大地降低步骤 14中的遍历运算量，提高预存储的滤波器系数� �� 的选择效率。

在上述步骤 14之后， 本实施倒还可以进一步包括以下步骤：

步骤 15， 针对每一量化场景， 通过遍历方式， 获得采用所述预存储的滤 波器系数组中的每一滤波器系数组对应的滤波 器， 对该量化场景分别进行滤 波所能够获得的最小滤波性能损失， 进而建立该最小性能损失对应的滤波器 系数组与该量化场景之间的映射关系。

也就是说， 对于一个量化场景， 该量化场景所映射的滤波器系数组， 是 在预存储的滤波器系数组中对于该量化场景具 有最小滤波性能损失的滤波器 系数组。这样， 通过以上步骤 15， 本实施例给出了 Μ组滤波器系数组和所有 量化场景之间的映射关系， 进而可以保存这些映射关系， 以便于在实际滤波 中据此选择合适的滤波器系数组。

为了帮助理解上述步骤， 以下对上述步骤做进一步的说明。

上述步骤 11具体包括：

步骤 Π 1 ,根据所述时变系统中定义时变场景的每个参� �的取值范围，对 每个参数分别进行量化， 得到量化后的量化参数；

步骤 112， 根据量化后的量化参数， 确定由量化参数组所定义的量化场 景， 其中， 每个量化参数组唯一对应于一个量化场景。

本实施例中， 可以根据每个量化场景的量化参数组， 生成该量化场景对 应的滤波器系数组， 确定该量化场景对应的滤波器。

对于时变场景来说， 通过 P个参数可以定义无限多个时变场景， 然而在 实际的预存储滤波器系数组选择过程中， 只能对有限个场景进行处理。 本实 施例通过以上步骤 11 , 对可能的时变场景进行量化， 得到量化场景， 具体来 说， 是对用于定义时变场景的各个参数进行量化。 即对于参数 ft， _e [l,P]， 将 其可能值从小到大依次量化为 ， . ·ν ， 为其最小可能值， v¾为 其最大可能值， ^表示参数 ft量化后能够得到的量化值的个数。 每一个量化 场景由 _e [l,Ay, l,U唯一确定。 籍此， 可将所有可能场景 量化为 A^I¾ 个场景。 例如， 以 UTE系统中的频域信道估计为例，在给定滤波器 结构的情况下， 最优的滤波器系数取决于当前信道的时延扩展 和信噪比， 也即， 其时变场景 由时延扩展和信噪比这两个参数共同决定。 每个参数都有一个取值范围， 例 如， 小区半径决定了时延扩展的取值范围， 基于现有技术中的小区半径， 时 延扩展的取值范围大概在 0〜5000ns (该时延扩展对应的 RMS delay大约在 0〜1000ns)。 再例如， 对于当前信道的信噪比， 则可以根据系统正常工作时 要求的最低信噪比， 作为信噪比的下限 （最低信噪比可以根据相关技术标准 或设备参数查找的）； 再根据没有噪声的理想情况下（可以考虑电路 板的热噪 声） 信道的信噪比， 来确定最高信噪比， 作为信噪比的上限。 对于其它各种 参数， 本领域技术人员都能够根据影响该参数的因素 或相关技术标准， 确定 该参数的一个可能的取值范围， 此处不再一一赘述。

对于 ΕΤΈ系统中的频域信道估 % 如果： 0 按照时延扩展 RMS Delay 量化为 50:50： (表示以 50ns为起始点， 50ns为步长进行量化）， 则 0〜1000 ns 的时延扩展范围量化后得到一共 21个值。 ii) 信噪比 SNR量化为 5:2.5 : (表 示以- 5dB为起始点， 2.5为步长进行量化）， 贝^ 5dB〜25dB的信噪比的范 量 化后得到一共 15个值。每个量化场景由任一时延扩展的量化� ��数和任一信噪 比的量化参数组成的组合所定义， 这样可以定义得到 2】 X 15= 315个量化场 景， 每个量化场景唯一对应于一组量化参数组。

上述步骤 12具体包括：

步骤 121， 对于预定的滤波器结构， 根据每个量化场景对应的量化参数 组， 确定每个量化场景对应的滤波器系数组。 这里， 根据量化参数组， 确定 对应的滤波器系数组， 可以按照现有技术的根据参数确定对应的滤波 器系数 的方式进行确定。

步骤 i22， 根据每个量化场景对应的滤波器系数组， 确定每个量化场景 对应的滤波器。

步骤 123， 针对每一量化场景， 确定采用任一量化场景对应的滤波器对 该每一量化场景进行滤波， 相对于采用该每一量化场景对应的滤波器对该 每 一量化场景进行滤波所引起的滤波性能损失。 这里， 滤波性能损失可以通过 滤波饬真或根据滤波器系数生成过程通过理论 计算得到的， 以下分别举例说 明：

假设 , , _e [l，N]表示用第 个场景对应的滤波器对第 个场景进行滤波所 引起的滤波性能损失。 所谓通过仿真得到性能损失矩阵， 可按照如下方式得 到：

1 ) 使 ffi第 个量化场景对应于滤波器系数对第 个量化场景进行滤波， 得到滤波输出的信噪比为 SV¾ _; ，，./ _e [l,N]； 使用第 个量化场景对应于滤波器 系数对第 J个量化场景进行滤波， 得到滤波输出的信噪比为 SW^, [1,N]。

2 ) 如果信噪比用 dB表示， 则 5Λ¾ Λ [1 ] ; 如果信噪比用 线性值表示， 贝 υ / _; [ , ΛΠ

所谓根据滤波器系数生成过程通过理论计算得 到， 请参考图 2 用的滤波器结构， 其中包括多个延迟单元 22和多个加法器 21， 用于对输入 信号 X(n)进行滤波处理， 得到滤波后的输出信号 Y(n)。 图 2 中， _ai 、 a ₂ ...a _p 为反馈环路的系数， bj、 b ₂ ...b _Q 为前向通路的系数， 这里， 中的下标 P表

2

：滤波器中横向结构网络的阶数， 即横向结构网络由 P个前向通路 m 要 P个延迟单元 22 , b _Q 中的下标 Q表示滤波器中反馈网络的阶数， 即反馈网 络由 Q个反馈环路构成， 需要 Q个延迟单元 22。 图 2表 ^ '以n下^的T滤^波 i^ffi枉， 可以^以下公式表示：

p

Yin) - Z h _p . X{n― p) - aj(n― q )

引起的误差

E h— h

E \ |/i -- W ⁽ · Χ

E 2 R

：:: 1 - 2 Re(w ⁱ⁰

0[1, Λ'·]

匕式中： /ί表示滤波样点 Υ(η) i理想值， 表示滤波样点的实际滤波 输出的值， 即图 2中的 Y(n); X表示输入信号 X(n); W«表示用于第 个量 化场景对应的滤波器系数组； (" ₂ ,--' ,^,'- ]， 这里表示为一个行 向量； R 表示第 j个量化场景下 ffi于 ii 输出 Y ( n ) 的输入样点的自相关 矩 阵 ， 即

I [X(n), X(n - 1), - - - X(n -- P), Y(n - 1), Υ(η - 2), - - - Υ(η -- 0) Ri 表示第 /个量化场景下向量 I与 h的互相关向量， 是一个列向：

2 ) i

上述步骤 13具体包括：

步骤 131， 生成一滤波性能损失的损失矩阵 L, 失矩阵的元素 /^ f，./ _e [l，N]表示位于损失矩阵 L第 行第. /列的元素，其含义为第 个量化场 景对应的滤波器对第 ./个量化场景进行滤波，相对于第. /个量化场景对应的滤 波器对第. /个量化场景进行滤波所引起的性能损失。

步骤 132， 对选择基础进行一次以上的选择处理， 直至选择基础中符合 预定的性能损失容忍门限要求的元素的数目为 零， 其中， 每次选择处理具体 包括：

从损失矩阵中剔除列号集合所涉及的所有列上 的元素， 以及剔除与所述 列号集合的值相等的所有行上的元素， 得到本次选择处理的选择基础；

在本次选择处理的选择基础上， 选择出具有最多的、 符合预定的性能损 失容忍门限要求的元素的行， 将该行的行号增加到行号集合中， 同 B寸将该行 上符合所述性能损失容忍门限要求的所有元素 的列号增加到列号集合中。

步骤 133， 在所述行号集合中行号的数量大于预存储的滤 波器系数组的 组数 （M ) 确定所述行号集合中所有行号所对应的量化场 景的集合， 进 而根据所述量化场景的集合中所有量化场景的 量化参数组， 确定对应的滤波 器系数组， 得到一组待选择的滤波器系数组。

步骤 134， 在所述行号集合中行号的数量小于或等于预存 储的滤波器系 数组的组数 （M ) 时， 减小所述预定的性能损失容忍门限后， 返回步骤 132。

步骤 13根据性能损失矩阵 L和容忍门限 来确定出初始的 Μ' ( Μ'通常 会大于 Μ ) 组滤波器系数组。 这里的容忍门限 ρ是指， 用第 /个量化场景对 应的滤波器对第_/个量化场景进行滤波时， 所能容忍的的最大的滤波性能损 失。该门限是 ^来保证通过初始选择得到的 M'组滤波器系数在覆盖所有量化 场景^， 所引入的最大性能损失不会超过 。 这里， 上述步骤 131〜134具体 可以通过以下的矩阵操作来实现：

步骤 I , 定义一个标志矩阵 F , 其大小为 Ν Χ Ν， 标志矩阵中的元素 0,./ _e [l, N]表示位于标志矩阵 F第 i行第 列的元素。 根据损失矩阵 L和容 忍门限 p对标志矩阵中的元素进行初始置位， 当 小于容忍门限 P时 , 即用 第 个量化场景对应的滤波器对第 个量化场景进行滤波所引起的性能损失小 于容忍门限 P时， 将标志矩阵中的相应位置的 . 置 1 , 反之则置 0。 具体可用 如下伪代码表示： for i：:: 1： N

forj = 1 : N

./;, 0

if < p

Λ i

end

end

步骤 2 ,选择初始的 M'组的待选择的滤波器系数组，步骤 2具悻又可以包 括以下步骤：

步骤 A，判断标志矩阵 F中是否存在不为 0的元素，若是，则进入步骤 B; 否则进入步骤 D;

步骤 B, 从标志矩阵 F中选择出具有最大的元素和值 的一行， 将该行的 行号记录到行号集合 中， 并将该行上所有不为 0 的元素的列号 记录到列号集合^ Afo^ 中， 然后进入步骤 C。

步骤 C,根据列号集合记录的列号，对标志矩阵 F中的元素进行置位处理， 所述置位处理包括： 将标志矩阵 F中列号等于列号集合记录的列号的列上的 所有元素置 0， 以及将标志矩阵 F中行号等于列号集合记录的列号的行上的 所有元素置 0, 然后返回步骤 A。

步骤 D，判断行号集合中记录的行号的数量是否大� �预存储的滤波器系数 组的组数 （M)， 若是， 则进入步骤 E， 否则， 减小所述容忍门限 ^后返回步 骤 。

这里， 减小容忍门限 p， 可以使得行号集合中记录的行号的数量增加。 由 于行号集合中记录的行号的数量即是待选择的 滤波器系数组的数量 M' ， M' 越大， 后续的运算量也就越大， 即容忍门限 越大， 滤波性能要求就越低， 而后续计算量越小； 容忍门限 p越小， 滤波性能要求就越高， 而后续 算量 越大。 因此， 可以综合考虑系统的滤波性能需求和后续的运 算量来设定。

步骤 E， 确定行号集合中记录的行号对应的量化场景的 集合， 进而根据所 述量化场景的集合中所有量化场景的量化参数 组，确定对应的滤波器系数组， 得到一组待选择的滤波器系数组。 例如， 在行号集合中记录的行号记录的行 号包括 3、 5、 27、 122等行号时， 则这些行号对应的量化场景包括第 3、 5、 27、 122 等量化场景， 于是可以根据这些量化场景的量化参数组， 确定这些 量化场景对应的滤波器系数组， 得到待选择的滤波器系数组。

上述步骤 2具体可用如下伪代码表示：

m ― 0;

N

while ： ,. > 0

= ^ar g ^χ (β:) Ι ^? ' "' ^Λ

记录系数^ 覆盖的场景为 Afapp 即 Fd^ , :)中不为 0的场景序号 inildxArrayim ) ^ f _ifia /7save the index f j -------- 0, i G idxMapped , j -------- 1,… Λ' ^? — ;

f] _; = 0, / = L ' * - N, j G idxA4apped;

m ― m + 1

end

M ^f - m

通过滤波器系数组初始选择器， M'个待选择的滤波器系数组的序号记录 在数组 中。 如果 _Μ '小于 _M 可以将容忍门限换成一个更小的值， 重 复步骤 1和 2， 直至 M'大于或等于 M

上述步骤 14具体包括：

步骤 141， 根据预存储的滤波器系数组的组数， 丛所述待选择滤波器系数 组选择出预存储的滤波器系数组的所有可能的 组合。

步骤 142, 针对每一组合， 通过遍历方式， 确定该每一组合中的每一滤波 器系数组对应的滤波器对每一量化场景分别进 行滤波所能够获得的最小性能 损失， 计算该每一组合针对所有量化场景所能获得的 最小性能损失的和值。

步骤 143， 选择出最小的所述和值所对应的组合， 作为预存储的滤波器系 数组。

上述步骤 14是从 M'个滤波器系数组中选出 M组用于预存储， 因此共有 ί/ ( 种可能的组合方式。 这里通过对所有可能的组合进行遍历式检验， 找 到最佳的 Μ个用于预存储的滤波器系数组。具体过程可� ��如下伪代码表示如 C possible combinations of M

// candidates from initial candidates set

// C is with size U x M

for / =^ 1 : U

C ^f = C(/,:);/7 get ith possible combination in C

E( = 0

for j：::： 1： N

:E(0 E(0 -f _miR

end

end

i _best = arg min[E(i)]， / ^ 1,2，， . . U

idxArrayBest ^ C(U;

这样， 通过以上步骤 14进行滤波器系数组的二次筛选， 最佳的 M个预存 储的滤波器系数组就被选择出来了。

上述歩骤 15用以生成量化场景与预存储的滤波器系数组� ��的滤波器系数 组之间的映射关系， 该映射关系可以根据性能损失矩阵 L确定下来， 也即在 预存储的滤波器系数组中， 找到对于每一个量化场景的滤波性能损失最小 的 一组的滤波器系数组， 作为该量化场景所映射的滤波器系数组。 这里通过映 射表 T的形式表示出来。 具体可用如下伪代码表示如下：

for j - 1： N

/c _m ；.,™ argmin(' _! ) ,Α' e i bcArrayBest

τ(./Χ

end

以上介绍了本实施例所述的时变系统中预存储 滤波器系数组的选择方法， 基于以上选择方法， 本发明实施例还提供了一种时变系统中预存储 滤波器系 数组的选择装置， 请参照图 3 , 所述选择装置具体包括：

场景量化器， ^于对所述时变系统中定义时变场景的每个参� �进行量化， 获得由量化参数组所定义的量化场景， 每个量化参数组包括有所述每个参数 量化后得到的量化参数；

性能损失生成器， 用于根据所获得的量化场景和预定的滤波器结 构， 确定 每一量化场景对应的滤波器系数组， 以及采用具有该滤波器系数组的滤波器 对所有量化场景分别进行滤波所引起的滤波性 能损失； 滤波器系数组初始选择器， ffi于根据所述滤波性能损失， 选择出数量大于 预存储的滤波器系数组的第一组数的待选择的 量化场景， 将该待选择的量化 场景对应的滤波器系数组作为待选择的滤波器 系数组；

滤波器系数组二次选择器， 用于通过遍历方式， 从待选择的滤波器系数组 中选择出针对所有量化场景具有最佳滤波性能 的、 组数等于所述第一组数的 滤波器系数组， 作为所述预存储的滤波器系数组。

优选地， 所述选择装置还可以包括： 映射表生成器， 用于针对每一量化场 景， 通过遍历方式， 获得采用所述预存储的滤波器系数组中的每一 滤波器系 数组对应的滤波器， 对该量化场景分别进行滤波所能够获得的最小 滤波性能 损失， 进而建立该最小性能损失对应的滤波器系数组 与该量化场景之间的映 射关系。

优选地， 所述场景量化器具体用于：

根据所述^变系统中定义 ^变场景的每个参数的取值范围，对每个参数� � 别进行量化， 得到量化后的量化参数；

根据量化后的量化参数， 确定由量化参数组所定义的量化场景， 其中， 每 个量化参数组唯一对应于一个量化场景。 对于预定的滤波器结构， 根据每个量化场景对应的量化参数组， 确定每个 量化场景对应的滤波器系数组；

根据每个量化场景对应的滤波器系数组， 确定每个量化场景对应的滤波 ；

针对每一量化场景，确定采用任一量化场景对 应的滤波器对该每一量化场 景进行滤波， 相对于采用该每一量化场景对应的滤波器对该 每一量化场景进 行滤波所引起的滤波性能损失。

优选地，所述性能损失生成器进一歩用于通过 滤波仿真或根据滤波器系数 生成过程通过理论计算， 得到所述滤波性能损失。

优选地， 所述滤波器系数组初始选择器具体用于：

生成一滤波性能损失的损失矩阵， 所述损失矩阵的元素^表示第 i个量 化场景对应的滤波器对第 个量化场景进行滤波，相对于第„/个量化场� �对应 的滤波器对第 个量化场景进行滤波所引起的性能损失；

对选择基础迸行一次以上的选择处理，直至选 择基础中符合预定的性能损 失容忍门限要求的元素的数目为零， 其中， 每次选择处理具体包括： 从损失 矩阵中剔除列号集合所涉及的所有列上的元素 ， 以及剔除与所述列号集合的 值相等的所有行上的元素， 得到本次选择处理的选择基础； 在本次选择处理 的选择基础上， 选择出具有最多的、 符合预定的性能损失容忍门限要求的元 素的行， 将该行的行号增加到行号集合中， 同时将该行上符合所述性能损失 容忍门限要求的所有元素的列号增加到列号集 合中；

在所述行号集合中行号的数量大于预存储的滤 波器系数组的组数时，确定 所述行号集合中所有行号所对应的量化场景的 集合， 进而根据所述量化场景 的集合中所有量化场景的量化参数组， 确定对应的滤波器系数组， 得到一组 待选择的滤波器系数组。

优选地，所述滤波器系数组初始选择器进一步 用于在所述行号集合中行号 的数量小于或等于预存储的滤波器系数组的组 数时， 减小所述预定的性能损 失容忍门限后， 返回所述对选择基础进行一次以上的选择处理 ， 直至选择基 础中符合预定的性能损失容忍门限要求的元素 的数目为零的步骤。

优选地， 所述滤波器系数组二次选择器具体用于：

根据预存储的滤波器系数组的组数，从所述待 选择滤波器系数组选择出预 存储的滤波器系数组的所有可能的组合；

针对每一组合， 通过遍历方式， 确定该每一组合中的每一滤波器系数组对 应的滤波器对每一量化场景分别进行滤波所能 够获得的最小性能损失， if算 该每一组合针对所有量化场景所能获得的最小 性能损失的和值；

选择出最小的所述和值所对应的组合， 作为预存储的滤波器系数组。 以上所述仅是本发明的实施方式， 应当指出， 对于本技术领域的普通技术 人员来说， 在不脱离本发明原理的前提下， 还可以作出若干改进和润饰， 这 些改进和润饰也应视为本发明的保护范围。