b樣條曲線演算法

㈠ 求;B-樣條曲線的deboor演算法和C++程序

voidCGraph151View::seekPoint(intk)//表示k階2{3CPoint *interim;//存儲B樣數組指針4CDC *pDC =GetDC();5CGraph151Doc *pDoc = GetDocument();//6intj = pDoc->Length-1;7intlength=1.0/0.001;//B樣曲線的點的個數8interim =newCPoint[length];//保存B樣曲線點9intindex_length=0;//記錄interim數組下一個下標10//畫控制點圖11pDC->MoveTo(pDoc->m_point[0]);12for(intn=0;n<pDoc->Length;n++)13{14pDC->LineTo(pDoc->m_point[n]);15pDC->MoveTo(pDoc->m_point[n]);16}17///結束1819for(floatt=pDoc->T[0];t<pDoc->T[pDoc->T_Length-1];t +=0.001)20{21inti;22chushihua();//還原m_point 和 im_point點數組使與初始化相同23for(intr=1;r<k;r++)24{25guodian();//更新一下m_point點數組 的數據26for( i=r;i<=j;i++)27{28floatpara1;29floatpara2;3031if( (pDoc->T[i+k-r] - pDoc->T[i]) !=0.0)32{33para1 = (t - pDoc->T[i]) / (pDoc->T[i+k-r] - pDoc->T[i]);34para2 = (pDoc->T[i+k-r] - t) /(pDoc->T[i+k-r] - pDoc->T[i]);35}36else37{38para1=0.0;39para2=0.0;40}4142intx =int( (para1 * pDoc->m_point[i].x) + (para2 * pDoc->m_point[i-1].x) );43inty =int( (para1 * pDoc->m_point[i].y) + (para2 * pDoc->m_point[i-1].y) );44pDoc->im_point[i].SetPoint(x,y);45}4647}4849interim[index_length++].SetPoint( pDoc->im_point[3].x,pDoc->im_point[3].y);//存儲計算出的點50}5152//畫圖 B樣曲線53CPen pen(PS_SOLID,1,RGB(255,0,0));54pDC->SelectObject(pen);55pDC->MoveTo(interim[0]);56for(intn=0;n<length;n++)57{58pDC->LineTo(interim[n]);59pDC->MoveTo(interim[n]);60}61//B樣曲線畫圖結束62}

㈡ 三次均勻B樣條插值演算法

三次樣條插值（Cubic Spline Interpolation）簡稱Spline插值，是通過一系列形值點的一條光滑曲線，數學上通過求解三彎矩方程組得出曲線函數組的過程。
實際計算時還需要引入邊界條件才能完成計算。一般的計算方法書上都沒有說明非扭結邊界的定義，但數值計算軟體如Matlab都把非扭結邊界條件作為默認的邊界條件。