毫米波雷达图解算法原理（基于TI雷达）_毫米波雷达算法-程序员宅基地

技术标签： 算法  

毫米波雷达数据处理原理

前言

原始数据bin文件包可以由结束语下第一个链接获取

由于网上可以找到很多毫米波雷达处理的流程都是基于抽象原理讲解。所以本文章将会基于实际例子进行数据转换的说明。
小编也是刚入门毫米波雷达，如若有讲解不对的地方还请指出。
在未来的几个月里，小编将会把此次处理过程制作成为动画发布，届时将把链接附上，所以文章中如有不对的地方还请大神帮助指出。
11.17日笔记：（小编太懒了，就先不出动画了）
本篇文章基于TI AWR1642雷达生成的.bin文件数据进行说明，同样适用于其他TI系列雷达。

基础

我们将对毫米波雷达的基本参数进行说明，在后续的使用中将会涉及到它们：
（以下名词前4个由下包含至上的结构进行排列，5、6为硬件基础）

  1. IQ正交采样 复数形式，采样方式决定基本单位大小，正交采样为32bit即4个字节（实数16bit，2个字节，虚数16bit，2个字节）（所以读取第xx帧的时候要多乘个2）；
  2. 采样点数 sampls ，每次脉冲的采样点数；
  3. 脉冲 chirps 每次脉冲即是一次发射，雷达需要发射n次脉冲以达到1帧；
  4. 帧数 每帧包含n次脉冲。
  5. 发射天线 TX发射天线数量；
  6. 接收天线 RX 接收天线数量；

此次文件设定采样的参数：
2 字节 IQ正交采样
256 sampls 采样点数
128 chirps 脉冲数
1 TX 发射天线
4 RX 接收天线

bin文件解读

1642SRR2m.bin文件为采样回读的文件，大小为125 MB (131,072,000 字节)。
其包含了250帧数据。
1帧解析：4字节x256samplsx128chirpsx4RXx1TX=524288字节
即524288*250=131，072，000字节

读取出的结构（实际是根据操作过程猜测的数据，咱也不知道为啥是这样的？？）（目前认为极有可能是TI采用的硬件加速的独用ram有关）：
实数1
实数2
虚数1
虚数2
实数3
实数4
虚数3
虚数4
…
…
…

将bin文件读入到matlab中可以看到一共存在65536000行数据
将bin文件读入到matlab中可以看到一共存在65536000行数据
1帧：524288字节（由于实数虚数的长度为16bit而不是8bit，所以每行数据的单位是16bit而不是8bit）
1帧：262144行

以下我们取1帧进行操作：

在读取出bin原始数据后和对其进行计算前需要将其转换成合适的矩阵格式以方便计算。

对数据矩阵进行操作前的转换——

由 262144行1列数据经过实数虚数结合后：
变为 1行131072列数据：
在这里插入图片描述

对数据进行矩阵排列，变为：chirps/TX行，samples·RX·TX列：

再变换为以发送通道作为行数的阵列（RX*TX）：

最终变换（三维数组，samples行*chirp列）：

开始计算结果矩阵

排列好的数据需要经过一维FFT（距离维FFT），二维FFT（速度维FFT），以及角度维FFT计算。
FFT:快速傅里叶变换

下面将以matlab语言进行展示程序

一维FFT（距离）

可以理解为对每组chirp进行FFT（个人理解）：
N为FFT点数，取N=sampls。
在这里插入图片描述

二维FFT（速度）

可以理解为对每组sampls进行FFT（个人理解）：
N为FFT点数，取N=chirp。
fft_shift：将零频点移动到频谱中央（因为速度存在正负）
在这里插入图片描述

角度维FFT

可以理解为对每组RX*TX通道进行FFT（个人理解）：
N为FFT点数，取N=180。
fft_shift：将零频点移动到频谱中央（因为角度存在正负）
在这里插入图片描述

（假设利用结果已经获取目标）对目标进行计算，得出距离速度角度信息

已知获得的目标位置[row,col,pag]（可以利用峰值计算单个目标，可以利用CFAR计算多个目标）

