疣有什么危害| neo什么意思| brush什么意思| alba手表是什么牌子| 摩羯座是什么星座| 北面属于什么档次| 粿是什么意思| 扁桃体看什么科室| 软肋是什么意思| cp感什么意思| 沙土地适合种什么农作物| 狗狗咳嗽吃什么药好得快| 女生下面长什么样| 始于初见止于终老是什么意思| 多此一举是什么生肖| 鹦鹉爱吃什么| 乙型肝炎表面抗体高是什么意思| 动次打次是什么意思| 血止不住是什么原因| 低密度脂蛋白偏低是什么意思| 心肌标志物是查什么的| 感染幽门螺杆菌吃什么药| 右手中指发麻是什么原因| 刘姥姥和贾府什么关系| 1994年出生属什么| 什么时间泡脚最好| 失恋什么意思| 这是什么字| 欣五行属什么| 梦见家里好多蛇是什么预兆| 幽门螺旋杆菌吃什么药| 闪购是什么意思| 左眼跳什么预兆| 小孩睡觉流鼻血是什么原因引起的| 乙醇是什么| cyl是什么意思| 做梦抓到很多鱼是什么征兆| 先入为主是什么意思| 诺贝尔为什么没有数学奖| 乌鱼子是什么意思| 杉菜是什么意思| 手术后为什么不让睡觉| 消融手术是什么意思| 吕洞宾属什么生肖| 有样学样是什么意思| 麒麟飞到北极会变成什么| 依靠是什么意思| 复方对乙酰氨基酚片是什么药| 滴虫病女性有什么症状| 屁股大什么原因| 高血糖吃什么水果最好| 嗓子有异物感吃什么药| b1是什么意思| 秋分节气的含义是什么| 湖南有什么好玩的| 白羊座前面是什么星座| 天然气是什么味道| 悬壶济世是什么意思| 胃糜烂吃什么药效果好| 医学ns是什么意思| 丙字五行属什么| 霾是什么意思| 高温中暑吃什么药| 丰富多腔的腔是什么意思| 花千骨什么时候上映的| 青钱柳有什么功效与作用| 指南针是什么时候发明的| 赶集什么意思| 霉菌性中耳炎用什么药| 女性胆囊炎有什么症状| 蒲公英有什么作用| 泥鳅不能和什么一起吃| 军长相当于地方什么官| 便秘是什么症状| 什么季节补钙效果最好| 全身皮肤瘙痒是什么原因| 红烧排骨用什么排骨比较好| 阿米替林片是治什么病的| 什么品牌奶粉最好| 七月八号是什么星座| 女性腰肌劳损吃什么药| 什么水果含糖量低| 肝早期硬化身体有什么症状| 说风就是雨什么意思| 体检吃早餐有什么影响| 非你莫属是什么意思| 5月3日是什么星座| 胃痛看什么科| 大浪淘沙下一句是什么| 荣耀是什么品牌| 兔死狗烹是什么生肖| 胸口闷痛挂什么科| 什么家庭不宜挂八骏图| 梦见自己头发白了是什么意思| 二甲双胍缓释片什么时候吃| 五戒十善是什么| 尿酸吃什么药最有效果| 为什么下巴经常长痘痘| 四级什么时候报名| 风向是什么意思| 脾虚生痰吃什么中成药| 腠理是什么意思| 复印病历需要什么证件| 胃痛胃胀什么原因引起的| 大脚趾发黑是什么原因| 蟑螂屎长什么样| 八月十号是什么星座| 不拘小节是什么意思| 属鸡和什么属相相克| 什么叫翡翠| 985是什么学校| 动物都有什么| 拉烂屎是什么原因| 黄喉是什么动物身上的| 暄字五行属什么| 冬瓜和什么不能一起吃| 姓陈取什么名字好听| 闻鸡起舞是什么意思| 头不由自主的轻微晃动是什么病| aldo是什么牌子| 脚腕筋疼是什么原因| 梦到和死人说话是什么意思| 可乐喝多了有什么危害| 右肺中叶少许纤维灶是什么意思| 巧克力的原料是什么| 总感觉自己有病是什么心理病| 哈达是什么意思| 1130是什么星座| 孩子不愿意吃饭是什么原因| 黑茶金花是什么菌| 护理专业是做什么的| 湿气重是什么意思| 西瓜坏了是什么样| 氨曲南是什么药| 转氨酶高是什么原因| 想法是什么意思| 半月板后角变性什么意思| 正值当年什么意思| 婴儿感冒吃什么药| 阴历六月十五是什么日子| 台甫是什么意思| 喝苦丁茶有什么好处| 蓦然回首什么意思| 石头五行属什么| 过路车是什么意思| 为什么辰不能见亥| 爷爷的爸爸叫什么| clarks是什么牌子| 老想喝水是什么原因| 老农民韩美丽结局是什么| 内科是什么| 防中暑喝什么水| 什么的眼光| 吃什么食物对眼睛好| 雷蒙欣氨麻美敏片是什么药| 三书六礼指的是什么| 半夜流鼻血是什么原因| 金枝玉叶什么生肖| 接站是什么意思| 道听途说是什么意思| 标本是什么意思| 黄瓜敷脸有什么作用与功效| 五年存活率是什么意思| 什么不及| 肝脏钙化灶是什么意思| 凌晨是什么时辰| 白子画什么时候爱上花千骨的| 狼吞虎咽的意思是什么| 2026是什么年| 吃南瓜子有什么好处| 胎盘粘连是什么原因造成的| 一月十二号是什么星座| 子宫肌瘤吃什么中药可以消除掉| 五蕴皆空是什么意思| 进字五行属什么| 丹参长什么样子图片| 冬至有什么忌讳| 手脚出汗是什么原因| 属鼠的是什么命| 心脏跳得快是什么原因| 春天什么花开| 隐是什么意思| 活化部分凝血活酶时间偏高是什么意思| 查处是什么意思| 好事将近是什么意思| md是什么学位| 黑今读什么| la帽子是什么牌子| 雪莲菌泡牛奶有什么功效| 被口什么感觉| 肌酸有什么用| 什么字寓意好| 乙肝表面抗原大于250是什么意思| 鸽子不能和什么一起吃| 澳门什么时候回归的| 雌激素是什么东西| 缓释是什么意思| taco是什么| 起伏不定是什么意思| 俗不可耐是什么意思| 山人是什么意思| 白色的猫是什么品种| 皮肤溃烂化脓用什么药| 老年人全身无力是什么原因| 鱼吃什么| 煮牛肉放什么料| 高铁为什么会晚点| 网盘是什么东西| 血氧饱和度低于90有什么危害| 死去活来是什么生肖| 排卵期出血是什么颜色| 手指疣初期什么样子| 什么水果最好吃| 开水烫伤用什么方法好的最快| 养儿防老下一句是什么| 学习是什么| 化妆棉是干什么用的| 2333是什么意思| 子宫内膜单纯性增生是什么意思| 人得猫癣用什么药| 子宫内膜9mm意味着什么| 长智齿是什么原因引起的| 预激综合征是什么病| 美容师都要学什么| 下九流指的是什么| 里急后重吃什么药最好| 治疗呼吸道感染用什么药最好| 甲状腺结节忌口什么| 剁椒能做什么菜| 秋葵有什么营养| 棉花糖是什么做的| hkc是什么牌子| 上钟什么意思| 985是什么| 脾胃虚弱吃什么药| 目瞪口呆是什么生肖| 消化不良吃什么食物好| 汉语拼音什么时候发明的| 12.18是什么星座| 固执什么意思| 脓毒症是什么引起的| 沉鱼落雁闭月羞花是什么意思| 什么症状| 猪和什么属相不合| 什么牌子助听器好| 桃花眼的女人什么命| 皮下脂肪瘤挂什么科| 菜场附近开什么店好| 验血肝功能看什么指标| 与时俱进是什么意思| 十月八号是什么星座| 宝宝支气管炎吃什么药| 回眸一笑百媚生什么意思| 痛风什么东西不可以吃| 8月10日是什么星座| 细菌感染发烧吃什么药| 头部ct挂什么科| 小孩为什么经常流鼻血| 为什么长不胖一直很瘦| apm是什么牌子| 中央委员是什么级别| 女生来大姨妈要注意什么| 蛋白糖是什么糖| 脑供血不足用什么药好| 百度

