udsheft.m

这是我们组参加2008年全国大学生数学建模时A题(数码相机定位)所用到得程序

%通过（x1,y1）(x2,y2)建立直线方程，求得有圆的下切线上的切点
function udsheft(bw,x1,y1,x2,y2)  
Y=[];upsheft=[];
for y=1:1024;   
    x=(y-y1)*(x2-x1)/(y2-y1)+x1; 
    x=floor(x);        
    if x>2;             
      if bw(x,y)==0     
          upsheft=[upsheft;x,y];  
      end
    end
end
upsheft