版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、用matlab实现三次 NURBS插值曲线作者:大漠孤狼 发表于matlab乐园()作者:这是我在大学时做大学生研究计划时写的,当时刚学会matlab ,编写了这个程序,用了很多循环,效率不高.当时我并不清楚循环是matlab的弱点,等明白了,也不做这方面的工作了,也就懒的去改写了.如果谁需要用,就自己改吧.算法也有一些问题,我就不多说了,自己看吧 一 、三次NURBS插值算法 给定平面控制顶点di(i=1,2,n)及对应的权因子(i=1,2,n),可定义一条三次NURBS曲线。先对控制顶点进行参数化,得一矢量: U=u0,u1,u2,un+4 则三次NURBS曲线的分段方程形
2、式为: uuk+3,uk+4 k=0,1,2,n-3 (1) 首先证明一条性质: 若三点dk,dk+1,dk+2共线,且满足 则三次NURBS曲线插值点dk+1. 证明:由(1)式可得: 由以知可得 故性质得证。 下面构造三次 NURBS插值曲线; 取一型值点列V0,V1,V2,Vn,用下列方法决定各点Vi处的切矢Ti: 设顺序五个数据点Vi-2,Vi-1,Vi+2,在中间数据点Vi的切矢ti被局部定义为 ti=(1-a)pi-1+api 其中,pi=pi+1-pi a=, A= B= 在非周期情况下,首末各两数据点处的切矢可采用贝塞尔条件确定,即用如下差分矢量 p-1=2p0p1 p-2=2
3、p-1p0 pn=2pn1pn2 pn+1=2pnpn1 仍按前述公式确定t0,t1,tn-1,tn。 下面,我们令d3i=pi (i=1n) 然后对数据点d3i进行参数化 u0=u1=u2=u3=0 (i=2,3,,n) u3n+1=u3n+2=u3n+3=u3n+4=1 在区间u3i,u3(i+1)插入节点u3i+1 ,u3i+2 u3i+1= 即得节点矢量。 下面在点d3i ,d3(i+1) (i=0 , 1 , ,n-1) 之间插入两个辅助控制顶点 d3i+1 , d3i+2 ,构成新控制顶点组,得分段三次NURBS方程。 我们按下列原则判断两型值点 pi (d3i) , pi+1(
4、d3(i+1) ) 之间的曲线段是否有拐点出现:如果矢量乘积ai-1ai 与aiai+1同向,则pi 与pi+1之间无拐点,此时两切线Li , Li+1之间有一交点 Vi*=Vi+ 否则,Vi 与Vi+1 之间有一拐点,记 Vi*=(Vi+Vi+1)2 令 T0 = Ti = min , (i = 1 , 2 , 3 , n-1) Tn = 我们按下列方法决定 Vi 与 Vi+1 之间的辅助点 d3i+1 , d3i+2 讨论在区间 u3i+2 , u3i+3上的曲线段 uu3i+2 , u3i+3 令 d3i = Vi , i = 0 , 1, 2 , n d3i+1 = Vi + (0 i
5、 Ti ) d3i+2 = Vi - (0 i+1 Ti+1 ) 由此三式得 由性质可得 P( u3i+2 ) = d3i , 即此曲线插值点 d3i ( Vi ) .二、基函数计算及整体表示1, 基函数根据伯恩斯坦多项式可计算Ni , 3( u )=+ + + 整理后得:Ni+1,0 的系数(-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-u(1+i)-1/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-u(1+i)-1/(u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-u(1+i)*u3+(1/(u(4+i)-u(1+
6、i)/(u(3+i)-u(i)/(u(2+i)-u(1+i)*u(4+i)+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-u(1+i)*u(i)+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-u(1+i)*u(1+i)+2/(u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-u(1+i)*u(i)+1/(u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-u(1+i)*u(2+i)+2/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-u(1+i)*u(1+i)+1/(u(3+i)-
7、u(i)/(u(3+i)-u(1+i)/(u(2+i)-u(1+i)*u(3+i)*u2+(-1/(u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-u(1+i)*u(i)2-2/(u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-u(1+i)*u(i)*u(2+i)-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-u(1+i)*u(i)*u(1+i)-1/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-u(1+i)*u(1+i)2-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-u
8、(1+i)*u(i)*u(3+i)-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-u(1+i)*u(3+i)*u(1+i)-2/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-u(1+i)*u(4+i)*u(1+i)*u+1/(u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-u(1+i)*u(i)2*u(2+i)+1/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-u(1+i)*u(4+i)*u(1+i)2+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-u(1+i)*u
9、(i)*u(3+i)*u(1+i)Ni+2,0(u)的系数(1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)+1/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(3+i)-u(2+i)+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u3+(-1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(4+i)-1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(3+i)-1/(u(4+i)-u(1+i)/(
10、u(4+i)-u(2+i)/(u(3+i)-u(2+i)*u(2+i)-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(i)-2/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(3+i)-2/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(3+i)-u(2+i)*u(4+i)-1/( u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(1+i)*u2+(1/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(3+i)-u(2+i)*u(
11、4+i)2+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(3+i)2+1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(4+i)*u(1+i)+1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(4+i)*u(3+i)+2/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(i)*u(3+i)+2/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(3+i)-u(2+i)*u(4+i)*u(2+i)+1
12、/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(1+i)*u(3+i)*u-1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(4+i)*u(1+i)*u(3+i)-1/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(3+i)-u(2+i)*u(4+i)2*u(2+i)-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(i)*u(3+i)2Ni+3,0(u) 的系数1/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(4+i
13、)-u(3+i)*u(4+i)3-3/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(4+i)-u(3+i)*u(4+i)2*u+3/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(4+i)-u(3+i)*u(4+i)*u2-1/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(4+i)-u(3+i)*u31, 整体表示为了便于编程,采用矩阵表示NURBS曲线。三次NURBS的分段表示为 uuk+3,uk+4 k=0,1,2,n-3 设D=idi T i= 0 , 1 , 2 , , n = 1 , 2 , , n i= 0 , 1 , 2
14、, , n Qk =N0 , 3(u) , N1 , 3(u) , Nn , 3(u) i= 0 , 1 , 2 , , n U =1 u u2 u3 Qk 表示uuk+3 ,uk+4 内的基函数矢量,则NURBS的矩阵表示为 由此可见,NURBS曲线可由矩阵 D Q 表示。进一步,可由分析矩阵而得到曲线的性质。三、插值源程序%给定点列V( i ),随机给定权因子,构造三次NURBS插值点列V( i )%载入数据 注意:每一行必为矢量(三个分量),否则会出错%参数化V=0.5 2 0;1 3 0;2 8 0;3 10 0;5 11 0;8 30 0;10 10 0;s=size(V);n=s(
15、1,1);u(1)=0;u(2)=0;u(3)=0;u(4)=0;u(3*n-1)=1;u(3*n)=1;u(3*n+1)=1;u(3*n+2)=1;sum=0;for i=1:(n-1) %计算边矢及求和 a(i+2,1:3)=V(i+1,1:3)-V(i,1:3); A=sqrt(dot(a(i+2,1:3),a(i+2,1:3); sum=sum+A; %求和endfor i=3:n sumb=0; for j=2:(i-1) %? b=V(j,1:3)-V(j-1,1:3); B=sqrt(dot(b,b); sumb=sumb+B; end u(3*i-2)=sumb/sum;end
16、for i=1:(n-1) u(3*i-1)=2/3*u(3*i-2)+1/3*u(3*i+1); u(3*i)=1/3*u(3*i-2)+2/3*u(3*i+1);end%计算切矢a(2,1:3)=2*a(3,1:3)-a(4,1:3);a(1,1:3)=2*a(2,1:3)-a(3,1:3);a(n+2,1:3)=2*a(n+1,1:3)-a(n,1:3);a(n+3,1:3)=2*a(n+2,1:3)-a(n+1,1:3);for i=1:n+2 b(i,1:3)=cross(a(i,1:3),a(i+1,1:3); A(i)=sqrt(dot(b(i,1:3),b(i,1:3);end
17、for i=1:n M=A(i)/(A(i)+A(i+2); t(i,1:3)=(1-M)*a(i+1,1:3)+M*a(i+2,1:3);end%随机产生权因子w=rand(3*n-2,1); %w(i)在之间随机产生 %建立三维数组Qfor k=1:(3*n-5) for i=1:3*n-2 if i=k Q(1:4,i,k)=1/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(4+i)-u(3+i)*u(4+i)3;. -3/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(4+i)-u(3+i)*u(4+i)2;3/(u(4+i)-u(1+i)/. (
18、u(4+i)-u(2+i)/(u(4+i)-u(3+i)*u(4+i);-1/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/. (u(4+i)-u(3+i); elseif i=(k+1) Q(1:4,i,k)=-1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-. u(2+i)*u(4+i)*u(1+i)*u(3+i)-1/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/. (u(3+i)-u(2+i)*u(4+i)2*u(2+i)-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/. (u(3+i)-u(2+i)*u(i)*
19、u(3+i)2;(1/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(3+i)-. u(2+i)*u(4+i)2+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(3+i)2+1/. (u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(4+i)*u(1+i)+1/(u(4+i)-. u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(4+i)*u(3+i)+2/(u(3+i)-u(i)/. (u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(i)*u(3
20、+i)+2/(u(4+i)-u(1+i)/(u(4+i)-. u(2+i)/(u(3+i)-u(2+i)*u(4+i)*u(2+i)+1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/. (u(3+i)-u(2+i)*u(1+i)*u(3+i);(-1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-. u(2+i)*u(4+i)-1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(3+i)-1/. (u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(3+i)-u(2+i)*u(2+i)-1/
21、(u(3+i)-u(i)/. (u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(i)-2/(u(3+i)-u(i)/(u(3+i)-u(1+i)/. (u(3+i)-u(2+i)*u(3+i)-2/(u(4+i)-u(1+i)/(u(4+i)-u(2+i)/(u(3+i)-. u(2+i)*u(4+i)-1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)*u(1+i);.(1/(u(4+i)-u(1+i)/(u(3+i)-u(1+i)/(u(3+i)-u(2+i)+1/(u(4+i)-u(1+i)/.(u(4+i)-u(2+i)/(u(3
22、+i)-u(2+i)+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(3+i)-.u(2+i); elseif i=(k+2) Q(1:4,i,k)=1/(u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-. u(1+i)*u(i)2*u(2+i)+1/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-. u(1+i)*u(4+i)*u(1+i)2+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-. u(1+i)*u(i)*u(3+i)*u(1+i);(-1/(u(3+i)-u(i)/(u(2+i)-u(i)/
23、. (u(2+i)-u(1+i)*u(i)2-2/(u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-. u(1+i)*u(i)*u(2+i)-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-. u(1+i)*u(i)*u(1+i)-1/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-. u(1+i)*u(1+i)2-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-. u(1+i)*u(i)*u(3+i)-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-. u(1+i)*
24、u(3+i)*u(1+i)-2/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-. u(1+i)*u(4+i)*u(1+i);(1/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-. u(1+i)*u(4+i)+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-. u(1+i)*u(i)+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-. u(1+i)*u(1+i)+2/(u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-. u(1+i)*u(i)+1/(u(3+i)-u(i)
25、/(u(2+i)-u(i)/(u(2+i)-. u(1+i)*u(2+i)+2/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-. u(1+i)*u(1+i)+1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-. u(1+i)*u(3+i);(-1/(u(3+i)-u(i)/(u(3+i)-u(1+i)/(u(2+i)-. u(1+i)-1/(u(4+i)-u(1+i)/(u(3+i)-u(i)/(u(2+i)-u(1+i)-1/. (u(3+i)-u(i)/(u(2+i)-u(i)/(u(2+i)-u(1+i); elseif i=(k+3) Q(1:4,i,k)=-u(i)3/(u(i+3)-u(i)/(u(i+2)-u(i)/. (u(i+1)-u(i);3*u(i)2/(u(i+3)-u(i)/(u(i+2)-u(i)/(u(i+1)-u(i);. -3*u(i)/(u(i+3)-u(i)/(u(i+2)-u(i)/(u(i+1)-u(i);1/(u(i+3)-u(i)/. (u(i+2)-u(i)/(u(i+1)-u(i); else Q(1:4,i,k)=0;
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 4.7 图形的位似 浙教版数学九年级上册课件
- 体育《弯道跑》说课稿
- 5年中考3年模拟试卷初中道德与法治七年级下册02第2课时我与集体共成长
- 1-抖音商业模式分析
- 小学音乐三年级下册教案
- 多需求类存货配给策略综述
- 多彩多姿的现代雕塑 课件 2024-2025学年赣美版初中美术八年级上册
- (统考版)2023版高考化学一轮复习第十章化学实验第1讲常见仪器的使用和实验基本操作学生用书
- 汽车涂装技术(彩色版配实训工单)课件 任务三 施工和修补底漆
- 书店装修协议
- 车位租赁合同电子版
- 化妆品行业标准操作程序《玻璃瓶检验标准》
- 《煤矿窄轨铁道质量标准及检查评级办法》(83)煤生字第892号
- 关于幼儿数学教育与区域游戏相融合的思考
- 金相显微镜点检卡
- 短信营销培训文档
- 天之至私,用之至公,禽之制在气熊厚音
- 《人类成长与社会环境》形考作业1-4答案
- 担保公司担保业务审查委员会工作规程
- 活动执行时间推进表
- 关于学习考察气象灾害防御立法工作的调研报告
评论
0/150
提交评论