B样条曲线正算反算贝塞尔曲线拼接曲面车身CAD作业答案知识讲解.doc

《B样条曲线正算反算贝塞尔曲线拼接曲面车身CAD作业答案知识讲解.doc》由会员分享，可在线阅读，更多相关《B样条曲线正算反算贝塞尔曲线拼接曲面车身CAD作业答案知识讲解.doc（20页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Good is good, but better carries it.精益求精，善益求善。B样条曲线正算反算贝塞尔曲线拼接曲面车身CAD作业答案-1、 贝塞尔曲线的拼接用matlab画代码如下：%Bylyqmathclc;clearall;closeall;p=12;48;615;918;p=p;t=linspace(0,1,200);n=size(p,2)-1;r=0;fork=0:nr=r+prod(1:n)/(prod(1:k)*prod(1:n-k)*p(:,k+1)*(t.k.*(1-t).(n-k);endplot(r(1,:),r(2,:),p(1,:),p(2,:),-or)

2、2、 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.3+3*t.2-3*t+1)+p1(1)*(1/6)*(3*t.3-6*t.2+4).+p2(1)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p3(1)*(1/6)*t.3;y=p0(2)*(1/6)*(-t.3+3*t.2-3*t+1)+p1(2)*(1/6)*(3*t.3-6*t.2+4).+p2(2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p3(2)*(

3、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.3+3*t.2-3*t+1)+p2(1)*(1/6)*(3*t.3-6*t.2+4).+p3(1)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p4(1)*(1/6)*t.3;y=p1(2)*(1/6)*(-t.3+3*t.2-3*t+1)+p2(2)*(1/6)*(3*t.3-6*t.2+4).+p3(2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p4(2)*(1/6)*t.3;%p

4、lot(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*t.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(2)*(1/6)*(-t.3+3*t.2-3*t+1)+p3(2)*(1/6)*(3*t.3-6*t.2+4).+p4(2)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p5(2)*(1/6)*t.3;%plot(p0(1)p1

5、(1)p2(1)p3(1),p0(2)p1(2)p2(2)p3(2);holdon;plot(x,y,r);x=p3(1)*(1/6)*(-t.3+3*t.2-3*t+1)+p4(1)*(1/6)*(3*t.3-6*t.2+4).+p5(1)*(1/6)*(-3*t.3+3*t.2+3*t+1)+p6(1)*(1/6)*t.3;y=p3(2)*(1/6)*(-t.3+3*t.2-3*t+1)+p4(2)*(1/6)*(3*t.3-6*t.2+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(

6、1),p0(2)p1(2)p2(2)p3(2);holdon;plot(x,y,r);x=p4(1)*(1/6)*(-t.3+3*t.2-3*t+1)+p5(1)*(1/6)*(3*t.3-6*t.2+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.3+3*t.2-3*t+1)+p5(2)*(1/6)*(3*t.3-6*t.2+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(

7、2)p2(2)p3(2);holdon;plot(x,y,r);plot(p0(1)p1(1)p2(1)p3(1)p4(1)p5(1)p6(1)p7(1),p0(2)p1(2)p2(2)p3(2)p4(2)p5(2)p6(2)p7(2);执行：3、 Byt8(0,0,1,4,3,9,5,7,6,2,7,6,9,5,11,3)B样条曲线的反算functionByangtiao8(p)t=0:0.005:1;holdonfori=1:5x=p(1,i)*(1/6)*(-t.3+3*t.2-3*t+1)+p(1,i+1)*(1/6)*(3*t.3-6*t.2+4).+p(1,i+2)*(1/6)*(

8、-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);end4、 plot(p(1,1)p(1,2)p(1,3)p(1,4)p(1,5)p(1,6)p(1,7)p(1,8),p(2,1)p(2,2)p(2,3)p(2,4)p(2,5)p(2,6)p(2,7)p(2,8);双三次B样条曲面的算法/TestVie

9、w.cpp:implementationoftheCTestViewclass/#includestdafx.h#includeTest.h#defineROUND(a)int(a+0.5)/四舍五入#includemath.h/数学头文件#includeTestDoc.h#includeTestView.h#ifdef_DEBUG#definenewDEBUG_NEW#undefTHIS_FILEstaticcharTHIS_FILE=_FILE_;#endif/CTestViewIMPLEMENT_DYNCREATE(CTestView,CView)BEGIN_MESSAGE_MAP(CT

10、estView,CView)/AFX_MSG_MAP(CTestView)ON_COMMAND(ID_MENUDrawHermite,OnMENUDrawHermite)/AFX_MSG_MAP/StandardprintingcommandsON_COMMAND(ID_FILE_PRINT,CView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_DIRECT,CView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_PREVIEW,CView:OnFilePrintPreview)END_MESSAGE_MAP()/CTestView

11、construction/destructionCTestView:CTestView()/TODO:addconstructioncodehere/9个型值点的初始化p100=410;p101=532;p110=210;p111=386;p120=410;p121=316;p130=310;p131=162;p140=510;p141=50;p150=710;p151=162;p160=610;p161=316;p170=810;p171=386;p180=610;p181=532;/9个型值点导数的初始化p200=100;p201=-100;p210=100;p211=-100;p220=

12、100;p221=-100;p230=100;p231=-100;p240=-100;p241=-100;p250=-100;p251=100;p260=-100;p261=100;p270=-100;p271=100;p280=-100;p281=-100;CTestView:CTestView()BOOLCTestView:PreCreateWindow(CREATESTRUCT&cs)/TODO:ModifytheWindowclassorstylesherebymodifying/theCREATESTRUCTcsreturnCView:PreCreateWindow(cs);/CTe

13、stViewdrawingvoidCTestView:OnDraw(CDC*pDC)CTestDoc*pDoc=GetDocument();ASSERT_VALID(pDoc);/TODO:adddrawcodefornativedatahereCPenMyPen,*pOldPen;MyPen.CreatePen(PS_SOLID,3,RGB(0,0,255);/蓝色笔绘制特征多边形pOldPen=pDC-SelectObject(&MyPen);pDC-MoveTo(p100,p101);pDC-Ellipse(p100-2,p101-2,p100+2,p101+2);/绘制特征多边形顶点f

14、or(inti=1;iLineTo(p1i0,p1i1);pDC-Ellipse(p1i0-2,p1i1-2,p1i0+2,p1i1+2);pDC-SelectObject(pOldPen);MyPen.DeleteObject();/CTestViewprintingBOOLCTestView:OnPreparePrinting(CPrintInfo*pInfo)/defaultpreparationreturnDoPreparePrinting(pInfo);voidCTestView:OnBeginPrinting(CDC*/*pDC*/,CPrintInfo*/*pInfo*/)/TO

15、DO:addextrainitializationbeforeprintingvoidCTestView:OnEndPrinting(CDC*/*pDC*/,CPrintInfo*/*pInfo*/)/TODO:addcleanupafterprinting/CTestViewdiagnostics#ifdef_DEBUGvoidCTestView:AssertValid()constCView:AssertValid();voidCTestView:Dump(CDumpContext&dc)constCView:Dump(dc);CTestDoc*CTestView:GetDocument(

16、)/non-debugversionisinlineASSERT(m_pDocument-IsKindOf(RUNTIME_CLASS(CTestDoc);return(CTestDoc*)m_pDocument;#endif/_DEBUG/CTestViewmessagehandlersvoidCTestView:Hermite(CDC*pDC)/绘制Hermite三次插值样条inta44=2,-2,1,1,-3,3,-2,-1,0,0,1,0,1,0,0,0;/Mh矩阵系数intb42;/边界点for(inti=0;iSelectObject(&MyPen);pDC-MoveTo(p1i0,p1i1);for(doublet=0.0;tLineTo(x,y);pDC-SelectObject(pOldPen);MyPen.DeleteObject();voidCTestView:MultiMatrix(inta44,intb42)/矩阵相乘inti,j,k;for(i=0;i4;i+)for(j=0;j2;j+)resultij=0;/矩阵清零for(i=0;i2;i+)for(j=0;j4;j+)for(k=0;kSetWindowText(案例17:三次Hermit样条曲线);Hermite(&dc);-