数码相机的定位2008年全国大学生数学建模程序运行

wy722526

model:
!五个像点的坐标;
A1y=72/3.78;
A1z=68/3.78;
B1y=69/3.78;
B1z=32/3.78;
C1y=63/3.78;
C1z=-47/3.78;
D1y=-43/3.78;
D1z=81/3.78;
E1y=-43/3.78;
E1z=-29/3.78;
v=-1577/3.78;
@free(A1y);@free(A1z);@free(B1y);@free(B1z);@free(C1y);@free(C1z);@free(D1y);@free(D1z);@free(E1y);@free(E1z);
@free(v);

init:
Ax=1200;Cx=1000;
endinit
!以四条边的长度尽量接近100为目标,含义是尽量保持标靶不边形;

min=(((ax-cx)^2+(ay-cy)^2+(az-cz)^2)^0.5-100)^2+
(((cx-dx)^2+(cy-dy)^2+(cz-dz)^2)^0.5-100)^2+
(((ex-dx)^2+(ey-dy)^2+(ez-dz)^2)^0.5-100)^2+
(((ex-ax)^2+(ey-ay)^2+(ez-az)^2)^0.5-100)^2;
!(((ax-dx)^2+(ay-dy)^2+(az-dz)^2)^0.5-1.414*100)^2+
(((ex-cx)^2+(ey-cy)^2+(ez-cz)^2)^0.5-1.414*100)^2+
(((cx-ax)*(cx-dx)+(cy-ay)*(cy-dy)+(cz-az)*(cz-dz)))^2+
(((cx-dx)*(ex-dx)+(cy-dy)*(ey-dy)+(cz-dz)*(ez-dz)))^2+
(((ex-ax)*(ex-dx)+(ey-ay)*(ey-dy)+(ez-az)*(ez-dz)))^2;
!3(cx-ax)=10(bx-ax);
!3(cy-ay)=10(by-ay);
!3(cz-az)=10(bz-az);
!(((ax-bx)^2+(ay-by)^2+(az-bz)^2)^0.5-30)^2+
!边互相垂直;
!((Cx-ax)*(dx-cx)+(Cy-ay)*(dy-cy)+(Cz-az)*(dz-cz))=0;
((cx-ax)*(ex-ax)+(cy-ay)*(ey-cy)+(cz-az)*(ez-az))=0;
!((ex-dx)*(ax-ex)+(ey-dy)*(ay-ey)+(ez-dz)*(az-ez))=0;
!((ex-dx)*(cx-dx)+(ey-dy)*(cy-dy)+(ez-dz)*(cz-dz))=0;
!点A与像A的对应关系;
!ax^2*(A1y^2+A1z^2)=v^2*(ay^2+az^2);
ay^2*(A1y^2+v^2)=A1y^2*(ay^2+ax^2);
az^2*(A1y^2+v^2)=A1z^2*(ay^2+ax^2);
!点C与像C的对应关系;
cy^2*(C1y^2+v^2)=C1y^2*(cy^2+cx^2);
cz^2*(C1y^2+v^2)=A1z^2*(cy^2+cx^2);
!点D与像D的对应关系;
dy^2*(D1y^2+v^2)=D1y^2*(dy^2+dx^2);
dz^2*(D1y^2+v^2)=D1z^2*(dy^2+dx^2);
!点E与像E的对应关系;
ey^2*(E1y^2+v^2)=E1y^2*(ey^2+ex^2);
ez^2*(E1y^2+v^2)=E1z^2*(ey^2+ex^2);

@free(ax);@free(ay);@free(az);!@free(bx);!@free(by);!@free(bz);@free(cx);@free(cy);@free(cz);
@free(dx);@free(dy);@free(dz);@free(ex);@free(ey);@free(ez);
end
没有结果，请高手指点，是新手~！

wy722526

最好大家给个结果！自己做了一下修改！