t=5:1:99;
x=-0.207806-0.610303*t+0.206455*t.^2-0.0144935*t.^3+...
    0.000517774*t.^4-8.394241977754047*(10^(-6))*t.^5+...
    6.133353112035975*(10^(-8))*t.^6-1.6673218267444805*(10^(-10))*t.^7;
y=158.211+1.86595*t-0.266798*t.^2+0.0141407*t.^3+...
    -0.000325412*t.^4+3.043275597680807*(10^(-6))*t.^5+...
    -9.899171274615063*(10^(-9))*t.^6;
z=t;
R=30;
plot3(x,y,z);
hold on 

%開始沿軸線方向繪制其法向的圓。
[M N]=size(t);
theta=0:0.6:2*pi;
[MM NN]=size(theta);
pp=repmat(0,[NN 3 N]);
%橫向網(wǎng)格的繪制
for i=1:N
    tt=t(i);
    
    T=rr(tt);
    T=T/(T(1)^2+T(2)^2+T(3)^2)^0.5;
    
    B=cross(rr(tt),rrr(tt));
    B=B/(B(1)^2+B(2)^2+B(3)^2)^0.5;
    
    M=cross(B,T);
    
    for j=1:NN
        pp(j,:,i)=0*R*T+R*cos(theta(j))*B+R*sin(theta(j))*M+[x(i) y(i) z(i)];
    end
    for j=1:NN-1
        plot3([pp(j,1,i) pp(j+1,1,i)],[pp(j,2,i) pp(j+1,2,i)],[pp(j,3,i) pp(j+1,3,i)],'r');
    end
    plot3([pp(1,1,i) pp(NN,1,i)],[pp(1,2,i) pp(NN,2,i)],[pp(1,3,i) pp(NN,3,i)],'r');
end
%縱向網(wǎng)格的繪制
for i=1:NN
    for j=1:N-1
       plot3([pp(i,1,j) pp(i,1,j+1)],[pp(i,2,j) pp(i,2,j+1)],[pp(i,3,j) pp(i,3,j+1)],'r'); 
   end
end

axis equal
grid on

figure
%面狀繪圖
hold on
for i=1:NN-1
    for j=1:N-1
        
        q1=[pp(i,1,j) pp(i,2,j) pp(i,3,j)];
        q2=[pp(i+1,1,j) pp(i+1,2,j) pp(i+1,3,j)];
        q3=[pp(i+1,1,j+1) pp(i+1,2,j+1) pp(i+1,3,j+1)];
        q4=[pp(i,1,j+1) pp(i,2,j+1) pp(i,3,j+1)];    
        patch([q1(1) q2(1) q3(1) q4(1)],[q1(2) q2(2) q3(2) q4(2)],[q1(3) q2(3) q3(3) q4(3)],'r'); 
    end
end
for j=1:N-1
    q1=[pp(NN,1,j) pp(NN,2,j) pp(NN,3,j)];
    q2=[pp(1,1,j) pp(1,2,j) pp(1,3,j)];
    q3=[pp(1,1,j+1) pp(1,2,j+1) pp(1,3,j+1)];
    q4=[pp(NN,1,j+1) pp(NN,2,j+1) pp(NN,3,j+1)];
    patch([q1(1) q2(1) q3(1) q4(1)],[q1(2) q2(2) q3(2) q4(2)],[q1(3) q2(3) q3(3) q4(3)],'r');
end

axis equal
grid on
axis([-35 210 -250 200 0 140])