【正文】 39。*V/2)。 fai0=fai0+XX(4)。 a26=x(h)。 O=a1*(X(h)Xs0)+b1*(Y(h)Ys0)+c1*(Z(h)Zs0)。for v=1:+inf。s=0。⑧解求外方位元素。 實習原理利用至少三個已知地面控制點的坐標、與其影像上對應的三個像點的影像坐標、根據共線方程，反求該像片的外方位元素。將所求得的外方位元素的改正數與規(guī)定的限差比較，通常對φ，ω，κ的改正數△φ，△ω，△κ給予限差，迭代結束。 j=i+1。 a3=sin(fai0)*cos(omig0)。 x1=f*O/Q。 lh=[lx ly]39。 R=[a1,a2,a3。m=[m(1),m(2),m(3),m(4)/pi*180*3600,m(5)/pi*180*3600,m(6)/pi*180*3600]。\n旋轉矩陣R：\n39。Ys=%g\t39。,m0)。fprintf(fp,39。,mZs)。κ=%g\n39。fprintf(fp,39。 fprintf(fp,39。m0=sqrt(V39。 Zs0=Zs0+XX(3)。 a25=f*cos(ka0)y(h)/f*(x(h)*sin(ka0)+y(h)*cos(ka0))。 for h=1:4。ka0=0。Z=I(:,5)39。⑦計算法方程的系數矩陣ATA與常數項ATL，組成法方程式。共線條件方程式： 將共線式線性化并取一次小值項得： 解算過程①獲取已知數據。否則用新的近似值重復④～⑧步驟的計算，直到滿足要求為止。 sij=sqrt((x(i)x(j))^2+(y(i)y(j))^2)。 b1=cos(om