以下参数为用户雷达设定的参数，全部已知（部分参数请参 照参考文章1 代码部分）
fb = ((row-1) × fs)/N; %差拍频率 fs:采样率/N：一维FFT点数
fd = (col-M/2-1)/(M × Tc); %多普勒频率 M：二维FFT点数/Tc：chirp总周期
fw = (pag-Q/2-1)/Q; %空间频率 Q：角度FFT点数
R = c × (fb-fd)/2/K; %距离公式 ****
v = lambda × fd/2; %速度公式 lambda：波长
theta = asin(fw × lambda/d); %角度公式 d：天线阵列间距
angle = theta × 180/pi;
得出：
目标距离：R m
目标速度：v m/s
目标角度：angle °

结束

参考文章：
1.回波3DFFT处理（测距、测速、测角）：
此博主包含原始bin文件，需要用网盘下载。请支持博主一下
2.调皮连续波：TI毫米波雷达IWR6843ISK室内人员检测算法（附MATLAB程序和数据）
此博主有公众号，大家可以关注。

如有不对之处还请大家帮忙指正
1.11.17-感谢某位研究生纠正前方TX/RX搞反的情况
2.11.18-之前一维二维FFT的对象搞反了，这里请注意冒号的用法。
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45830602/article/details/126966809

智能推荐

攻防世界_难度8_happy_puzzle_攻防世界困难模式攻略图文-程序员宅基地

文章浏览阅读645次。这个肯定是末尾的IDAT了，因为IDAT必须要满了才会开始一下个IDAT，这个明显就是末尾的IDAT了。，对应下面的create_head()代码。，对应下面的create_tail()代码。不要考虑爆破，我已经试了一下，太多情况了。题目来源：UNCTF。_攻防世界困难模式攻略图文

达梦数据库的导出（备份）、导入_达梦数据库导入导出-程序员宅基地

文章浏览阅读2.9k次，点赞3次，收藏10次。偶尔会用到，记录、分享。1. 数据库导出1.1 切换到dmdba用户su - dmdba1.2 进入达梦数据库安装路径的bin目录，执行导库操作　　导出语句：./dexp cwy_init/[email protected]:5236 file=cwy_init.dmp log=cwy_init_exp.log　注释：　　 cwy_init/init_123..._达梦数据库导入导出

js引入kindeditor富文本编辑器的使用_kindeditor.js-程序员宅基地

文章浏览阅读1.9k次。1. 在官网上下载KindEditor文件，可以删掉不需要要到的jsp，asp，asp.net和php文件夹。接着把文件夹放到项目文件目录下。2. 修改html文件，在页面引入js文件：<script type="text/javascript" src="./kindeditor/kindeditor-all.js"></script><script type="text/javascript" src="./kindeditor/lang/zh-CN.js"_kindeditor.js

STM32学习过程记录11——基于STM32G431CBU6硬件SPI+DMA的高效WS2812B控制方法-程序员宅基地

文章浏览阅读2.3k次，点赞6次，收藏14次。SPI的详情简介不必赘述。假设我们通过SPI发送0xAA，我们的数据线就会变为10101010，通过修改不同的内容，即可修改SPI中0和1的持续时间。比如0xF0即为前半周期为高电平，后半周期为低电平的状态。在SPI的通信模式中，CPHA配置会影响该实验，下图展示了不同采样位置的SPI时序图[1]。CPOL = 0，CPHA = 1：CLK空闲状态 = 低电平，数据在下降沿采样，并在上升沿移出CPOL = 0，CPHA = 0：CLK空闲状态 = 低电平，数据在上升沿采样，并在下降沿移出。_stm32g431cbu6

计算机网络-数据链路层_接收方收到链路层数据后,使用crc检验后,余数为0,说明链路层的传输时可靠传输-程序员宅基地

