各种不同几何形状布局布阵下的GDOP相对值图

%布阵方式
clc; clear all; close all;
%%（1）直线形布站
s(:,:,1) = [-90 -60 -30 0 30 60 90;...
0 0 0 0 0 0 0]*1e3;
%%（2）楔形布站
s(:,:,2) = [-45 -30 -15 0 15 30 45;...
-90*cos(pi/6) -60*cos(pi/6) -30*cos(pi/6) 0 -30*cos(pi/6) -60*cos(pi/6) -90*cos(pi/6)]*1e3;
%%（3）工字形布站
s(:,:,3) = [-30 0 30 0 -30 0 30;...
30 30 30 0 -30 -30 -30]*1e3;
%%（4）梯形布站
s(:,:,4) = [-30 0 30 -45 -15 15 45;...
0 0 0 -30*cos(pi/6) -30*cos(pi/6) -30*cos(pi/6) -30*cos(pi/6)]*1e3;
%%（5）圆形布站
r = 3.45714e4;
theta0=2*pi/7;
for i = 1:7%七部
s(:,i,5)=[r*cos(theta0*i);r*sin(theta0*i)];
end
%%（6）十字形布站
s(:,:,6) = [-60 -30 0 30 60 0 0;...
0 0 0 0 0 30 -30]*1e3;
index_s = 2;%步战形式编号

sd_angle = 2*pi/180; %方位角标准差
sd_zhanzhi = 10; %站址标准差
Qa = eye(size(s,2))*(sd_angle)^2; %方位角方差
% Qb = eye(size(s,2)*2)*sd_zhanzhi^2;
x = -400e3:10e3:400e3; y = -400e3:10e3:400e3;
for i = 1:81
for j=1:81
u = [x(i);y(j)];
[rms1(j,i), ~] = crlb_AOA(u,s(:,:,index_s),Qa);%仅有测角误差时的CRLB的绝对值
rms2(j,i) = rms1(j,i)/norm(u);
% [rms1(j,i), ~] = crlb_AOA(u,s(:,:,index_s),Qa,Qb);
% rms2(j,i) = rms1(j,i)/norm(u);
% diff_rms(j,i) = crlb_AOA1(u,s(:,:,1),Qa,Qb);
end
end
figure(1)
switch index_s
case 1
title('直线形布阵时的GDOP图');
Q=[5 10 20 30 50 70 90 120 300]; %等位线标识数值
case 2
title('楔形布阵时的GDOP图');
Q=[5 10 20 30 50 70 90 120 300];
case 3
title('工字形布阵时的GDOP图');
Q=[5 10 20 30 50 70 90 120 300];
case 4
title('梯形布阵时的GDOP图');
Q=[5 10 20 30 50 70 90 120 300];
case 5
% title('不存在平台导航误差且圆形布阵时的GDOP图');
% Q=[1 3 5 7 9 12 15 17];
title('圆形布阵时的GDOP图');
Q=[5 10 20 30 50 70 90 120 300];
case 6
title('十字形布阵时的GDOP图');
Q=[5 10 20 30 50 70 90 120 300];
end
hold on;
[c,handle]=contour(x*1e-3,y*1e-3,rms1*1e-3,Q);
clabel(c,handle);hold on; %等位线标识
plot(s(1,:,index_s)*1e-3,s(2,:,index_s)*1e-3,'rd')
xlabel('x/km'); ylabel('y/km');

figure(2)
switch index_s
case 1
title('直线形布阵时的GDOP相对值图');
Q=[0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.12 0.15 0.17 0.2 0.3 0.5];
case 2
title('楔形布阵时的GDOP相对值图');
Q=[0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.12 0.15 0.17 0.2 0.3 0.5];
case 3
title('工字形布阵时的GDOP相对值图');
Q=[0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.12 0.15 0.17 0.2 0.3 0.5];
case 4
title('梯形布阵时的GDOP相对值图');
Q=[0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.12 0.15 0.17 0.2 0.3 0.5];
case 5
% title('不存在平台导航误差且圆形布阵时的GDOP相对值图');
% Q=[0.005 0.01 0.018 0.025 0.045];
title('圆形布阵时的GDOP相对值图');
Q=[0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.12 0.15 0.17 0.2 0.3 0.5];
case 6
title('十字形布阵时的GDOP相对值图');
Q=[0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.12 0.15 0.17 0.2 0.3 0.5];
end
hold on;[c,handle]=contour(x*1e-3,y*1e-3,rms2,Q);
clabel(c,handle);hold on;
plot(s(1,:,index_s)*1e-3,s(2,:,index_s)*1e-3,'rd')
xlabel('x/km'); ylabel('y/km');

D107

各种不同几何形状布局布阵下的GDOP相对值图相关推荐

各种不同几何形状布局布阵下的GDOP相对值图matlab仿真
目录 1.算法描述 2.仿真效果预览 3.MATLAB核心程序 4.完整MATLAB 1.算法描述全球四大导航卫星系统包括:美国GPS.俄罗斯GLONASS(格洛纳斯).中国北斗和欧洲Galileo ...
opencv 检测几何图形_使用OpenCV + ConvNets检测几何形状
opencv 检测几何图形 A simple yet powerful pipeline for detecting shapes in scanned documents 一个简单而强大的管道,用于 ...
create_metrology_model创建测量几何形状所需的数据结构（原理）
目录 create_metrology_model(算子) 描述二维计量的基本原理创建计量模型数据结构提供近似值修改模型参数修改对象参数对齐计量模型应用测量访问结果清理记忆注意参 ...
20款华丽的几何形状字体【免费下载】
这里手机的字体使用几何形状设计.流畅简洁的线条,完美的圆形的角度建立一个完整性的设计感.使用几何形状生成出每一个优雅而现代的字母.这些字体可以用于标题和正文.由于他们的设计适合任何干净简约设计,因此很 ...
Geometry（几何对象定义空间位置和关联几何形状）
Geometry 描述几何对象定义空间位置和关联几何形状. 讨论在许多地理处理工作流中,您可能需要使用坐标和几何信息运行特定操作,但不一定想经历创建新(临时)要素类.使用光标填充要素类.使用要素类 ...
电池的几何形状和标准尺寸
除了可重复使用性和电池化学性能,电池的包装方式也是一个重要因素,因为它应该能够在不损害安全标准的情况下安装在设备中.电池的几何形状和大小在其应用于电路中也起着至关重要的作用. 电池从18世纪就开始使用 ...
2020年最值得关注的5个logo设计趋势：简单的几何形状logo
LOGO设计是图形设计中一块非常重要的领域,logo设计趋势相对于图形设计会更加稳定,每年变化不大.主要原因是logo的使用时间很长,如果不是企业做重大革新,一般不会重新设计logo,在2019年那些 ...
漫步线性代数二——线性方程的几何形状
漫步线性代数二--线性方程的几何形状 2016年08月15日 23:10:10 会敲键盘的猩猩阅读数:1818 几何形状理解这个主题的方法是举例说明.我们以两个极其简单的方程开始,可以说大家在没有 ...
Opencv-Python图像像素均值、方差，绘制几何形状
图像像素均值.方差计算数组元素的均值和方差函数cv.meanStdDev().传入数组数据,返回数组数据的均值和方差.在图像数组中,可以计算出图像数组的像素均值和方差. 使用np.min和np.ma ...

各种不同几何形状布局布阵下的GDOP相对值图

各种不同几何形状布局布阵下的GDOP相对值图相关推荐

最新文章

热门文章