計算機輔助設計研究專題報告-數(shù)據(jù)的光順與擬合

1、整體 clear alla=load(C:UsersthinkpadDesktopclf.txt);b=a;x=b(:,1);y=b(:,2);z=b(:,3);plot3(x,y,z,+); X & Yclear alla=load(C:UsersthinkpadDesktopclf.txt);b=a;x=b(:,1);y=b(:,2);z=b(:,3);plot(x,y,+); X&Y 去噪聲點并修正clear alla=load(C:UsersthinkpadDesktopclf.txt);b=a;x=b(:,1);y=b(;,2);plot(x,y,o);hold onp=polyfi

2、t(x,y,3);for i=1:length(x)m(i)=polyval(p,x(i)-y(i);sum=0; sum=sum+m(i)*m(i);endk=sqrt(sum/length(x);t=0;for i=1:length(x) if abs(m(i)=k t=t+1; y(i)=polyval(p,x(i); endendplot(x,y,*r);hold on插值xi=4680.00:1:5040.00;yi=spline(x,y,xi);plot(xi,yi,.b) 1:15 f(x) = p1*x2 + p2*x + p3 p1 = 0.007571 (0.0075, 0

3、.007643) p2 = -70.86 (-71.53, -70.19) p3 = 1.661e+005 (1.645e+005, 1.676e+005)16:25 f(x) = p1*x + p2 p1 = 0.5478 (0.5449, 0.5507) p2 = -2289 (-2302, -2275)26: 44 f(x) = p1*x3 + p2*x2 + p3*x + p4 p1 = 0.0002025 (0.0001959, 0.0002091) p2 = -2.893 (-2.987, -2.8) p3 = 1.378e+004 (1.333e+004, 1.422e+004)

4、 p4 = -2.187e+007 (-2.257e+007, -2.117e+007)45:65f(x) = p1*x + p2 p1 = -0.2425 (-0.2445, -0.2404) p2 = 1457 (1447, 1466)66:78 f(x) = p1*x + p2 p1 = -0.333 (-0.3365, -0.3296) p2 = 1892 (1875, 1908)79:100f(x) = p1*x3 + p2*x2 + p3*x + p4 p1 = -4.327e-005 (-4.609e-005, -4.046e-005) p2 = 0.6354 (0.5944,

5、0.6765) p3 = -3110 (-3309, -2911) p4 = 5.074e+006 (4.751e+006, 5.396e+006)101:118 f(x) = p1*x + p2 p1 = 0.2357 (0.2332, 0.2383) p2 = -867.2 (-879.8, -854.7)119:135 f(x) = p1*x3 + p2*x2 + p3*x + p4 p1 = 0.0001109 (0.000106, 0.0001159) p2 = -1.642 (-1.715, -1.568) p3 = 8100 (7738, 8462) p4 = -1.332e+0

6、07 (-1.392e+007, -1.273e+007)136:164 f(x) = p1*x3 + p2*x2 + p3*x + p4 p1 = -5.707e-005 (-5.857e-005, -5.558e-005) p2 = 0.847 (0.8247, 0.8693) p3 = -4190 (-4301, -4079) p4 = 6.909e+006 (6.724e+006, 7.093e+006)165:181f(x) = p1*x3 + p2*x2 + p3*x + p4 p1 = -8.58e-006 (-1.18e-005, -5.363e-006) p2 = 0.134

7、8 (0.08631, 0.1833) p3 = -704.7 (-948.3, -461.2) p4 = 1.227e+006 (8.191e+005, 1.635e+006)X&z 去噪聲點并修正clear alla=load(C:UsersthinkpadDesktopclf.txt);b=a;x=b(:,1);z=b(;,3);plot(x,z,o);hold onp=polyfit(x,z,3);for i=1:length(x)m(i)=polyval(p,x(i)-z(i);endsum=0;for i=1:length(x) sum=sum+m(i)*m(i);endk=sqr

8、t(sum/length(x);t=0;for i=1:length(x) if abs(m(i)=k t=t+1; z(i)=polyval(p,x(i); endendplot(x,z,*r);hold on插值xi=4680.00:1:5040.00;zi=spline(x,z,xi);plot(xi,zi,.b)1:30f(x) = p1*x3 + p2*x2 + p3*x + p4 p1 = 2.262e-006 (-9.388e-007, 5.464e-006) p2 = -0.0321 (-0.07732, 0.01313) p3 = 151.8 (-61.15, 364.8)

9、p4 = -2.394e+005 (-5.736e+005, 9.492e+004)31:60f(x) = p1*x3 + p2*x2 + p3*x + p4 p1 = 9.025e-007 (5.064e-007, 1.299e-006) p2 = -0.0131 (-0.01876, -0.007428) p3 = 63.33 (36.3, 90.36) p4 = -1.021e+005 (-1.45e+005, -5.91e+004)61:90 f(x) = p1*x2 + p2*x + p3 p1 = -0.0001227 (-0.0001265, -0.0001189) p2 = 1.166 (1.129, 1.202) p3 = -2767 (-2856, -2678)91:120f(x) = p1*x + p2 p1 = -0.01897 (-0.01911, -0.01883) p2 = 92.14 (91.45, 92.83)121:150 f(x) = p1*x2 + p2*x + p3 p1 = 0.0002289 (0.0001707, 0.0002871) p2 = -2.265 (-2.84