新闻中心

EEPW首页 > 模拟技术 > 设计应用 > 一文搞懂《扩展卡尔曼滤波融合IMU和ODOM》

国家信访局通报一起追责案例

作者: 时间:2025-08-03 来源:硬十 收藏
1、什么是
百度     纽约南区的联邦检察官杰弗里·伯曼23日在记者会上说,这种“无耻的大规模网络攻击”是“司法部起诉过的最大一起由国家发起的黑客行动”。

(Kalman filtering)是一种利用线性系统状态方程,通过系统输入输出观测数据,对系统状态进行最优估计的算法。由于观测数据中包括系统中的噪声和干扰的影响,所以最优估计也可看作是滤波过程。
就像移动靶在运动中,每一刻都有 “实际在哪”“跑多快”“往哪个方向拐” 这些真实状态(对应系统状态向量x)。但狙击手看不到这些 “真实值”,只能通过瞄准镜(传感器)得到带误差的观测(比如 “好像在 100 米处,可能偏左半米”)。

观测数据里的 “噪声和干扰”,就像瞄准镜前的树叶、雾气、树木遮挡

当树叶挡住视线时(观测数据突然不准),不会完全相信这一帧的观测,而是更多依赖 “基于目标运动规律的预判”(比如 “刚才目标在匀速跑,不可能突然跳 10 米,这帧数据肯定是被树叶挡了,忽略它”)。

