該樓層疑似違規已被系統折疊 隐藏此樓檢視此樓
A=[1790,3.9,
1800,5.3,
1810,7.2,
1820,9.6,
1830,12.9,
1840,17.1,
1850,23.2,
1860,31.4,
1870,38.6,
1880,50.2,
1890,62.9,
1900,76.0,
1910,92.0,
1920,106.5,
1930,123.2,
1940,131.7,
1950,150.7,
1960,179.3,
1970,204.0,
1980,226.5,
1990,251.4]
t0=A(1,1);x0=A(1,2);
t=A(:,1);x=A(:,2);
fun=@(a,t) a(1)*x0*exp(a(2)*(t-t0))./(x0*a(3)*exp(a(2)*(t-t0))+a(1)-x0*a(3))
a=nlinfit(t,x,fun,[400 0.1 2]);
xm=a(1),r=a(2),c=a(3)
tt=1790:1990;
plot(t,x,'o',tt,fun(a,tt))
format short g--
A=[1790,3.9,
1800,5.3,
1810,7.2,
1820,9.6,
1830,12.9,
1840,17.1,
1850,23.2,
1860,31.4,
1870,38.6,
1880,50.2,
1890,62.9,
1900,76.0,
1910,92.0,
1920,106.5,
1930,123.2,
1940,131.7,
1950,150.7,
1960,179.3,
1970,204.0,
1980,226.5,
1990,251.4]
t0=A(1,1);x0=A(1,2);
t=A(:,1);x=A(:,2);
fun=@(a,t) a(1)*x0*exp(a(2)*(t-t0))./(x0*exp(a(2)*(t-t0))+a(1)-x0)
a=nlinfit(t,x,fun,[400 0.1]);
xm=a(1),r=a(2)
tt=1790:1990;
plot(t,x,'o',tt,fun(a,tt))
format short g--
怎麼把上面兩個圖像放在一個裡面呢?謝謝了……