人口问题 matlab分析,人口问题数据拟合的matlab程序.doc

人口问题数据拟合的matlab程序.doc

人口问题数据拟合的MATLAB程序拟合拟合数据人口问题X1949195419591964196919741979198419891994Y546067708191981031131181线性模型用一阶多项式BPOLYFITX,Y,1ZB2B1XPLOTX,Y, R ,X,Z,XLABEL X 用矩阵运算AONESSIZEX ,X BA\Y ZB1B2XPLOTX,Y, R ,X,Z,XLABEL X 用线性回归AONESSIZEX ,X B,C,R,J,RREGRESSY ,AB回归系数C回归系数的置信区间R残差J拟合数据的置信区间R相关系数F值、P值ZB1B2XZ1ZJ,1 Z2ZJ,2 PLOTX,Y, R ,X,Z,X,Z1,X,Z2,XLABEL X ESQRTSUMYZ2/8ZZ1Z196EZZ2Z196EPLOTX,Y, R ,X,Z,X,ZZ1,X,ZZ22非线性模型YB2EXPB1X转化为线性函数AONESSIZEX ,X Y1LOGYB1,R,J,RREGRESSY1 ,ABEXPB11B12ZB1EXPB2XESQRTSUMYZ2/8Z1Z196EZ2Z196EPLOTX,Y, R ,X,Z,X,Z1,X,Z2用非线性函数拟合(缺点初值不合适，就得不到解)X49545964697479848994Y54606770819198103113118FUNINLINE B1EXPB2X , B , X B02001B,R,JNLINFITX,Y,FUN,B0ZB1EXPB2XPLOTX,Y, R ,X,ZNLINTOOLX,Y,FUN,B0拟合曲线图。预报X1949195419591964196919741979198419891994Y54606770819198103113118AONESSIZEX ,X B,C,R,J,RREGRESSY ,AZB1B2XESQRTSUMZY2/8X1X19992005ZZB1B2X1Y99ZZ11196EZZ11196EY05ZZ12196EZZ12196EZ1ZJ,1 Y991Y051Z2ZJ,2 Y992Y052PLOTX,Y, R ,X1,ZZ,X1,Z1,X1,Z2,XLABEL X 人员疏散问题X2550100200500Y1934495661B023参数初值FUNINLINE B1X/B2X , B , X 拟合函数B,R,JNLINFITX,Y,FUN,B0拟合函数的系数、残差ZB1X/B2XZ1ZJ,1 Z2ZJ,2 PLOTX,Y, R ,X,Z,X,Z1,X,Z2ESQRTSUMYZ2/3ZZ1Z196EZZ2Z196EPLOTX,Y, R ,X,Z, R ,X,ZZ1,X,ZZ2NLINTOOLX,Y,FUN,B0拟合曲线图。CLFX1949195419591964196919741979198419891994Y546677819198103111118PLOTX,Y, R , MARKERSIZE ,20,XLABEL X 解代数方程X1,X2SIZEXAONESX2,1,X BA\Y 线性回归BREGRESSY ,AY1B1B2XN1一阶多项式模型PPOLYFITX,Y,NY2POLYVALP,X指数模型转化为线性模型YLLOGYPLPOLYFITX,YL,1Y3EXPPL2EXPPL1XB2EXPPL2指数模型非线性拟合失败FUNINLINE B1EXPB2X , B , X B000001,30,DETA,R,JNLINFITX ,Y ,FUN,B0非线性拟合XX2550100200500YY1934495661BETA02,3FUN1INLINE B1XX/B2XX , B , XX BETA,R,JNLINFITXX,YY,FUN1,BETA0PLOTX,Y1,X,Y2,X,Y3,X,Y, R , MARKERSIZE ,20,XLABEL X 解方程零点的MATLAB语言零点解方程513005/WSQRT052W2用内联函数求零点FUNINLINE 513005/WSQRT052W2 , W W0FZEROFUN,03用函数文件求零点FUNCTIONYFUNWY513005/WSQRT052W2W0FZEROFUN,03用循环语句搜寻零点FORI1100W030001IY3005/WSQRT052W2IFY51W0W,BREAKENDENDW0一元函数直接赋值做图W02000104Y513005/WSQRT052W2XZEROSSIZEWPLOTW,Y,W,X,GRID一元函数用内联函数赋值做图YINLINE 513005/WSQRT052W2 , W W0100105Y_CHARVECTORIZEYYFY_CHAR,WPLOTW,Y求解微分方程的MATLAB语言程序微分方程一阶微分方程求表达式(符号运算)SDSOLVE DXRX1X/K S1DSOLVE DXRX1X/K , X0001 转换为函数值R03K8C1001SSUBSST0100SSSUBSS, T ,TPLOTT,SS,GRID求数值解建立函数文件FUNMFUNCTIONYFUNX,TY03X1X/8调用ODE23二、三阶龙格库塔法解微分方程T,XODE23FUN,0,20,001PLOTT,X高阶微分方程初值问题求数值解要化为一阶方程组建立M文件FUNSMFUNCTIONFYYFUNST,YFYYY11Y22Y2Y1T,YODE23FUNS ,0,25,025,0PLOTT,YPLOTT,Y,1,T,Y,2, R 高阶微分方程边值问题ZDSOLVE XD2Y3DYX2 , Y10 , Y50 , X X10015SSUBSZPLOTX,SLORENTZ方程组周期轨线FUNCTIONDOTODONET,YDOT10Y2Y128Y1Y2Y1Y3Y1Y28/3Y3TT,YYODE45ODONE,0200,123PLOT3YY,1,YY,2,YY,3奇怪吸引子FUNCTIONDOT1ODONEQT,YDOT110Y2Y1328Y1Y2Y1Y328Y1Y2Y1Y3Y1Y28/3Y3T1,Y1ODE45ODONEQ,0200,123PLOT3Y1,1,Y1,2,Y1,3, R 试卷合理的均衡分配方案的计算机程序MATLAB原题在大学生数学建模竞赛的评卷工作中，M个评委(M个评委来自不同的学校)要完成N份试卷的打分，竞赛试卷来自K个学校，第J个学校有竞赛试卷LJ份，N∑LJ(J＝1K)为节省人力，每份试卷只要由其中P(P32HH32ENDENDMH,JIHH1IFH33H1JJ1ENDENDABIFN40BREAKENDBBLN1ENDMM即所求矩阵求数列2/13/25/38/513/821/13。。。。的前N(＝100)项和求数列2/1，3/2，5/3，8/5，13/8，21/13。。。的前N(＝100)项和A12A23B11B22N100SUM0FORI1NAI2AIAI1BI2BIBI1SUMSUMAI/BIENDSUM结果SUM1621030差分方程的MATLAB解法差分方程的一般形式为A(N1＝RA(N＋B计算程序A1A0