l3_1.m

利用幂法反幂法等进行程序运算,可以获得各种适用结果.

function l3_1
t=[1790:10:1980];
x=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76 92 106.5 123.2 131.7 150.7 179.3 204 226.5];
x1=log(x./x(1));
t1=t-t(1);
r=x1/t1;
tt=[min(t):max(t)];
xx=tt*r;
xx=x(1)*exp(r*(tt-t(1)));
plot(tt,xx)
hold on
plot(t,x,'.')