function?[X0Y0Z0uvwk]=bursa7para（oldXYZnewXYZ）
%布爾沙模型求七參數
%輸入參數：oldXYZnewXYZ均為空間直角坐標
%輸出7參數：3個平移參數，3個旋轉參數，1個尺度參數
X1=oldXYZ（:1）;
Y1=oldXYZ（:2）;
Z1=oldXYZ（:3）;
X2=newXYZ（:1）;
Y2=newXYZ（:2）;
Z2=newXYZ（:3）;
m=length（oldXYZ）;%重合點個數
if（m<3）
????return;
end
P=eye（3*m）;
B=[ones（m1）?zeros（m1）?zeros（m1）?X1?zeros（m1）?-Z1?Y1;
???zeros（m1）?ones（m1）?zeros（m1）?Y1?Z1?zeros（m1）?-X1;
???zeros（m1）?zeros（m1）?ones（m1）?Z1?-Y1?X1?zeros（m1）];
L=[X2;Y2;Z2];
N=B‘*P*B;
%M=B‘*P*L;
dx=（inv（N））*B‘*P*L;
X0=dx（11）;
Y0=dx（21）;
Z0=dx（31）;
u=dx（51）/dx（41）;
v=dx（61）/dx（41）;
w=dx（71）/dx（41）;
k=dx（41）-1;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????687??2017-03-31?19:55??布爾沙求七參數\bursa7para.asv

?????文件????????706??2017-03-31?19:56??布爾沙求七參數\bursa7para.m

?????文件????????532??2017-03-31?20:00??布爾沙求七參數\my7para.asv

?????文件????????545??2017-03-31?20:04??布爾沙求七參數\my7para.m

?????文件????????150??2017-03-31?19:03??布爾沙求七參數\newXYZ.txt

?????文件????????150??2017-03-31?18:59??布爾沙求七參數\oldXYZ.txt

?????文件????????164??2017-03-31?20:05??布爾沙求七參數\result.txt

?????文件?????153857??2017-03-31?12:19??布爾沙求七參數\基于MATLAB的布爾莎模型七參數解算實現_朱小美.caj

?????文件????????134??2017-03-31?21:38??布爾沙求七參數\說明.txt

?????目錄??????????0??2017-09-17?18:54??布爾沙求七參數

-----------?---------??----------?-----??----

???????????????156925????????????????????10