Ch4 程序设计初步.pptVIP

Fibonacci[1170-1250] 程序文件ffib.m 程序YY21 第一类文件又称为脚本文件. 程序sum1.m 程序p3.m 程序p2.m 程序mutiif.m 程序for1.m 程序for2.m 程序for4.m 程序for3.m 程序while1.m 程序while2.m 程序yy1.m 程序for4.m 程序yy6.m 程序yy3.m 关于零点问题的更进一步讨论，有兴趣的读者可参阅同济大学计算数学教研室编写的《现代数值数学和计算》。 程序yy6.m 程序yy5.m 程序yy2.m 运行结果: 程序如下: 例 有一球以 的速度从水平线以 向斜上方抛 去, 该球着地后仍保持原方向弹跳, 但速度是原来的 以后每着地一次, 速度都比上次减少 求此球 在第三次着地前的轨迹. 例 用曲边梯形方法近似计算积分 其中 定积分的思想方法: 为求积分 首先将区间进行 等分, 在每个小区间上, 近似以梯形 面积 代替代替曲边梯形面积. 由此: 程序如下 积分值为: 可以通过数值积分方法的到相应的积分值: 程序如下 声称为整体变量 记录迭代次数 积分值为: 例 数列与递归函数 问题的提出 提出了这样一个问题: 一对 13个月份中的兔子的对数如下表 小兔子二个月后可以生兔子, 而成熟兔子每月可生一对 小兔子. 假如去年12月底养一对小兔子, 问到今年年底 共有多少对兔子. 所示: 233 144 89 55 33 21 12 11 10 9 8 7 13 8 5 3 2 1 1 6 5 4 3 2 1 12 13个月份中的兔子的对数 以 表示各月份中兔子的对数, 则 定义 称数列为 为 数列, 若数列有关系: 法国数学家奇拉特在 死后400年时证明了 有关系 从而证明了 相应程序如下: 在MatLab的命令窗口中输入 计算结果: 四、Monte Carlo 方法 Monte Carlo方法是用概率统计的方法, 来研究一些具 体的数学问题. 例 求圆周率 问题描述 在矩形 中任取一个点, 则该点可能落在圆内, 其中 为 的面积. 也有可能落在圆外. 由几何概率知道: 落在区域 内的 概率为 为估计概率, 今产生随机数: 其中: 且随机变量 均服从区间 由此得到问题的解法. 上的均匀分布. 则 下面这段程序给出问题的求解方法. 记录有多少个点在圆内. 例 用Monte Carlo方法估计定积分 相应程序为: Monte Carlo方法综述 1.确定问题的概率模型; 2.从已知概率分布中抽样; 3.建立所需要的统计量. 若计算积分 该如何进行? 例 用Monte Carlo方法估计定积分 相应程序如下: 在MatLab下进行数值积 分, 得积分值: 例 利用Monte Carlo方法计算积分: 其中: 问题分析: 因在所给区域上, 有 记空间区域为 而积分所 所表示的空间区域为: 由此得相应程序为: 五、练习 1.编写程序, 实现功能: 由用户从键盘上输入3个数, 然 后按从大到小进行输出. 2.编写程序, 要求产生20个随机整数, 输出其中小于平均 值的偶数. 3.级数 求 使误差小于 4.编写程序, 求矩阵的乘积及矩阵的点积, 当阶数不同时, 显示信息:“矩阵阶数不同”并返回用户界面. 5.找出 中个位数为1的素数. 并记录其个数和它 们的和. 6.生成一个 的矩阵, 再将其中的元素按先行后列的 规则, 进行从大到小的排列. 7.考虑下面程序执行的结果: 六、实验 1.由键盘输入不超过100个正整数, 若不到100个, 以0结 结束. 输入结束后, 求出这些数中小于平均值的个数, 以 及这些数的平均数. 2.用Mont Carlo方法, 计算由两球体 所围成的体积. 并用数值积分进行验证. 3.由计算机产生