最小二乘法的多項(xiàng)式擬合matlab實(shí)現(xiàn)

1、用最小二乘法進(jìn)行多項(xiàng)式擬合(matlab實(shí)現(xiàn))西安交通大學(xué)徐彬華算法分析：,1,2,3,.,m), 一共m+1個(gè)數(shù)據(jù)點(diǎn)，取多項(xiàng)式 P(x),對(duì)給定數(shù)據(jù)(X. . (i=0 使ffim2 心工 1p（巧f F=imnM j-0使得函數(shù)P(x)稱為擬合函數(shù)或最小二乘解，令似的心(工)=£/工=mm其中，aO, a1,a2,an為待求未知數(shù)，n為多項(xiàng)式的最高次幕，由此，該問題化 為求/二/(詢馬尸q的極值問題。由多元函數(shù)求極值的必要條件:得至上C1加 71二j«0 匕0一梟用=。j=0,1，,nj=0,1,n這是一個(gè)關(guān)于aO, a1,a2,an的線性方程組，用矩陣表示如下:Z福1

2、QmN巧i0mj0m/=0my* r】i - 0rrlZnz=0mz”r =0Z -r1 一Qm2廣f-o因此，只要給出數(shù)據(jù)點(diǎn)（j）及其個(gè)數(shù)m,再給出所要擬合的參數(shù) n,則即 可求出未知數(shù)矩陣（a0, a1,a2,an）試驗(yàn)題1編制以函數(shù)xk n n為基的多項(xiàng)式最小二乘擬合程序，并用于對(duì) k 0F列數(shù)據(jù)作三次多項(xiàng)式最小二乘擬合（取權(quán)函數(shù)wi=1）xi-1.0-0.50.00.51.01.52.0yi-4.447-0.4520.5510.048-0.4470.5494.552總共有7個(gè)數(shù)據(jù)點(diǎn)，令m=6第一步：畫出已知數(shù)據(jù)的的散點(diǎn)圖，確定擬合參數(shù)n;x=-1.0:0.520;y=-4.447,-0

3、.452,0.551,0.048,-0.447,0.549,4.552; plot(x,y,'*')xlabel 'x 軸ylabel 'y 軸title '散點(diǎn)圖'hold on321啜0-5-1-0.500.51152x軸散點(diǎn)圖5因此將擬合參數(shù)n設(shè)為3.第二步：計(jì)算矩陣A=f=0z-0m? (J削5注意到該矩陣為(n+1)*(n+1)矩陣,Vi-O的關(guān)系為i+j-2,由此可建立循環(huán)來求矩陣的各個(gè)元多項(xiàng)式的幕跟行、列坐標(biāo) 素，程序如下：m=6;n=3;A=zeros(n+1);for j=1:n+1for i=1:n+1for k=1:m+1A

4、(j,i)=A(j,i)+x(k)A(j+i-2) end end end;再來求矩陣B=B=0 0 0 0;for j=1:n+1 for i=1:m+1B(j)=B(j)+y(i)*x(i)A(j-1) end end第三步：寫出正規(guī)方程，求出a0,a1,an.B=B'a=inv(A)*B;第四步：畫出擬合曲線x=-1.0:0.000120;z=a(1)+a(2)*x+a(3)*xA2+a(4)*xA3;plot(x,z)legend(離散點(diǎn)','y=a(1)+a(2)*x+a(3)*x.A2+a(4)*xA3')title('擬合圖')總程

5、序附下：x=-1.0:0.520;y=-4.447,-0.452,0.551,0.048,-0.447,0.549,4.552; plot(x,y,'*')xlabel 'x 軸ylabel 'y 軸title '散點(diǎn)圖'hold onm=6;n=3;A=zeros(n+1);for j=1:n+1for i=1:n+1for k=1:m+1A(j,i)=A(j,i)+x(k)A(j+i-2) endend end;B=0 0 0 0;for j=1:n+1for i=1:m+1B(j)=B(j)+y(i)*x(i)A(j-1)endendB=B'a=inv(A)*B;x=-1.0:0.0001:2.0;z=a(1)+a(2)*x+a(3)*x.A2+a(4)*x.A3;plot(x