2008年全国大学生数学建模竞赛A题优秀论文

format long;
[p1txt1raw1]=xlsread（‘d:/1.xls‘）;
[p2txt2raw2]=xlsread（‘d:/2.xls‘）;
[p3txt3raw3]=xlsread（‘d:/3.xls‘）;
[p4txt4raw4]=xlsread（‘d:/4.xls‘）;
[p5txt5raw5]=xlsread（‘d:/5.xls‘）;

qd=ones（64）;
P=imread（‘picture.bmp‘）;
imshow（P）;
hold on;
%第1组 圆13
[min1 i]=min（p1（:1））;
qd（11）=i;
[min3 i]=min（p3（:1））;
qd（12）=i;
size1=size（p1）;
size3=size（p3）;
flag=1;
while flag>0
    flag=0;
    for i=1:size1（1）
         if p1（i1）<（p1（qd（11）1）-p3（qd（12）1））/（p1（qd（11）2）-p3（qd（12）2））*（p1（i2）-（p1（qd（11）2）））+p1（qd（11）1）
             qd（11）=i;
             flag=1;
         end
    end
    for i=1:size3（1）
         if p3（i1）<（p1（qd（11）1）-p3（qd（12）1））/（p1（qd（11）2）-p3（qd（12）2））*（p3（i2）-（p1（qd（11）2）））+p1（qd（11）1）
             qd（12）=i;
             flag=1;
         end
    end
end

for i=1:1024
    y（i）=（p1（qd（11）1）-p3（qd（12）1））/（p1（qd（11）2）-p3（qd（12）2））*（i-（p1（qd（11）2）））+p1（qd（11）1）;
end
plot（[1:1024]y）;
hold on;

%另一条切线
[max1 i]=min（p1（:1））;
qd（13）=i;
[max3 i]=min（p3（:1））;
qd（14）=i;
size1=size（p1）;
size3=size（p3）;
flag=1;
while flag>0
    flag=0;
    for i=1:size1（1）
         if p1（i1）>（p1（qd（13）1）-p3（qd（14）1））/（p1（qd（13）2）-p3（qd（14）2））*（p1（i2）-（p1（qd（13）2）））+p1（qd（13）1）
             qd（13）=i;
             flag=1;
         end
    end
    for i=1:size3（1）
         if p3（i1）>（p1（qd（13）1）-p3（qd（14）1））/（p1（qd（13）2）-p3（qd（14）2））*（p3（i2）-（p1（qd（13）2）））+p1（qd（13）1）
             qd（14）=i;
             flag=1;
         end
    end
end

for i=1:1024
    y（i）=（p1（qd（13）1）-p3（qd（14）1））/（p1（qd（13）2）-p3（qd（14）2））*（i-（p1（qd（13）2）））+p1（qd（13）1）;
end
plot（[1:1024]y）;
hold on;

%第2组 圆45
[min4 i]=min（p4（:1））;
qd（21）=i;
[min5 i]=min（p5（:1））;
qd（22）=i;
size4=size（p4）;
size5=size（p5）;
flag=1;
while flag>0
    flag=0;
    for i=1:size4（1）
         if p4（i1）<（p4（qd（21）1）-p5（qd（22）1））/（p4（qd（21）2）-p5（qd（22）2））*（p4（i2）-（p4（qd（21）2）））+p4（qd（21）1）
             qd（21）=i;
             flag=1;
         end
    end
    for i=1:size5（1）
         if p5（i1）<（p4（qd（21）1）-p5（qd（22）1））/（p4（qd（21）2）-p5（qd（22）2））*（p5（i2）-（p4（qd（21）2）））+p4（qd（21）1）
             qd（22）=i;
             flag=1;
         end
    end
end

for i=1:1024
    y（i）=（p4（qd（21）1）-p5（qd（22）1））/（p4（qd（21）2）-p5（qd（22）2））*（i-（p4（qd（21）2）））+p4（qd（21）1）;
end
plot（[1:1024]y）;
hold on;

%另一条切线
[max4 i]=min（p4（:1））;
qd（23）=i;
[max5 i]=min（p5（:1））;
qd（24）=i;
size4=size（p4）;
size5=size（p5）;
flag=1;
while flag>0
    flag=0;
    for i=1:size4（1）
         if p4（i1）>（p4（qd（23）1）-p5（qd（24）1））/（p4（qd（23）2）-p5（qd（24）2））*（p4（i2）-（p4（qd（23）2）））+p4（qd（23）1）
             qd（23）=i;
             flag=1;
         end
    end
    for i=1:size5（1）
         if p5（i1）>（p4（qd（23）1）-p5（qd（24）1））/（p4（qd（23）2）-p5（qd（24）2））*（p5（i2）-（p4（qd（23）2）））+p4（qd（23）1）
          

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     334756  2008-12-13 15:15  A2105373论文数码相机定位.pdf

     文件      25088  2008-09-20 10:20  A2105373数据边缘点集1.xls

     文件      24576  2008-09-20 10:20  A2105373数据边缘点集2.xls

     文件      23552  2008-09-20 10:20  A2105373数据边缘点集3.xls

     文件      23552  2008-09-20 10:20  A2105373数据边缘点集4.xls

     文件      22528  2008-09-20 10:20  A2105373数据边缘点集5.xls

     文件      10558  2008-09-22 02:53  A2105373源代码gongqiexian.m

     文件        151  2008-09-22 03:25  A2105373源代码jiaodian.m

     文件        301  2008-09-22 01:17  A2105373源代码jiaodian2.m

     文件       2214  2008-09-22 02:50  A2105373源代码m1.m

     文件       1287  2008-09-21 15:36  A2105373源代码moni.m

     文件       1296  2008-09-22 03:15  A2105373源代码moni2.m

     文件       1293  2008-09-22 03:16  A2105373源代码moni3.m

     文件    3152954  2008-09-19 13:41  A2105373源代码picture.bmp

     文件        867  2008-09-22 03:45  A2105373源代码代码说明.txt

     文件    1895708  2009-04-21 21:56  A透视投影中圆心投影定位和双目立体系统标定8高教社杯.pdf

     文件     972428  2009-04-21 15:58  A透视投影中圆心投影定位和双目立体系统标定01.jpg

     文件     871496  2009-04-21 16:04  A透视投影中圆心投影定位和双目立体系统标定02.jpg

     文件    1003428  2009-04-21 16:10  A透视投影中圆心投影定位和双目立体系统标定03.jpg

     文件     816805  2009-04-21 16:12  A透视投影中圆心投影定位和双目立体系统标定04.jpg

     文件     888022  2009-04-21 16:15  A透视投影中圆心投影定位和双目立体系统标定05.jpg

     文件     918886  2009-04-21 16:17  A透视投影中圆心投影定位和双目立体系统标定06.jpg

     文件     962007  2009-04-21 16:22  A透视投影中圆心投影定位和双目立体系统标定07.jpg

     文件     913950  2009-04-21 16:24  A透视投影中圆心投影定位和双目立体系统标定08.jpg

     文件    1044872  2009-04-21 16:26  A透视投影中圆心投影定位和双目立体系统标定09.jpg

     文件     876787  2009-04-21 16:28  A透视投影中圆心投影定位和双目立体系统标定10.jpg

     文件     919152  2009-04-21 16:30  A透视投影中圆心投影定位和双目立体系统标定11.jpg

     文件     414212  2009-02-24 11:27  A26_170256_.rar

     目录          0  2009-05-16 07:55  A2105373数据边缘点集

     目录          0  2009-05-16 07:55  A2105373论文

............此处省略8个文件信息