【雷达原理】雷达杂波抑制方法

目录

一、杂波及其特点

1.1 什么是杂波?

1.2 杂波的频谱特性

二、动目标显示(MTI)技术

2.1 对消原理

2.2 数字对消器设计

三、MATLAB仿真

3.1 对消效果验证

3.2 代码


一、杂波及其特点

1.1 什么是杂波?

杂波是相对目标回波而言的,它是指对雷达正常工作时的目标回波有干扰的回波信号。对于工作在地面的雷达,杂波主要来源于固定建筑物、围栏、树木等静止目标或慢速运动目标产生的回波以及地表反射的回波。

由于单个杂波分量具有随机的相位和幅度,因此合成的总杂波也是随机的,具有类似热噪声的统计特性,并且杂波信号的强度要远大于接收机内热噪声的强度。因此,雷达在强杂波背景下检测目标的能力主要取决于信杂比,而不是信噪比。

杂波抑制的目的就是采用数字信号处理技术滤除或者抑制固定目标杂波和慢速目标杂波,从而保留运动目标回波。

1.2 杂波的频谱特性

雷达工作时可能会碰到的杂波包括地物、海浪、云雨及敌人施放的箔条等。除了孤立的建筑物可认为是固定的点目标外,大多数杂波均属于分布杂波且包含内部运动。当天线不扫描时,固定杂波的功率谱是位于nf_{r}(n=0,±1,±2,...) 位置上的离散谱线,其中f_{r}为脉冲重复频率。当天线扫描时,由于回波脉冲数有限,将引起谱线的展宽。此外,由于天线扫描引起双程天线方向图对回波信号调幅,也会导致功率谱展宽。

一般杂波谱可用高斯函数表示为:

G(f)=G_{0}e^{-f^2/2\sigma _{s}^2}

其中,\sigma _{s}=0.265/n,n为单程天线方向图3dB宽度内的脉冲数。

二、动目标显示(MTI)技术

2.1 对消原理

当固定目标、地杂波等与运动目标处于同一距离单元时,前者的回波通常较强,以至于运动目标的回波被淹没其中,故必须设法对二者进行区分。

因固定目标回波的多普勒频率为零,慢速运动的杂波中所含的多普勒频移也集中在零频附近,它们的回波经相位检波后,输出信号的相位将不随时间变化或者不随时间做缓慢变化,反映在幅度上则为其幅度不随时间变化或者不随时间做缓慢变化。相反运动目标回波经相检输出后,因其相位随时间变化较大,反映在幅度上也是幅度随时间变化较快。

因此若将同一距离单元在相邻重复周期内的脉冲做相减运算,则固定目标回波将被完全对消,慢速杂波也将得到很大程度的衰减,只有运动目标回波得以保留。这样就可以将固定目标、慢速目标与运动目标区别开来。

2.2 数字对消器设计

数字对消器本质上就是一种数字滤波器,根据所设计滤波器的幅频响应特性,可以滤除相应频率的信号分量。常用的MTI滤波器为延迟线对消器,根据对消次数的不同,可分为单延迟线对消器、双延迟线对消器和多延迟线对消器。

1、单延迟线对消器

1-1 单延迟线对消器结构图

 单延迟线对消器的结构如图1-1所示,滤波器的输出为输入的两个脉冲进行运算得到,该滤波器也称为双脉冲对消器或者一次对消器,其用数学公式可表示为:

y(t)=x(t)-x(t-T)

其中,T为脉冲重复周期PRT。

对上式求Z变换,可得:Y(z)=X(z)-z^{-1}X(z)

由此可知其系统函数为:H(z)=1-z^{-1}

由傅里叶变换与Z变换的关系,z=e^{jwT}=e^{j2\pi fT},得到对消器的功率增益为:

\left | H(w) \right |^2=H(w)H^{*}(w)=(1-e^{-jwT})(1-e^{jwT})=4(sin(\pi fT))^2

进而得到其幅频响应为:

\left | H(w) \right |=2\left | sin(\pi fT) \right |

