验孕棒什么时候测最准确| 沦落什么意思| 为什么总是做梦| 滥竽充数的充是什么意思| 平起平坐代表什么生肖| 指甲紫色是什么原因| 人为什么会长痔疮| 喝紫苏水有什么功效| 木姜子是什么| 心脏搭桥是什么意思| 三叉神经痛用什么药| 乙巳年是什么命| 意识是什么意思| 虞是什么意思| 确立是什么意思| 牙龈一直肿不消什么原因| 检查怀孕挂什么科| 信手拈来是什么意思| 什么病不能吃竹笋| 颈椎不舒服挂什么科| 梦见蛇是什么意思啊| itp是什么意思| 男人性功能太强是什么原因| 头皮痛是什么原因| 试纸什么时候用最准确| 8五行属什么| 什么品种的狗最粘人| 米干是什么| 吃什么对心脏有好处| 什么是马甲线| 首发是什么意思| 榴莲不能和什么一起吃| idh是什么意思| 手麻是什么病| sassy是什么意思| 胆结石挂什么科室| 砚是什么东西| 卵巢早衰吃什么药| 做梦梦到蛆是什么意思| 王景读什么| 医保断了一个月有什么影响| 男人吃蚂蚱有什么好处| 微波炉不热是什么原因| 阴瑜伽是什么意思| 法界是什么意思| 花甲之年是什么意思| 汇总压缩是什么意思| 肠粉是用什么粉做的| 为什么外阴老是长疖子| 慢慢地什么填词语| 晧字五行属什么| 肾阴虚有什么症状| 牛的尾巴有什么作用| 膝盖咔咔响吃什么药| 国历是什么意思| 血栓是什么| 口坐读什么| 女性为什么不适合喝茉莉花茶| 拉红尿是什么原因| 财五行属什么| 宫颈纳氏囊肿是什么意思严重吗| 白头翁是什么意思| 梦见人头是什么征兆| 女性潮红是什么意思| 青稞是什么东西| 床上放什么可以驱蟑螂| 胃胀吃什么| 诸事皆宜是什么意思| 床垫什么材质的好| 舌头生疮是什么原因引起的| 什么叫做流年| 肾有问题挂什么科| cab是什么意思| 霜降是什么季节| 末梢神经炎吃什么药| 外露什么意思| 什么动物怕热| 什么样的人容易垂体瘤| 双肾泥沙样结石是什么意思| 宫颈短是什么意思| 农历5月20日是什么星座| 月经期间喝什么比较好| 孕妇d2聚体高是什么原因| 女人吃什么排卵最快| 四大菩萨分别保佑什么| 吃完桃子不能吃什么| 愚痴是什么意思| 眼睛五行属什么| 打疫苗前后要注意什么| 天荒地老是什么生肖| y是什么元素| 腺样体肥大吃什么药| 氯化钠敷脸有什么作用| 七宗罪是什么| 梦见别人拉屎是什么意思| 奥美拉唑主要治什么| 为什么会干呕| 导语是什么| 左心室高电压什么意思| 变卖是什么意思| 宝宝吃什么增强抵抗力| 今年50岁属什么| chanel是什么牌子| 甲亢不能吃什么食物| 美丽的动物是什么生肖| 便秘吃什么蔬菜| 尿蛋白定量是什么意思| 不撞南墙不回头是什么意思| 尿里带血是什么原因| 舌头上有红点是什么原因| 胖大海是什么东西| 高岗为什么自杀| 亚麻跌是什么意思| 两个夫一个车是什么字| 1月26是什么星座| 容易手麻脚麻是什么原因| 肝浸润是什么意思| 吃榴莲对身体有什么好处| 吃了榴莲不可以吃什么| 生肖龙和什么生肖相冲| ml什么单位| 哪是什么意思| 示字旁与什么有关| 为什么射精是流出来的| 皮肤湿疹用什么药膏| 什么是棱长| 内科包括什么| 口腔溃疡挂什么科| 华侨是什么| 异丙醇是什么| 微博id是什么| 女装什么牌子好| 胆固醇是什么东西| pppd是什么意思| 喉咙挂什么科室| 青岛为什么叫青岛| 肠鸣吃什么药| 他达拉非是什么药| 45岁属什么| 黄埔军校现在是什么学校| 玩是什么意思| 维生素b4又叫什么| 脑梗挂号挂什么科室| 水手是干什么的| 为什么大拇指只有两节| 高血糖吃什么降得快| 内痔疮吃什么药最好| 52年属什么生肖| 什么水果补血效果最好| 急性肠胃炎可以吃什么| 运六月有什么说法| 下联是什么| 湿热体质吃什么药| le是什么元素| 蛾子吃什么| 风口浪尖是什么意思| 男士睾丸疼是什么原因| 嘴角烂是什么原因| 文爱什么意思| 突然消瘦是什么原因| 好嘛是什么意思| 联通查话费打什么号码| 乙肝表面抗原阴性是什么意思| 懿字五行属什么| 什么时间量血压最准确| 小孩睡觉出很多汗是什么原因| 元宵节的习俗是什么| 樱花什么时候开花| 什么是红肉| 县团级是什么级别| 小王子讲了什么故事| 脖子短适合什么发型| 牙膏什么牌子好| 接风是什么意思| 头孢过敏用什么药代替| 橄榄菜长什么样子图片| 天天喝可乐有什么危害| 什么叫银屑病| april什么意思| 月经期间可以吃什么水果| 胃窦在胃的什么位置| 看诊是什么意思| 腰间盘突出用什么药好| 肌无力是什么症状| 补气血吃什么中成药最好| lgbtq是什么意思| 霍光和卫子夫什么关系| 微喇裤配什么鞋子好看| 鬼压床是什么原因造成的| 左耳疼痛什么原因引起| 大拇指旁边是什么指| 血压为什么会高| 灯光什么| 丁胺卡那又叫什么药名| 170是什么尺码| 为什么腿老是抽筋| 拉肚子是什么原因导致的| 什么下奶最快最多| 女孩子为什么会痛经| 北京有什么特产| 西安有什么好吃的特产| 玄武岩属于什么岩石| 红豆和什么搭配最好| 什么的愿望| 肺部微结节是什么意思| 为什么端午节要吃粽子| 熬中药用什么锅最好| 久负盛名的负是什么意思| 肠胃感冒是什么症状| 性格开朗是什么意思| 锑对人体有什么危害| 检查hpv需要注意什么提前注意什么| 月半是什么意思| ck是什么牌子的包包| thirty什么意思| 描述是什么意思| 吃什么可以变白| drg是什么意思| 尿液发黄是什么病| 什么的山野| 羊得布病什么症状| 时光荏苒是什么意思| 蟑螂的天敌是什么| 早晨起来口干舌燥是什么原因| 做梦吃酒席什么预兆| 爱马仕是什么| 背疼什么原因| 脾大是什么原因引起的| 暴力熊是什么牌子| 夏天适合穿什么衣服| 双肺纹理增多增粗是什么病| 送什么礼物给女老师| 吃什么降血糖最快| 做什么生意好挣钱| 甲状腺什么原因引起的| 嘴唇干燥是什么原因| 背上长痘痘是什么原因| 九月15是什么星座| 思诺思是什么药| 肾阴虚吃什么食物补| 什么心丧气| 小孩老咳嗽是什么原因| 病字旁加且念什么| 受精卵着床是什么意思| 防蓝光眼镜有什么好处| 梦见小猪仔什么意思| 什么是优质碳水| 抗hbc阳性是什么意思| 什么是普拉提| 脾虚湿盛吃什么中药| 神采奕奕是什么意思| 开小差是什么意思| 宫颈欠光滑是什么意思| 大便不规律是什么原因| 什么奔什么走| 鼻涕由清变黄说明什么| 昱字五行属什么| 宥怎么读什么意思| 李小龙和丁佩什么关系| 舟状腹见于什么疾病| 茉莉什么时候开花| 男性夜间盗汗是什么原因| 五谷杂粮是什么| 胎膜早破是什么原因引起的| 百度

