%参数r,Nm估计（Matlb代码）
t=0:1:38;%令1981年为0，2018年为38，间隔为1年
n=[9.94,10.09,10.23,10.37,10.51,10.67,10.84,11.02,11.19,...
11.35,11.51,11.65,11.78,11.92,12.05,12.18,12.30,12.42,...
12.53,12.63,12.72,12.80,12.88,12.96,13.04,13.11,13.18,...
13.25,13.31,13.38,13.44,13.51,13.57,13.64,13.71,13.79,...
13.86,13.93];%1981年到2018年的人口数据
n1=[9.94,10.09,10.23,10.37,10.51,10.67,10.84,11.02,11.19,...
11.35,11.51,11.65,11.78,11.92,12.05,12.18,12.30,12.42,...
12.53,12.63,12.72,12.80,12.88,12.96,13.04,13.11,13.18,...
13.25,13.31,13.38,13.44,13.51,13.57,13.64,13.71,13.79,...
13.86];%1981年到2017年的人口数据
n2=[10.09,10.23,10.37,10.51,10.67,10.84,11.02,11.19,11.35,...
11.51,11.65,11.78,11.92,12.05,12.18,12.30,12.42,12.53,...
12.63,12.72,12.80,12.88,12.96,13.04,13.11,13.18,13.25,...
13.31,13.38,13.44,13.51,13.57,13.64,13.71,13.79,13.86,...
13.93];%1982年到2018年的人口数据
dn=(n2-n1)./n2;%每一年的人口增长率
a=polyfit(n2,dn,1);%最小二乘法的拟合公式
r=a(2),Nm=-r/a(1) %r,Nm的值
我改成这样就不会改了
from numpy import polyfit
n = [9.94, 10.09, 10.23, 10.37, 10.51, 10.67, 10.84, 11.02, 11.19,
11.35, 11.51, 11.65, 11.78, 11.92, 12.05, 12.18, 12.30, 12.42,
12.53, 12.63, 12.72, 12.80, 12.88, 12.96, 13.04, 13.11, 13.18,
13.25, 13.31, 13.38, 13.44, 13.51, 13.57, 13.64, 13.71, 13.79,
13.86, 13.93] # 1981年到2018年的人口数据
n1 = [9.94, 10.09, 10.23, 10.37, 10.51, 10.67, 10.84, 11.02, 11.19,
11.35, 11.51, 11.65, 11.78, 11.92, 12.05, 12.18, 12.30, 12.42,
12.53, 12.63, 12.72, 12.80, 12.88, 12.96, 13.04, 13.11, 13.18,
13.25, 13.31, 13.38, 13.44, 13.51, 13.57, 13.64, 13.71, 13.79,
13.86] # 1981年到2017年的人口数据
n2 = [10.09, 10.23, 10.37, 10.51, 10.67, 10.84, 11.02, 11.19, 11.35,
11.51, 11.65, 11.78, 11.92, 12.05, 12.18, 12.30, 12.42, 12.53,
12.63, 12.72, 12.80, 12.88, 12.96, 13.04, 13.11, 13.18, 13.25,
13.31, 13.38, 13.44, 13.51, 13.57, 13.64, 13.71, 13.79, 13.86,
13.93] # 1982年到2018年的人口数据
dn = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
for t in range(0, 36): # 令1981年为0，2018年为38，间隔为1年
dn[t] = (n2[t] - n1[t]) / n2[t] # 每一年的人口增长率
a = polyfit(n2, dn, 1) # 最小二乘法的拟合公式
r = a[2]
Nm = -r / a[1] # r,Nm的值