文章浏览阅读1.2k次，点赞2次，收藏8次。数据链路层习题自测问题1.数据链路(即逻辑链路)与链路(即物理链路)有何区别?“电路接通了”与”数据链路接通了”的区别何在?2.数据链路层中的链路控制包括哪些功能?试讨论数据链路层做成可靠的链路层有哪些优点和缺点。3.网络适配器的作用是什么?网络适配器工作在哪一层?4.数据链路层的三个基本问题(帧定界、透明传输和差错检测)为什么都必须加以解决？5.如果在数据链路层不进行帧定界，会发生什么问题？6.PPP协议的主要特点是什么？为什么PPP不使用帧的编号？PPP适用于什么情况？为什么PPP协议不_接收方收到链路层数据后,使用crc检验后,余数为0,说明链路层的传输时可靠传输

软件测试工程师移民加拿大_无证移民，未受过软件工程师的教育（第1部分）-程序员宅基地

文章浏览阅读587次。软件测试工程师移民加拿大 无证移民，未受过软件工程师的教育(第1部分) (Undocumented Immigrant With No Education to Software Engineer(Part 1))Before I start, I want you to please bear with me on the way I write, I have very little gen...

随便推点

Thinkpad X250 secure boot failed 启动失败问题解决_安装完系统提示secureboot failure-程序员宅基地

文章浏览阅读304次。Thinkpad X250笔记本电脑，装的是FreeBSD，进入BIOS修改虚拟化配置（其后可能是误设置了安全开机），保存退出后系统无法启动，显示：secure boot failed ，把自己惊出一身冷汗，因为这台笔记本刚好还没开始做备份.....根据错误提示，到bios里面去找相关配置，在Security里面找到了Secure Boot选项，发现果然被设置为Enabled，将其修改为Disabled ，再开机，终于正常启动了。_安装完系统提示secureboot failure

C++如何做字符串分割（5种方法）_c++ 字符串分割-程序员宅基地

文章浏览阅读10w+次，点赞93次，收藏352次。1、用strtok函数进行字符串分割原型： char *strtok(char *str, const char *delim);功能：分解字符串为一组字符串。参数说明：str为要分解的字符串，delim为分隔符字符串。返回值：从str开头开始的一个个被分割的串。当没有被分割的串时则返回NULL。其它：strtok函数线程不安全，可以使用strtok_r替代。示例：//借助strtok实现split#include <string.h>#include <stdio.h&_c++ 字符串分割

2013第四届蓝桥杯 C/C++本科A组 真题答案解析_2013年第四届c a组蓝桥杯省赛真题解答-程序员宅基地

文章浏览阅读2.3k次。1 .高斯日记 大数学家高斯有个好习惯：无论如何都要记日记。他的日记有个与众不同的地方，他从不注明年月日，而是用一个整数代替，比如：4210后来人们知道，那个整数就是日期，它表示那一天是高斯出生后的第几天。这或许也是个好习惯，它时时刻刻提醒着主人：日子又过去一天，还有多少时光可以用于浪费呢？高斯出生于：1777年4月30日。在高斯发现的一个重要定理的日记_2013年第四届c a组蓝桥杯省赛真题解答

基于供需算法优化的核极限学习机(KELM)分类算法-程序员宅基地

文章浏览阅读851次，点赞17次，收藏22次。摘要：本文利用供需算法对核极限学习机(KELM)进行优化，并用于分类。

metasploitable2渗透测试_metasploitable2怎么进入-程序员宅基地

文章浏览阅读1.1k次。一、系统弱密码登录1、在kali上执行命令行telnet 192.168.26.1292、Login和password都输入msfadmin3、登录成功，进入系统4、测试如下：二、MySQL弱密码登录：1、在kali上执行mysql –h 192.168.26.129 –u root2、登录成功，进入MySQL系统3、测试效果：三、PostgreSQL弱密码登录1、在Kali上执行psql -h 192.168.26.129 –U post..._metasploitable2怎么进入

Python学习之路：从入门到精通的指南_python人工智能开发从入门到精通pdf-程序员宅基地

文章浏览阅读257次。本文将为初学者提供Python学习的详细指南，从Python的历史、基础语法和数据类型到面向对象编程、模块和库的使用。通过本文，您将能够掌握Python编程的核心概念，为今后的编程学习和实践打下坚实基础。_python人工智能开发从入门到精通pdf

推荐文章

热门文章

相关标签