數(shù)學實驗擬合與插值

數(shù)學實驗擬合與插值第一頁，共二十三頁，2022年，8月28日北京科技大學數(shù)學實驗曲線擬合插值第二頁，共二十三頁，2022年，8月28日已知平面上若干個點（xi,yi)i=1,…n,

求一個較簡單的函數(shù)（曲線）y=f(x),

使f(x)盡可能的靠近數(shù)據(jù)點，在某種意義下達到最優(yōu)。f(x)稱為擬合函數(shù)。+++++++++xyy=f(x)(xi,yi)i------最小二乘法曲線擬合第三頁，共二十三頁，2022年，8月28日

p=polyfit（x，y，n）多項式數(shù)據(jù)擬合功能：將給定向量x，y對應的（x[i]，y[i]）作為數(shù)據(jù)點，擬合成n次多項式；向量x，y具有相同的維數(shù)；p為多項式的系數(shù)向量。多項式的擬合y=a1xn+a2xn-1+……...+anx+an+1p=[a1,a2,...,an,an+1]第四頁，共二十三頁，2022年，8月28日poly2str(p,‘x’)將多項式表示成習慣的形式p是多項式系數(shù)，字符‘x’為自變量polyval（p，X）

按數(shù)組規(guī)則計算X處多項式的值功能：計算多項式p的變量在點陣x處的值；X可以為向量或矩陣，計算結果是與X同維的向量或矩陣。第五頁，共二十三頁，2022年，8月28日x=1:0.1:2;y=[2.1,3.2,2.1,2.5,3.2,3.5,3.4,4.1,4.7,5.0,4.8];p2=polyfit(x,y,2)%多項式擬合,次數(shù)是2,p2為擬合多項式的系數(shù)p3=polyfit(x,y,3);p7=polyfit(x,y,7);disp(‘二次擬合函數(shù)'),f2=poly2str(p2,'x')disp(‘三次擬合函數(shù)'),f3=poly2str(p3,'x');disp(‘七次擬合函數(shù)'),f7=poly2str(p7,'x');p2=1.3869-1.26082.141二次擬合函數(shù)f2=1.3869x^2-1.2608x+2.141第六頁，共二十三頁，2022年，8月28日7x1=1:.01:2;y2=polyval(p2,x1);%多項式p2在x1處的值y3=polyval(p3,x1);y7=polyval(p7,x1);plot(x,y,'rp',x1,y2,'--',x1,y3,'k-.',x1,y7);第七頁，共二十三頁，2022年，8月28日8練習：用多項式擬合余弦函數(shù)y=cosx。第八頁，共二十三頁，2022年，8月28日插值對函數(shù)f(x),其函數(shù)形式可能很復雜，假如可以獲得f(x)在區(qū)間[a,b]上的一組n+1個不同的點

上的函數(shù)值求一個簡單函數(shù)p(x)，使得：并且用p(x)近似代替f(x),這就是插值問題。函數(shù)p(x)為函數(shù)f(x)的插值函數(shù)。(1)式稱為插值條件。xi為插值結點，點x稱為插值點。點x在插值區(qū)間內叫內插，否則叫外插?！?.（1）第九頁，共二十三頁，2022年，8月28日p(x)f(x)x0x1x2x0y1y0y2第十頁，共二十三頁，2022年，8月28日對于被插函數(shù)f(x)和插值函數(shù)p(x)在節(jié)點xi處得函數(shù)值必然相等但在節(jié)點外p(x)的值可能就會偏離f(x)因此p(x)近似代替f(x)必然存在著誤差整體誤差的大小反映了插值函數(shù)的好壞為了使插值函數(shù)方便在計算機上運算，一般插值函數(shù)都使用多項式和有理函數(shù)。第十一頁，共二十三頁，2022年，8月28日一維多項式插值

yi=interp1（x，y，xi，method）功能：x,y是已知數(shù)據(jù)點；yi是插值點xi處的值；當輸入的x是等間距時，可在插值方法method前加*，以提高處理速度。method表示不同的插值方法的字符串，有下面四種插值方法可選:第十二頁，共二十三頁，2022年，8月28日‘nearest’：最近點插值，插值點處的值取與該插值點距離

最近的數(shù)據(jù)點函數(shù)值；‘linear’：分段線性插值，用直線連接數(shù)據(jù)點，插值點的值

取對應直線上的值；‘spline’：三次樣條函數(shù)插值，該方法用三次樣條曲線通過

數(shù)據(jù)點，插值點處的值取對應曲線上的值；‘cubic’：分段三次Hermite插值，確定三次Hermite函數(shù)，

根據(jù)該函數(shù)確定插值點的函數(shù)值。缺省時表示分段線性插值。第十三頁，共二十三頁，2022年，8月28日x=0:6;y=cos(x);xi=0:.25:6;yi1=interp1(x,y,xi,'*nearest');yi2=interp1(x,y,xi,'*linear');yi3=interp1(x,y,xi,'*spline');yi4=interp1(x,y,xi,'*cubic');plot(x,y,'ro',xi,yi1,'--',xi,yi2,'-',xi,yi3,'k.-',xi,yi4,'m:')legend(‘原始數(shù)據(jù)’,‘最近點插值’,‘線性插值’,’樣條插值‘,’立方插值’)例6-19：用以上4種方法對y=cosx在[0,6]上的一維插值效果進行比較。第十四頁，共二十三頁，2022年，8月28日第十五頁，共二十三頁，2022年，8月28日二維多項式插值

zi=interp2（x，y，z，xi，yi，method）功能：已知數(shù)據(jù)點（x,y,z），運用method指定的方法，計算插值點（xi,yi）處得函數(shù)值zi。method指定的方法同一維多項式插值。第十六頁，共二十三頁，2022年，8月28日t=-2:0.5:2[x,y]=meshgrid(t);z=x.*exp(-x.^2-y.^2);

[x1,y1]=meshgrid(-2:0.1:2);z1=x1.*exp(-x1.^2-y1.^2);

figure(1)subplot(1,2,1),mesh(x,y,z),title(‘數(shù)據(jù)點')subplot(1,2,2),mesh(x1,y1,z1),title(‘函數(shù)圖象')例6-21：用以上4種方法對在[-2,2]上的二維多項式插值效果進行比較。第十七頁，共二十三頁，2022年，8月28日第十八頁，共二十三頁，2022年，8月28日[xi,yi]=meshgrid(-2:.125:2);zi1=interp2(x,y,z,xi,yi,'*nearest');zi2=interp2(x,y,z,xi,yi,'*linear');zi3=interp2(x,y,z,xi,yi,'*spline');zi4=interp2(x,y,z,xi,yi,'*cubic');figure(2)subplot(2,2,1),mesh(xi,yi,zi1),title(‘最近點插值')subplot(222),mesh(xi,yi,zi2),title(‘線性插值')subplot(223),mesh(xi,yi,zi3),title(‘樣條插值')subplot(224),mesh(xi,yi,zi4),title(‘立方插值')第十九頁，共二十三頁，2022年，8月28日第二十頁，共二十三頁，2022年，8月28日定義：對于結點,若函數(shù)s(x)滿足1)都在區(qū)間[a,b]上連續(xù)2)在每個小區(qū)間上是三次多項式3),i=0,1,2,…n則稱s(x)為三次樣條插值函數(shù)。第二十一頁，共二十三頁，2022年，8月28日定義：若不僅已知函數(shù)f(x)在結點上的函數(shù)值,還已知一階導數(shù)值yi’．若函數(shù)Ih(x)滿足1)都在區(qū)間[a,b]上連續(xù)2)Ih(x)在每個小區(qū)間上是次數(shù)不大于3的多項式3)（i=0,1,2,…n）則稱Ih(x)為分段三次Hermite插值多項式。第二十二頁，共二十三頁，2022年，8月28日北京科技大學數(shù)學實驗

下表是1971年到1990年我國總人口的統(tǒng)計數(shù)字，試根據(jù)1971年到1985年這15年人口的統(tǒng)計數(shù)字用多種方法預測未來20年的人口數(shù)字，并比較1986年到1990年間預測人口數(shù)字與實際統(tǒng)計數(shù)字的差異，在你所使用的幾種預測方法中找出一種較為合理的預測方法。年份人口統(tǒng)計數(shù)字年份人口統(tǒng)計數(shù)字19718.5229198110.007219728.7177198210.1654197