单延迟线对消器的幅频响应曲线如图1-2所示,取脉冲重复周期T=10^{-3},则脉冲重复频率f_{r}=1000,由滤波器的幅频特性曲线可知,在频率f=nf_{r}(n=0,\pm 1,\pm 2,...)处有零点,因而该滤波器起到了抑制固定目标回波和慢速运动目标回波的作用。

1-2 单延迟线对消器的幅频响应

2、双延迟线对消器

2-1 双延迟线对消器结构图

 双延迟线对消器的结构如图2-1所示,它由两个单延迟线对消器级联而成。双延迟线对消器经常被称为三脉冲对消器或者二次对消器,其数学表达式如下:

y(t)=x(t)-2*x(t-T)+x(t-2T)

对上式求Z变换,可得:Y(z)=X(z)-2z^{-1}X(z)+z^{-2}X(z)

由此可知其系统函数为:H(z)=1-2z^{-1}+z^{-2}=(1-z^{-1})^2,从系统函数的表达式可以看出其等效于两个一次对消器级联。

由此得到二次对消器的功率增益为:

\left | H(w) \right |^2=\left | H_{1}(w) \right |^2\cdot \left | H_{1}(w) \right |^2=16(sin(\pi fT))^4

式中,\left | H_{1}(w) \right |^2为单延迟线对消器的功率增益。

进而得到其幅频响应为:

\left | H(w) \right |=4\left | sin(\pi fT) \right |^2

双延迟线对消器的幅频响应曲线如图2-2所示,可以看出,双延迟线对消器比单延迟线对消器具有更好的响应(更深的凹口和更平坦的通带响应)。单延迟线对消器的频率响应较差,原因在于其阻带没有宽的凹口。而双延迟线对消器无论在阻带还是通带上都比单延迟线对消器有更好的频率响应,因此得到了更广泛的应用。

2-2 延迟线对消器的幅频响应

3、多延迟线对消器

3-1 多延迟线对消器

多延迟线对消器的结构如图3-1所示,它是由多个单延迟线对消器级联而成,N延迟线对消器的数学表达式为:

y(t)=\sum_{n=0}^{N}w_{n}x(t-nT)

式中,N为对消器的次数,对消器的系数w_{n}为二项式系数,用下式计算:

w_{n}=(-1)^nC_{N}^{n}=(-1)^n\frac{N!}{(N-n)!n!},n=0,1,...,N

N次对消器的系统函数可表示为:

H(z)=(1-z^{-1})^N=\sum_{n=0}^{N}w_{n}z^{-n}

其幅频响应和相频响应分别为:

\left | H(w) \right |=\left | 2sin(wT/2) \right |^N=\left | 2sin(\pi fT) \right |^N

\varphi (w)=N(\frac{\pi }{2}-wT/2)=N(\frac{\pi }{2}-\pi fT)

可见,相位响应φ(w)与w是线性关系,所以对消器是一种线性相位滤波器,回波信号通过它后,相位关系不产生非线性变化。

三、MATLAB仿真

3.1 对消效果验证

由第二节的内容可知,MTI技术实质上就是设计一个滤波器,通过对多组脉冲回波的同一个距离单元加权求和即可实现。对于静止目标在该距离单元频谱固定,加权求和后将固定的频谱相互抵消;对于运动目标在该距离单元频谱存在多普勒频移,加权求和后目标频谱不会完全抵消。

通过MATLAB仿真多组具有静止目标和运动目标的回波信号,对比不同次数的对消器对杂波的抑制效果。
案例及对应的代码会在博文【雷达测速原理及实现方法】写完后更新。

3.2 代码

clc;
clear;
close all;

%% 单延迟线对消器
PRT = 1e-3;
PRF = 1/PRT;
f = (-2000:1:2000);
H1 = 2*abs(sin(pi*f*PRT));
H1 = H1/max(H1);
H1_d = mag2db(H1);
figure;
subplot(2,1,1);plot(f/PRF,H1);
xlabel('归一化频率f/fr');ylabel('幅度');title('延迟线对消器的幅频响应');
subplot(2,1,2);plot(f/PRF,H1_d);
xlabel('归一化频率f/fr');ylabel('幅度(dB)');ylim([-50 0]);

