矩阵与数值分析实习题.docx

矩阵与数值分析上机实习

设,其精确值为.

编制按从大到小的顺序,计算的通用程序

编制按从小到大的顺序,计算的通用程序

按两种顺序分别计算并指出有效位数〔编制程序时用单精度〕

通过本上机题，你明白了什么

秦九韶算法。n次多项式，用秦九韶算法编写通用的程序计算函数在x0点的值，并计算在23点的值.

(提示：编写程序时，输入系数向量和点x0，输出结果，多项式的次数可以通过向量的长度来判断

分别用Gauss消元法和列主元消去法编程求解方程组Ax=b，其中

，

.

编程求解题3中矩阵A的LU分解及列主元的LU分解〔求出L,U和P〕，并用LU分解的方法求A的逆矩阵及A的行列式

编制程序求解矩阵A的cholesky分解，并用程序求解方程组Ax=b，其中

,

用追赶法编制程序求解方程组Ax=b，其中

，

.

编程求解矩阵A的QR分解；

分别应用Jacobi迭代法和Gauss-Seidel迭代法求解如下方程组

分别应用Newton迭代法和割线法计算〔1〕非线性方程在[1,2]上的一个根；〔2〕在[-4,-3]上的一个根。

采用二分法计算非线性方程，查找区间为[-4,-2]

函数fx=11+x2，在[-5,5]上分别取2，1，

用三次样条插值上题中的插值节点，并画图比拟插值效果。

〔提示：原函数在两个端点-5,5的导数值可作为边界条件〕

fx=x2sinx，分别用复合梯形法和复合Simpson公式计算积分-22f(x)dx，区间分为20，

用2点，3点和5点Gauss积分法分别计算定积分，并与真值作比拟。

S

常微分方程

dudx

分别用Euler法，改良的Euler法，Runge-Kutta法去求解该方程，步长选为0.1,0.05,0.01.画图观察求解效果。

要求：1.考试前提交作业〔word形式提交，包括代码和实验结果〕，主题写“学号+姓名”.

发送至邮箱

张宏伟老师：

孟兆良老师：

程明松老师：

董波老师：

2.可用任何一种语言编程