数值分析 练习题.doc

已知函数，给出该函数在点0， 0.5， 1， 1.5， 2上的二次最小二乘拟合多项式. 给出公式中的参数，使其代数精度达到最高，并指明代数精度，用它计算积分 . 写出下面线性代数方程组的Jacobi迭代格式， 并说明它是否收敛 用Euler公式求解下面的常微分方程初值问题，取步长: 利用Cholesky分解计算方程组的解 用牛顿插值法计算满足下面条件的三次多项式 . 下面题目为编程题，请勿手算，请只用Matlab语言. 七、(10分) 已知切比雪夫函数定义如下： ， 编程画出到在区间[-1，1]上的图像. 已知有两个Matlab函数文件定义如下： 文件1 function v = f(x) v = x – cos(x); 文件2 function z = g(x) z = x*log(x)-1; 请编写一个二分法求根的程序，给定初始区间，给出求解这两个函数根的恰当的调用方式. 解： 给出一个收敛的迭代方法计算非线性方程的根，要求得到的根的近似值至少有5位有效数字. 利用LU分解求解下面线性代数方程组 . 求函数在区间[-1，1]上的二次最佳平方逼近多项式. 已知函数有函数值表如下 利用三次插值方法近似计算. 确定下面积分公式中的参数，使其代数精度达到最高，并指明代数精度： 下面题目为编程题，请勿手算，请只用Matlab语言. 六、(16分) 已知微分方程初值问题如下： , 利用Matlab中的函数求出该微分方程初值问题的数值解，并画出其函数图像； 利用上面小题的计算结果，求函数在区间[1，3]上的复合梯形积分的积分值. 编写一个Gauss-Seidel迭代方法求解下面线性代数方程组的解， 矩阵的维数由命令行输入： ，， 要求程序设置有最大迭代步数以及控制精度. 计算如下矩阵的LU分解. 已知函数在节点,,上的函数值以及一阶导数值如下表： 0 1 2 0 0 4 0 2 0 求一个5次多项式,满足 已知勒让德递推公式 推导三点高斯勒让德求积公式 并用该公式计算. 考虑线性方程组Ax=b，其中（为实数）,. （1）写出该方程的高斯-赛德尔迭代格式. （2）证明：当时，高斯-赛德尔迭代发散. 用欧拉公式计算如下微分方程在步长时的近似解. 以下三题为Matlab编程题： 六、(12分)已知函数的实测数据组如下表，试用最小二乘法求二次多项式曲线与此数据拟合. 用牛顿法计算方程在附近的根，要求满足精度. 用乘幂法和反幂法计算对称正定矩阵A的条件数，其中和分别是矩阵A的最大和最小特征值. 用高斯消去法或三角分解法求解线性代数方程组 已知函数在节点处的函数值及一阶导数值如下： i 0 1 2 3 1 2 3 5 8 3 4 0 0 求四次插值多项式 求函数在上的线性最佳一致逼近. 考虑线性方程组Ax=b，其中,. 试写出该方程的逐次超松弛迭代格式（松弛因子为0.9），并讨论该格式的敛散性. 用改进欧拉法 计算常微分方程初值问题 的数值解，取步长. 已知数据 1 3 5 7 9 2 2 1 1 2 试用最小二乘法求形如的拟合曲线. 以下两题为Matlab编程题： 七、(10分) 用牛顿法计算方程在在附近的根，要求满足精度. 八、（12分） 用乘幂法计算矩阵 的按模最大特征值与对应的特征向量. 设，，将进行Cholisky分解，其中：下三角矩阵，并由此求解线性方程组． 对线性方程组 设法导出使Jacobi迭代和 Gauss—Seidel迭代收敛的格式，要求分别写出格式，并说明收敛的理由。 确定系数和，使得下列函数 是一个三次样条，且最小。 4、（15分）求拟合下列数据的线性最小二乘多项式解。 0 0.2 0.4 0.6 0.8 1.0 1.2 0.9 1.9 2.8 3.3 4.0 5.7 6.5 确定下列带权求积公式中的系数使其代数精度最高，指明代数精度： 并用上述公式近似计算积分。