神马品牌网隔壁的女孩完整版在线:求一段代码的解释,有关matlab的
来源：百度文库 编辑：神马品牌网 时间：2024/05/03 01:03:11
对于以下代码,我不能翻译成数学公式,希望哪位有兴趣的朋友帮个忙!
>> T=[0.25,0.5,0.75,1,1.5,2,2.5,3,3.5,4,4.5,5,6,7,8,9,10,11,12,13,14,15,16];
>> C=[30,68,75,82,82,77,68,68,58,51,50,41,38,35,28,25,18,15,12,10,7,7,4];
>> plot(T,log(C))
>> b=polyfit(T(6:end),log(C(6:end)),1)
b =
-0.1954 4.7912
>> a=polyfit(T(1:5),log(exp(b(2)+b(1).*T(1:5))-C(1:5)),1)
a =
-1.8558 4.7692
>> fun=inline('e(1).*(exp(e(2).*T)-exp(e(3).*T))','e','T');
>> ee=nlinfit(T,C,fun,[exp(4.8),-0.20,-1.86])
ee =
114.4331
-0.1855
-2.0079
>> C2=ee(1).*(exp(ee(2).*T)-exp(ee(3).*T));
>> plot(T,C,'*',T,C,'r',T,C2)
>> T=[0.25,0.5,0.75,1,1.5,2,2.5,3,3.5,4,4.5,5,6,7,8,9,10,11,12,13,14,15,16];
>> C=[30,68,75,82,82,77,68,68,58,51,50,41,38,35,28,25,18,15,12,10,7,7,4];
>> plot(T,log(C))
>> b=polyfit(T(6:end),log(C(6:end)),1)
b =
-0.1954 4.7912
>> a=polyfit(T(1:5),log(exp(b(2)+b(1).*T(1:5))-C(1:5)),1)
a =
-1.8558 4.7692
>> fun=inline('e(1).*(exp(e(2).*T)-exp(e(3).*T))','e','T');
>> ee=nlinfit(T,C,fun,[exp(4.8),-0.20,-1.86])
ee =
114.4331
-0.1855
-2.0079
>> C2=ee(1).*(exp(ee(2).*T)-exp(ee(3).*T));
>> plot(T,C,'*',T,C,'r',T,C2)
嗯……大概是多边形拟合和线性拟合？貌似没有什么公式可以翻译的，你看看matlab的帮助，里面对polyfit和nlinfit函数的解释
plot是绘图，不用管它，行首没有>>的行是matlab的输出结果，也不用管
以上。
手头没有matlab的老狼