一类可控制的三次参数样条曲线的拐点、奇点的讨论

本文是在文[3]基础上的继续与发展,内容包括: (1)构造了一类三次参数曲线。 (2)证明了一个定理:在非退化的情形下,由曲线方程所确定的平面三次参数样条曲线段:(ⅰ)既没有奇点又没有多余拐点;(ⅱ)没有拐点的充分必要条件是: 1/2<λv_i/(1-λ)<2。     

对某些奇异积分的逐次分半加速法

对于f(x)∈C~(2p 1)[0,1],尤拉-麦克洛林求和公式可写为:(h_k=1/2~k) (?)(1)对此,在逐次减半加速法中计算序列: T_(mi)=(γ_mT_(m-1,i 1)-T_(m-1),:)/(γ_m-1)(m=1,2,3,……) (2)其中γ_m=4~m,形成γ序列(4,4~2,4~3,……)。 对于f(x)∈C~(2p 1)(0,1],我们建立公式: (?)(3)其中(?) (4)对于f(x)=x~βln~nx(β>-1,n=0,1,2,…)则有(?) (5)(3),(5)右边的前半部是由于f(x)在x=0上的奇异性而引起的误差。如将α=2~(β 1)以(n 1)次加入γ序列,就可怕除此项误差。选出十个计算实例作为说明。     

二维电磁衍射场的静场近似

设通过保角变换: ζ=x+jy=Aζ_1+A_1ζ_1~(-1)+A_2ζ_2~(-2)+……使无限长导体柱的正截面外部变成ζ_1平面上的单位园外部。由二维的Helmholtz公式出发,求得当波长远较柱截面尺寸为大的平面电磁波以垂直于柱轴的方向投射时: (1)E_1平行于柱轴,E_1=exp[jb(ycosα-xsinα)],则远区衍射场 (2)H_1平行于柱轴,H_1=exp[jk(ycosα-xsina)],则远区衍射场 其中:S=截面积, u=cosθ+jsinθ=(x+jy/γ), p=2π(∈μ)(1/2)A(e~(jα)A_1-e~(jα)A),p=P_α+jP_y所相应的矢量P=i_xP_x+i_yP_y就是导体柱在入射波的电场下所感应的等效电矩。 在椭柱(长短半径各为a,b)的情形中: A=(a+b)/2,A_1=(a-b)/2,S=πaba=b就是圆柱的情形;b=0就是薄片的情形,利用Babinet原理,可推得平面上无限长开槽的情形——此二情形都已有准确解,与本文结果相比较,当ka→0时,只差高阶无限小。     

与坐标系无关的双圆弧样条

设通过n 2个节点P_0,P_1……P_n,P_(n 1)的拟合曲线是由2(n 1)段光滑连接的圆弧(P_0,Q_0,P_1,Q_1……,P_n,Q_n,P_(n 1)构成,且在边缘节点P_0,P_(n 1)上给定切向。则由判据 integral from P_0 to P_(n 1) k~2ds=min (k:曲率;s:弧长)可以近似得到与解三次样条相似的线性方程组。由此确定各节点上的切向。选定相邻圆弧的其他切点Q_0,Q_1,……Q_n使(?)(i=0,1,……n),因而决定每段圆弧的半径、圆心。这样一组圆弧,称为“双圆弧样条”,与坐标系无关,而且有效地降低了曲率起伏,使整个拟合曲线的曲率变化得均匀而缓和。     

一类二维样条插值函数的最佳误差估计

本文主要结果为下述定理。 定理:设x(uw)是矩形域上关于该矩形上均匀分割的二维双三次样条插值函数,且x(uw)满足条件(5),则x(uw)在矩形域R边界上的节点处的四阶混合偏导数有估计式: |S_(i,0)|≦|A[i,n—1]||ε_(n,0)| |B[i,n—2]||ε_(0,0)|=[0,-4,(-1)~2 4,…(-1)~i 4]/[0,-4,(-1)~2…(-1)~n 4]|ε_(n,0)| sum from h=i to n-2 (-1)~(k(k-2)-(i 1)(i-2))[0,-4,(-j)~2 4…(-1)~i 4]/[0,-4,(-1)~2 4,…(-1)~(k 1) 4][0,-4,(-1)~2 4,…,(-1)~(k 2)4] (-1)~(i(i 1)/2)/[0,-4,(-1)~2 4,…(-1)~n 4]|ε_(0,0)|其中等号成立的条件分别为: A[i,n—1] B[i,n—2] ε_(n0),ε_(00)>0 A[i,n—1] B[i,n—2] ε_(nm),ε_(0m)>0 其中 i=1,2,…,n—1. j=1,2 …,m—1.     

R~n空间中初值问题的广义拟线性化方法

本文用广义拟线性化方法在Rn 空间中研究了初值问题x′=f(t ,x) ,x(0 ) =x0 的解 ,首先建立了一个比较定理 ,本文的主要结论是定理 2 ,该定理推广了文〔1〕的结果     

R^n空间中初值问题的广义拟线性化方法

本文用广义拟线性化方法在R^n空间中研究了初值问题x′=f(t,x),x(0)=x0的解，首先建立了一个比较定理，本文的主要结论是定理2，该定理推广了文[1]的结果。     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

加力燃烧室内火焰扩张的研究

本文用数值方法初步研究气流参数(速度、温度、余气系数)对管内紊流火焰扩张的影响。所用的计算模型有两个:一个是修改的EBU模型,假设系数C_(EBU)随来流参数而变化。另外一个是混合长度模型,假设当紊流强度很低时混合长度系数取一固定值,而当紊流强度较大时混合长度系数在一定范围内变化(随来流温度和油气比而变)。 在这个假设下计算结果表明:ⅰ)当来流紊流强度较大时,火焰扩张随来流参数而变,这与文献[2]实验数据基本相符;ⅱ)当来流紊流强度较小时,气流参数对火焰扩张基本无影响,这与文献[1]基本相符。     

