設(shè)置m文件
n&DRh.@ function A=nihe(x,y,n)%定義多項(xiàng)式擬合系數(shù)求解
函數(shù),% x、y 為輸入數(shù)據(jù)量,n 為擬合次數(shù)
!$ItBn/_ m=length(x); %測(cè)量數(shù)據(jù)長(zhǎng)度
p~1!O]qLt X1=zeros(1,2*n); %生成 X 矩陣
;4tmnC>OnA for i=1:2*n
BU
|]4 X1(i)=sum(x.^i); end
A2LqBirkl X2=[m,X1(1:n)];
8VxjC1v+ X3=zeros(n,n+1);
'mx_]b^O for j=1:nX3(j,:)=X1(j:j+n);
wcDHx#~ end X=[X2;X3];
^LO=&Cq Y=zeros(1,n); %生成 Y 向量
vs]#?3+ for k=1:nY(k)=sum(x.^k.*y); end
?nN3K Y=[sum(y),Y]; Y=Y';
A2.[P== A=X\Y; %求得擬合系數(shù)向量 A
x