[求助]新人求助.......

kk12

<P>麻烦大家看看这个程序，出错的地方用红色high light了<BR>这个程序是想求出线形方程组的系数的<BR>clear;clc<BR>x=input('Please input the set of data points x:');<BR>y=input('Please input the set of data points y:');<BR>a=zeros(0,0)<BR>b=zeros(0,0)<BR>c=zeros(0,0)<BR>A=zeros(0,0)</P>
<P>for i=1:6<BR>    a=[a;1 x(i) x(i)^2 x(i)^3]<BR>end</P>
<P>for i=1:5 <BR>    A=[A;zeros(1,4*(i-1)) a(i,:) zeros(1,4*(5-i))]<BR>    A=[A;zeros(1,4*(i-1)) a(i+1,:) zeros(1,4*(5-i))]<BR>end</P>
<P>for i=2:5<BR>    b=[b;0 1 2*x(i) 3*x(i)^2] <BR>    c=[c;0 0 2 6*x(i)]<BR>end</P>
<P>for i=1:4<BR>    A=[A;zeros(1,4*(i-1)) b(i,:) -b(i,:) zeros(1,4*(4-i))]<BR>    A=[A;zeros(1,4*(i-1)) c(i,:) -c(i,:) zeros(1,4*(4-i))]<BR>end</P>
<P>A=[A;0 0 2 6*x(1) zeros(1,16)]<BR>A=[A;zeros(1,16) 0 0 2 6*x(6)]</P>
<P>for i=1:14<BR>    <FONT color=#ff0000>y=[y;zeros(1,1)];<BR></FONT>end</P>
<P>X=A\y<BR></P>

suffer

<FONT color=#ff0000>y=[y zeros(1,1)];<BR></FONT>话说回来，你这个程序写的不太好，就算这次调通了下次很可能还是出错的，另外通用性很差<BR>

kk12

<P>謝謝幫忙，我再修改一下</P>[em04]