线性代数机算的速度和精度6.pptVIP

线性代数机算的速度和精度 机算中精度和速度的重要性 在用笔算时，通常都用整数矩阵来演示和做题，几乎不涉及误差，也就没有精度问题。至于矩阵计算速度的低下，即使在二、三阶的低价矩阵中，就已暴露无遗，但那从来不是纯粹进行理论探讨的数学家所关心的内容，甚至很怕读者触及这个敏感问题，这是用笔算解矩阵方程的根本弱点。承认这个弱点必然导致计算机的引入和课程的改造。 线性代数进入应用领域，必定要使用计算机，速度和精度两个现实问题就不可避免地摆在我们面前。在这里，我们将注重精度问题，并只着重于MATLAB命令中包含的、或会对计算精度作出提示的问题做一介绍。 1．双精度浮点数的精度 按照IEEE标准，表达一个数需要8个字节，也就是64个二进制位。双精度浮点数η用下式表示 其中M是一个小于一大于1/2的二进制分数，称为尾数，占用52个二进制位表示； 而指数E是一个带符号的二进制整数。占11个二进制位，总共可表示2048个整数，即可以表示从-1023到+1024的数集，它决定了数的动态范围。数的正负号反映在S上，它只占一位。总计1+52+11=64位。 MATLAB中数的精度 浮点数的量化步长可以代表它的相对精度。它是由M的位数决定的。52位二进制数的量化步长是2 -52=2.2204×10 –16。该数的动态范围取决于指数部分。因为2 –1023≈10 –307及21024≈10