原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
人口问题数据拟合的MATLAB程序
人口问题数据拟合的MATLAB程序拟合
%拟合数据? 人口问题 x=[1949? 1954? 1959?? 1964? 1969? 1974? 1979? 1984? 1989? 1994]; y=[5.4? 6.0? 6.7? 7.0? 8.1? 9.1?? 9.8?? 10.3? 11.3?? 11.8]; %? 1 线性模型 %用一阶多项式 b=polyfit(x,y,1) z=b(2)+b(1).*x; plot(x,y,r*,x,z),xlabel(x) %用矩阵运算 A=[ones(size(x)), x]; b=A\y z=b(1)+b(2).*x; plot(x,y,r*,x,z),xlabel(x) %用线性回归 A=[ones(size(x)), x]; [b,c,r,j,R] =regress(y,A) % b 回归系数? c 回归系数的置信区间? r 残差? j 拟合数据的置信区间? R 相关系数 F值、p值 z=b(1)+b(2).*x; z1=z+j(:,1); z2=z+j(:,2); plot(x,y,r*,x,z,x,z1,x,z2),xlabel(x) e=sqrt(sum((y-z).^2)/8) zz1=z-1.96*e;? zz2=z+1.96*e; plot(x,y,r*,x,z,x,zz1,x,zz2)
% 2 非线性模型 y=b(2)exp(b(1)x) %转化为线性函数 A=[ones(size(x)), x]; y1=log(y); [b1,r,j,R]=regress(y1,A) b=[exp(b1(1)) b1(2)] z=b(1).*exp(b(2).*x); e=sqrt(sum((y-z).^2)/8) z1=z-1.96*e;? z2=z+1.96*e plot(x,y,r*,x,z,x,z1,x,z2) %用非线性函数拟合（缺点 初值不合适，就得不到解） x=[49? 54 59?? 64? 69 74 79? 84? 89? 94]; y=[5.4? 6.0? 6.7? 7.0? 8.1? 9.1?? 9.8?? 10.3? 11.3?? 11.8]; fun=inline(b(1).*exp(b(2).*x),b,x); b0=[2? 0.01]; [b,r,j]=nlinfit(x,y,fun,b0) z=b(1).*exp(b(2).*x); plot(x,y,r*,x,z) nlintool(x,y,fun,b0) %拟合曲线图。 %预报 x=[1949? 1954? 1959?? 1964? 1969? 1974? 1979? 1984? 1989? 1994]; y=[5.4? 6.0? 6.7? 7.0? 8.1? 9.1?? 9.8?? 10.3? 11.3?? 11.8]; A=[ones(size(x)), x]; [b,c,r,j,R] =regress(y,A) z=b(1)+b(2).*x; e=sqrt(sum((z-y).^2)/8) x1=[x 1999 2005] zz=b(1)+b(2).*x1; y99=[zz(11)-1.96*e? zz(11)+1.96*e] y05=[zz(12)-1.96*e zz(12)+1.96*e] z1=[z+j(:,1) y99(1) y05(1)]; z2=[z+j(:,2) y99(2) y05(2)]; plot(x,y,r*,x1,zz,x1,z1,x1,z2),xlabel(x)
%人员疏散问题 x=[25 50 100 200 500]; y=[1.9 3.4 4.9 5.6 6.1]; b0=[2 3]; %参数初值 fun=inline( b(1).*x./(b(2)+x),b,x); %拟合函数 [b,r,j]=nlinfit(x,y,fun,b0) %拟合函数的系数、残差 z=b(1).*x./(b(2)+x);? z1=z+j(:,1); z2=z+j(:,2); plot(x,y,*r,x,z,x,z1,x,z2) e=sqrt(sum((y-z).^2)/3) zz1=z-1.96*e; zz2=z+1.96*e; plot(x,y,*r,x,z,r,x,zz1,x,zz2) ? nlintool(x,y,fun,b0) %拟合曲线图。
clf x=[1949 1954 1959 1964 1969 1974 1979 1984 1989 1994]; y=[5.4 6 6.7 7 8.1 9.1 9.8 10.3 11.1 11.8]; plot(x,y,.r,MarkerSize,20),xlabel(x)
?
%解代数方程 [x1,
文档评论（0）