当视线清晰时(观测可靠),滤波会用新观测修正之前的预判(比如 “果然在 105 米,和我想的差不多,稍微调一下瞄准点”)。
卡尔曼滤波就像一个经验丰富的狙击手:用目标的运动规律(系统模型)做预判,用带噪声的观测(传感器数据)做修正,最终在干扰中锁定真实位置。这个过程既不盲目相信单次观测(过滤噪声),也不固执依赖历史判断(结合新信息),从而实现 “最优估计”。
2、扩展卡尔曼滤波(EKF)和 卡尔曼滤波(KF)
卡尔曼最初提出的滤波理论只适用于线性系统,Bucy,Sunahara等人提出并研究了扩展卡尔曼滤波(Extended Kalman Filter,简称EKF),将卡尔曼滤波理论进一步应用到非线性领域。EKF的基本思想是将非线性系统线性化,然后进行卡尔曼滤波,因此EKF是一种次优滤波。其后,多种二阶广义卡尔曼滤波方法的提出及应用进一步提高了卡尔曼滤波对非线性系统的估计性能。二阶滤波方法考虑了Taylor级数展开的二次项,因此减少了由于线性化所引起的估计误差,但大大增加了运算量,因此在实际中反而没有一阶EKF应用广泛。
EKF全称ExtendedKalmanFilter,即扩展卡尔曼滤波器,一种高效率的递归滤波器(自回归滤波器)。

扩展卡尔曼滤波(EKF)是卡尔曼滤波(KF)的 “升级版”,专门用来解决卡尔曼滤波搞不定的场景。我们可以继续用 “狙击枪瞄准移动靶” 的比喻来理解它们的关系:

1. 卡尔曼滤波(KF):只能对付 “规矩的目标”

卡尔曼滤波的核心前提是:目标的运动规律是 “线性的”
比如移动靶在匀速直线运动(速度不变、方向不变),或者匀加速运动(每秒多跑 1m/s)—— 这些运动可以用简单的线性公式描述(比如 “下一秒位置 = 现在位置 + 速度 × 时间”)。
这时候,狙击手只要用 “匀速 / 匀加速模型” 就能精准预判,卡尔曼滤波的公式也能直接套用,计算简单又高效。

2. 扩展卡尔曼滤波(EKF):能搞定 “调皮的目标”

但现实中,很多目标的运动是 “非线性的”:

  • 比如移动靶突然转弯、绕圈,甚至做 S 形走位(运动方向和速度的关系不是简单的线性公式);

  • 或者传感器的观测方式是非线性的(比如用雷达测角度,角度和位置的关系是三角函数,不是直线关系)。

这时候,卡尔曼滤波就失效了 —— 因为它的公式只能处理线性关系,面对转弯的目标,用 “匀速直线” 模型预判会偏差越来越大。
EKF 的作用就是 “把非线性问题掰弯了变成线性的”
它会在 “当前估计的状态点” 附近,用一条直线(切线)近似目标的非线性运动轨迹(就像用很多小线段拼接成曲线)。这样一来,就能继续套用卡尔曼滤波的线性公式来计算,只不过每次都要重新算这条 “切线”(专业上叫 “线性化”,用雅克比矩阵实现)。

