空间误差分析：统一的应用导向处理（Matlab代码实现）

 👨‍🎓个人主页：研学社的博客 

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

本文关于一维、二维和三维空间误差测量的研究。

📚2 运行结果

部分代码：

clear
ic=2;   % skip step 1

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if ic==1
% part I
%        use y=sqrt(x/n), x has chi-square distribution with deg n
clear global n p
global  n p

q=[0.5 .8 .9 .95];
l=[(0.1:0.025:0.2) (0.3:.1:.9) (1:9) (10:10:100)]';
NL = length(l);
NQ = length(q);

for jj=1:NQ
  p=q(jj);
  x =1;      %initial guess for iteration
% x =0.01;   %initial guess for fsolve
for ii=1:NL
  n=l(ii);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% use fsolve from optimization toolbox %%
%%  xs=fsolve('x2disc2',x);
%%  xo= xs;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%% iteration %%%%%%%%%%%%

      p0=x2cdf(x);
   for j=1:2:9
      while (p0 > p)
      x=x-(0.1)^j;
      p0=x2cdf(x);
      end

      while (p0 < p)
      x=x+(0.1)^j;
      p0=x2cdf(x);
      end
   end
     xo=x;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  r(ii,jj)=sqrt(xo/n)
  pause
end
end

DAT =[l  r];
save \mfile\chi.mat  l r

fname='chi.dat'
fid=fopen(fname,'wt')
fprintf(fid,'%8.2f %8.4f %8.4f %8.4f %8.4f\n',DAT')
fclose(fid)
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if ic==2
% part II
load chi.mat
NQ =4
adj= [ 0.7  0.5  0.15  -0.1];

semilogx(l,r)
xlabel('n')
ylabel(' R(p)/RMS ')
% title(' R(p)/RMS vs n,  p=0.5, 0.8, 0.9, 0.95')
grid

for kk=1:NQ
if kk==1   p=0.5
   elseif kk==2 p=0.8
   elseif kk==3 p=0.9
   elseif kk==4 p=0.95
end

   text(.5,r(5,kk)+adj(kk),['p=' num2str(100*p) '%'])
end
axis('square')
figure(1)
pause
pf3('f62')
end
 

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]David Hsu (2023). Spatial Error Analysis: A Unified Application-Oriented Treatment 

🌈4 Matlab代码实现