Matlab画二次及三次Bezier曲线,8控制点的B样条曲线

1、二次Bezier曲线代码：m-文件函数：functionbezier2(p0,p1,p2)t=0:0.001:1;x=(p2(1)-2*p1(1)+p0(1)*t.2+2*(p1(1)-p0(1)*t+p0(1);y=(p2(2)-2*p1(2)+p0(2)*t.2+2*(p1(2)-p0(2)*t+p0(2);plot(p0(1)p1(1)p2(1),p0(2)p1(2)p2(2),b),holdonplot(x,y,r);执行：bezier2(1,3,4,18,7,6)1234567三次Bezier曲线代码：functionbezier3(p0,p1,p2,p3)t=0:0.001:1;x

2、=(1-t).3*p0(1)+3壮.*(1-t).2*p1(1)+3*t.2.*(1-t)*p2(1)+t.3*p3(1);y=(1-t).3*p0(2)+3壮.*(1-t).2*p1(2)+3*t.2.*(1-t)*p2(2)+t.3*p3(2);plot(p0(1)p1(1)p2(1)p3(1),p0(2)p1(2)p2(2)p3(2),b);holdon;plot(x,y,r);执行：bezier3(0,3,5,20,7,2,9,1)0123456789B样条曲线代码（方法一）：functionByangtiao8(p)t=0:0.005:1;holdonfori=1:5x=p(1,i)

3、*(1/6)*(-t.3+3*t.2-3*t+1)+p(1,i+1)*(1/6)*(3壮.3-6*t.2+4).+p(1,i+2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p(1,i+3)*(1/6)*t.3;y=p(2,i)*(1/6)*(-t.3+3*t.2-3*t+1)+p(2,i+1)*(1/6)*(3*t.3-6*t.2+4).+p(2,i+2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p(2,i+3)*(1/6)*t.3;plot(x,y,k);endplot(p(1,1)p(1,2)p(1,3)p(1,4)p(1,5)p(1,6)p(1,7)p(1,8

4、),p(2,1)p(2,2)p(2,3)p(2,4)p(2,5)p(2,6)p(2,7)p(2,8);执行：Byangtiao8(4,6,3,1,7,9,15,11;0,9,11,15,15,7,6,12)B样条曲线代码(方法二)：functionByt8(p0,p1,p2,p3,p4,p5,p6,p7)t=0:0.001:1;%m=-13-31;3-630;-3030;1410;x=p0(1)*(1/6)*(-t.A3+3*t.A2-3*t+1)+p1(1)*(1/6)*(3*t.A3-6*t.A2+4).+p2(1)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p3(1)*(1/

5、6)*t.3;y=p0(2)*(1/6)*(-t.A3+3*t.A2-3*t+1)+p1(2)*(1/6)*(3*t.A3-6*t.A2+4).+p2(2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p3(2)*(1/6)*t.3;%plot(p0(1)p1(1)p2(1)p3(1),p0(2)p1(2)p2(2)p3(2);holdon;plot(x,y,r);x=p1(1)*(1/6)*(-t.A3+3*t.A2-3*t+1)+p2(1)*(1/6)*(3*t.A3-6*t.A2+4).+p3(1)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p4(1)*(1/6)*

6、t.3;y=p1(2)*(1/6)*(-t.A3+3*t.A2-3*t+1)+p2(2)*(1/6)*(3*t.A3-6*t.A2+4).+p3(2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p4(2)*(1/6)*t.3;%plot(p0(1)p1(1)p2(1)p3(1),p0(2)p1(2)p2(2)p3(2);holdon;plot(x,y,r);x=p2(1)*(1/6)*(-t.3+3*t.2-3*t+1)+p3(1)*(1/6)*(3壮.3-6*t.2+4).+p4(1)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p5(1)*(1/6)*t.3;y=p2

7、(2)*(1/6)*(-t.A3+3*t.A2-3*t+1)+p3(2)*(1/6)*(3*t.A3-6*t.A2+4).+p4(2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p5(2)*(1/6)*t.3;%plot(p0(1)p1(1)p2(1)p3(1),p0(2)p1(2)p2(2)p3(2);holdon;plot(x,y,r);x=p3(1)*(1/6)*(-t.A3+3*t.A2-3*t+1)+p4(1)*(1/6)*(3*t.A3-6*t.A2+4).+p5(1)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p6(1)*(1/6)*t.3;y=p3(2)

8、*(1/6)*(-t.A3+3*t.A2-3*t+1)+p4(2)*(1/6)*(3*t.A3-6*t.A2+4).+p5(2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p6(2)*(1/6)*t.3;%plot(p0(1)p1(1)p2(1)p3(1),p0(2)p1(2)p2(2)p3(2);holdon;plot(x,y,r);x=p4(1)*(1/6)*(-t.A3+3*t.A2-3*t+1)+p5(1)*(1/6)*(3*t.A3-6*t.A2+4).+p6(1)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p7(1)*(1/6)*t.3;y=p4(2)*(1/6)*(-t.A3+3*t.A2-3*t+1)+p5(2)*(1/6)*(3*t.A3-6*t.A2+4).+p6(2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p7(2)*(1/6)*t.3;%plot(p0(1)p1(1)p2(1)p3(1),p0(2)p1(2)p2(2)p3(2);holdon;plot(