该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

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--

怎么把上面两个图像放在一个里面呢？谢谢了……