Weibull分布定时截尾寿命试验参数的近似联合置信域

研究了Weibull分布大样本定时截尾寿命试验，利用一种新的途径给出了试验参数的近似联合置信域。设产品寿命X服从Weibull分布，对n个受试样本x1,x2，…，xn进行定时（时间为t0）截尾试验，得到观察数据Si=min(x1,x0)，i=1,2,…，n。利用似然方程，分别得到参数λ，b的极大似然估计λ和b。进一步考察联合分布（λ，b）可以用二元正态分布来近似。再利用多元正态分布与X^2分布的关系，可以推导出未知参数u=（λ，b）的置信度为1－a的联合置信区域（x-u)′I0（x-y)≤xa^2。     

定数截尾试验下双参数指数寿命分布尺度参数的EB估计

本文讨论了在无替换定效截尾试验方案下,当产品寿命为双参数指数分布时,尺度参数(失效率)久的经验Bayes(简记EB)估计问题及其收敛速度。设在给定λ,μ下,产品寿命T服从双参数指数分布,其概率密度为 受试产品有n个,试验中前r个产品依次出现的失效时间为t_(1)≤t_(2)≤……≤t_(r)。令 则(x,y)为(μ,λ)的充分统计量。记(x,y)的联合边缘密度为f(x,y),若取二次损失函数,则λ的Bayes点估计为 利用密度函数及其偏导数的核估计,构造出λ的EB估计为 φ_(1n)(x,y)与φ_(1m)(x,y)的Bayes风险分别为 在一定的正则性条件下,我们证明了 这表明,λ的EB估计的收敛速度q可任意接近于1/2。     

关于二阶马尔可夫过程的样本序列是否存在无后效性的讨论

本文指出了工程界关于高阶马尔可夫过程的一个错误定义,证明了(p=2)满足这个定义的平稳高斯过程是不存在的。 本文还指出由二阶微分方程 x″(t) a_1x′(t) a_2x(t)=ε(t) (其中ε(t)是白高斯过程)描写的随机过程x(t)的任意均匀采样序列都不能是AR(2)序列,而由下面微分方程 x″(t) a_1x′(t) a_2x(t)=ε′(t) βε(t)描写的随机过程x(t),当β~2>[max(c_1~2,c_2~2)]时(c_1、c_2是特征方程z~2 a_1z a_2=0的根),至少存在一个采样间隔Δ_1,使相应的样本序列是AR(2)模型,因此是一个二阶广义马尔可夫序列。     

台劳公式拉格郎奇型余项中的“θ”与函数近似值计算

用台劳公式表达函数f(x),其拉格郎奇型余项f~((n))(x_0+θh)/n！h~n中的“θ”,除下列二点外,吾人所知甚少,即:(i)O<θ<1(本文中一切θ均如此),(ii)若f(x)在所论区间上更有f~((n+1))(x),又f~((n+1))(x)在x_0连续,且不为零;则当h→0时,θ→1/(n+1)。(参阅Franklin:Treatise On Advanced calculus 138—139页又格列本卡:数学分析教程第一卷第二分册353页)本文意在补充若干材料,以利于函数之近似计算。     

