clear all
%¶¨òå2Îêy
w0=778.06;
h=1.0545*10^(-34);
k=1.38*10^(-23);
X=[93
123
153
183
213
243
273
303
333
363
393
423
453
483
513
573
603
663]; %%%%% XÎaT
Y=[777.421
777.326
777.219
777.084
776.965
776.656
776.166
775.667
775.07
774.664
774.394
773.885
773.407
772.861
772.31
771.138
770.541
769.219]; %%%%%%%%%%% Y=w_T
X0=[14 5.52*10^14]';
%options=optimset('Largescale','off');
%x=lsqnonlin(@fit_simp,X0,[],[],options,X,Y); %%% êä3öBoíw1
[x,res]=lsqcurvefit(@fit_simp,X0,X,Y); %%% êä3öBoíw1
res
BB=x(1);
ww1=x(2);
w00=5.4*10^14
fun=inline('3.09*10^(-6)+2.63*10^(-9).*t-1.08*10^(-12)*t.^2+2*(3.21*10^(-6)+3.56*10^(-9).*t-1.62*10^(-12)*t.^2)','t');
%fun=inline('3.09*10^(-6)+2.63*10^(-9).*t-1.08*10^(-12)*t.^2','t');
for i=1:length(X)
delta(i)=quadl(fun,0,X(i));
end
delta0=-w0*1.56*delta';
ww2=5.52*10^14;
n1=1./(exp(h.*ww1./(k.*X))-1);
n2=1./(exp(h.*ww2./(k.*X))-1);
Y_new=w0+BB*(1+n1+n2)+delta0;
%plot(X,delta0,'+')
plot(X,Y,'*')
hold on
plot(X,Y_new)
% gamma1=1;
%%%μ÷óÃoˉêycaculate %Dèòa×¢òaμÄêÇ ¶¨òåμÄoˉêyÃû3ƱØDëè«2¿êÇó¢ÎÄ×Öĸ
% delta0=caculINT(100,w0,gamma1);
% w_result=w0+delta0+Y_new
这是我用的源代码,请大神帮忙看看,为啥拟合不出来。