第七章 最小平方差的方法(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) =