www.pudn.com > ceshitwoccd.zip > Untitled5.m, change:2015-05-08,size:669b


% clear 
% clc 
 
load('pppr.mat'); 
load('qqqr.mat'); 
 
 
ycr=GetYY(214,226,ppp,qqq); 
zcr=GetZZ(214,226,ppp,qqq); 
 
ycr1=GetYY(400,326,ppp,qqq); 
zcr1=GetZZ(400,326,ppp,qqq); 
 
 
theta=-atan((zcr1-zcr)/(ycr1-ycr)); 
 
i=92; 
    a=num2str(i); 
    str=strcat(a,'.bmp'); 
    p=ReadPICData(str); 
    pp=p; 
% end 
% load('pp.mat'); 
% for i=1:1 
    for j=1:800 
        datay(j)=GetYY(j,pp(j),ppp,qqq); 
        dataz(j)=GetZZ(j,pp(j),ppp,qqq); 
    end 
% end 
 
% for i=1:1 
    for j=1:800 
        ty(j)=(datay(j)-ycr)*cos(theta)+(dataz(j)-zcr)*sin(theta); 
        tz(j)=(dataz(j)-zcr)*cos(theta)-(datay(j)-ycr)*sin(theta); 
    end 
% end 
 
plot(ty,tz)