40;-40,40;—40,40/40,-40];

D=[40z-40;120,-120;120,-120;120z-40;8

0,80;80,80;-80,80；-80,80;80,-80];

Q=[2,-2,l,l;-3,3,-2z-l;0,0,1,0；1,0,0,0];

plot(A(:z1),A(：/2))；

[m,n]=size(A)；

holdon

fori=l:m-l

fort=0：0.001:1

A

T=[t3,t*2rt,l]；

Px=[A(i,l),A(i+lz1),B(i,l),B(i+1,1)];

py=[A(i,2),A(i+1,2),B(i,2),B(i+1,

2)];

x=T*Q*Px，;

y=T*Q*Py,；

Pxl=[A(i,l),A(i+1,1),C(i,l),C

(i+bl)l；

Py1=,A(i+1,2),C(i,2),C

<i+l,2)];

xl=T*Q*Px1*;

y1=T*Q*Py1'；

Px2=[A(irl),A(i+1,1),D(i,l),D(i+

1,1)1;

Py2=[A(i,2),A(i+1,2),D(i,2),D(i+l,2)];

X2=T*Q*PX2'；

y2=T*Q*Py2*;

plot(x,y,;

plot(xl,yl,'r1);

pl0t(x2,y2,'r1);

end

end

二、繪制Bezier曲線

1、二次Bernstein基函數性質分析

三次Bernstein基函數性質分析

end

三、繪制B樣條曲線

1、繪制均勻二次三點B樣條曲線

x=[l001502002503003504004505001;

y=[10050100501505010050100];

plot(x,y,*k\'LineWidth',1);k=1ength(x)；

B=l/2;

fori=1：k-2

fort=0:0.005：1

Bx=B*(tA2-2*t+1)*x(i)+B*(-2*t*2+2*t+

1)*x(i+1)+…

B*(tA2)*x(i+2);

By=B*(L2—2*t+1)*y(i)+B*(-2*t入2

+2*t+l)*y(i+1)+...

B*(tA2)*y(i+2);

holdon

plot(Bx,By,1r1,'LineWidth',2);

end

end

2、繪制三次B樣條曲線

x=[100120220270370420420];

y=[300200200100100200300];

plot(x,y,'k','LineWidth\l);k=1ength(x)；

B=1/6；

fori=1:k-3

fort=0:0.005:1

Bx=B*(-tA3+3*tA2-3*t+l)*x(i)

+B*(3*tA3-6*tA2-0*t+4)*x(i+l)+...

B*(-3*t7+3*L2+3*t+l)*x(i+2)+

B*(l*t*3-0*tA2-0*t+0)*x(i+3);

By=B*(-tr+3*t人2-3*t+1)*y(i)+B*(3

*tA3-6*tA2-0*t+4)*y(i+1)+...

AA

B*(-3*t34-3*t2+3*t+1)*y(i+2)+B*(1

*tA3-0*t"2-0*t+0)*y(i+3);

holdon

;

p1Ot(Bx,Byz*r*LineWidth*,2);

end

end

3、在同?控制頂點下(四個控制定點)繪制三次Bezier曲

線和三次B樣條曲線

hoidon

x=[220270300320340360370420]；

y=[200100150150150150100200];

plot(x,y?,—b1,'LineWidth',2);pI

ot(x,y,**r*z*LineWidth,,3)；

k=length(x);B=1/6;

fori=l:k-3

fort=0:0.001:1

Bx=B*(-tA3+3*tA2—3*t+1)*x(i)+B

*(3*t^3-6*tA2—0*t+4)*x(i+1)+...

B*(-3*13+3*tA2+3*t+1)*x(i+

2)+B*(l*t"3-0*t"2-0*t+0)*x(i+3);

By=B-+3*tA2-3*t+l)*y(i)+B*(3*t

-3-6*t*2—0*t+4)*y(i+l)+...

B*(-3*tA3+3*t、2+3*t+l)*y(i+2)+

B*(l*tA3-0*t"2—0*t+0)*y(i+3);

plot(Bx,By,*r'z'LineWidLh'z2);

end

end

四、有理曲線繪制各種參數圖形

繪制有理二次Bezier曲線

x=[100150200250300350400450500];

y=[l0050100501505010050100];

p1ot(x,y,'k*,*LineWidth',1)；k=1eng

th(x);B=1/2;

fori=1:k-2

fort=0:0.005：1

Bx=(B*(t*2-2*t+1)*x(i)+B*(-2*t

A2+2*t+l)*x(i+1)+...

B*(tA2)*x(i+2))/(B*(tA2-2*t+l)+B*

(-2*t-2+2*t+l)+...

B*(tA2));

By=(B*(tA2-2*t+1)*y(i)+B*(-2*t^2+2

*t+l)*y(i+1)+...

B*(tA2)*y(i+2))/(B*(t22*t+l)+

B*(-2*t人2+2*t4-1)+...

B*(tA2));

hoIdon

p1ot(Bx,By,*r*>'LineWidth',2);

end

end

均勻有理B樣條曲線各種參數圖形

繪制有理三次B樣條曲線

h0Idon

x=[220270300320340370420];

y=[200100140150130100200];

p1ot(x,y,'-b*>'LineWidth',2);

k=length(x);B=l/6；

fori=l:k-3

fort=0:0.001：1

Bx=(B*(-tA3+3*tA2-3*t+1)*x(i)+B*

(3*t"3-6*tA2-0*t+4)*x(i+l)+...

B*(-3*tA34-3*t*2+3*t+l)*x(i+2)+B

*(1*t八3-0*t人2-0*t+0)*x(i+3))/(B*+

tA2-3*t+l)+B*(3*tA3—6*tA2-0*t+4)+.

B*(-3*tA