第五部分智能仪器的典型数据处理算法.pptVIP

第五章 智能仪器的典型数据处理算法 5. 2 常用函数算法 智能仪器的数值计算程序包括： 常用数值计算 四则运算 指数函数 对数函数 常用函数计算 三角函数 乘方 开方 等。 应用举例：如计算标准偏差需要平方、开方运算，计算分贝需要求对数，计算坐标需要三角函数等。 一、三角函数的计算 利用三角函数的诱导公式，不难把任意角的三角函数转化为锐角三角函数。 由于tan(x)=sin(x)/cos(x)，cot(x)=1/tan(x)，所以只要解决了sin(x)和cos(x)的精确的数值计算问题，其它的三角函数的计算问题也就迎刃而解了。 在智能仪器中应用，可以利用sin(x)和cos(x)的泰勒展式来求相应的近似值，需要结合机器计算采用相应的技巧。 课堂上只介绍求sin(x)近似值的计算方法，其余部分大家可以查阅相关文献。 1 正弦函数泰勒展式的处理方法 首先考虑正弦函数y=sin(x)在x∈[0,π]的函数值计算问题。利用y=sin(x)在x=0处的泰勒展式 不难得到 对于x∈[0,π]而言，上面两式的收敛性都不成问题。 对于求sin(x)的近似值来说，无论是计算量，还是误差控制，利用后面一个公式计算的优越性更多一些。 2．计算sin(x)近似值的数学形式 记SINTNV(x,n)表示利用sin(x)/x的泰勒展式取前n+1项之和作为sin(x)/x的近似值，亦即 可以立即得到 sin(x)≈x[SINTNV(x,n)] 不难理解，在一定条件下，n值取得愈大，利用上面两式得到的sin(x)的近似值的精度愈高。 所以，接下来的问题是SINTNV(x,n)数值计算方法以及如何根据精度要求确定n值。 3．SINTNV(x,n)数值计算方法 SINTNV(x,n)可以改写为 可以立即得到如下的递推格式： 计算源代码略。 4．根据精度要求确定项数n 对于给定的n值，利用SINTNR(x,n) 作为sin(x)/x的近似值，根据莱布尼兹级数理论，截断误差Rn(x)满足关系式 Rn(x)≤x2n/(2n+1)!≤x2n/(2n)! 为此，记 则对于给定的x∈[0,π]以及ε0，可以通过求解下面的不等式 SINTNR(x,n)ε 来寻找所需要的n。记 SINTRN(x,ε)=min{n|SINTNR(x,n)ε} 不难编写一个小程序来计算SINTRN(x,ε)。 教材中图5-1为一计算正弦函数框图（与该方法略有不同）。 二、对数函数值计算方法 求对数函数值最关键的问题实际上是计算[1,2]内的数的自然对数。 对于计算大于2的数x的对数来说总可以找到整数K,使得2Kx 2K+1，记u=x/2K,则有1u2, ln(x) =ln(2Ku)=Kln(2)+ln(u) 此时还有1u2。 对于计算(0,1)内的数的对数来说，也可以利用 ln(x)=-ln(1/x) 把它化为求大于1的数的对数问题。 求一般的对数来说，可以利用换底公式转化为自然对数的计算问题。 1.对数函数泰勒展式的处理方法 对数函数ln(1+x)在x=0处的泰勒展式 只是在(-1,1)内处处收敛，且而收敛速度太慢，所以对于数值计算来说，实际意义不大。为此，在写出ln(1-x)在x=0处的泰勒展式 它的收敛区域为[-1,1]。 将上面两式相减，再经数学处理即得 而且收敛区间为(-1,1)。 2．简单的数学处理方法 对适当给定的x，记LOGTNV(x,n)表示利用上面的(*)式右边取前n+1项和所得到的结果，亦即 假如要计算某个正数a(不妨假设a1)的对数的近似值，有 而且对任意a1，均有0x1，从而能保证(*)式右边收敛，所以对于适当的n值，有 ln a =2·.x·LOGTNV(x,n) 3．LOGTNV(x,n)的数值计算方法 模仿秦九韶算法的处理方法，可以先把LOGTNV(x,n)表示为 不难得到下面的递推关系 而且有y1 =LOGTNV(x,n)，从而形成了完整的算法。 4．误差估计公式的推导 对任意a1，记x=(a-1)/(a+1)，记 由于x1，所以，如果不计一个常数因子，可以把Rn(x)作为 ln a =2·x·LOGTNV(x,n-1) 的截断误差。 在上面Rn(x)的表达式中，我们可以把分母中的k值替换为较小的n值，并进一步略去常数1，从而得到有限形式的误差估计公式 5. 根据精度确定项数n 对任意a1（实际上可以约定a3）以及eps0,记 那么LOGTRN(x，eps)