数值分析实验报告.doc

重 庆 交 通 大 学 学 生 实 验 报 告 实验课程名称 数值分析 开课实验室 数学实验室 实验内容 P238 第二题 ，P240 第二题 ，P238 第七题 学院 数学与统计2013年级信息与计算科学专业班 2班 学 生 姓 名 高腾 学 号 631322020213 开 课 时 间 2015 至 2016 学年 第 1 学期 一、P238 第二题 为求方程在附近的一个根，设将方程改写成下列等价形式，并建立相应的迭代公式。 ，迭代公式； ，迭代公式； ，迭代公式。 （一）程序： 第一小问： zyoye1.m function y = zuoye1() %UNTITLED1 Summary of this function goes here % Detailed explanation goes here x = 1.5; while 1 == 1 y = 1 + 1/x^2; if abs(y-x) 10^(-5) x = y; else break end end 第二小问： zuoye2.m function y = zuoye2() %UNTITLED1 Summary of this function goes here % Detailed explanation goes here x = 1.5; while 1 == 1 y = (1 + x^2)^(1/3); if abs(y-x) 10^(-5) x = y; else break end end 第三小问： zuoye3.m function y = zuoye2() %UNTITLED1 Summary of this function goes here % Detailed explanation goes here x = 1.5; i = 0; while i 10 i = i + 1; y = 1 / (x-1)^(1/2) if abs(y-x) 10^(-5) x = y; else break end end 结果分析与检验： 第一小问：x = 1.4656 第二小问：x = 1.4656 第三小问：迭代前十次结果为y = 1.4142，y = 1.5538，y = 1.3438，y = 1.7055，y = 1.1906，y = 2.2907，y = 0.8802，y = 0.0000 - 2.8892i，y = 0.3317 + 0.4659i，y = 0.3321 - 1.0570i 由上述实验结果可知：第三个迭代公式是发散的，无法逐渐逼近跟。 二、P240 第二题 多项式求根是一个病态问题，考虑多项式 求解扰动方程。 产生系数。 取用MATLAB求根函数计算扰动方程的根，分析对根的影响。 程序： zuoye.m function [ output_args ] = zuoye( e ) %ZUOYE Summary of this function goes here % Detailed explanation goes here syms x p = 1; for i = 1:10 p = p * (x-i); end expand(p) t = p + e*x^9; sort(solve(t)) end 结果分析与检验 当时： 当时： 当时： 由上述实验结果可知： ，，，，，，，，， 当越来越小时，求得的根越逼近方程的根。 二、P238 第七题 用下列方法求在附近的根。根的准确值，要求计算结果准确到四位有效数字。 用牛顿法； 用弦截法，取； 程序： nd.m function y = nd( x ) %ND Summary of this function goes here % Detailed explanation goes here y = x^3 - 3*x - 1; end ndd.m function y = ndd( x ) %NDD Summary of this function goes here % Detailed explanation goes here y = 3*x^2 - 3; end 第一小问： zuoye_nd.m function y = zuoy