2010-6方程求根实验.ppt

螺旋线与平面的交点 (6) 西安交通大学 魏平 实验6 螺旋线与平面的交点 6.1 实验目的 本次实验目的在于使学生能够借助于计算机,利用数学 软件解决一些具体的非线性方程的求根问题。 6.2 实验内容 螺旋线与平面相交的情况多种多样，根据螺旋线与平面 方程的不同可以相交，也可以不相交。在相交的情况下，可 以交于一点，也可以交于好多点。对于各种相交的情况，要 求其交点的坐标并不是一件容易的事。本次实验就以此为背 景讨论几种方程求根方法。以下面的具体问题为例：已知螺 旋线的参数方程为 平面的方程为 求该螺旋线与平面的交点。 等价变形得 wp61.m theta=0:0.01:8*pi; y1=4*(cos(theta)+sin(theta)); y2=2-0.5*theta; plot(theta,y1,theta,y2) 从图形可见在 内直线与三角曲线有4个交点。 6.3 问题求解 将采用多种方法求螺旋线与平面的交点坐标，包括二分法、 迭代法和弦截法等。将螺旋线的参数方程代入平面方程后可得： (1) 区间二分法 定理：若函数 在区间[a, b]上连续, 且 则 在区间内至少有一个根。 [ ] a b . ] [ [ [ ] ] ] . . . ] [ ] [ ] [ [ 取每次二分小区间的中点可得一个收敛序列 该实数序列的极限就是方程的根。即 若要求精度 输入初值： 输出结果：STOP No No No Yes Yes Yes ★ 计算 内的交点 wp62.m a=5;b=6;x1=a;x2=b; while abs(x1-x2)0.000001 x3=(x1+x2)/2; f3=4*cos(x3)+4*sin(x3)+0.5*x3-2; if f3~=0 f1=4*cos(x1)+4*sin(x1)+0.5*x1-2; f2=4*cos(x2)+4*sin(x2)+0.5*x2-2; if f1*f30 x2=x3; else x1=x3; end end end x3 back Matlab (2) 迭代法 设方程 有实数根，若能够将方程等价地转 化成 取一个初始值 代入它的右端算得 依次再计算 类推可得序列 称此序列为由迭代函数 产生的迭代序列， 为 若该迭代序列收敛，则它的极限就是方程 的一个根 。 称为方程根的k次近似值. 称 收敛的初始值的取值范围为迭代收敛域 迭代初始值。 使得迭代法 back 方程可以等价地转化为 此时若取迭代初始值 用迭代法 若继续取迭代初始值 x=-5:0.1:3;y1=x;y2=sqrt(3-x); plot(x,y1,x,y2) hold on x0=0;x1=sqrt(3-x0);s=[];ss=[]; while abs(x1-x0)0.000001 x0=x1;s=[s,x0];x1=sqrt(3-x0);ss=[ss,x1]; end plot(s,ss,’r*’) 输入MATLAB软件命令： solve(‘x^2+x-3=0’) 例1 解方程 back 关于迭代法的收敛性有著名的李普希兹定理 若李普希兹条件成立，则用第k次近似值代替方程的 精确解 同时由 ,其误差有下面的估计式： 迭代法 定理 算法： 迭代格式： 迭代函数： continue 关于牛顿迭代法有下面的收敛性定理 : 定理 对于方程 ，若存在区间 ，使得 (1) (2) 在区间 上连续且符号一致； (3) 对于任意的 则当迭代初始值 ,且 与 保持同号时， 牛顿迭代法产生的迭代序列 收敛于方程在区间 上的唯一解 (3) 牛顿迭代法 ★ 计算 内的交点 可以证明在区间 内函数 选取初始迭代值9.5，编写程序：wp64 x0=9.5; x1=x0-(4*(cos(x0)+sin(x0))+0.5*x0-2)/(4*(cos(x0)-sin(x0))+0.5); while abs(x1-x0)0.000001 x0=x1; x1=x1-(4*(cos(x1)+sin(x1))+0.5*x1-2)/ (4*(cos(x1)-sin(x1))+0.5); end x1 Matlab (4) 弦截法 下面介绍的弦截法则是避开了求导数运算， 具体又分为单点弦截法和双点弦截法。 ★ 单点