ChMatLab基础.ppt

* 程序Ch1_2_3.m * 程序mutiif.m * 程序Ch1_3.m * 程序for1.m * 程序for2.m * 程序for3.m * 程序for4.m * * 程序while1.m * * 程序while2.m * 程序Numtest.m * * 程序yy1.m * 程序yy2.m * 程序yy2.m * * * 程序YY21 * 关于素数的进一步讨论, 可参阅李尚志编写的数学实验》 * 程序yy4.m * * 程序yy5.m * * * 程序yy6.m * * * Fibonacci[1170-1250] * 程序文件ffib.m Eratosthenes筛法 古希腊的另一名学者给出了解决这一问题的方法, 即 后人称之为Eratosthenes筛法的方法. Eratosthenes筛法的基本思想方法是: 将自然数从2开 始排列至某一正整数 然后从中划去所有2的整数倍, （不包括2）再划去3的整数倍（不包括3）, …… , 将该 过程延续下去, 一直划到 剩下的数即是不超过 的 所有素数. 例 求 中素数, 并将这些素数写入一个矩阵中. 注意函数 的使用方法. 程序如下: 筛法 使用方法 ⑴ 找出向量 中所有非零元的位置. 结果为 ⑵ 按表达式去寻找向量中满足条件的元的位置. 在上例中若执行: 有没有比筛法更好的方法? 若已知素数 则 一定是偶数. 故只需 检查 是否为素数. 若它是素数, 则一定不能被它 前面的素数整除. 故只需用前面的素数加以判定即可. 该方法称为试除法. 当 很大时, 该方法将极为有效. 程序如下: 跳过所有的偶数 进一步的讨论: 若以 表示不超过 的素数个数, 表示区 间 中素数的个数, 通过计算: 观察计算结果, 你能有什么样的结论? 例 若两个连续自然数的乘积 是素数, 则称这两个连 程序如下: 续自然数是亲密数对. 该数是亲密素数. 求 中亲 密数对的对数及对应的亲密素数的和. 运行结果: 程序如下: 例 有一球以 的速度从水平线以 向斜上方抛 去, 该球着地后仍保持原方向弹跳, 但速度是原来的 以后每着地一次, 速度都比上次减少 求此球 在第三次着地前的轨迹. 例 用曲边梯形方法近似计算积分 其中 定积分的思想方法: 为求积分 首先将区间进行 等分, 在每个小区间上, 近似以梯形 面积 代替代替曲边梯形面积. 由此: 程序如下 积分值为: 可以通过数值积分方法的到相应的积分值: 程序如下 声称为整体变量 记录迭代次数 积分值为: 例 数列与递归函数 问题的提出 提出了这样一个问题: 一对 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的命令窗口中输入 计算结果: * 程序p1.m 程序：sort1.m * 程序gama.m * 程序gama.m * 程序sum1.m * 程序csum.m * 程序Ch1_1.m * 程序Ch1_2_1.m * * 程序Ch1_2_2.m 3.控制结构 脚本文件中的三种基本控制结构. 1.顺序结构 顺序结构表示是按照程序中语句的排列顺序依次执行, 直到程序的最后一句. 这是程序中最简单的形式. 例 从键盘上输入两数, 求出其平方和再将其输