matlab 画xos函数,振荡积分的数值计算与Matlab实现

第 29卷 第 12期(下) 2013年 12月 赤 峰 学 院 学 报 (自 然 科 学 版 ) Journal of Chifeng University(Natural Science Edition) Vo1．29N0．12 Dec．2013 振荡积分的数值计算与 Matlab实现 陈卫红 (安 徽 财 贸职 业 学院 ，安徽 合 肥 230601) r D r D 摘 要：本文提出利用样条函数计算 f f(x)sintxcLx及 l f(x)costxdx类型的振荡积分，在每个比较小的子区间采用分部积 a J 分法，避 免了整体利用分部积分 需要计算函数在 区间端点处的高阶导数 ．能提 高计算的精确度． 关键词 ：振荡积分 ；样条插值 ；Matlab实现 中图分类号：0246 文献标识码：A 文章编号：1673—260X(2013)12—0004—01 1 背 景 导 引 在信号与信 息处理中，常常要计算如下形式的振荡积分： I=f R(x)f(x)dx，其中R(x)为冲击函数或振荡函数，称为 积分的核．例如：当t一。。时 l f(x)sintxdx与 『f(x)costxdx都 是振荡积分，利用插值计算时，为了使插值函数能很好地逼 近它们，就要求插值多项式的次数足够高，但是，高次插值 有 Rung现象，严重影响积分的计算效果．而分段低次插值光 滑性差，不能满足实际要求．样条函数光滑性好，有很好的磨 光性 ，且消除了 Rung现象 ，能有效地计算振荡积分 ． 2 积分计算原理分析 设 f(X)的样条插值函数为 s(x)，则 r b f b Ii(I)=J f(x)costxdx f s(x)costxdx f h f h I2(1)=f f(x)sintxdx-~-J s(x)sintxdx 首先 ，将 积分 区间[a,bln等 分 ，h=(b—a)／n，结 点 坐标 为 xi=a+ih(i=0，1，⋯，n) 令 I1 (f)=J s(x)costxdx 因此I．(I) ∑I。。(I)：∑f (x)COSt d I = s(x)c。stxdx= 1 s(x)sintx +}s’(x)c。stx 一 1()()sintx 一}s．tl(x)c。s 由于在每一个区间上的s(x)为三次多项式 ，因此 s⋯(x)为 一 常数，且在子区间『xI，x 11上 s-，．(x)： = 令 s”㈦ =T 则有 s’㈥ =s’(a)=f(a)，S’㈥ =s’ )=f(h) 于是有 I。。(I)={-[f(x。+。)sintx。+广f(xOsintxi]+}[s’(x。+．)c。s投i+lSI )c。stxi] 一 i 卜 【costxi+lCOStXi] ：r ∑r。 f(x)costxdx= s(x)costxdx=∑ 故I。(f)=i J I (I) {-[f(h)sintb—fl(a)sintal+} (h)c。stb—f(a)costa] 一 一 1[T ~sintb- 2sinth ， n I 删 利用类似方法得到 I =l f(x)sintxdx一一}[f(b)c0stb—f(a)costa] +} (h)sintb—ft(a)sinta]+t13[T．c。stb—rI 。sta] 一 2sinth 。 n-I 如果积分区问为『0，2盯]，则上述两个公式就变为 I (I) J 0 fl(x)costxdx 一 (2竹)一f(0)] + — nsi nA- 一 ~r ∑n-I 一 Ti) i +—{一 ∑ s nI[2 ± "fill i： 0