总结:关系就像 “基础款计算器” 和 “科学计算器”

  • 卡尔曼滤波是 “基础款”:只能算加减乘除(线性问题),简单直接,速度快。
  • 扩展卡尔曼滤波是 “科学计算器”:能处理三角函数、指数等复杂运算(非线性问题),通过 “近似线性化” 扩展了适用范围,但计算更复杂一点。

  • 两者的核心逻辑(“预判 + 修正” 的滤波思想)完全一致,只是 EKF 能应对更贴近现实的非线性场景。

3、EKF融合IMU和ODOM的计算过程

咱们可以把 EKF 融合 IMU 和 ODOM 的过程,想象成一个 “猜位置” 的游戏,你需要结合两种 “线索” 来不断修正自己的猜测,让结果越来越准。

IMU(Inertial Measurement Unit)和 ODOM(Odometry)在机器人领域以及其他涉及运动感知和定位的系统中是非常重要的概念,以下是它们的详细介绍:

IMU(惯性测量单元)

  • 定义

    IMU 是一种能够测量物体三轴加速度和三轴角速度的组合式传感器,有的 IMU 还集成了磁力计 ,可以测量物体的磁场信息,辅助确定物体的方向。它是一个小型化、高度集成的系统,能够实时提供物体的运动状态信息。
  • 工作原理

    • 加速度计

      :基于牛顿第二定律,通过检测内部质量块在加速度作用下产生的力,来测量加速度。比如,当物体加速时,质量块会相对于传感器内部的固定结构产生位移,通过检测这种位移并将其转换为电信号,就可以得到加速度的大小和方向。
    • 陀螺仪

      :利用科里奥利力原理工作。当陀螺仪的壳体发生旋转时,内部的旋转质量块会受到科里奥利力的作用,从而产生微小的振动或偏移,通过检测这些变化就可以测量出角速度。


ODOM(里程计)

  • 定义

    里程计是一种用于测量物体在运动过程中移动距离和方向的装置或算法,通过对运动过程中的一些参数进行累计和计算,来确定物体的位置变化。
  • 工作原理

    • 轮式里程计

      :常见于轮式机器人,通过安装在轮子上的编码器来记录轮子的转动圈数。结合轮子的半径,就可以计算出机器人在平面上移动的距离。同时,通过一些额外的传感器(如安装在差速驱动机器人两个轮子上的编码器,对比两个轮子的转动圈数差异),可以推算出机器人的转向角度和方向变化。


先给两个 “线索源” 做个简单画像:

  • IMU:像你身上的 “运动传感器”,能实时感觉到自己在怎么动 —— 比如走得快还是慢、有没有转弯(用角速度测转弯,用加速度测加速)。但它有个毛病:猜多了会 “跑偏”(比如站着不动,它可能会慢慢觉得你在挪位置,这叫 “漂移”)。

  • ODOM:像你手里的 “路标导航”,每隔一会儿告诉你当前的位置和速度(比如轮式里程计通过轮子转了多少圈算位置,或者视觉里程计通过摄像头看环境算位置)。它相对稳定,不会轻易跑偏,但更新没那么快(比如轮子转半圈才出一次数据),而且偶尔会 “骗你”(比如轮子打滑时,它算的位置就不准了)。


EKF 做的事,就是用 IMU 的 “实时感觉” 先“猜”,再用 ODOM 的 “路标信息” 来修正,最后得到一个又快又准的结果。具体分两步:

第一步:用 IMU “猜”(预测阶段)

假设你现在要走路去超市,EKF 先根据 IMU 的信息,实时猜你每一秒的位置和速度。比如:

  • 上一秒你在 “小区门口”(上一时刻的位置),速度是 “1 米 / 秒向前”(上一时刻的速度)。

  • 这一秒 IMU 告诉你:“你没转弯(角速度为 0),也没加速(加速度为 0)”(IMU 的测量值)。

  • 那 EKF 就会猜:“这一秒你应该在小区门口往前 1 米的地方,速度还是 1 米 / 秒”(预测的当前位置和速度)。


但 EKF 知道 IMU 会 “跑偏”,所以猜的时候会留个 “误差范围”:比如 “位置大概在 1 米左右,误差可能有 0.5 米”(这个误差会随着 IMU 用得越久慢慢变大,就像猜多了越来越没谱)。

第二步:用 ODOM “修正”(更新阶段)

