myfun.asv

来自「能够非线性函数的曲线拟合」· ASV 代码 · 共 32 行

function [A,a,E]= myfun(T,K)
logT=log(T);
logK=log(K);
daoT=T.^(-1);
lenT=length(T);
A=ones(3);
x=[];
for r=1:lenT-2
    A(:,2)=logT(r:r+2);
    A(:,3)=daoT(r:r+2);
    b=logK(r:r+2);
    xt=A\b;
    xt=xt';
    x=[x;xt];
end
logA=mean(x(:,1));
A=exp(logA);
a=mean(x(:,2));
E=mean(x(:,3));

h1=stem(T,K,'bo');
    set(h1,'MarkerFaceColor','green');
    set(h1,'LineStyle','none');
    set(get(h1,'BaseLine'),'LineStyle','none');
    hold on;

for t=200:10:1300
    k=A*(t^a)*exp(E/t);
plot(t,subfun(t,A,a,E));
    
function k = subfun(t,A,a,E)
k=A*(t^a)*exp(E/t);