求解第一类积分方程的正则化—小波方法及其数值试验

1 方法的描述 第一类(Fredholm)积分方程是指形如 (1.1)的积分方程,其中核k(x,y)和右端函数f(x)给定,u(x)是未知函数.许多物理、化学、力学和工程应用问题都能导致第一类积分方程.求解第一类积分方程的一个本质性困难是方程的不适定性,即解的存在性、唯一性和稳定性遭到破坏.常用的数值方法有奇异值分解(SVD)方法、Tikhonov正则化方法、投影方法、正则化-样条方法、再生核方法等.本文提出一种新的正则化-小波方法,在第一类积分方程有多个解时,可以求出具有最小范数的数值解;如果原积分方程有唯一解,则所得的数值解收敛于准确解.数值试验表明,该方法是可行的. 我们在L~2[a,b]中考虑第一类(Fredholm)积分方程,即假设方程(1.1)中积分算子K∈L~2([a,b]×[a,b])及右端f(x)∈L~2[a,b]给定.为保证数值求解算法的稳定性,我们先用正则化方法处理该方程,将不适定问题化为泛函极值问题来求解,然后利用多重正交样条小波基构造求解格式.由于我们给出了直接计算低阶的多重正交样条小波基函数的一般公式,使得解法可以在计算机迅速实现.     

争鸣

问题问题109已知函数f(x)满足:f(x y) f(x-y)=2f(x)·f(y),且f(0)≠f(π2)=0,求f(π)及f(2π)的值.解法1令x=y=0,得f(0)=1.令x=y=π2,得f(π)=-1.令x=y=π,得f(2π)=1.解法2令x=y=0,得f(0)=1.令x=32π,y=π2,得f(2π)=-f(π).再令x=y=π,得f(2π) 1=2f2(π),∴2f2(π) f(π)-1=0.∴f(π)=12或f(π)=-1,从而f(2π)=-12或f(2π)=1.问题出在哪里?问题110人教版高一数学(上)P8,有下面一段话:容易知道,对于集体A,B,C,如果A B,B C,那么A C.事实上,设x是集合A的任意一个元素,因为A B,所以x∈B,又因为B C,所以x∈C,从而A C.这个证明严格吗?…     

关于求解第二类Volterra积分方程的Picard迭代收敛性的一点注记

考虑第二类Volterra积分方程: φ(x)+integral from n=0 to x(K(x,y)φ(y)dy)=f(x),x∈[0,L],(1)其中f(x)∈C([0,L]),核函数 K(x,y)对y可积,且     

西北工业大学·高等数学(B卷)

一、填空题 (本题共 5小题 ,每小题 4分 ,满分 2 0分 )1 .设 z =e- ( yx xy) ,则 dz| ( 1,2 ) =2 .由曲面 z =4-12 (x2 y2 )与平面 z =2所围成的立体的体积等于3.设Σ是平面 x y z =6被圆柱 x2 y2 =1所载下的部分取上侧 ,则 Σzdxdy =4.设 f (x)是以 2π为周期的周期函数 ,在区间 (-π,π]上有 f (x) =1 -x,　 -π 相似文献   

一类正线性算子的渐近展开(续)

<正> 设 f(x)∈C_(2π),U_n(f,x)=1/π(?),其中(?)设 m 是正整数,以 X_m 记 f(x)∈C_(2n),且 f~((2m))(x)∈C_(2π)的 f(x)的全体,并设(?)     

二重福里哀级数的线性求和法

§1.引言设C(Ω)为定义在上的,对x,y均以2π为周期的连续函数空间。对于任意的f(x,y)∈c(Q),借助数组确定如下的二元三角多项式:此处,a_(kl),b_(kl),c_(kl),d_(kl)为f(x,y)的福里哀系数。现在,考虑下面问题: 当A={λ_(kl)~(mn)}满足何种条件时,对于任意的f(x,y)∈C(Ω),关系式     

如何探讨抽象函数的性质

