- 1、本文档共28页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数 学 实 验
报 告
实验五 迭代(一)——方程求解
学院:数学与信息科学学院
班级:09级数学(4)班
姓名:***
学号:***
实验五 素数
实验名称 迭代(一)——方程求解 实验目的 用迭代法求解方程及线性方程组 实验环境 Mathematica4.0系统 实验的基本理论与方法 方程求根
给定实数域上光滑的实值函数f(x)以及初值定义数列
(1)
称为f(x)的一个迭代序列。
给定迭代函数f(x)以及一个初值利用(1)迭代得到数列如果数列收敛于一个,则有
. (2)
即是方程x=f(x)的解。由此启发我们用如下的方法球方程g(x)=0的近似解。
将方程g(x)=0改写为等价的方程x=f(x), (3)
然后选取一初值利用(1)做迭代。迭代数列收敛的极限就是方程g(x)=0的解。
线性方程组的迭代求解
给定一个n元线性方程组
(6)
或写成距阵的形式Ax=b, (7)
其中是n阶方程, 及均为n维列向量。熟知,当距阵A的行列式非零时方程(7)有唯一的解。迭代法是求解这些问题的有效方法之一。用迭代的方法求解线性方程组的思想与上节介绍的方程求解的方法是类似的。假设我们可以将方程组(7)改写成x=Mx+f
其中是n 阶距阵,是n维列向量。任意给定初始向量,由迭代 (9)
确定向量序列如果收敛到向量,则有
即是方程组(7)的解。 实验的内容与步骤 一、方程求根
1、考察用迭代函数f(x)=2sin(x)求解方程g(x)=2sin(x)-x=0的解的情况。
(1)在同一直角坐标系中,画出y=f(x)及y=x的图像。从图上观察,方程x=2sin(x)有几个解?
在计算机中打开Mathematica4.0系统;
点击鼠标进入工作区后,输入以下语句
按Shift和Enter键运行。
运行结果:
3、素数的判别
(1)对n=2,3,… ,100,观察2n-1被n整除所得的余数;再取其他的整数m(如3,4,5),观察mn-1被n整除的情况:观察当n为素数时的结果。
m=2时
输入程序:
运行结果:
从运行结果可以发现:当n为素数时,2^n-1被n整除所得的余数都是1。
m=3时
输入程序:
运行结果:
m=4时
输入程序:
运行结果:
(2)对n=2,3,..,300,判断哪些Mersenne数Mn=2n-1是素数?
输入程序:
运行结果:
二、素数的分布
1、素数沿数轴的分布
将素数从小到大排序.用表示相邻的素数间的间隔.计算并画图
n=1000时画出(pn,dn)的图:
t:=
Table[{Prime[i],Prime[i+1]-Prime[i]},{i,1,1000}]
ListPlot[t,PlotStyle?RGBColor[1,0,0]]
n=10000时画出(pn,dn)的图:
t:=
Table[{Prime[i],Prime[i+1]-Prime[i]},{i,1,10000}]
ListPlot[t,PlotStyle?RGBColor[1,0,0]]
2、在二维坐标平面上标出点列观察趋于无穷的趋势,并和函数,以及比较.类似观察,以及, ,和
对任意整数计算, , 看看哪个更接近
文档评论(0)