信计102陈有超-实验报告.docVIP

数学建模实验?一、 实验目的与意义： 1、练习初等问题的建模过程；熟悉数学建模步骤 2、练习Matlab基本编程命令； ?二、 实验要求： 、较能熟练应用Matlab基本命令和函数； 、注重问题分析与模型建立，了解建模小论文的写作过程； 、提高Matlab的编程应用技能。 三、 实验学时数：学时 四、 实验类别：综合性 五、 实验内容与步骤： 练习： 基本命令 ：循环、绘图、方程（组）求解 作业： 1、某大学青年教师从31岁开始建立自己的养老基金，他把已有的积蓄10000元也一次性地存入，已知月利率为0.001（以复利计），每月存入700元，试问当他60岁退休时，他的退休基金有多少？又若他退休后每月要从银行提取1000元，试问多少年后他的基金将用完？（以复利计），第二月的养老基金为，然后利用递推程序计算出第个月的养老金，然后再换算成对应的年份。 递推公式： 程序：即建立m文件fun1（求退休时的退休基金）和fun2（求退休基金用完时间） （1）程序fun1： function y fun1 n n n-30; n n*12; x 10700*1.001; for i 2:n x x+700; x x*1.001; end y x; （2）程序fun2: function y fun2 n % n 317780; m 0; while n 0 n n*1.001; n n-1000; m m+1; end y m/12; 运行结果： 程序fun1运行结果： y fun1 31 y 1.8575e+004 注：31岁年终养老基金数，即18575元。 y fun1 60 y 3.1778e+005 注：60岁退休那年年终养老基金总数，即317780元。 程序fun2运行结果： y fun2 317780 y 31.9167 注：该青年老师的养老金在31.9167年后用完，即能用31年零11个月。 2、试对公平席位分配问题进行编程求解。 模型：表示按比例法分的人数，然后剩余的要分的人数按Q值法：， 其中表示第方的总人数，表示第方已占有的席位数。 程序：建立m文件fun3 function [x1,x2,x3] fun n1,n2,n3,x x1 floor x*n1/ n1+n2+n3 ; x2 floor x*n2/ n1+n2+n3 ; x3 floor x*n3/ n1+n2+n3 ; for i x1+x2+x3+1:x Q1 n1.^2/ x1* x1+1 ; Q2 n2.^2/ x2* x2+1 ; Q3 n3.^2/ x3* x3+1 ; if Q1 Q2Q1 Q3 x1 x1+1; end if Q2 Q1Q2 Q3 x2 x2+1; end if Q3 Q1Q3 Q2 x3 x3+1; end end 运行结果： [x1,x2,x3] fun3 103,63,34,21 x1 11 x2 6 x3 4 3、编程求解差分方程的阻滞增长模型，分别令从1.8逐渐增加，考察序列收敛、2倍周期收敛、4倍周期收敛……，直至一片混乱的情况，试以b为横坐标，收敛点为纵坐标作图。（与7.3节图8比较）。 模型：，其中的值可以给定，即初值也可以给定。 程序：建立M文件fun4和fun5 程序fun4： function y fun4 x,b y b*x* 1-x ; 程序fun5: function fun5 fun4,x0,b,n m 0; for i b 1 :b 3 :b 2 m m+1; y m,1 feval fun4,x0,i for j 2:n 2 y m,j feval fun4,y m,j-1 ,i ; end end plot [b 1 :b 3 :b 2 ],y :,n 1 +1:n 2 ,b- ; 运行结果： fun4 @fun5,0.2,[1.8,4.0,0.01],[200,300] ;