新闻中心

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

江西团开放日  创新引领高质量发展,坚决打赢脱贫攻坚战

作者: 时间:2025-08-04 来源:硬十 收藏
1、什么是
百度   科技创新人才  中国专利金奖获奖专利的发明人、获得3项以上(含)发明专利的独立完成人、以第二作者及以上身份获得6项以上(含)发明专利的主要完成人,其专利在京落地转化并取得显著经济社会效益的。

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



评论


相关推荐

技术专区

关闭
酒酿蛋什么时候吃效果最好 肚子疼用什么药好 mps是什么意思 喝酒不能吃什么东西 儿茶酚胺是什么
为什么天气热皮肤就痒 脚气是什么症状 语五行属什么 质子泵抑制剂是什么药 神经衰弱是什么意思
lp是什么的简称 兔儿爷是什么意思 三焦不通吃什么中成药 塔克是什么食物 什么果
午餐肉是什么肉 伊朗说什么语言 体育总局局长什么级别 耘字五行属什么 saucony是什么牌子
明信片是什么hcv7jop5ns3r.cn 泛醇是什么hcv8jop3ns9r.cn 老花眼有什么症状0735v.com 心肾不交吃什么中成药hcv8jop6ns6r.cn 脱发吃什么维生素wuhaiwuya.com
阿托伐他汀治什么病hcv7jop6ns1r.cn 肌酐偏低是什么意思hcv9jop6ns2r.cn 心律不齐吃什么药hcv7jop6ns3r.cn 什么是肛漏hcv7jop6ns2r.cn 三七泡酒有什么功效hcv8jop6ns8r.cn
生肠是什么hcv8jop8ns8r.cn 每天喝豆浆有什么好处hcv8jop1ns6r.cn 如果怀孕了会有什么预兆hcv8jop3ns7r.cn 星座之王是什么座hcv9jop3ns2r.cn 卵泡是什么东西zhiyanzhang.com
睾丸癌是由什么引起的hcv8jop2ns4r.cn 身心疲惫是什么意思hcv8jop3ns6r.cn ganni是什么牌子hcv8jop5ns2r.cn 洗牙后要注意什么hcv7jop6ns1r.cn 男性尿道炎是什么原因引起的hcv9jop1ns8r.cn
百度