在近几年的高考试卷中出现过不少有关抽象函数的题目,要求研究抽象函数的定义域和值域、反函数、奇偶性、单调性、周期性等,下面逐一加以例析.一、定义域这类问题一般是给出y=f(x)和g(x)的定义域,求解复合函数y=f(g(x))的定义域.解决的关键是将g(x)看成一个整体,来替代y=f(x)中的x,从而转化为求解不等式.例1函数y=f(x)的定义域为[-12,21],求函数y=f(cosx)的定义域.分析与简解:因为函数y=f[g(x)]中的g(x)相当于f(x)中的自变量x.所以?21≤cosx≤12,解三角不等式得kπ 3π≤x≤kπ 2π3(k∈Z).解题的关键是始终要明白定义域是自变量的取值范围…     

On the Best Degree of Approximation by Euler (E,q)-Means

Let f(x)∈L_(2π) and its Fourier series by f(x)～α_0/2+sum from n=1 to ∞(α_ncosnx+b_nsinx)≡sum from n=0 to ∞(A_n(x)). Denote by S_n (f,x) its partial sums and by E_n~q(f,x) its Euler (E, q)-means, i. e. E_n~q(f,x)=1/(1+q)~π sum from m=0 to n((?)q~(n-m)S_m(f,x)), with q≥0 (E_n~0≡S_n). In [1] Holland and Sahney proved the following theorem. THEOREM A Ifω(f,t) is the modulus of continuity of f∈C_(2π), then the degree of approximation of f by the (E,q)-means of f is givens by##特殊公式未编改     

新题征展(36)

A　题组新编1 . ( 1 )函数 y =π - x2 -x2 -π(　　 ) ;( 2 )设 e为自然对数的底 ,则函数y = eπ - x2| 4- x| - 4(　　 ) ;( 3)函数 y =12 sin(πx) .( 1ax - 1 12 ) 3 3(　　 ) ;( 4 ) f ( x)不是常函数 ,且 f( x)满足f ( 8 x) =f( 8- x) ,f ( x 2 ) =f( x - 2 ) ,则 f ( x) (　　 ) .( A)是奇函数 ,不是偶函数( B)是偶函数 ,不是奇函数( C)是奇函数 ,也是偶函数( D)既不是奇函数 ,也不是偶函数2 .( 1 ) f( x)为奇函数是 f ( 0 ) =0的(　　 ) ;( 2 ) sinθ <0是θ在第三或第四象限的(　　 ) ;( 3) p为假或 q为假是 p为真且 q为真…     

新二元三角插值多项式的构造

设Ω=[-πxπ,-πyπ],C(Ω)表示关于x,y均以2π为周期的连续函数空间.若f(x,y)∈C(Ω),取结点组为(xk,yl)=(2k+2n 1)π,(2l 2+m 1)πk=0,1,2,…,2n,l=0,1,2,…,2m,则我们获得一个二元三角插值多项式Cn,m(f;x,y)=M1N∑k=2n0∑l=2m0f(xk,yl).1+2∑nα=1cosα(x-xk)+2∑mβ=1cosβ(y-yl)+4∑nα=1∑mβ=1cosα(x-xk)cosβ(y-yl)其中M=2m+1,N=2n+1.为改进其收敛性,本文构造一个新的因子ρα,β,使得带有该因子ρα,β的二元三角插值多项式Ln,m(f;x,y)可以在全平面上一致地收敛到每个连续的f(x,y),且具有最佳逼近阶.     

Finite sections for singular integral operators by weighted Chebyshev polynomials

A finite section method for the approximate solution of singular integral equations with piecewise continuous coefficients on intervals is considered. The problem is transformed in such a way that results which were previously obtained for singular integral equations on the unit circle using localization methods in Banach algebras are applicable to it. Thus, necessary and sufficient conditions for the stability of the approximation method can be proved.     

Application of the Fourier Transform to the Solution of Singular Integral Equations