负二项分布可靠性增长的经典分析方法

本文用经典分析方法研究了负二项分布的可靠性增长问題。研制规划由m个阶段组成,每个阶段的失败概率(即不可靠性)满足条件:q_1≥q_2≥…≥q_m,在此条件下,我们得到了qi,i=1,…,m的约束极大似然估计以及最末阶段不可靠性q_m的经典置信上限。我们也讨论了在Bayes分析方法中,如何选取参数的验前分布问题并给出了趋势检验方法。几何分布可靠性增长的相应结果可从负二项分布可靠性增长的结果直接给出。最后用数值例说明了这些方法。     

非稳态四阶椭圆方程的Galerkin有限元法

针对平面有界凸多角形区域上一类非稳态四阶椭圆方程 u t+Δ2 u =f (x,t) (x,t)∈Ω× [0 ,T]u = u n=0 (x,t)∈ Ω× [0 ,T]u(x,0 ) =u0 (x) x∈Ω的初边值问题进行了 Galerkin有限元方法分析。首先给出了所讨论问题的 Galerkin有限元方法的离散格式 ,其次对所讨论问题的解与其离散问题的解之间的误差进行了分析研究 ,最后利用椭圆投影算子的性质 ,获得了一定模意义下的一些误差估计     

不等式Sum(1/a_i)(from i=1 to i=n)≥n~2/Sum(a_i)(from i=1 to i=n)的两个推广

从n个正数的调和平均值 ,不大于其算术平均值出发 ,导出了一个分式不等式 ,又应用排序的方法 ,从两个方面将这个不等式加以推广 ,得到了两个更具一般形式的分式不等式和指数不等式 :∑ni=1biai ≥n∑ni =1 bi∑ni=1 ai,∑ni =1 aim≥(∑ni=1 ai) mnm - 1 (n∈N ,m∈Z)。     

南昌航空工业学院学报(自然科学版) 2003年总目录(1～4期)

基于卡尔曼滤技术的相对导航技术之研究丁兆平 (1- 1)……………………………………………………………………………最佳平方逼近多项式病态改善刘建国　黄　华 (1- 7)………………………………………………………………………………椭圆外切 2n +1边形中切顶线三角形有向面积的定值定理及其应用喻德生 (1- 10 )……………………………………………总量控制与排污权交易的理论分析张秋根　魏立安　何钱昌 (1- 13)……………………………………………………………自动功率因数调整器中的计数式模糊控制方法赵文龙　周　洪等 (1- 18)……………     

关于A.M.Ostrowski的切线—平行线程序的一类收条件

本文给出非线性泛函方程P(x)=0用A. M. Ostrowski的切线——平行线程序: (?)_n=X_n-Γ_nP(Xn),X_n+1=X_n+Γ_nP(Xn)(n=0,1,2…)此处Γ_n=[P'(Xn)]~(+1),可以求解的一类收饮条件。它推广了M. A. [3]和徐利治[5]的结果。     

一类Jacobi矩阵特征值反问题

给定三个互异实数α，β，γ及三个不同的非零定向量x=(x1,x2,……xn)^T,y=(y1,y2,…,yn)^T,z=(z1,z2,…zn)^T，构造n阶Jacobi矩阵J使（α,x),(β,y),(γ,z)是J的第p,q,r个特征对，给出了这一类Jacobi矩阵特征值反问题有解和有惟一解的充分必要条件及求解这类问题的数值算法。     

EXTENDED SELF SIM ILARITY OF PASSIVE SCALAR IN RAYLEIGH-BNARD CONVECTION FLOW BASED ON WAVELET TRANSFORM

INTRODUCTIONTo fully developed turbulence flow,Kol-mogorov ( 1 941 ) define structure function of qorder,i.e.Dq( r) =〈| u( x+ r) - u( x) | q〉∝ rζ( q) ,and foundζ( q) =q3.Similar to velocityfield,for temperature T( x) by definition struc-ture function,i.e.Rq( r) =〈| T( x+ r) - T( x) | q〉∝ rξ ( q) ,Obukhov also gave out a theoreticalscaling rule ofξ( q) =q3.After the 60 's,becauseof the discovery of intermittence in turbulent ex-periment,it is well known thatζ( q) ,ξ( q)≠ q3,…