【精品PPT】powerpoint 演示文稿 - 螺旋线与平面的交点.ppt

* 实验3 螺旋线与平面的交点 3.1 实验目的 本次实验目的在于使学生能够借助于计算机,利用数学软件解决 一些具体的非线性方程的求根问题。 3.2 实验内容 螺旋线与平面相交的情况多种多样，根据螺旋线与平面方程的不 同可以相交，也可以不相交。在相交的情况下，可以交于一点，也可 以交于好多点。对于各种相交的情况，要求其交点的坐标并不是一件 容易的事。本次实验就以此为背景讨论几种方程求根方法。以下面的 具体问题为例：已知螺旋线的参数方程为 , , , , 平面的方程为 , 求该 螺旋线与平面的交点。 3.3 问题求解 将采用多种方法求螺旋线与平面的交点坐标，包括二分法、迭代 法和弦截法等。 将螺旋线的参数方程代入平面方程后可得： 等价变形得 ， luoxuanxianyuzhixian.m theta=0:0.01:8*pi; y1=4*(cos(theta)+sin(theta)); y2=2-0.5*theta; plot(theta,y1,theta,y2) 从图形可见在 内直线 与三角曲线有4个交点。 Matlab (1) 区间二分法 . (1) 区间二分法 定理：若函数 在区间[a, b]上连续, 且 , 则 在区间内至少有一个根。 (1) 区间二分法 定理：若函数 在区间[a, b]上连续, 且 , 则 在区间内至少有一个根。 [ ] a b . ] [ [ [ ] ] ] . . . ] [ ] [ ] [ [ 取每次二分小区间的中点可得一个收敛序列 ， 该实数序列的极 就是方程的根。即 。 若要求精度 ， 则要做等分区间的次数至少是： ★ 计算 内的交点 Matlab 实际上： erfengfaqiujie.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 执行结果为：x3 =5.37587451934814 该点的近似坐标为： （2.46346668293075，-3.15140157740808，5.37587451934814） 查看图形 continue back Matlab (2) 牛顿迭代法 设方程 有实数根， 若能够将方程等价地转化成 ， 取一个初始值 代入 的右端算得 ， 依次再计算 ， 类推可得序列 ， ， 称此序列为由迭代函数 产生的迭代序列， 为迭代初始值。 若该迭代序列收敛，则它的极限就是方程 的一个根 。 称为方程根的k次近似值. 称使得迭代法收敛的初始值的取值范围为迭 代收敛域。 图示点击 例3.1 容易验证方程 有两个不相等的实数根。 输入MATLAB软件命令 solve(‘x^2+x-3=0’) 执行结果为： ans = [ -1/2+1/2*13^(1/2)] [ -1/2-1/2*13^(1/2)] 该方程可等价地转化为 ， 由此而决定的迭代格式为 ， continue back ， ，显然该迭代序列不可能收敛到方程的根。 若取迭代初始值 但是，方程也可以等价地转化为 ， 此时若取迭代初始值 ， 用迭代法 ， 则经过34次迭代后就可得到方程的一个正的近似解 ， 则经过22次迭代后就可得到方程的一个负的近似解 若继续取迭代初始值 ， Run-iteration x=-5:0.1:3;y1=x;y2=sqrt(3-x); plot(x,y1,x,y2) hold on x0=1;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*) continue figure back 关于迭代法的收敛性有著名的李普希兹定理 若李普希兹条件成立，则用第k次近似值代替方程的精确解 其误差 , 有下面的估计式： 同时由 可得要满足精度至少要做的迭代次数： 即： ★ 牛顿迭代法 迭代格式：