实验2非线性方程f[x]0解法.docVIP

《数值计算方法》实验报告  PAGE \* MERGEFORMAT 20 《数值计算方法》实验报告  PAGE \* MERGEFORMAT 21 Solution of Nonlinear Equations f(x)=0 1.实验描述 1．P40.1:参照程序2.1求解出单调收敛的不动点。 2．P49.1:已知初值，时间和末值，求解汇率。 3. P54.3：已知函数求函数在区间内的极值和根。 4. P69.1：已知运动方程求解运动时间和距离。 2.实验内容 1．使用程序2.1求解下面每个函数的不动点（尽肯能多）近似值，答案精确到小数点后12位。 同时，构造每个函数的图和直线y=x来显示所有不动点。 (a) (b) (c) (d) 2．如果在240个月内每月付款300美元，求解满足全部年金A为500000美元的汇率I的近似值（精确到小数点后10位）。 3.一个计算机程序使用点（x0,y0），（x1,y1），（xn,yn）可画出函数y=f(x)的图形，通常还标记出图形的纵向高度，而且必须写出一个子程序来确定函数f在区间[a,b]内的最大值和最小值。 （a）构造一个算法寻找值Ymax=maxk{yk}和Ymin=mink{yk}. (b)写一个MATLAB程序寻找函数f（x）在区间[a,b]内根的近似值位置和极值。 （c）使用（b）中的程序寻找第1题和第2题中根的位置和极值，并用真实值进行比较。 4.设投射方的运动程为 求当撞击地面时经过的时间，精确到小数点后10位。 求水平飞行行程，精确到小数点后10位。 3.实验结果及分析 1. P40.1: 算法： （1）输入函数g，p0,tol,max1,令k=2。 （2）判断kmax1是否成立，如成立输出结果，如不成立，执行（3）。 （3）令p(k)=g(p(k-1)),err=|p(k)-p(k-1)|。 (4)判断errtol是否成立，如成立，输出结果，如不成立，令k=k+1,跳转到（2）。 start Input g,p0,tol;k=2 kmax1 Y N Let p(k)=g(p(k-1));err=|p(k)-p(k-1)| k=k+1 errto1 output Y N output end 对（a）令，简单画出与的图像，可令x=0,,. 图1. 与 运行程序后，当输入 [k,p,err,P]=fixpt(@g,-1,1e-12,100) 输出k=3 p=2 err=0 P =-1 2 2 通过图形我们知道不动点有3个，X=2 是吸引不动点，其余2个为排斥不动点。 对（b）令，简单画出与的图像，可令x=1,,. 图2.与 运行程序后，当输入[k,p,err]=fixpt(@g,1,1e-12,100) 输??：k =37 p =0.768169156736961 err =5.226929999935237e-013 通过图形我们知道不动点有1个，X=0.768169156736961 是吸引不动点。因而函数只有一个吸引不动点 对（c）令，简单画出与的图像，可令x=1,, 图3.与 运行程序后，当输入[k,p,err]=fixpt(@g,0,1e-12,100) 输出k=38 p=0.768169156736961 err=5.226929999935237e-013 运行程序后，当输入[k,p,err]=fixpt(@g,1.5,1e-12,100) 输出k=38 p=0.768169156736633 err=5.176969963827105e-013 通过图形我们知道不动点有2个，X=0.768169156736961,X=0.768169156736633是吸引不动点。因而函数有两个吸引迭代点。 令（4）令，简单画出与的图像，可令x=1,,。 图4. 与 当输入[k,p,err,P]=fixpt(@g,1,1e-12,100) 输出k =2 p =1 err =0 P = 1 1 通过图形我们知道不动点有2个，X=1是吸引不动点.因而函数只有一个吸引迭代点。 2.p49.1 算法： 输入f,a,b,delta,max1,令k= 1。 f(a)=feval(f,a);f(b)=feval(f,b);k=1. （3）判断f(a)*f(b)0是否成立，若不成立输出结果，若成立执行步骤（4）。 （4）判断kmax1是否成立,若成立输出结果，若不成立执行步骤（5）。 （5）令dx=yb*(b-a)/(yb-ya),c=b-x,yc=feval(f,c). (6) 判断yc=0是否成立，若成立输出结果，若不成立