哈尔滨工程大学数值大作业2014-附fortran程序研讨.docxVIP

B班大作业要求： 1. 使用统一封皮； 2. 上交大作业内容包含： 一 摘要 二 数学原理 三 程序设计（必须对输入变量、输出变量进行说明；编程无语言要求，但程序要求通过） 四 结果分析和讨论 五 完成题目的体会与收获 3. 提交大作业的时间：本学期最后一次课，或考前答疑；过期不计入成绩； 4. 提交方式：打印版一份；或手写大作业，但必须使用A4纸。 5. 撰写的程序需打印出来作为附录。 课 程 设 计 课程名称： 设计题目： 学 号： 姓 名： 完成时间： 题目一：非线性方程求根 一 摘要 非线性方程的解析解通常很难给出，因此非线性方程的数值解就尤为重要。本实验通过使用常用的求解方法二分法和Newton法及改进的Newton法处理几个题目，分析并总结不同方法处理问题的优缺点。观察迭代次数，收敛速度及初值选取对迭代的影响。 用Newton法计算下列方程 （1） ， 初值分别为，，； （2） 其三个根分别为。当选择初值时给出结果并分析现象,当，迭代停止。 二 数学原理 对于方程f(x)=0，如果f(x)是线性函数，则它的求根是很容易的。牛顿迭代法实质上是一种线性化方法，其基本思想是将非线性方程f(x)=0逐步归结为某种线性方程来求解。 设已知方程f(x)=0有近似根xk(假定) ，将函数f(x)在点xk进行泰勒展开，有 于是方程f(x)=0可近似的表示为 这是个线性方程，记其根为xk+1，则xk+1的计算公式为 ，k=0,1,2,… 这就是牛顿迭代法或简称牛顿法。 三 程序设计（本程序由Fortran语言编制） （1）对于，按照上述数学原理，编制的程序如下 program newton implicit none real :: x(0:50),fx(0:50),f1x(0:50)!分别为自变量x，函数f(x)和一阶导数f1(x) integer :: k write(*,*) x(0)= read(*,*) x(0) !输入变量：初始值x(0) open(10,file=1.txt) do k=1,50,1 fx(k)=x(k-1)**3-x(k-1)-1 f1x(k)=3*x(k-1)**2-1 x(k)=x(k-1)-fx(k)/f1x(k) !牛顿法 write(*,(I3,1x,f11.6)) k,x(k) !输出变量：迭代次数k及x的值 write(10,(I3,1x,f11.6)) k,x(k) if(abs(x(k)-x(k-1))1e-6) exit !终止迭代条件 end do stop end （2）对于，按照上述数学原理，编制的程序如下 program newton implicit none real :: x(0:50),fx(0:50),f1x(0:50)!分别为自变量x，函数f(x)和一阶导数f1(x) integer :: k write(*,*) x(0)= read(*,*) x(0) !输入变量：初始值x(0) open(10,file=1.txt) do k=1,50,1 fx(k)=x(k-1)**3+94*x(k-1)**2-389*x(k-1)+294 f1x(k)=3*x(k-1)**2+188*x(k-1)-389 x(k)=x(k-1)-fx(k)/f1x(k) !牛顿法 write(*,(I3,1x,f11.6)) k,x(k) !输出变量：迭代次数k及x的值 write(10,(I3,1x,f11.6)) k,x(k) if(abs(x(k)-x(k-1))5e-6) exit !终止迭代条件 end do stop end 四 结果分析和讨论 （1）对于方程，当初值分别为，，时，所得结果如下 k xk 初始值1 初始值2 初始值3 0 x0 1 0.45 0.65 1 x1 1.500000 -3.012102 5.791591 2 x2 1.347826 -2.046517 3.909853 3 x3 1.325200 -1.395849 2.686963 4 x4 1.324718 -0.916236 1.926420