震荡函数积分方法论述.docVIP

震荡函数积分 震荡函数作为一大类函数，在工程领域都有重要的应用，但由于这些函数的震荡性，使 得对这类函数在某一区间上的数值积分积分变得异常困难。其主要是积分的误差的控制问 题，或是收敛性的讨论，例如对于一简单的震荡函数/(x) = sin(-),在任一包含原点0附 x 近的区间上的积分，经典的数值积分公式，例如梯形公式，辛普森法则，牛顿一科茨公式等, 对它的收敛性难以保证，或是误差难以控制。这类函数的特点是函数值在积分区间上急剧震 荡，无法找到一个确且的点使得这类函数的函数值在这个值附近的波动都比较小，或是这个 这个点太难找了，或是即就是找到了，最终的误差也难以把握，可能我的会得到一个不收敛 的函数值。 因此，我们有必要对这类函数另行讨论，针对其中的一部分找到较好的方法。 在此我们仅限于对连续的震荡函数的进行讨论其数值积分。 直观上，函数f(x)在某一区间上震荡，就是说，函数没有断点，在这个区间上，函数多 次回到y=b这条直线上，又多次远离y=b这条直线，但不一定函数回到这条直线上与y=b 任意相邻的两交点距离相同，同吋也不一定每次函数远离y二b的这条直线的最远点的纵坐 标都相同，即我们可以看到所谓的震荡函数只是函数值在急剧波动，对函数本身的要求较广, 1 ￡ 如/(x) = -sin(x),这一函数；或是类似于这类fCx) = sin(p), k不等于0。如此看来， 任一连续的周期函数也是震荡函数,II任一震荡函数可以视为振幅和周期不断随位置变化的 周期函数，又由魏尔斯特拉斯第二逼近定理，我们知道任一区间上的周期函数可以由一个三 角多项式函数进行逼近，由此我们在此只讨论三角性或是由其经过有限次的运算构成的函 数，便可以解决一大类问题了。由上述的讨论我们很自然地给出了震荡函数的定义如下； 1.1.1定义：震荡函数是如下形式的一大类函数 fM = [g(sin(m(x)))/z(x)/ g.m.he C[a,b]} 其中我们称函数g(r)为调频函数，加⑴称为调频函数g(/)的内含函数，加兀)为f(x)伸 缩系数，其屮g(/),加兀)均为连续函数，gCv)是一致有界函数，称此类震荡函数为基本 的震荡函数,且以上的函数加劝,m(x)都是在积分区I可上除去至多有限个点n+1阶可导。 1.1.2定义：XOG (c, J)称为g(X)的诡点，当且仅当，满足 lim m(x) = * x—x0 且|g(m(x))| M0,Vxe [c,d] 明显以上定义的震荡函数f(x)的结构较为复杂，一下我们先以一大类较为简单的震荡函数为 例来说明这类函数的求数值积分方法。 1」.3定义：一下函数称为0型震荡函数 fQ(x)三{比 sin(ox + b) + k2]/z(x)/a,b为常数，h(x)e C[c,d]} 即调频函数为线性函数，我们在此扔掉调频函数的o次项，因为这一项为非震荡项，且同吋 丢掉了线性函数的系数，这是因为由积分的线性性得到的。 对上述的0型震荡函数，我们再进一步做如下形式的化简： 由 三 角 函 数 的 和 角 公 式， 扎(x) = sin(ax) cos(Z?)/z(x) + cos(?x) sin(Z?)7z(x) 又因为sin(x)和cos(Q仅仅相差龙/2的平移变换，因此而sin(/?)和cos(b)均为常数，由积 分的线性性上述办⑴由可剔除为./；(%) = sin(^)/?(x),x€ [c,d]的最简单形式。 1.1.4计算兀(兀)的数值积分 TOC \o 1-5 \h \z d d 记：A) = J = j sin(6L￥)/i(x)6k (1) c c 通过多项式插值的积分公式，现在重多项式插值开始假设要求(1)中人的值，可以选取[c,d] 中的结点，％0，兀1，无2，，，，兀「对伸缩系数进行插值 n X — X ? ￡(x) = n (0 in) 7=0 兀 _ X ； z J 这些是基本的插值多项式，在节点上伸缩系数h的次数最多是n次的多项式是： n p(x) = Xh(xi)li(x) i=0 然后我们可以简单的得出以上(1)中的人;的积分值： d n d j sin(?x) p(x)dx = E 力(xjJ sin(cu:)￡(jr)/r 用这种方法我们可以得到一个用于任何伸缩系数h的公式: n zXA/(^) z=0 d 其屮 Az = f sin( ax )lj(x)dx 以上公式的误差为: i=0 d d -j sin(tzx)p{x)dx - j sin(ax)h(x)dx c c 所以可得: d =J sin(or)[p(x) - 7z(x)] j |sin(6?x)[/?(%) - /?(x)]|dr 乂因为：sin(or) 1,以及由插值误差不超过: “ 1 z=o 右上我们可得