matlab矩阵相加结果不对,用matlab算结果不对,求大神给看看。

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

这是我用的源代码,请大神帮忙看看,为啥拟合不出来。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值