晕死,重装MATLAB后还是一样,出错....我的程序是这样写的:能帮忙看看吗<BR>
<P 0cm 0cm 0pt"><B normal">clear;clc<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">x=[12.5 12.625 14.125 14.5 17.25 17.75]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">y=[17 16.5 23 26.5 41 49]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">a=zeros(0,0)<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">b=zeros(0,0)<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">c=zeros(0,0)<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">A=zeros(0,0)<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">for i=1:6<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"> a=[a;1 x(i) x(i)^2 x(i)^3]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">end<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">for i=1:5<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"> A=[A;zeros(1,4*(i-1)) a(i,:) zeros(1,4*(5-i))]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"> A=[A;zeros(1,4*(i-1)) a(i+1,:) zeros(1,4*(5-i))]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">end<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">for i=2:5<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"> b=[b;0 1 2*x(i) 3*x(i)^2] <p></p></B></P>
<P 0cm 0cm 0pt"><B normal"> c=[c;0 0 2 6*x(i)]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">end<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">for i=1:4<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"> A=[A;zeros(1,4*(i-1)) b(i,:) -b(i,:) zeros(1,4*(4-i))]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"> A=[A;zeros(1,4*(i-1)) c(i,:) -c(i,:) zeros(1,4*(4-i))]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">end<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">A=[A;0 0 2 6*x(1) zeros(1,4*(4))]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">A=[A;zeros(1,4*(4)) 0 0 2 6*x(6)]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">d=y(1)<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">for i=2:5<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"> d=[d;y(i);y(i)]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">end<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">d=[d; y(6)]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">for i=1:10<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"> d=[d;zeros(1,1)]<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">end<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">X=A\d<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">X=X'<p></p></B></P>
<P 0cm 0cm 0pt"><B normal"><p> </p></B></P>
<P 0cm 0cm 0pt"><B normal">%picture<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">xi=[12.5:0.1:17.75];<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">yi=interp1(x,y,xi,'spline');<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">plot(xi,yi);<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">hold on<p></p></B></P>
<P 0cm 0cm 0pt"><B normal">scatter(x,y);<p></p></B></P>[em04] |