等高线建模与函数拟合建模.docVIP

班级：油工10902 姓名：潘建生 序号：69 题目：1. 用给定的多项式，如y = x3 ? 6x2 + 5x ? 3 ，产生一组数据( , ) i i x y ， i = 1,2,L,m ，再在i y 上添加随机干扰（可用rand 产生（0,1）均匀分布随机数，或 用randn 产生N(0,1)分布随机数），然后用i x 和添加了随机干扰的i y 作3 次多项式拟 合，与原系数比较。如果作2 或4 次多项式拟合，结果如何？ 分析： 有题目可知，首先我们得必须产生一组若干（x ,y）的书然后再产生干扰数y,根据所得数组用多项式拟合函数polyfit可求的多项式系数，而由系数可用polyval函数求得一系列的点，从而话出图像进行比较。 解析代码如下： x0=1:2:40; %x0=linspace(1,60,30); n=length(x0); for i=1:n y0(i)=x0(i)^3-6*x0(i)^2+5*x0(i)-3+rand(1); end a0=polyfit(x0,y0,2) a1=polyfit(x0,y0,3) a2=polyfit(x0,y0,4) for i=1:n y1(i)=polyval(a1,x0(i)); y2(i)=polyval(a2,x0(i)); y3(i)=x0(i)^3-6*x0(i)^2+5*x0(i)-3; end subplot(2,2,1),plot(x0,y0,*) ,title(二次拟合图像) subplot(2,2,2),plot(x0,y1,+) ,title(三次拟合图像) subplot(2,2,3),plot(x0,y2,o) ,title(四次拟合图像) subplot(2,2,4),plot(x0,y3,o) ,title(原函数图像) 求解结果截图如下： 1.系数结果图： 2.图像比较图： 结论：有求出的系数与所画的图像比较可知，拟合的次数过高与过低多不是太理想，太低所得的函数相差较远，不能达到解题的目的，太高有相当于做无用功，因为前面的系数已经无限接近零，或者说等于零了，我们拟合再高的次数也没多大帮助，所以对于拟合函数时最好先估算一下，然后找出比较适合的一个区间，得到一个最优解。 表9： 4800 4400 4000 3600 3200 2800 2400 2000 1600 1200 800 400 0 1350 1370 1390 1400 1410 960 940 880 800 690 570 430 290 210 150 1370 1390 1410 1430 1440 1140 1110 1050 950 820 690 540 380 300 210 1380 1410 1430 1450 1470 1320 1280 1200 1080 940 780 620 460 370 350 1420 1430 1450 1480 1500 1550 1510 1430 1300 1200 980 850 750 550 500 1430 1450 1460 1500 1550 1600 1550 1600 1600 1600 1550 1500 1500 1550 1550 950 1190 1370 1500 1200 1100 1550 1600 1550 1380 1070 900 1050 1150 1200 910 1090 1270 1500 1200 1100 1350 1450 1200 1150 1010 880 1000 1050 1100 880 1060 1230 1390 1500 1500 1400 900 1100 1060 950 870 900 936 950 830 980 1180 1320 1450 1420 400 1300 700 900 850 810 380 780 750 740 880 1080 1130 1250 1280 1230 1040 900 500 700 780 750 650 550 650 760 880 970 1020 1050 1020 830 800 700 300 500 550 480 350 510 620 730 800 850 870 850 780 720 650 500 200 300 350 320 370 470 550 600 670 690 670 620 580 450 400 300 100 150 250 X/Y 0 400 800 1200 1600 20