[ÇóÖú] matlab±Æ½ü·½³ÌʽÌáÈ¡

clear;clc;format long


load('data03v_0');
load('data03v_30');
load('data03v_60');

x0_1=data_03v_0(:,1);
y0_1=data_03v_0(:,2);
z0_1=data_03v_0(:,3);
a0_1=[x0_1 y0_1 z0_1];

x0_2=data_03v_0(:,4);
y0_2=data_03v_0(:,5);
z0_2=data_03v_0(:,6);
a0_2=[x0_2 y0_2 z0_2];

x0_3=data_03v_0(:,7);
y0_3=data_03v_0(:,8);
z0_3=data_03v_0(:,9);
a0_3=[x0_3 y0_3 z0_3];

x30_1=data_03v_30(:,1);
y30_1=data_03v_30(:,2);
z30_1=data_03v_30(:,3);
a30_1=[x30_1 y30_1 z30_1];

x30_2=data_03v_30(:,4);
y30_2=data_03v_30(:,5);
z30_2=data_03v_30(:,6);
a30_2=[x30_2 y30_2 z30_2];

x30_3=data_03v_30(:,7);
y30_3=data_03v_30(:,8);
z30_3=data_03v_30(:,9);
a30_3=[x30_3 y30_3 z30_3];

x60_1=data_03v_60(:,1);
y60_1=data_03v_60(:,2);
z60_1=data_03v_60(:,3);
a60_1=[x60_1 y60_1 z60_1];


x60_2=data_03v_60(:,4);
y60_2=data_03v_60(:,5);
z60_2=data_03v_60(:,6);
a60_2=[x60_2 y60_2 z60_2];

x60_3=data_03v_60(:,7);
y60_3=data_03v_60(:,8);
z60_3=data_03v_60(:,9);
a60_3=[x60_3 y60_3 z60_3];

X0=[a0_1(:,1) a0_2(:,1) a0_3(:,1)];
Y0=[a0_1(:,2) a0_2(:,2) a0_3(:,2)];
Z0=[a0_1(:,3) a0_2(:,3) a0_3(:,3)];

X30=[a30_1(:,1) a30_2(:,1) a30_3(:,1)];
Y30=[a30_1(:,2) a30_2(:,2) a30_3(:,2)];
Z30=[a30_1(:,3) a30_2(:,3) a30_3(:,3)];

X60=[a60_1(:,1) a60_2(:,1) a60_3(:,1)];
Y60=[a60_1(:,2) a60_2(:,2) a60_3(:,2)];
Z60=[a60_1(:,3) a60_2(:,3) a60_3(:,3)];

n=0.1;
xA=[];
yA=[];
zA=[];
centerline = [9,15,6;
              8,8.8,6];
for w=[26:61 61:72 1:26]
aa=[X0(w,:) X30(w,:) X60(w,:) centerline(1,1) centerline(2,1);
    Y0(w,:) Y30(w,:) Y60(w,:) centerline(1,2) centerline(2,2);
    Z0(w,:) Z30(w,:) Z60(w,:) centerline(1,3) centerline(2,3)];

% x=[aa(1,1:9)];
% y=[aa(2,1:9)];
% z=[aa(3,1:9)];

x=[aa(1,1:11)];
y=[aa(2,1:11)];
z=[aa(3,1:11)];

[X,Y]=meshgrid(min(x):n:max(x),min(y):n:max(y));
Z=griddata(x,y,z,X,Y,'cubic');
% mesh(X,Y,Z)
surf(X,Y,Z)

hold on
view(-60, 10)
axis([0 10 10 23 1 15])
xA=[xA x];
yA=[yA y];
zA=[zA z];
plot3(xA,yA,zA,'.')

plot3([centerline(1,1) centerline(2,1)],[centerline(1,2) centerline(2,2)],[centerline(1,3) centerline(2,3)])

m(w)=getframe
pause
hold off
end