A method for solving a certain system of singular integral equations with constant coefficients is proposed. It is based on a procedure for reducing singular equations to equations with continuous difference kernel; the solution of the latter is constructed by using the classical Fourier transform in the class of absolutely integrable functions. Explicit expressions for the solution of the singular integral equations under consideration are obtained.     

非线性积分微分方程组奇摄动边值问题

本文讨论含积分算子的非线性微分方程组Ｒｏｂｉｎ边值问题的奇摄动,在适当假设条件下通过对角化技巧,利用逐步逼近法证明了解的存在,并得到直到Ｏ（ε＾Ｎ＋１）的按范数界限的一致有效估计。     

Approximate solution of a class of singular integral equations of second kind

A simple method based on polynomial approximation of a function is employed to obtain approximate solution of a class of singular integral equations of the second kind. For a hypersingular integral equation of the second kind, this method avoids the complex function-theoretic method and produces the known exact solution to Prandtl's integral equation as a special case. For a particular singular integro-differential equation of the second kind, this also produces an approximate solution which compares favourably with numerical results obtained by various Galerkin methods. The convergence of the method for both the equations is also established.     

THE COLLOCATION METHODS FOR SINGULAR INTEGRAL EQUATIONS WITH CAUCHY KERNELS

1 Introduction Singular integral equations (SIEs) with Cauchy kernels Of the formoften arise in mathematical models of physical phenomena. Since closed-form solutions to SIEsare generally not available, much att.ntion has been focused on numerical methods of solution.In the past twenty years, various collocation methods for SIEs have been the topic of a greatmany of papers, most of which can be found in two surveys[213]. The early works in the fieldis to study tile numerical solutions for…     

高阶奇异积分的小波逼近及数值计算

该文所讨论的是在Ｈａｄａｍａｒｄ主值意义下，高阶奇异积分（犛犳）（狋）＝∫ 犳（狓）（狓－狋）狀＋１ｄ狓，　　狀≥１的小波逼近及数值计算．特别是当小波函数未知时，借助于方程（３．１），对高阶奇异积分作数值计算，建立了收敛性定理．     

On parameters influencing the stability of quadrature methods for singular integral equations with conjugation

In this paper, we study an approximation method for solving singular integral equations with conjugation on an open arc. The stability of the method depends on the invertibility of certain operators which belong to well-known algebras. We investigate properties of these operators and show how to choose the parameters of the approximation method so that the Fredholm indices of the operators mentioned become equal to zero.     

Volterra积分方程的蒙特卡罗数值求解方法

本文讨论了第一类、第二类以及具有奇异核的Volterra积分方程的数值解问题.利用重要抽样蒙特卡罗随机模拟方法获得积分方程解的近似计算结果.通过对文献中算例的实现表明文中所提方法扩展了Volterra型积分方程的数值求解方法,     

关于亚音速绕流问题的强奇性基本解法

本文对求解亚音速流的偶极子基本解法作了新的处理,导出了一个关于求解偶极子强度的强奇性积分方程,给出了强奇性积分有效主值的定义及计算公式.由此可以导出多种整体连续分布的数值基本解法.适用于亚音速气动力计算.     

Semi-analytic treatment of the three-dimensional Poisson equation via a Galerkin BIE method

A systematic treatment of the three-dimensional Poisson equation via singular and hypersingular boundary integral equation techniques is investigated in the context of a Galerkin approximation. Developed to conveniently deal with domain integrals without a volume-fitted mesh, the proposed method initially converts domain integrals featuring the Newton potential and its gradient into equivalent surface integrals. Then, the resulting boundary integrals are evaluated by means of well-established cubature methods. In this transformation, weakly-singular domain integrals, defined over simply- or multiply-connected domains with Lipschitz boundaries, are rigorously converted into weakly-singular surface integrals. Combined with the semi-analytic integration approach developed for potential problems to accurately calculate singular and hypersingular Galerkin surface integrals, this technique can be employed to effectively deal with mixed boundary-value problems without the need to partition the underlying domain into volume cells. Sample problems are included to validate the proposed approach.