走了一会儿,ODOM 这个 “路标” 说话了:“根据轮子转的圈数,你现在应该在‘离小区门口 1.2 米’的地方”(ODOM 的测量值)。

这时候 EKF 要做两件事:

  • 对比 “猜的位置(1 米)” 和 “ODOM 说的位置(1.2 米)”,算出差了 0.2 米。

  • 看看这两个线索谁更可信,再决定怎么修正。

比如:

  • 如果地面很平,轮子没打滑,ODOM 很靠谱(权重高),那就多听 ODOM 的,把猜测从 1 米调成 1.15 米(靠近 1.2 米)。

  • 如果地面坑坑洼洼,轮子可能打滑了,ODOM 不太准(权重低),就少听 ODOM 的,调成 1.05 米(靠近原来的猜测)。

修正后,误差范围也会缩小(比如从 0.5 米变成 0.2 米),因为两个线索交叉验证了,结果更靠谱了。

循环往复,越猜越准

之后就重复这两步:

  • 没 ODOM 数据时,用 IMU 实时猜(预测),误差慢慢变大;

  • 有 ODOM 数据时,用它修正猜测(更新),误差缩小。

这样一来,既利用了 IMU “反应快” 的优点(实时跟紧运动),又用 ODOM “不跑偏” 的优点(修正累积误差),最后得到的位置、速度就会又快又准。
简单说,EKF 就像个 “聪明的裁判”:

  • 让 IMU 当 “实时报信员”,随时喊 “现在往哪动、动多快”;
  • 让 ODOM 当 “定期校准员”,每隔一会儿喊 “实际位置在这”;
  • 裁判根据两者的靠谱程度,综合出一个最准的结果。


4、晦涩难懂的预测公式

这个方程是扩展卡尔曼滤波器(EKF)的核心预测方程,确实比较抽象。我来一步步解释,让它变得更容易理解。

状态预测方程的直观理解

这个方程描述了如何根据上一时刻的状态估计当前的控制输入,来预测当前时刻的系统状态

我们可以用一个生活中的例子来类比:

上一时刻的状态:就像你现在知道自己的位置和速度

控制输入:就像你知道自己正在向前走,或者正在转弯
预测的状态:基于现在的位置、速度和动作,预测出下一刻你会在哪里

数学符号详细解释

让我们拆解这个方程:

符号含义:

  • x?:带帽子的 x 表示 "估计值"(不是真实值,因为有误差)

  • k|k-1:竖线后面的数字表示 "基于哪个时刻的信息"

    • x?_k|k-1

      表示:基于 k-1 时刻的信息,对 k 时刻状态的预测
    • x?_k-1|k-1

      表示:基于 k-1 时刻的信息,对 k-1 时刻状态的估计(这是上一轮更新得到的结果)
  • f():这是一个函数,表示状态如何从一个时刻转移到下一个时刻

  • u_k:控制输入,例如 IMU 测量的加速度和角速度

  • w_k:过程噪声,表示系统的不确定性(例如,你的动作可能不完全按照计划执行)

更具体的例子:移动机器人

假设我们有一个在平面上移动的机器人,它的状态包括:

  • 位置:x, y

  • 朝向:θ(角度)
  • 速度:v_x, v_y(x 和 y 方向的速度)


总结

状态预测方程x?_k|k-1 = f(x?_k-1|k-1, u_k, w_k)的核心思想是:

  • 基于上一时刻的状态估计

  • 加上控制输入的影响
  • 考虑系统的不确定性
  • 预测当前时刻的状

通过这个方程,EKF 能够利用 IMU 等高频率的控制输入数据,不断预测系统的状态变化。
5、过程噪声协方差矩阵
过程噪声协方差矩阵

1. 矩阵结构与状态对应

矩阵的行和列严格对应 EKF 状态向量的顺序,通常是:
[x, y, z, roll, pitch, yaw, vx, vy, vz, vroll, vpitch, vyaw, ax, ay, az]

每个元素 process_noise_covariance[i][j] 的含义:

  •  i=j(对角线元素):表示单个状态量的预测噪声强度(不确定性大小)。
  •  i≠j(非对角线元素):表示两个状态量之间的噪声相关性(通常设为 0,假设各状态噪声独立)。

2. 逐行 / 列解读(以你的配置为例)

