原创力文档- 1、本文档共35页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
实验9-人口预测与数据拟合
;实验13 人口预测与数据拟合;据人口统计年鉴，知我国从1949年至1994年人口数据资料如下： (人口数单位为：百万);如何确定a,b?;1 、曲线拟合问题的提法: ;数据插值;２、 用什么样的曲线拟合已知数据?;３、 拟合函数组中系数的确定;MATLAB中的数据拟合命令;曲线拟合举例;实验过程：;MATLAB中的数据拟合命令;多项式拟合举例;实验过程：;二、人口预测线性模型 ; 英国人口学家Malthus根据百余年的人口统计资料，于1798年提出了著名的人口自然增长的指数增长模型。;1. 由前100年的数据求出美国的人口增长Malthus模型。;解：;％ sy09_3.m %
% This program is to predict the number of population %
format long
t1=[1790;1800;1810;1820;1830;1840;1850;1860;1870;1880];
t2=[1890;1900;1910;1920;1930;1940;1950;1960;1970;1980];
x1=[3.9;5.3;7.2;9.6;12.9;17.1;23.2;31.4;38.6;50.2];
x2=[62.9;76.0;92.0;106.5;123.2;131.7;150.7;179.3;204.0;226.5];
lnx1=log(x1); lnx2=log(x2);;a12=sum(t1);a11=10;a21=a12;a22=sum(t1.^2); d1=sum(lnx1);d2=sum(lnx1.*t1); A=[a11,a12;a21,a22];D=[d1;d2]; ab=inv(A)*D; disp(a=);disp(ab(1)); disp(b=);disp(ab(2)); for i=1:10 xx1(i)=exp(ab(1)+ab(2)*t1(i)); end for i=1:10 xx2(i)=exp(ab(1)+ab(2)*t2(i)); end plot(t1,x1,r*--,t1,xx1,b+-, t2,x2,g*--,t2,xx2,m+-); ;a= -49.79535457790735
b=0.02859807120038;四、人口预测的Logistic模型;假设;例2 中国人口预测;中国人口Logistic模型—zgrk_lgstic.m;中国人口预测的Logistic模型;多项式拟合人口模型——sy09_04.m;a =1.0e+006 *
-0.00000000000014
0.00000000107892
-0.00000304878595
0.00381927346813
-1.79012132225427;例3: 海底光缆线长度预测模型;;解： 用12次多项式函数拟合光缆走势的曲线图如下;% sy09_05.m This program is to fit the data by polynomial %
format long
t=linspace(0,20,21);
x=linspace(0,20,100);
P=[9.01,8.96,7.96,7.97,8.02,9.05,10.13,11.18,12.26,13.28,13.32,12.61,11.29,10.22,9.15,7.90,7.95,8.86,9.81,10.80,10.93];
[a,s]=polyfit(t,P,12);
yy=polyval(a,x);
disp(yy=);disp(yy);
plot(x,yy,r*--,t,P,b+-);
L=0;
for i=2:100
L=L+sqrt((x(i)-x(i-1))^2+(yy(i)-yy(i-1))^2);
end
disp(The length of the label is L=);disp(L);;上机实验任务
文档评论（0）