代码
clear;clc;
n=[1 2 4 8];
t=[ 7.21 6.88 6.32 5.84];
logt=log(t);logn=log(n);
p=polyfit(logn,logt,1);
beta=p(1)
alfa=exp(p(2))
x=[0:0.01:10]' ;
y=alfa*(n.^beta);
fprintf('模型:logt= log %.4f +%.4f logn\n',alfa,beta);
s='n 实际值 模型计算出的值\n';
fprintf(s);
for i=1:3
fprintf('%d %.2f %.4f\n',n(i),t(i),y(i));
end
y=alfa*(x.^beta);
plot(n,t,'r*',x,y,'b')
结果