注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

如何绘制GDOP图?  

2012-03-20 13:08:53|  分类: 绘图 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

clc;
clear;
x1=-30;x2=-30;x3=60;xt=-60;
y1=51.96;y2=-51.96;y3=0;yt=0;
z1=0;z2=0;z3=0;zt=0.1;
z=10;
y=-100:1:100;x=-100:1:100;
 
for i=1:201
for j=1:201
m=x(i);
n=y(j);
r1=((m-x1).^2+(n-y1).^2+(z-z1).^2).^(1/2);
r2=((m-x2).^2+(n-y2).^2+(z-z2).^2).^(1/2);
r3=((m-x3).^2+(n-y3).^2+(z-z3).^2).^(1/2);
rt=((m-xt).^2+(n-yt).^2+(z-zt).^2).^(1/2);
c11=(m-x1)/r1;c21=(m-x2)/r2;c31=(m-x3)/r3; ct1=(m-xt)/rt;    % 各个方向余弦
c12=(n-y1)/r1;c22=(n-y2)/r2;c32=(n-y3)/r3;   ct2=(n-yt)/rt;    
c13=(z-z1)/r1;c23=(z-z2)/r2;c33=(z-z3)/r3;     ct3=(z-zt)/rt;
c=[(ct1+c11) (ct2+c12) (ct3+c13);(ct1+c21) (ct2+c22) (ct3+c23);(ct1+c31) (ct2+c32) (ct3+c33)];   % 此处加上了ct项表示什么意义?
 
b=inv(c'*c)*(c');    % 最小二乘中的那个矩阵?
b11=b(1);
b12=b(4);
b13=b(7);
b21=b(2);
b22=b(5);
b23=b(8);
b31=b(3);
b32=b(6);
b33=b(9);
 
sigmap=0.1;
sigmas=0.005;
eta=0.3;
 
sigma11=sigmap.^2+sigmas.^2*2;
sigma22=sigmap.^2+sigmas.^2*2;
sigma33=sigmap.^2+sigmas.^2*2;
sigma12=eta*sigmap.^2+sigmas.^2;
sigma13=eta*sigmap.^2+sigmas.^2;
sigma21=eta*sigmap.^2+sigmas.^2;
sigma23=eta*sigmap.^2+sigmas.^2;
sigma31=eta*sigmap.^2+sigmas.^2;
sigma32=eta*sigmap.^2+sigmas.^2;
 
sigmax2=b11*b11*sigma11+b11*b12*sigma12+b11*b13*sigma13+b12*b11*sigma21+b12*b12*sigma22+b12*b13*sigma23+b13*b11*sigma31+b13*b12*sigma32+b13*b13*sigma33;
sigmay2=b21*b21*sigma11+b21*b22*sigma12+b21*b23*sigma13+b22*b21*sigma21+b22*b22*sigma22+b22*b23*sigma23+b23*b21*sigma31+b23*b22*sigma32+b23*b23*sigma33;
sigmaz2=b31*b31*sigma11+b31*b32*sigma12+b31*b33*sigma13+b32*b31*sigma21+b32*b32*sigma22+b32*b33*sigma23+b33*b31*sigma31+b33*b32*sigma32+b33*b33*sigma33;
 
gdop(j,i)=(sigmax2+sigmay2+sigmaz2).^(1/2);
 
end
end
figure(1);
[c,handle]=contour(gdop,25);%[c,handle]=contour(gdop,25);
clabel(c,handle);
xlabel('x方向(单位:km)')
ylabel('y方向(单位:km)')
title('GDOP图')
text(1,10,'{\T1}')如何绘制GDOP图? - 阿英 - Mr.Right

  评论这张
 
阅读(983)| 评论(3)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016