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

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  A\2105373\论文\数码相机定位.pdf

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

     文件     414212  2009-02-24 11:27  A\26_170256_.rar

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

     目录          0  2009-05-16 07:55  A\2105373\论文

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