c=[];d=[];fork=4:12H=hilb(k);B=eye(k,k);L=sqrtm(H)*B;F=(inv(L))*H*(inv(L));c=[c,cond(H)];d=[d,cond(F)];endlc=log(d/c);n=4:12;p=polyfit(n,lc,1);plot(n,lc,\'*\',n,polyval(p,n));residdure=c-exp(polyval(p,n))???Errorusing==polyfitat48XandYvectorsmustbethesamesize.怎么改????急!!!!
|