依題意就變成了以下的方程組
%a+x(1)*b+x(1)^2*c+x(1)^3*d=y(1);
%a+x(2)*b+x(2)^2*c+x(2)^3*d=y(2);
%a+x(3)*b+x(3)^2*c+x(3)^3*d=y(3);
%寫成矩陣的形式:xx(3*4)*aa(4*1)=y(3*1) 其中3*1表示3行1列
x=[1 2 3];y=[2; 1; 5];%依題意寫入知道的X和Y的三組值
c=[0 1 2 3];
xx=zeros(3,4);
for i=1:3
xx(i,:)=x(i).^c;
end
aa=xx\y %因?yàn)槭且粋€超定方程,即方程個數(shù)小于未知數(shù)個數(shù)只能得到一個特解
結(jié)果:
aa =
3.3636
0
-2.1364
0.7727