基于提升小波变换的雷达生命信号去噪技术
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第 34 卷 第 3 期 2014 年 3 月
光 学 学 报 ACTA OPTICA SINICA
Vol.34,No.3 March,2014
基于提升小波变换的雷达生命信号去噪技术
杨秀芳 张 伟 杨宇祥
(西安理工大学机械与精密仪器工程学院,陕西 西安 710048)
摘要 对于噪声干扰严重、非线性、非平稳、多奇异点 的 微 弱 雷 达 生 命 信 号 而 言,去 噪 是 对 有 用 信 号 进 行 分 析 前 的 必要手段。基于多普勒效应原理和雷达噪声的统计特 性,建 立 了 雷 达 生 命 信 号 的 模 型,分 别 用 传 统 小 波 变 换 和 提 升小波变换对强噪声干扰下的生命雷达信号进行了去噪处理,结果表明被强噪声污 染 的 雷 达 生 命 信 号 可 以 用 传 统 小波变换的方法和提升小波变换法对其有效去噪,提 升 小 波 变 换 去 噪 效 果 优 于 传 统 小 波 变 换 去 噪 效 果 ,其 信 噪 比 (SNR)和均方误差(MSE)两个性能指标均高于传统小波去噪。对雷达生命信号进行 去 噪 处 理 时,使 用 的 小 波 基 函 数 是 sym8,分 解 层 数 为 3 层 。 关 键 词 信 号 处 理 ;雷 达 生 命 信 号 去 噪 ;提 升 小 波 变 换 ;信 噪 比 ;均 方 误 差 中 图 分 类 号 TP312 文 献 标 识 码 A doi:10.3788/AOS201434.0328003
被人体表面反射,由 于 心 跳 和 呼 吸 运 动 引 起 的 多 普
勒效应将产生一种窄带调相反射波 S(t):
S(t)= Acos[2πf0t+θ0 +θ′+φ(t)], (2) 式中 A 为回波信号的幅 值,θ′ 为 由 被 测 对 象 与 雷 达
间的距离所产生 的 相 位,φ(t)为 由 多 普 勒 效 应 产 生 的 附 加 相 移 。由 多 普 勒 效 应 理 论 可 知
1 引 言
生 命 信 号 是 一 种 非 线 性 、非 平 稳 、多 奇 异 点 的 微 弱信号,一般正常人 的 呼 吸 和 心 跳 频 率 分 别 为 每 分 钟 16~20 次 和 60~100 次 。 性 别 的 差 异 、年 龄 的 不 同、运动加强与否等 都 会 引 起 心 跳 和 呼 吸 频 率 的 变 化 ,且 这 些 信 号 在 采 集 过 程 中 容 易 受 到 仪 器 、人 体 活
动 等 干 扰 因 素 的 影 响 而 引 入 噪 声[1],给 分 析 人 体 生 理状态参数带来很 大 困 难,因 此 信 号 去 噪 成 为 分 析 前的必要步骤,而 小 波 去 噪 是 最 常 用 的 方 法。 文 献 [2]介 绍 了 小 波 分 解 与 重 构 法 去 噪 ,指 出 该 方 法 对 于 已知噪声频率范围,且 噪 声 和 信 号 的 频 带 相 互 分 离 时 非 常 有 效 ,对 实 际 中 广 泛 存 在 的 白 噪 声 ,去 噪 效 果
12ArAcos[2π(2f0 +fd)t+θ0 +θ]+ 12ArAcos(2πfdt+θ′+φ0).
(5) (6)
经低通滤波器后获得测量信号 So(t): So(t)= 12ArAcos(2πfdt+θ′+φ0).
(7)
雷达回波中的噪声 信 号 成 分 复 杂,无 法 给 出 时 域 表 达式,综合考虑认为 其 服 从 高 斯 分 布 且 其 功 率 谱 密 度 函 数 也 呈 高 斯 状 ,则 雷 达 生 命 信 号 的 表 达 式 为
φ(t)= 2π×2 Xλ(t),
(3)
式中 X(t)为 人 体 体 动 位 移,其 值 围 绕 平 均 距 离
x(t)上下波动,以 Δx(t)表示人体 体 动 位 移 的 波 动
幅度,则 X(t)= x(t)+Δx(t),λ 为雷达波波长值。
于是ຫໍສະໝຸດ Baidu普勒附加相移可以表示为
φ(t)=2π×2xλ(t)+2π×2Δxλ(t)=φ0 +2πfdt,
E-mail:yxf5078@163.com
0328003-1
光 学 学 报
则 较 差 。 文 献 [3]提 出 模 极 大 值 法 去 噪 ,它 主 要 适 用 于 信 号 中 混 有 白 噪 声 ,且 含 有 较 多 奇 异 点 的 情 况 ,可 以有效地保留信号 奇 异 点 的 信 息,去 噪 后 的 信 号 没 有 多 余 震 荡 ,但 小 波 分 解 尺 度 的 选 择 对 其 非 常 重 要 , 尺度小,会产生伪 极 值 点,尺 度 大,会 丢 失 某 些 重 要 的局部奇异性。文献[4]用 小 波 变 换 的 方 法 对 光 纤 温度传感系统的测量信号实现了有效的降噪处理。 文 献 [5-6]指 出 相 关 性 去 噪 在 分 析 信 号 的 边 界 方 面 虽 然 去 噪 效 果 比 较 稳 定 ,但 其 计 算 量 较 大 ,并 且 需 要 估计噪声方差。文献[7]提 出 小 波 阈 值 去 噪 主 要 适 用于信号中混有白 噪 声 的 情 况,原 始 信 号 的 特 征 尖 峰点可得到很好的 保 留,其 中 的 软 阈 值 去 噪 可 使 估 计 信 号 实 现 最 大 均 方 误 差 最 小 化 ,计 算 速 度 很 快 ,但 其去噪效果依赖于信噪比的大小和合适的阈值。文 献 [8]用 平 移 不 变 法 去 噪 ,主 要 适 用 于 信 号 中 混 有 白 噪声且含有若干个 不 连 续 点 的 情 况,可 以 有 效 去 除 阈值法中产生的伪 吉 布 斯 现 象,得 到 比 阈 值 法 更 小 的 均 方 误 差 ,信 噪 比 也 得 到 一 定 的 提 高 ,但 计 算 速 度 没有阈值法快。上述的小波变换方法用于信号去噪 时在自身构造及对信号进行时频域变换时都需要耗 费大量系统资源[9],文 献 [10]用 自 适 应 提 升 小 波 变 换 对 心 音 信 号 进 行 了 降 噪 处 理 ,取 得 了 不 错 的 效 果 。 文 献 [11-13]用 提 升 小 波 变 换 分 别 对 电 极 法 测 量 的 弱噪声干扰的生 命 信 号、弱 噪 声 干 扰 的 Cusp 信 号、 弱噪声干 扰 的 Heavisine信 号、弱 噪 声 干 扰 的 多 普 勒雷达信号实现了去 噪 处 理,文 献 [14]分 别 用 小 波 变换和提升小波变换对受到较大噪声污染的 X 射 线脉冲 星 信 号 进 行 降 噪 处 理,结果提升小波的降 噪 效 果 优 于 传 统 的 小 波 变 换 。 参 考 以 上 的 文 献 ,考 虑 到 雷达生命信号本身的特点和各种小波去噪方法的计 算速度,本文提出将提升小波变换用于雷达生命信号 去噪处理,结果实现了强噪声干扰下的雷达生命信号 的去噪处理,而且方法简单,计算量小。
Abstract For the nonlinear,non-stationary and weak radar life signal with serious noise interference and multiple singularity points,denoising is the essential means before analyzing the useful signal.The radar life signal model is established based on the principle of Doppler effect and the statistical characteristics of noise.Wavelet transform and
lifting wavelet transform are applied respectively on the radar life signal under strong noise background for denoising. The results show that for radar life signal under strong noise pollution,it is an efficient denoising method to use wavelet transform or lifting wavelet transform.The lifting wavelet transform denoising effect is better than traditional denoising methods of wavelet transform.The signal-to-noise ratio (SNR)and mean square error (MSE) are higher than those of traditional wavelet denoising.In the radar life signal denoising,the wavelet basis function is
该过程其实就是低 速 运 动 目 标 的 检 测 过 程,而 低 速
运动目标的检测离不开多普勒效应。
图1为基于多普勒效应的雷达生命信号模型。
现假设雷达的发射信号为
Sr(t)= Arcos(2πf0t+θ0),
(1)
式中 Ar 为发射信号的幅值,f0 为发射信号的频率,
θ0 为发射信号的初相位。当发射信号照射到人体后,
Denoising Technology of Radar Life Signal Based on Lifting Wavelet Transform
Yang Xiufang Zhang Wei Yang Yuxiang
(Faculty of Mechanical and Precision Instrument Engineering,Xi′an University of Technology, Xi′an,Shaanxi 710048,China)
sym8 and the number of decomposition layers is 3. Key words signal processing;radar life signal denoising;lifting wavelet transform;signal-to-noise ratio;mean square error OCIS codes 280.1415;280.5600;100.7410
(4)
式中φ0 为 平 均 距 离 x(t)引 起 的 固 定 附 加 相 位,fd 为人体微动引起的多普勒效应频移。从而 S(t)表示
为
S(t)= Acos[2π(f0 +fd)t+θ0 +θ′+φ0]= Acos[2π(f0 +fd)t+θ].
则 Sr(t)与 S(t)在 混 频 器 中 混 频 : So(t)=Sr(t)×S(t)= Arcos(2πf0t+θ0)×Acos(2π(f0 +fd)t+θ)=
0328003-2
杨 秀 芳 等 : 基 于 提 升 小 波 变 换 的 雷 达 生 命 信 号 去 噪 技 术
S(t)= 12ArAcos(2πfdt+θ′+φ0)+c(t),(8)
收 稿 日 期 :2013-09-24;收 到 修 改 稿 日 期 :2013-11-11 基 金 项 目 :陕 西 省 教 育 厅 科 学 研 究 计 划 (自 然 科 学 专 项 )(12JK0527) 作 者 简 介 :杨 秀 芳 (1965— ),女 ,副 教 授 ,硕 士 生 导 师 ,主 要 从 事 光 电 测 试 与 信 息 提 取 方 面 的 研 究 。
2 雷达生命信号检测原理
雷达生命信号指生命探测雷达回波中包含的雷
达 发 射 载 频 信 号 、生 命 信 号 和 杂 波 信 号 等 ,对 雷 达 生
命信号的检测就是将经人体表面微动反射后的雷达
回波信号与雷达发 射 信 号 进 行 混 频 解 调,得 出 与 人
体生理运动幅度成比例的呼吸和心跳信号的过程。
光 学 学 报 ACTA OPTICA SINICA
Vol.34,No.3 March,2014
基于提升小波变换的雷达生命信号去噪技术
杨秀芳 张 伟 杨宇祥
(西安理工大学机械与精密仪器工程学院,陕西 西安 710048)
摘要 对于噪声干扰严重、非线性、非平稳、多奇异点 的 微 弱 雷 达 生 命 信 号 而 言,去 噪 是 对 有 用 信 号 进 行 分 析 前 的 必要手段。基于多普勒效应原理和雷达噪声的统计特 性,建 立 了 雷 达 生 命 信 号 的 模 型,分 别 用 传 统 小 波 变 换 和 提 升小波变换对强噪声干扰下的生命雷达信号进行了去噪处理,结果表明被强噪声污 染 的 雷 达 生 命 信 号 可 以 用 传 统 小波变换的方法和提升小波变换法对其有效去噪,提 升 小 波 变 换 去 噪 效 果 优 于 传 统 小 波 变 换 去 噪 效 果 ,其 信 噪 比 (SNR)和均方误差(MSE)两个性能指标均高于传统小波去噪。对雷达生命信号进行 去 噪 处 理 时,使 用 的 小 波 基 函 数 是 sym8,分 解 层 数 为 3 层 。 关 键 词 信 号 处 理 ;雷 达 生 命 信 号 去 噪 ;提 升 小 波 变 换 ;信 噪 比 ;均 方 误 差 中 图 分 类 号 TP312 文 献 标 识 码 A doi:10.3788/AOS201434.0328003
被人体表面反射,由 于 心 跳 和 呼 吸 运 动 引 起 的 多 普
勒效应将产生一种窄带调相反射波 S(t):
S(t)= Acos[2πf0t+θ0 +θ′+φ(t)], (2) 式中 A 为回波信号的幅 值,θ′ 为 由 被 测 对 象 与 雷 达
间的距离所产生 的 相 位,φ(t)为 由 多 普 勒 效 应 产 生 的 附 加 相 移 。由 多 普 勒 效 应 理 论 可 知
1 引 言
生 命 信 号 是 一 种 非 线 性 、非 平 稳 、多 奇 异 点 的 微 弱信号,一般正常人 的 呼 吸 和 心 跳 频 率 分 别 为 每 分 钟 16~20 次 和 60~100 次 。 性 别 的 差 异 、年 龄 的 不 同、运动加强与否等 都 会 引 起 心 跳 和 呼 吸 频 率 的 变 化 ,且 这 些 信 号 在 采 集 过 程 中 容 易 受 到 仪 器 、人 体 活
动 等 干 扰 因 素 的 影 响 而 引 入 噪 声[1],给 分 析 人 体 生 理状态参数带来很 大 困 难,因 此 信 号 去 噪 成 为 分 析 前的必要步骤,而 小 波 去 噪 是 最 常 用 的 方 法。 文 献 [2]介 绍 了 小 波 分 解 与 重 构 法 去 噪 ,指 出 该 方 法 对 于 已知噪声频率范围,且 噪 声 和 信 号 的 频 带 相 互 分 离 时 非 常 有 效 ,对 实 际 中 广 泛 存 在 的 白 噪 声 ,去 噪 效 果
12ArAcos[2π(2f0 +fd)t+θ0 +θ]+ 12ArAcos(2πfdt+θ′+φ0).
(5) (6)
经低通滤波器后获得测量信号 So(t): So(t)= 12ArAcos(2πfdt+θ′+φ0).
(7)
雷达回波中的噪声 信 号 成 分 复 杂,无 法 给 出 时 域 表 达式,综合考虑认为 其 服 从 高 斯 分 布 且 其 功 率 谱 密 度 函 数 也 呈 高 斯 状 ,则 雷 达 生 命 信 号 的 表 达 式 为
φ(t)= 2π×2 Xλ(t),
(3)
式中 X(t)为 人 体 体 动 位 移,其 值 围 绕 平 均 距 离
x(t)上下波动,以 Δx(t)表示人体 体 动 位 移 的 波 动
幅度,则 X(t)= x(t)+Δx(t),λ 为雷达波波长值。
于是ຫໍສະໝຸດ Baidu普勒附加相移可以表示为
φ(t)=2π×2xλ(t)+2π×2Δxλ(t)=φ0 +2πfdt,
E-mail:yxf5078@163.com
0328003-1
光 学 学 报
则 较 差 。 文 献 [3]提 出 模 极 大 值 法 去 噪 ,它 主 要 适 用 于 信 号 中 混 有 白 噪 声 ,且 含 有 较 多 奇 异 点 的 情 况 ,可 以有效地保留信号 奇 异 点 的 信 息,去 噪 后 的 信 号 没 有 多 余 震 荡 ,但 小 波 分 解 尺 度 的 选 择 对 其 非 常 重 要 , 尺度小,会产生伪 极 值 点,尺 度 大,会 丢 失 某 些 重 要 的局部奇异性。文献[4]用 小 波 变 换 的 方 法 对 光 纤 温度传感系统的测量信号实现了有效的降噪处理。 文 献 [5-6]指 出 相 关 性 去 噪 在 分 析 信 号 的 边 界 方 面 虽 然 去 噪 效 果 比 较 稳 定 ,但 其 计 算 量 较 大 ,并 且 需 要 估计噪声方差。文献[7]提 出 小 波 阈 值 去 噪 主 要 适 用于信号中混有白 噪 声 的 情 况,原 始 信 号 的 特 征 尖 峰点可得到很好的 保 留,其 中 的 软 阈 值 去 噪 可 使 估 计 信 号 实 现 最 大 均 方 误 差 最 小 化 ,计 算 速 度 很 快 ,但 其去噪效果依赖于信噪比的大小和合适的阈值。文 献 [8]用 平 移 不 变 法 去 噪 ,主 要 适 用 于 信 号 中 混 有 白 噪声且含有若干个 不 连 续 点 的 情 况,可 以 有 效 去 除 阈值法中产生的伪 吉 布 斯 现 象,得 到 比 阈 值 法 更 小 的 均 方 误 差 ,信 噪 比 也 得 到 一 定 的 提 高 ,但 计 算 速 度 没有阈值法快。上述的小波变换方法用于信号去噪 时在自身构造及对信号进行时频域变换时都需要耗 费大量系统资源[9],文 献 [10]用 自 适 应 提 升 小 波 变 换 对 心 音 信 号 进 行 了 降 噪 处 理 ,取 得 了 不 错 的 效 果 。 文 献 [11-13]用 提 升 小 波 变 换 分 别 对 电 极 法 测 量 的 弱噪声干扰的生 命 信 号、弱 噪 声 干 扰 的 Cusp 信 号、 弱噪声干 扰 的 Heavisine信 号、弱 噪 声 干 扰 的 多 普 勒雷达信号实现了去 噪 处 理,文 献 [14]分 别 用 小 波 变换和提升小波变换对受到较大噪声污染的 X 射 线脉冲 星 信 号 进 行 降 噪 处 理,结果提升小波的降 噪 效 果 优 于 传 统 的 小 波 变 换 。 参 考 以 上 的 文 献 ,考 虑 到 雷达生命信号本身的特点和各种小波去噪方法的计 算速度,本文提出将提升小波变换用于雷达生命信号 去噪处理,结果实现了强噪声干扰下的雷达生命信号 的去噪处理,而且方法简单,计算量小。
Abstract For the nonlinear,non-stationary and weak radar life signal with serious noise interference and multiple singularity points,denoising is the essential means before analyzing the useful signal.The radar life signal model is established based on the principle of Doppler effect and the statistical characteristics of noise.Wavelet transform and
lifting wavelet transform are applied respectively on the radar life signal under strong noise background for denoising. The results show that for radar life signal under strong noise pollution,it is an efficient denoising method to use wavelet transform or lifting wavelet transform.The lifting wavelet transform denoising effect is better than traditional denoising methods of wavelet transform.The signal-to-noise ratio (SNR)and mean square error (MSE) are higher than those of traditional wavelet denoising.In the radar life signal denoising,the wavelet basis function is
该过程其实就是低 速 运 动 目 标 的 检 测 过 程,而 低 速
运动目标的检测离不开多普勒效应。
图1为基于多普勒效应的雷达生命信号模型。
现假设雷达的发射信号为
Sr(t)= Arcos(2πf0t+θ0),
(1)
式中 Ar 为发射信号的幅值,f0 为发射信号的频率,
θ0 为发射信号的初相位。当发射信号照射到人体后,
Denoising Technology of Radar Life Signal Based on Lifting Wavelet Transform
Yang Xiufang Zhang Wei Yang Yuxiang
(Faculty of Mechanical and Precision Instrument Engineering,Xi′an University of Technology, Xi′an,Shaanxi 710048,China)
sym8 and the number of decomposition layers is 3. Key words signal processing;radar life signal denoising;lifting wavelet transform;signal-to-noise ratio;mean square error OCIS codes 280.1415;280.5600;100.7410
(4)
式中φ0 为 平 均 距 离 x(t)引 起 的 固 定 附 加 相 位,fd 为人体微动引起的多普勒效应频移。从而 S(t)表示
为
S(t)= Acos[2π(f0 +fd)t+θ0 +θ′+φ0]= Acos[2π(f0 +fd)t+θ].
则 Sr(t)与 S(t)在 混 频 器 中 混 频 : So(t)=Sr(t)×S(t)= Arcos(2πf0t+θ0)×Acos(2π(f0 +fd)t+θ)=
0328003-2
杨 秀 芳 等 : 基 于 提 升 小 波 变 换 的 雷 达 生 命 信 号 去 噪 技 术
S(t)= 12ArAcos(2πfdt+θ′+φ0)+c(t),(8)
收 稿 日 期 :2013-09-24;收 到 修 改 稿 日 期 :2013-11-11 基 金 项 目 :陕 西 省 教 育 厅 科 学 研 究 计 划 (自 然 科 学 专 项 )(12JK0527) 作 者 简 介 :杨 秀 芳 (1965— ),女 ,副 教 授 ,硕 士 生 导 师 ,主 要 从 事 光 电 测 试 与 信 息 提 取 方 面 的 研 究 。
2 雷达生命信号检测原理
雷达生命信号指生命探测雷达回波中包含的雷
达 发 射 载 频 信 号 、生 命 信 号 和 杂 波 信 号 等 ,对 雷 达 生
命信号的检测就是将经人体表面微动反射后的雷达
回波信号与雷达发 射 信 号 进 行 混 频 解 调,得 出 与 人
体生理运动幅度成比例的呼吸和心跳信号的过程。