function?[?Lx0y0]?=?Computer_L（?Points?）
%功能：計算L系數的近似值和精確值
%輸入：Points為控制點的像點坐標及物方坐標
??????%格式為：點號?x?y?X?Y?Z
%輸出：L系數的精確值，像主點坐標想x0y0

%計算L系數近似值
X=zeros（2*size（Points1）11）;
B=zeros（2*size（Points1）1）;
for?i=1:size（Points1）;
????X（2*i-11:3）=Points（i4:6）;
????X（2*i-14）=1;
????X（2*i-19:11）=Points（i2）*Points（i4:6）;
????X（2*i5:8）=X（2*i-11:4）;
????X（2*i9:11）=Points（i3）*Points（i4:6）;
????B（2*i-11）=Points（i2）;
????B（2*i1）=Points（i3）;
end
L=X（1:11:）\（-B（1:111））;
%計算像主點近似值x0y0
x0=-（L（1）*L（9）+L（2）*L（10）+L（3）*L（11））/（L（9）*L（9）+L（10）*L（10）+L（11）*L（11））;
y0=-（L（5）*L（9）+L（6）*L（10）+L（7）*L（11））/（L（9）*L（9）+L（10）*L（10）+L（11）*L（11））;
%計算L系數精確值
M=zeros（2*size（Points1）12）;
W=zeros（2*size（Points1）1）;
n=0;???????????????????????????%迭代次數統計
while?n<10;???%L迭代次數為10
????for?i=1:6;
????????A=L（9）*Points（i4）+L（10）*Points（i5）+L（11）*Points（i6）+1;
????????M（2*i-11:4）=[Points（i4:6）?1]/A;
????????M（2*i-19:11）=Points（i2）*?M（2*i-11:3）;
????????R=（Points（i2）-x0）^2+（Points（i3）-y0）^2;
????????M（2*i-112）=（Points（i2）-x0）*R;
????????M（2*i5:8）=M（2*i-11:4）;
????????M（2*i9:11）=Points（i3）*?M（2*i-11:3）;
????????M（2*i12）=（Points（i3）-y0）*R;
????????W（2*i-11）=Points（i2）/A;
????????W（2*i1）=Points（i3）/A;
????end
????M=-M;
????L=（M‘*M）\M‘*W;
????x0=-（L（1）*L（9）+L（2）*L（10）+L（3）*L（11））/（L（9）*L（9）+L（10）*L（10）+L（11）*L（11））;
????y0=-（L（5）*L（9）+L（6）*L（10）+L（7）*L（11））/（L（9）*L（9）+L（10）*L（10）+L（11）*L（11））;
????n=n+1;
end
end

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

?????文件???????1570??2018-06-12?19:25??DLT的MATLAB實現\Computer_L.m

?????文件???????3492??2018-06-12?19:37??DLT的MATLAB實現\Computer_UnknownPoint.m

?????文件????????311??2018-06-27?22:54??DLT的MATLAB實現\DLT.m

?????文件????????330??2018-06-12?21:59??DLT的MATLAB實現\left_ran2.txt

?????文件????????328??2018-06-12?18:30??DLT的MATLAB實現\right_ran2.txt

?????文件????????499??2018-06-12?12:44??DLT的MATLAB實現\test1.txt

?????文件???????1466??2018-06-12?22:51??DLT的MATLAB實現\test2.txt

?????文件???????1883??2018-06-12?19:02??DLT的MATLAB實現\test3.txt

?????目錄??????????0??2019-02-26?17:47??DLT的MATLAB實現

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

?????????????????9879????????????????????9