%% 双延迟线对消器
H2 = 4*abs(sin(pi*f*PRT)).^2;
H2 = H2/max(H2);
H2_d = mag2db(H2);
figure;plot(f/PRF,H1_d);hold on;plot(f/PRF,H2_d);ylim([-50 0]);
xlabel('归一化频率f/fr');ylabel('幅度(dB)');title('延迟线对消器的幅频响应');
legend('单延迟线对消器','双延迟线对消器');

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/115996.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

IDEA中如何移除未使用的import

👨🏻‍💻 热爱摄影的程序员 👨🏻‍🎨 喜欢编码的设计师 🧕🏻 擅长设计的剪辑师 🧑🏻‍🏫 一位高冷无情的编码爱好者 大家好,我是全栈工…

链栈的练习

链栈练习 相关内容&#xff1a;栈的链式存储结构&#xff08;链栈&#xff09; //链栈的初始化、判空、入栈、出栈、读取栈顶元素 //链栈的结点&#xff1a;数据域、指针域 #include<stdio.h> #include<stdlib.h> typedef int Status; #define OK 1 #define ERRO…

算法通关村第五关-白银挑战队列经典问题

大家好我是苏麟 , 今天带来几道经典小题 . 大纲 两数之和 两数之和 相信大家对这道题还是很眼熟的 , 打开LeetCode第一道题就是它 , 对它可真的又爱又恨 , 很多新手朋友们想刷LeetCode但又不知道从哪开始就打开了第一题 , 结果就对算法失去了信心 . 这道题找对方法还是很容易…

2022最新版-李宏毅机器学习深度学习课程-P32 Transformer

一、 seq2seq 1. 含义 输入一个序列&#xff0c;机器输出另一个序列&#xff0c;输出序列长度由机器决定。 文本翻译&#xff1a;文本至文本&#xff1b;  语音识别&#xff1a;语音至文本&#xff1b;  语音合成&#xff1a;文本至语音&#xff1b;  聊天机器人&#…

R语言_RColorBrewer包--全平台可用

R语言_RColorBrewer包–全平台可用

微服务架构——笔记(2)

微服务架构——笔记&#xff08;2&#xff09; 一、客户客户端模块 文章来源B站视频 尚硅谷SpringCloud框架开发教程(SpringCloudAlibaba微服务分布式架构丨Spring Cloud)教程 本次笔记内容为消费者订单Module模块 1.1 项目名称、目录结构 1.2 Pom.xml <?xml version&q…

【kafka】记一次kafka基于linux的原生命令的使用

环境是linux&#xff0c;4台机器&#xff0c;版本3.6&#xff0c;kafka安装在node 1 2 3 上&#xff0c;zookeeper安装在node2 3 4上。 安装好kafka&#xff0c;进入bin目录&#xff0c;可以看到有很多sh文件&#xff0c;是我们执行命令的基础。 启动kafka&#xff0c;下面的…

【使用Python编写游戏辅助工具】第三篇:鼠标连击器的实现

前言 这里是【使用Python编写游戏辅助工具】的第三篇&#xff1a;鼠标连击器的实现。本文主要介绍使用Python来实现鼠标连击功能。 鼠标连击是指在很短的时间内多次点击鼠标按钮&#xff0c;通常是鼠标左键。当触发鼠标连击时&#xff0c;鼠标按钮会迅速按下和释放多次&#xf…

自学黑客,一般人我劝你还是算了吧!

写在开篇 笔者本人 17 年就读于一所普通的本科学校&#xff0c;20 年 6 月在三年经验的时候顺利通过校招实习面试进入大厂&#xff0c;现就职于某大厂安全联合实验室。 我为啥说自学黑客&#xff0c;一般人我还是劝你算了吧&#xff01;因为我就是那个不一般的人。 首先我谈下…

Jenkins安装(Jenkins 2.429)及安装失败解决(Jenkins 2.222.4)

敏捷开发与持续集成 敏捷开发 敏捷开发以用户的需求进化为核心&#xff0c;采用迭代、循序渐进的方法进行软件开发。在敏捷开发中&#xff0c;软件项目在构建初期被切分成多个子项目&#xff0c;各个子项目的成果都经过测试&#xff0c;具备可视、可集成和可运行使用的特征。…

