基于到达角的超声多频信号手势识别(Matlab代码实现）

           目录

💥1 概述

📚2 运行结果

🎉3 参考文献

👨‍💻4 Matlab代码

---

💥1 概述

DoA ：波达方向   AoA：波达角  

波达方向是指空间信号的到达方向（各个信号到达阵列参考阵元的方向角，简称波达方向），波达角与波达方向是一回事。

对于3D空间来说，DoA（或AoA）包括（θ，φ）。对于一个2D平面来说，DoA即为方位角θ。目前基于DoA的定位技术都是在2D平面的基础上进行的，此时DoA即为方位角θ。

📚2 运行结果

 

🎉3 参考文献

[1]曹轶超. 无线电定位算法研究及目标运动分析[D].东华大学,2010.

👨‍💻4 Matlab代码

主函数部分代码：

close all;
clear all;
clc;
% load data before plotting
%% Fig 1
load fig_8_3_data_3.mat;
% close all;
plot_ind = 7:17;
linewidth = 1.5;
markersize = 6;
figure;
semilogy(snrdB(plot_ind),sqrt(rvec1(plot_ind)),'-ob','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec2(plot_ind)),'--^c','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec3(plot_ind)),'-*r','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec4(plot_ind)),'-+m','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec5(plot_ind)),'-dg','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec6(plot_ind)),'>m-.','LineWidth',linewidth,'MarkerSize',markersize*1.2);
hold on; semilogy(snrdB(plot_ind), sqrt(crlb(plot_ind)),'--k','LineWidth',2);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec1(plot_ind)),'-ob','LineWidth',linewidth,'MarkerSize',markersize);

grid on;set(gca,'FontSize',13);
xlabel('{SNR [dB]}','Interpreter','Latex');
ylabel('{RMSE [Deg]}','Interpreter','Latex');

% set(gca,'FontSize',12)

% Zoomed plot ****************************************
bounds = [32 33 10^-1.35 10^-1.15]; % area to be zoomed
pos = [0.16 0.15 0.17 0.2];      % projected area [x, y, x_width, y_width]
vertex = [2 3];                 % vertex to be linked
p = gca;
% Calculate x,y points of zoomPlot
x1 = (pos(1)-p.Position(1))/p.Position(3)*diff(p.XLim)+p.XLim(1);
x2 = (pos(1)+pos(3)-p.Position(1))/p.Position(3)*diff(p.XLim)+(p.XLim(1));
y1 = 10^((pos(2)-(p.Position(2)))/(p.Position(4))*diff(log10(p.YLim))+log10(p.YLim(1)));
y2 = 10^(((pos(2)+pos(4)-p.Position(2))/p.Position(4))*diff(log10(p.YLim))+log10(p.YLim(1)));

rectangle('Position',[bounds(1) bounds(3) (bounds(2)-bounds(1)) (bounds(4)-bounds(3))]);
hold on;
offset = [0 -0.025 -0.0025 0];
if any(vertex==1)
    plot([bounds(1) x1], [bounds(4) y2 + offset(1)], '-.k'); % Line to vertex 1
end
if any(vertex==2)
    plot([bounds(2) x2], [bounds(4) y2 + offset(2)], '-.k'); % Line to vertex 2
end
if any(vertex==3)
    plot([bounds(2) x2], [bounds(3) y1 + offset(3)], '-.k'); % Line to vertex 4
end
if any(vertex==4)
    plot([bounds(1) x1], [bounds(3) y1 + offset(4)], '-.k'); % Line to vertex 3
end
% legend hide
z = axes('position',pos);

box on; % put box around new pair of axes
semilogy(snrdB(plot_ind),sqrt(rvec1(plot_ind)),'-ob','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec2(plot_ind)),'--^c','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec3(plot_ind)),'-*r','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec4(plot_ind)),'-+m','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec5(plot_ind)),'-dg','LineWidth',linewidth,'MarkerSize',markersize);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec6(plot_ind)),'>m-.','LineWidth',linewidth,'MarkerSize',markersize*1.2);
hold on;semilogy(snrdB(plot_ind),sqrt(crlb(plot_ind)),'--k','LineWidth',2);
hold on;semilogy(snrdB(plot_ind),sqrt(rvec1(plot_ind)),'-ob','LineWidth',linewidth,'MarkerSize',markersize);

axis(bounds);

set(gca,'xtick',[])
set(gca,'ytick',[])
set(gca,'xticklabel',{[]})
set(gca,'yticklabel',{[]})

% set(gcf,'position',[100,100,420*1.2,400*1.2])
% print -dpng -r600 sim1-3sensors.png

legend('PDP','MUSIC ($0.01^\circ$)','MLE ($0.01^\circ$)','2Q-order','Two-step','EM-ESPRIT ($2^\circ$)','CRLB','Interpreter','Latex');
grid on; set(gca,'FontSize',13)

set(gcf,'position', [100,100,350*1.5,350])