exm040406.m

来自「matlab电子书籍」· M 代码 · 共 13 行

x0=(0:0.1:1)';
y0=[-.447,1.978,3.11,5.25,5.02,4.66,4.01,4.58,3.45,5.35,9.22]';
m=length(x0);
n=3;				

X=zeros(m,n+1);
for k=1:n
	X(:,n-k+1)=(x0.^k);
end
X(:,n+1)=ones(m,1);

aT=(X\y0)'