【Solidity】Remix在线环境及钱包申请

好久没有学习区块链方面的知识了&#xff0c;目前通过自学大致掌握了Fabric联盟链的搭建&#xff0c;链码编写、部署&#xff0c;api调用&#xff0c;可以独立开发出一些基于fabric的应用&#xff0c;感觉开发出去中心化的应用还是很有意思的&#xff0c;因为他与之前开发的ssm…

el-table表格设置——动态修改表头

(1) 首先是form表单写表单设置按钮&#xff1a; &#xff08;1.1&#xff09;使用el-popover&#xff0c;你需要修改的是this.colOptions&#xff0c;colSelect: <el-popover id"popover" popper-class"planProver" placement"bottom" width&…

USB PD v1.0快速充电通信原理

1 原理 本篇文章讲的快速充电是指USB论坛所发布的USB Power Delivery快速充电规范&#xff08;通过VBUS直流电平上耦合FSK信号来请求充电器调整输出电压和电流的过程&#xff09;&#xff0c;不同于本人发布的另一篇文章所讲的高通Quick Charger 2.0规范&#xff0c;因为高通QC…

【sql注入】sql关卡1~4

前言&#xff1a; 靶场自取 level-1 测试注入点 POC: 1,1,1,1"",1/1,1/0 》存在注入点 爆破 POC: id-1andextractvalue(1,concat(0x7e,user(),0x7e))-- level-2 尝试注入点 POC1:admin POC2:admin POC3:adminandsleep(3)-- POC4: adminandif(1,1,0)0-- POC…

Adobe acrobat 11.0版本 pdf阅读器修改背景颜色方法

打开菜单栏&#xff0c;编辑&#xff0c;首选项&#xff0c;选择辅助工具项&#xff0c;页面中 勾选 替换文档颜色&#xff0c;页面背景自己选择一个颜色&#xff0c;然后确定&#xff0c;即可&#xff01;

产品手册应该如何组织内容,以便用户能够快速找到所需信息?

产品手册应该如何组织内容&#xff0c;以便用户能够快速找到所需信息&#xff1f;这是一个关乎用户体验和产品文档效力的重要问题。当用户需要了解产品的功能、操作指南或故障排除时&#xff0c;他们希望能够轻松地找到准确、清晰的信息&#xff0c;而不是在冗长的手册中迷失方…

Leetcode—421.数组中两个数的最大异或值【中等】明天写一下字典树做法!!!

2023每日刷题&#xff08;十九&#xff09; Leetcode—421.数组中两个数的最大异或值 算法思想 参考自灵茶山艾府 实现代码 class Solution { public:int findMaximumXOR(vector<int>& nums) {int maxValue *max_element(nums.begin(), nums.end());int highId…

Java CAS是什么,它的底层原理?

文章目录 前言一、CAS是什么二、CAS底层原理1、UnSafe类&#xff08;Native方法&#xff09;2、CAS思想&#xff08;自旋锁&#xff09;3、为什么使用CAS&#xff0c;不用synchronized&#xff1f;4、CAS缺点5、ABA问题&#xff0c;原子引用更新&#xff1f; 前言 对于CAS部分…

『精』Vue 组件如何模块化抽离Props

『精』Vue 组件如何模块化抽离Props 文章目录 『精』Vue 组件如何模块化抽离Props一、为什么要抽离Props二、选项式API方式抽离三、组合式API方式抽离3.1 TypeScript类型方式3.2 文件分离方式3.3 对文件分离方式优化 参考资料&#x1f498;推荐博文&#x1f357; 一、为什么要抽…

系列六、Mybatis的一级缓存

一、概述 Mybatis一级缓存的作用域是同一个SqlSession&#xff0c;在同一个SqlSession中执行两次相同的查询&#xff0c;第一次执行完毕后&#xff0c;Mybatis会将查询到的数据缓存起来&#xff08;缓存到内存中&#xff09;&#xff0c; 第二次执行相同的查询时&#xff0c;会…