新修改建模实验指导书模板.doc

资料内容仅供您学习参考，如有不当或者侵权，请联系改正或者删除。实验目的作为实践性非常强的课程, 安排上机实验的目的, 不但是为了验证教材和授课内容, 更重要的是, 要经过实验深入理解方法的设计原理与处理问题的技巧, 培养自行处理常规数学模型的能力和综合运用知识分析、 解决问题的能力。1、 经过上机实验加深课堂内容的理解。计算机的应用在数学建模的教学中占有重要地位, 在为解决实际问题而建立数学模型的过程中、 对所建模型的检验以及大量的数值计算中, 都必须用到计算机。数学建模的实验课的目的和任务是经过实验培养并提高学生的数学建模能力和计算机应用能力。2、 学会对模型计算结果的分析和处理。数学建模实验不只是编写程序得到一个数值结果, 我们应在掌握数学模型的基本原理和思想的同时, 注意方法处理的技巧及其与计算机的密切结合, 重视对结果的分析与讨论。最终数值结果的正确性或合理性是第一位的, 当结果不正确、 不合理、 或误差大时, 我们要能够分析原因, 对算法、 计算方法、 或模型进行修正、 改进。3、 培养学生解决实际问题的能力。经过对实际问题的分析, 抓住问题本质, 培养学生将实际问题转化为数学问题能力, 要求经过数学实验的学习, 初步掌握将实际问题转化为数学问题的方法, 能够建立简单的实际问题的数学模型。同时要求学生经过查阅文献, 撰写符合要求的数学建模论文形式, 使学生论文写作能力等得到培养。实验基本要求一、 上机前的准备工作1、 复习和掌握与本次实验有关的教学内容。2、 根据本次实验要求, 根据本次实验要求, 按教材和任课老师介绍的方法完成数学建模实验任务, 对数学建模的各种基本类型和方法都作适度的练习, 并对学过的计算机编程语言在实验过程中进行全面的实践和提高。二、 上机实验步骤1、 启动开发环境; 2、 建立源程序文件, 输入源程序; 3、 编译产生目标程序, 连接生成可执行程序, 运行程序, 输出结果; 4、 对数值计算结果进行分析, 讨论其合理性与正确性; 5、 整理实验报告。三、 实验报告实验报告是记录实验工作全过程的技术文档, 实验报告的撰写是科学技术工作的一个组成部分。实验中, 学生要对问题进行分析, 计算, 编程, 处理在实验时记录的相关实验数据, 课后完成实验报告上交。实验项目实验一、 Dynamical SystemsDynamical systems can be used to describe almost all real-world systems. When time is discrete, the dynamical systems can be represented by difference equations, which usually are difficult to analyze mathematically but easy to do computations using computer softwares. The problem now is to figure out an appropriate model for a specific real-world system, i.e., a difference equation, that is consistent with the real-world data and that encapsulates the important aspects of the real-world system.一、 实验目的与要求掌握运用软件求解动态系统模型, 经过研究散点图得到动态系统的内在性质和长期趋势。经过对数据进行处理, 归纳出动态系统模型。1、 用Excel对数据进行处理, 建立动态系统模型而且进行验证; 2、 用Excel画散点图, 对动态系统模型解的长期趋势进行分析; 3、 用Excel求解动态系统模型并估计均衡点; 4、 用Excel分析多元动态系统模型。二、 算法实例Example 1.1 Six years ago your parents purchased a home by financing $80,000 for 20 years paying monthly payments of $880.87 with a monthly interest of 1%.(a) They have made 72 payments and wish to know how much they owe on the mortgage.(b) They are considering paying off the remaining mortgage in 10 years, find out the new monthly payments.解答如下, 因每个月还款而减少: bn=bn+1-bn=0.01bn-880.87求解下列动力系统模型: bn+1= bn+0.01bn880.87b0=80000Where bn 表示n个月后的欠款。因而有b1=80000+0.01(80000)-880.87=79919.13b2=79919.13+0.01(79919.13)-880.87=79837.45用序列表示B=(80000, 79919.13, 79837.45,)Excel操作步骤: 1打开excel表格, 输入如下表格。2用智能标识把月份拉到231。3.在B4 中输入= B3 +0.01B3-880.87, 回车后下拉即可可到序列B=(80000, 79919.13, 79837.45,)。4 在插入图表XY散点图, 选中数据格就可得出下表: Example 1.2You wish to buy a new car and narrow your choices to a Saturn, Cavalier, and Hyundai. Each company offers you its best deal: Saturn $13,990.00$1000 down3.5% interest for up to 60 monthsCavalier$13,550.00$1500 down4.5% interset for up to 60 monthsHyundai$12,400.00$500 dwon6.5% interset for up to 48 monthsYou are able to spend at most $475 a month on a car payment. Use a dynamical system to determine which car to buy: 解答如下, 对三家公司分别建立动力系统模型: Saturn: bn=bn+1-bn=0.035bn-475 bn+1= bn+0.035bn-475b0=12990Cavalier: bn=bn+1-bn=0.045bn-475bn+1= bn+0.045bn-475b0=12050Hyundai: bn=bn+1-bn=0.065bn-475bn+1= bn+0.065bn-475b0=11900SaturnCavalierHyundai月份 余额 月份 余额 月份 余额 012990.00 012050.00 011900.00 112969.65 112117.25 112198.50 212948.59 212187.53 212516.40 312926.79 312260.96 312854.97 412904.23 412337.71 413215.54 512880.87 512417.91 513599.55 612856.70 612501.71 614008.52 712831.69 712589.29 714444.08 812805.80 812680.81 814907.94 912779.00 912776.44 915401.96 1012751.27 1012876.38 1015928.09 1112722.56 1112980.82 1116488.41 1212692.85 1213089.96 1217085.16 1312662.10 1313204.00 1317720.69 1412630.27 1413323.18 1418397.54 1512597.33 1513447.73 1519118.38 1612563.24 1613577.88 1619886.07 1712527.95 1713713.88 1720703.67 1812491.43 1813856.00 1821574.41 1912453.63 1914004.52 1922501.74 2012414.51 2014159.73 2023489.35 2112374.02 2114321.92 2124541.16 2212332.11 2214491.40 2225661.34 2312288.73 2314668.52 2326854.33 2412243.84 2414853.60 2428124.86 2512197.37 2515047.01 2529477.97 2612149.28 2615249.13 2630919.04 2712099.50 2715460.34 2732453.78 2812047.99 2815681.05 2834088.27 2911994.67 2915911.70 2935829.01 3011939.48 3016152.73 3037682.90 3111882.36 3116404.60 3139657.29 3211823.24 3216667.80 3241760.01 3311762.06 3316942.86 3343999.41 3411698.73 3417230.28 3446384.37 3511633.18 3517530.65 3548924.36 3611565.35 3617844.53 3651629.44 3711495.13 3718172.53 3754510.35 3811422.46 3818515.29 3857578.52 3911347.25 3918873.48 3960846.13 4011269.40 4019247.79 4064326.13 4111188.83 4119638.94 4211105.44 42 7.69 4311019.13 4320474.84 4410929.80 4420921.21 4510837.34 4521387.66 4610741.65 4621875.10 4710642.61 4722384.48 4810540.10 4822916.79 4910434.00 4923473.04 5010324.19 5024054.33 5110210.54 5124661.77 5210092.91 5225296.55 539971.16 5325959.90 549845.15 5426653.09 559714.73 5527377.48 569579.75 5628134.47 579440.04 5728925.52 589295.44 5829752.17 599145.78 5930616.02 608990.88 6031518.74 Excel操作步骤和例1相似, 故我们能够得出应当购买Saturn公司的汽车。Example 1.3你的父母正在考虑月息0.5%、 数额为100000美元的抵押贷款。试建立一个用每月还款表示的模型, 使得在360次还款后就能还清贷款。经过计算数值解的实验来确保360月( 30年) 还清贷款的值。解答如下: 由题意, , 则有建立动力系统模型: 由, 可得下面我们取进行实验分析。省略中间的表格, 从最后的数据能够看出到最后一个月可还清。练习: 1、 你的信用卡上有月付利息1.5%的欠款500美元。你每月偿还50美元而且不再有新的贷款。试写出能对此问题的变化确切模型的动力系统公式。2、 下列数据表示从1790到 的美国人口数据yearPopulationyearPopulation17903,929,000.00190075,995,000.0018005,308,000.00191091,972,000.0018107,240,000.001920105,711,000.0018209,638,000.001930122,755,000.00183012,866,000.001940131,669,000.00184017,069,000.001950150,697,000.00185023,192,000.001960179,323,000.00186031,443,000.001970203,212,000.00187038,558,000.001980226,505,000.00188050,156,000.001990248,710,000.00189062,948,000.00 281,416,000.00 求出能够相当好地拟合该数据的动力模型, 经过画出模型的预测值和数据值来测试你的模型。实验二、 Testing Proportionality比例性关系验证是计算机数据处理的重要内容, 对模型结果具有重大的影响, 它在工程技术中有着广泛的应用。对实际问题而言, 求出比例性常数是一个极其重要而又比较困难的问题。一、 实验目的掌握运用软件验证线性关系, 经过研究散点图与趋势线得到线性模型。经过对数据进行处理, 归纳出合适的模型。1、 用Excel画散点图与趋势线, 对模型的线性关系进行分析; 2、 用Excel对数据进行处理, 归纳出合适的模型; 3、 用Excel计算不同模型的偏差, 选择合适的模型。二、 算法实例Example 2.1Determine whether the following data support a proportionality argument for yz1/2. If so, estimate the slope. y3.55678Z3691215解答如下: Excel操作步骤如下: 1. 打开excel表格, 输入X、 Y数据。2. 在在插入函数算中Power函数, 令Power=1/2, 然后导入Z数据, 得3.算出斜率We estimate the model to be 4.画出y与z1/2的线性图, 如下: Example 2.1木材切割者希望利用容易得到的测量数据来估计木材的板英尺数。她们测量树木腰高处的直径, 以构建一个将预测板英尺数与直径相联系起来的函数的模型。利用下列数据供检验之用。x17192023252832383941y19253257711131232522592941. 考虑两个不同的假设, 每个假设生成一个模型, 充分分析每个模型。（1） 假设所有的树木都是正圆柱体而且高度大致相同。（2） 假设所有的树木都是正圆柱体而且高度与直径成比列2. 哪个模型看起来好些? 为什么? 解答如下: 1、 当假设所有的树木都是正圆柱体而且高度大致相同时, 则成比列关系, 因此我们根据已知数据能够模拟出中的系数。当假设所有的树木都是正圆柱体而且高度与直径成比列时, 则成比列系, 因此我们根据已知数据能够模拟出中的系数。 下面进行excel数据处理, 找出比列关系, 如下结果经过线性拟合可得如下即我们可得到两种不同假设下的拟合方程: 与2、 为比较两个模型, 我们做如下处理其中分别为拟合误差, 做出误差分布图, 蓝为, 红为的结果。经过比较可知第二种模型较好。练习: 1.下面的数据中, 是美国黄松在树身中部测得的直径, 是体积的度量, 即用10除后的板英寸数, 变换数据画图, 检验模型, 并估计模型中参数和。实验三、 Least Square Fitting曲线拟合的最小二乘法是计算机数据处理的重要内容, 也是函数逼近的另一种重要方法, 它在工程技术中有着广泛的应用。对实际问题而言, 拟合曲线的选择是一个极其重要而又比较困难的问题, 必要时可由草图观察选取几种不同类型的拟合曲线, 再以其偏差小者为优, 经检验后再决定最后的取舍。一、 实验目的掌握运用软件进行最小二乘的参数估计, 经过合适的数据变换选择适当的线性回归模型。经过对模型偏差的计算, 选择合适的模型。1、 用Excel对数据进行最小二乘参数估计, 掌握不同的操作工具; 2、 用Excel对数据进行变换, 选择合适的线性回归模型; 3、 用Excel计算不同模型的偏差, 选择合适的模型。一、 算法实例Example 3.1Use the follow equations to estimate the coefficients of the line such that the sum of the squared deviations between the line and the following data points is minimized., the slope, the intercept(a) x1.02.33.74.26.17.0y3.63.03.25.15.36.8(b)x29.148.272.792.0118140165199y0.04930.08210.1230.1540.1970.2340.2740.328For each problem, compute and to bound 步骤如下: 1.打开excel, 输入X、 Y数据2.利用excel的各种求和函数, 计算出下列值3.根据公式计算=0.564, =2.215 4.利用工具数据分析( 如果没有, 可先加载宏) 回归, 做出X,Y的图像。对于( b) , 可完全类似( a) 的步骤, 得出下结果。请同学们自己验证。=864 ,=746496, =1.4414 , =2.0776,=1245.37, =117269.3, =194.99=0.0016, =0.0029Example 3.2求使一数据点集与二次形模型间偏差平方和极小化的方法。使用这些方法对下列数据集找出三个参数的估计。x0.10.20.30.40.5y0.060.120.360.650.95解答如下: 能够直接利用课本关于最小二乘估计公式计算三个参数, 下面我们用excel进行实验模拟, 如下图从而三个参数的大小可直接得出练习: 1.为拟合模型做一个适当的变换并估计参数和。71421283542841133250280297实验四Data Experiments 当我们构造一个预测模型时, 总需要细心分析收集到的数据, 看数据存在什么样的倾向, 然后将数据变换得到新的函数, 最好写出其偏差, 这对拟合数学模型非常重要。一、 实验目的掌握运用软件进行数据分析, 选择合适的数据变换, 建立数学模型。经过对模型偏差的计算, 比较不同的模型。1、 用Excel对数据进行分析, 选择合适的数据变换; 2、 用Excel构造差商表进行数据分析, 选择合适的低阶多项式模型;3、 用Excel计算模型的偏差, 比较不同的模型二、 算法实例Example 4.1The following data measure two characteristics of a ponderosa pine. The variable X is the diameter of the tree,in inches,measured at breast height;Y is a measure of volume-number of board feet divided by 10.fit a model to the data.The express Y in terms of X.解答如下: 1.打开excel, 输入X、 Y数据, 并做出散点图2.对X、 Y数据进行变换, 得到x2、 ln y、 sqrt y。3.分别已Y和x2, ln y和X, sqrt y和X做散点图, 如下4.从上述可知模型1的误差比较小。Example 4.2The following data represent the population of the United States from 1790 to . Construct a scatterplot of the given data. Is there a trend in the data? Are any of the data points outliers? Construct a divided difference table. Is smoothing with a low-order polynomial appropriate? If so, choose an appropriate polynomial and fit using the least-squares criterion of best fit. Anylyze the goodness of fit by examining appropriate indicators and graphing the model, the data points, and the deviations.YearObserved populationYearObserved population17903,929,000190075,995,00018005,308,000191091,972,00018107,240,0001920105,711,00018209,638,0001930122,755,000183012,866,0001940131,669,000184017,069,0001950150,697,000185023,192,0001960179,323,000186031,443,0001970203,212,000187038,558,0001980226,505,000188050,156,0001990248,709,873189062,948,000 281,416,000解答如下: 首先均差计算公式可得下列差分表divided difference table均差YearObserved population1790392,900 1800530,800 13,790 18107,240,000 670,920 32856.518209,638,000 239,800 -21556-1813.75183012,866,000 322,800 4150856.866766.76542184017,069,000 420,300 487524.16667-20.8175185023,192,000 612,300 9600157.53.333333186031,443,000 825,100 1064034.66667-3.07083187038,558,000 711,500 -5680-544-14.4667188050,156,000 1,159,800 22415936.537.0125189062,948,000 1,279,200 5970-548.167-37.1167190075,995,000 1,304,700 1275-156.59.791667191091,972,000 1,597,700 14650445.833315.058331920105,711,000 1,373,900 -11190-861.333-32.67921930122,755,000 1,704,400 16525923.833344.629171940131,669,000 891,400 -40650-1905.83-70.74171950150,697,000 1,902,800 505703040.667123.66251960179,323,000 2,862,600 47990-86-78.16671970203,212,000 2,388,900 -23685-2389.17-57.57921980226,505,000 2,329,300 -2980690.166776.983331990248,709,873 2,220,487 -5440.64-82.0212-19.3047 281,416,000 3,270,613 52506.271931.56450.33962根据excel中”工具数据分析回归”, 可得如下图像 y = 670127x2 - 3E+06x + 8E+06实验五、 Simulation计算机模拟与预测方法在材料科学研究中得到愈来愈广泛的应用。蒙特卡罗模拟有时能相对容易地近似很复杂的随机系统, 而且与分析模型的应用范围常常受到限制相比, 蒙特卡罗模拟能够很容易在更广泛的条件下估计候选方案的性能, 在建立模型时能够排除许多繁琐的工作。一、 实验目的掌握运用软件进行Monte Carlo方法模拟确定型现象和概率型现象, 掌握随机数的生成, 理解Monte Carlo模拟法在存贮模型和排队模型中的应用。1、 用Matlab进行Monte Carlo模拟, 编写程序计算面积与体积; 2、 用Matlab进行Monte Carlo模拟, 编写程序模拟抛硬币与掷骰子; 3、 用Matlab编写程序模拟存贮模型, 选择合理的进货量与进货周期; 4、 用Matlab编写程序模拟排队模型, 分析计算结果。二、 算法实例Example 5.1Using Monte Carlo simulation, write an algorithm to calculate an approximation to by considering the number of random points selected inside the quarter circlewhere the quarter circle is taken to be inside the squareUse the equation that =area /area.步骤如下: 1.打开matlab7.0, 创立一个新的 .m文件2.在文件中输入如下程序3.调入程序, 运行4. 最后得Pi=3.1359 Example 5.2Using Monte Carlo simulation, write an algorithm to calculate that part of the volume of an ellipsoid that lies in the fist octant, .步骤如下: 输入如下程序, 运行后得V= 107.298Example 5.3Modify the ship harbor system algorithm to keep track of the number os ships waiting in the queue.步骤如下: 输入如下程序得 HT =95.8300, WT =28.5700, IDT =0.1695, MHT =197, MWT =121练习: 1.将一个正规骰字连续郑100次, 用蒙特卡洛方法模拟总和的结果。2.修正存储算法以了解未满足的需求, 以及加油站短缺汽油的天数。