马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
任意给定点坐标,可以拟合出B样条曲线
- % B样条曲线生成程序
- % 说明:给定8个控制顶点{(3 5),(2 4),(3 2),(6 1),(5 8),(10 6),(8 1),(6 0)}
- load c;x=800:25:1800;
- a(1,:)=x;a(2,:)=c;
- % a=[10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90;20 16 26 26 26 26 24 24 24 22 22 16 28 26 26 26 24]; % 控制顶点数组;
- %a=[4 3 2 3 6 5 10 8 6 4;6 5 4 2 1 8 6 1 0 -1];%增加两个控制顶点(4 6)(4 -1)从而使曲线通过原来的首末点;
- plot(a(1,:),a(2,:),':'); % 绘制控制多边形;
- hold on; % 默认为hold off,此命令用来保留控制多边形的图形;
- for i=1:38; % 由于4个顶点可以确定一条B样条,共有8个控制顶点,可确定5条B样条曲线;
- %for i=1:7;
- for u=0:0.001:1; % 参数u;
-
- b0=1.0./6.*(1-u).^3; % 基函数b0;
- b1=1.0./6.*(3.*u.^3-6.*u.^2+4); % 基函数b1;
- b2=1.0./6.*(-3.*u.^3+3.*u.^2+3.*u+1); % 基函数b2;
- b3=1.0./6.*u.^3; % 基函数b3;
-
- x=b0.*a(1,i)+b1.*a(1,i+1)+b2.*a(1,i+2)+b3.*a(1,i+3); % 确定曲线的横坐标x;
- y=b0.*a(2,i)+b1.*a(2,i+1)+b2.*a(2,i+2)+b3.*a(2,i+3); % 确定曲线的纵坐标y;
-
- line(x,y); % 绘制曲线坐标点;
- end
- end
- hold off;
复制代码 |