最小二乘法matlab多项式拟合

eval(['d',oo,'=diff(w^2,name_v(i).f1)']) end

for k=1:13

kk=name_v(k).f1; for j=1:length(x);

eval(['save_',num2str(j),'=subs(',['d',kk],',[t,y],[x(j),z(j)])']); end

eval(['sd',kk,'=save_1+save_2+save_3+save_4']); end

strings='';

strings=[strings,'sd',kk,'=0'];

[j1 j2 j3 j4 j5 j6 j7 j8 j9 j10 j11 j12 j13]=solve([char(sda0),'=0'], [char(sda1),'=0'],[char(sda2),'=0'],[char(sda3),'=0'],[char(sda4),'=0'],[char(sdw1),'=0'],[char(sdw2),'=0'],[char(sdw3),'=0'],[char(sdw4),'=0'],[char(sds1),'=0'],[char(sds2),'=0'],[char(sds3),'=0'],[char(sds4),'=0'],'a0','a1','a2','a3','a4','w1','w2','w3','w4','s1','s2','s3','s4');

% ['a0' 'a1' 'a2' 'a3' 'a4' 'w1' 'w2' 'w3' 'w4' 's1' 's2' 's3' 's4']=solve('sda0','sda1','sda2','sda3','sda4','sdw1','sdw2','sdw3','sdw4','sds1','sds2','sds3','sds4','a0' 'a1' 'a2' 'a3' 'a4' 'w1' 'w2' 'w3' 'w4' 's1' 's2' 's3' 's4'];