MATLAB语言程序设计训练题.docVIP

2013年数学素质拓展选修暨数学建模培训班 MATLAB语言程序设计训练题 1．编写脚本式M0时结束输入，求这些数的平均值和它们之和。 2．编写脚本式M文件实现：输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。 3．一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。编写输出全部水仙花数的脚本式M编写救援问题的MATLAB语言程序。 题目描述：假设洪水把一个村庄淹没了，只有那些比较高的房屋的屋顶露出水面，人们已经爬上屋顶等待求援。如图，在图的原点是大本营，救生船每次从大本营出发，救了人之后将人送回大本营。图中的点代表屋顶，每个屋顶由其位置坐标和其上的人数表示。救生船每次从大本营出发，以速度50米/分钟驶向下一个屋顶，达到一个屋顶后，救下其上的所有人，每人上船1分钟，船原路返回，达到大本营，每人下船0.5分钟。假设原点与任意一个屋顶的连线不穿过其它屋顶。假设屋顶数为50，随机给定每个屋顶的坐标和人数，求出所有人都到达大本营并登陆所用的时间。并输出所有屋顶的坐标和人数。 题目分析：这是一个求和问题，单独计算出每个屋顶上的人到达大本营的时间，然后将它们相加即可。计算公式如下： ， 其中，totalTime是要求的总时间，N是屋顶数（本题为50），speed是船行驶的速度（本题为50米/分钟），(xi,yi)是第i个屋顶的坐标，pi是第i个屋顶上的人数，因为要在计算后输出所有屋顶位置和人数所以需要保存所屋顶和人数的信息。 程序设计： S1：定义N行4列的数组roofs贮存屋顶信息，其中，第1列为屋顶的横坐标，第2列为屋顶的纵坐标，第3列为屋顶上的人数，第4列为营救完该屋顶的人所需的时间，初始值均为0。 S2：确定输入参数和输出参数。数组roofs的第1，2，3列数据、船行驶的速度speed、每个被救人员上船的时间up和下船的时间down作为输入参数。完成数组roofs第4列的计算后，数组roofs和营救所有屋顶的人所需的时间totalTime作为输出参数。 S3：计算营救每个屋顶的人所需的时间。 S4：计算营救所有屋顶的人所需的时间totalTime。 MATLAB语言程序： S2：把村民编号，按编号报上出生年份、出生月份、出生日期，并计算其出生年份加月份加日期的和。 S3：求出生年份加月份加日期的和与M相同的村民的编号，即幸运者编号。 S4：求获奖者人数。 S5：最后，如果有人获奖，则输出获奖者编号和获得的奖金数额；如果没有人获奖，就输出没有人获奖的信息。 MATLAB语言程序： 与矩阵的积的函数式M文件(不能直接使用运算“”)。 7．描述刘徽割圆术计算π的原理、方法和计算步骤，并编写实现计算的函数式M文件。采取不同的分割计算π的近似值，并将计算的结果与较准确的π值进行比较，对算法进行分析。 刘徽割圆术计算π的原理：从正六边形开始，逐步求正边形的边长与面积，并用它的面积近似其外接圆面积，从而得到π的近似值。 递推法：设边数为的正多边形边长为，如图，，所以 相应正多边形面积 . 从正六边形开始，逐步求正边形的边长与面积即可。 MATLAB语言程序：π值. digits(100) PI=vpa(pi) %输入n的值，调用函数calpi1计算π的近似值. n=input(‘请输入割圆的次数n：’); y=calpi1(n) 8．编写采用级数展开式 和 计算π的函数式M文件。调用编写的程序，尝试n取不同整数以求π的近似值，并将计算的结果与较准确的π值进行比较，对算法进行分析和比较。 9．描述Monte Carlo方法计算π的原理、方法和计算步骤，并编写实现计算的函数式M文件。调用编写的程序计算π的近似值，并将计算的结果与较准确的π值进行比较，并对算法进行分析。 Monte Carlo方法计算π的原理：在正方形0 x 1, 0 y1上随机的投大量的点，那么落在四分之一圆内的点数m与落在正方形内的点数n之比m/n 应为这两部分图形面积之比π/4，故π＝4 m/n。 计算方法和步骤：采用计算机模拟，产生区间[0, 1]上数目为n的一组随机数(x, y)，计算满足的点数m，用4 m/n的值作为π的近似值。 MATLAB语言程序：π值. digits(100) PI=vpa(pi) %输入n的值，调用函数calpi1计算π的近似值. n=input(‘请输入割圆的次数n：’); y=calpi4(n) 10．给出两种计算超越数e的方法，并通过数值计算实验进行算法分析。 方法一：函数在0点处的台劳展式为 . 取得, 因而。 方法二：当n取较大的值，用近似。 完成第6至10，请参考乐经良、向隆万、李世栋等编著，高等教育出版社出版