matlab 平方差,第七章 最小平方差的方法（The Method of Least Square ）.doc

第七章 最小平方差的方法(The Method of Least Square ).doc

第七章 最小平方差的方法(The Method of Least Square )

本章說明利用最小平方差法, 找出適合一組資料的曲線.

在本章中包含Matlab 的m-file

cheby.m

genpoly.m

orthpoly.m

將須要的m-file之檔案夾加入搜尋路徑中

path('c:\numerical', path)

註: 如果你有安裝Matlab Notebook 要執行下列input cells (綠色指令敘述) 之前必須先執行上面的cell – [path (…) ]

藍色的內容是Matlab [output cells]

1. cheby.m – 利用Chebyshev Polynomials 做為polynomial fit的基底函數. 計算normal equations的係數矩陣A與向量b Aij = Σk=0m Ti(zk) Tj(zk) , bi = Σk=0m yk Ti(zk)

顯示cheby.m的內容

type cheby.

function [A ,b]=cheby(m,n,z, y)

%compute values of (A)ij, and (b)i

%

A=zeros(n,n) ; b=zeros(n,1);

T=chebyval(m,n,z) ; %get Tj(z)

for i=1:n

s=0;

for k=1:m

s = s+y(k)*T(i,k); %compute (b)i

end ;

b(i) =