雷达图形仿真 | 如何把距离-方位矩阵谱转为扇形谱？

本文编辑：调皮哥的小助理

大家好，最近看到有读者提了一个问题，就是用MATLAB如何绘制雷达数据距离-方位的扇形图，如下所示：

通常距离-方位谱，我们都习惯绘制为矩形图，而不是扇形图，但是呢，有的时候需要绘制扇形图，比如论文或者给客户演示的界面中。那么针对这个问题如何解决呢？

调皮哥为大家专门做了程序，主要就是为了实用性，后续大家也不必为这种图形犯难，下面是效果图：

（图1：距离-方位扇形图）

这个界面只是初步效果，如果大家看起来不太满意可以随时修改后面我给出的程序，直到满足你自己的需求。程序很简单，就几行代码，核心代码如下：

%极坐标转为直角坐标形式
figure;%公众号：调皮连续波
X = R’*cosd(ang_ax); Y = R’*sind(ang_ax); %
clf;
pcolor(X,Y,pow2db(abs(range_az).^2));
axis equal tight %x轴单位刻度与 y轴单位刻度长度相等 , 最能体现出实际的曲线 ;
shading interp %着色，使色彩平滑过渡
axis off
initialAz = -90; endAz = 90; %标注文字
text((max®+10)*cosd(initialAz),(max®)*sind(initialAz),…
[num2str(initialAz) ‘^o’]);
text((max®+10)*cosd(endAz),(max®)*sind(endAz),…
[num2str(endAz) ‘^o’]);
text((max®+10)*cosd(0),(max®)*sind(0),[num2str(0) ‘^o’]);
colorbar; %热力值

这个程序片段更新在了《开源代码 | FMCW-MIMO雷达仿真MATLAB（2022年的最后一篇）》这篇文章的程序之中了，下载后运行即可看见，感兴趣的读者赶紧试试吧。

采用Pyhton绘制扇形图的代码可以参考这篇文章：https://github.com/moodoki/radical_sdk，效果如下图所示：

核心代码如下所示：

from radicalsdk.geometry import PolarToCartesianWarp
p2c = PolarToCartesianWarp() #投影到笛卡尔坐标系
cartesian_radar = p2c(np.abs(rf.range_azimuth_capon)[np.newaxis, …, np.newaxis])
plt.figure()
with np.errstate(divide=‘ignore’):
plt.imshow(np.log(cartesian_radar[0, …]))
plt.show()