下面结合常见机器人运动场景(如差速底盘、机械臂),解释对角线元素(非对角线多为 0,代表噪声独立)的物理意义:

状态维度
索引
对角线值
含义(预测时的不确定性)
位置
x
0.05
预测 x 方向位置时,假设的噪声强度(比如机器人在 x 轴运动,模型误差导致的位置偏差 “大概这么大”)

y
0.05
y 方向位置的预测噪声

z
0.06
z 方向位置(高度)的预测噪声(若机器人是地面移动设备,z 变化少,噪声可更小;若涉及飞行则需调大)
姿态(角)
roll
0.03
横滚角(roll)的预测噪声(比如模型假设 “姿态不变”,但实际因振动等导致的误差)

pitch
0.03
俯仰角(pitch)的预测噪声

yaw
0.1
偏航角(yaw)的预测噪声(地面机器人转向时,yaw 易受干扰,噪声通常更大)
线速度
vx
0.025
x 方向线速度的预测噪声(模型假设 “匀速”,但实际加减速、打滑会导致误差)

vy
0.025
y 方向线速度的预测噪声

vz
0.04
z 方向线速度(高度变化率)的预测噪声
角速度
vroll
0.01
横滚角速度的预测噪声

vpitch
0.01
俯仰角速度的预测噪声

vyaw
0.02
偏航角速度的预测噪声(转向时角速度易波动,噪声通常高于其他轴)
线加速度
ax
0.01
x 方向线加速度的预测噪声(加速度计本身有噪声,模型外推也会引入误差)

ay
0.01
y 方向线加速度的预测噪声

az
0.015
z 方向线加速度的预测噪声(受重力影响,若用于地面设备,az 需覆盖 “静止时 1g” 的误差)

3. 如何调参?核心逻辑

EKF 的本质是 **“信任模型预测” vs “信任传感器观测”** 的平衡,process_noise_covariance 控制对 “模型预测” 的信任度

  • 对角线值越小

    :认为 “模型预测的状态越准”,滤波会更依赖预测(适合模型精确的场景,比如高精度机械臂)。
  • 对角线值越大

    :认为 “模型预测的不确定性越高”,滤波会更依赖传感器观测(适合干扰多、模型粗糙的场景,比如地面机器人打滑)。


一句话总结

process_noise_covariance 是 EKF 的 “模型信任度旋钮”—— 值越小,越相信模型预测;值越大,越依赖传感器。调参时,先保持默认值跑场景,遇到 “收敛慢”(比如转向后角度回不来)就增大对应维度的噪声,遇到 “太敏感”(比如传感器一飘就跟丢)就减小噪声,反复迭代即可。



评论


相关推荐

技术专区

关闭
孺子可教也什么意思 75属什么生肖 聂的拼音是什么 达菲是什么药 体感是什么意思
为什么狐臭女很漂亮 胎位左枕前是什么意思 结膜炎是什么原因引起的 偷窥是什么意思 睡觉开风扇有什么危害
胃癌低分化是什么意思 纸醉金迷什么意思 元辰是什么意思 为什么医生不推荐特立帕肽呢 晚的反义词是什么
外交部长是什么级别 驾驶证扣6分有什么影响 8月15号什么星座 吃什么能胖起来 2021年是属什么年
睚眦是什么意思hcv8jop6ns4r.cn adh医学上是什么意思clwhiglsz.com 白血病有什么征兆hcv8jop7ns6r.cn 韬字五行属什么onlinewuye.com 甲亢属于什么科室hanqikai.com
破处是什么意思hcv9jop2ns1r.cn 霜降是什么意思hcv8jop5ns9r.cn 泾渭分明是什么意思hcv7jop7ns2r.cn 咽口水喉咙痛吃什么药hcv8jop7ns3r.cn 梦见死蛇是什么预兆hcv9jop2ns0r.cn
apart是什么意思hcv9jop7ns2r.cn 莫逆之交是什么意思hcv7jop5ns5r.cn 什么动物跑得快hcv8jop0ns7r.cn 老年人腿脚无力是什么原因weuuu.com 引产什么意思jasonfriends.com
2月22日什么星座hcv8jop1ns6r.cn 白内障有什么症状表现hcv7jop6ns6r.cn 空窗期什么意思hcv8jop8ns0r.cn 2月15日什么星座hcv8jop5ns8r.cn 低血压高是什么原因造成的hcv8jop7ns2r.cn
百度