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

用台劳公式表达函数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页)本文意在补充若干材料,以利于函数之近似计算。     

关于单叶解析函数一类子族的增长、掩盖定理

本文研究了单叶解析函数的一类子族βα族,令f(z)∈βα且z=0是f(z)-z的k+1阶零点,我们得到了f(z)的较为精细的增长、掩盖定理,所得结论推广了某些已知结果。     

Linear Codes with Two Weights or Three Weights from Two Types of Quadratic forms

Let F_qbe afinite field with q=pmelements,where pis an odd prime and mis apositive integer.Here,let D_0={(x_1,x_2)∈F_q~2\{(0,0)}:Tr(x_1~(pk1+1)+x_2~(pk2+1))=c},where c∈F_q,Tr is the trace function fromFF_qtoFpand m/(m,k_1)is odd,m/(m,k_2)is even.Define ap-ary linear code C_D =c(a_1,a_2):(a_1,a_2)∈F_q~2},where c(a_1,a_2)=(Tr(a_1x_1+a_2x_2))_((x1,x2)∈D).At most three-weight distributions of two classes of linear codes are settled.     

无穷区间二阶微分方程三点边值问题的可解性

运用非线性抉择原理讨论了无穷区间上二阶微分方程三点边值问题((1+t2)x′(t))′+f(t,x(t),x′(t))=0,t∈I=[0,+∞)x(0)-βx′(0)=αx(η),x(+∞)=0正解的存在性,其中α≥0,β≥0,0η+∞,f:I×I×R→I是一个L1-Carathēdory函数。     

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

本文主要结果为下述定理。 定理:设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.     

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

本文指出了工程界关于高阶马尔可夫过程的一个错误定义,证明了(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)模型,因此是一个二阶广义马尔可夫序列。     

关于Westergaard应力函数的可加性

当n个加载机构分别作用于相同的带裂纹结构上,并且对于每个加载机构情况下的Westergaard函数Z_(Ik)(z)(k=1,2,…n)都已知,则对应于n个加载机构同时作用的Westergaard函数Z_I(z)=ΣZ_(Ik)(z)。并且总的应力场(包括裂纹尖端的应力强度因子)可以按Z_I(z)来计算。     

平方阻尼在非线性缓冲系统中的作用

本文讨论平方阻尼(以参数a表征)在图(1)所示的缓冲系统中的作用,结果表明:平方阻尼对于降低相对位移峰值d_m总有良好的效果,但不一定降低加速度峰值a_m。在一定速度冲击V下,系统的平方阻尼有两种最佳值,使a_m最小的a_(am)和使缓冲效率γ最高的a_(cm),而平方阻尼最佳条件依赖于系统的缓冲弹簧的刚度变化类型: (1)线性弹簧情况:a_(am)=a_(cm),最佳条件为系统在首程中均匀减速; (2)渐硬弹簧情况:a_(am)=a_(cm),最佳条件为系统在首程中初始和终末加速度相等即; aV~2=ω_0~2q(dm) (3)渐软弹簧情况:a_(cm)>a_(am)。 本文的基础是所谓“首程显峰”,假定即系统在速度冲击V下的响应峰值d_m和a_m出现在第一行程中。当平方阻尼值不太小,缓冲弹簧刚度特性接近反对称时,该假定成立。     

非齐次Poisson过程的若干新性质

本文在一些文献的基础上,进一步讨论非齐次Poisson过程的某些性质,给出了若干新性质。 设{N(t),t≥0}是累积强度为的非齐次Poisson过程。迄今的文献(如[1][2]等)指出,(?)n≥1,前n个到达时刻τ_1,τ_2,…,τ_n的联合概率密度为 本文定理1指出,上式不仅是必要的,而且是充分的,并给出了充分性的证明。从而,得到了描写过程统计规律的一个刻画。然后,在N(t_i)=n_i(i=1,…,k,0相似文献   

常用的振动公式

正弦、位移、速度、频率 如果知道泣侈(D)’l{J峰一峰或涪隔胜,频率川赫竺,(11寸/秒)一」’川: v=汀fD 例:i艾D二一叫·,f二二一0112,则v=3 1 .11!、」/秒。方程l变换,如:则计算0一峰矢徽或峰速伙v(1)也按止个变以,},已知的_几个以D二v/汀fla)f二 V汀D比如D用毫米,t 11。) 如果用国际单位制或米制,也一样。前例中V=0.8米/秒.把结果除l,000,v就是米l秒。 正弦位移、加速度、频率 设D用时,频率用赫芝,现在要计算。一峰矢墩或峰加速度,加速度要)lj重力单位,则计算公式为: A=0 .05llf2D(2) 例:设D=1,f二10Hz,则人=5.泣19。方程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)次加入γ序列,就可怕除此项误差。选出十个计算实例作为说明。     

沿海近地层湍流特征分析

通过对杭州湾南北两侧共四个观测点近地层湍流特征分析发现:①湍流强度随风速增大而减小。低风速时,几乎与风速呈反比,高风速时基本保持为常数。可描述为:U≤U_0时,I_k=mk/U;U>U_0时,I_k=n_K。其中n_k=a_k/ln(z/z_0),mk=a_kU_0/ln(z/z_0)。a_x=1.0,ay=0.8,a_z=0.5。海岸附近,U_0=4m/s,深入内陆1.5千米以后,U_0=3m/s。②u_*/U=k/ψ(z/L)的 Businger 模式在海岸附近与观测结果吻合很好,在内陆,-1.0≤z/L≤0.3时二者吻合较好,在此范围以外,偏离较大。③不稳定条件下,水平风速标准差可用(?)描述,垂直风速标准差可用(?)描述。稳定条件下,可用(?)(?)描述。此外,稳定条件下,有σ_w/σ_u=0.48,σ_w/σ_v=0.59;不稳定条件下有σ_w/σ_u=0.53,σ_w/σ_v=O.67。     

关于共形紧致流形的一个注记

对由一个分裂定理确定的共形紧致流形的结构,给出了一个注记,并且证明:若(M,g)是一个n维共形紧致流形且Ric_M≥-(n-1)和λ_0(M)=n-2,则在H~1(L~2(M))中不存在任何一个k≥2正交调和形式组。     

直径≤d的对称本原有向图的广义本原指数集

给出了含有自环的直径≤d的全体n阶对称本原有向图的第k(1≤k≤n)个广义本原指数的上确界,并证明了这类有向图的第k个广义本原指数集为:当1≤k≤d时,E0d(n,k)=1,2,…,d-1+k,当d+1≤k≤n时,E0(n,k)=1,2,…,2d。     

一类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矩阵特征值反问题有解和有惟一解的充分必要条件及求解这类问题的数值算法。     

可修系统试验数据的比较

非齐次Poisson过程,特别是Weibull过程,在可靠性增长试验和加速寿命试验领域中起着非常重要的作用。本文考虑了多个Weibull过程强度函数是否相等的假设检验问题。在假定所有的形状参数相等的条件下,本文在Ⅰ型多台定数截尾、多台定时截尾和Ⅱ型多台定数截尾三种情况下,给出了检验假设H_0:α_1=α_2=…=a_k的小样本和大样本的统计方法。     

关于Davidson—Lanczos方法的收敛率

本文对文[1]提出的求解大型对称矩阵A的极端(几个最大或最小)特征值及相应特征向量的Davidson—Lanczos方法,用Rayleigh—Ritz逼近理论,研究了该方法的收敛率。证明了由该方法产生的规范正交向量{v_i}_i~m=1是Krylov子空间K_m≡Span(v_1,Av_1,…,A~(m-1)v_1)的一组基。设A的k个最大特征值为又,λ_1>λ_2>…>λ_k,相应的近似特征值为λ_i~(m)(i=1,…,k),得到 这里γ_i(γ_i>1),W_i和W_i~(m)是常数。     

带比例矩阵逆特征值问题

研究讨论了一类带比例矩阵的特征值反问题:任意给定2n-1(n≥2)个实数λ(n)1…λ(2)1λ(1)1λ(2)2…λ(n)n,求一个带比例矩阵A,使得λ(j)1和λ(j)j分别是其顺序主子阵Aj(1!j!n)的最小和最大特征值。文中给出了此问题有唯一解的充要条件以及有解的充分条件,并给出了解的表达式,最后用数值算例验证了结论的正确性。     

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算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

对数正态分布下基于截断数据的统计分析

加速寿命试验是在短时间内获得高可靠产品的失效数据的一种有效办法,因而得到普遍应用。但在低应力下,要在短时间内得到失效数据较为困难,常常进行到一段时间被迫停止,因而得到的是截断数据。本文考虑在有截断数据情形下,对数正态分布未知参数的估计。具体地讲,试验一共有k个加速应力S_1,S_2,…,S_k,有R个产品进行试验,在应力S_i下有n_i个产品。进行加速寿命试验,但是我们得到的不全是寿命失效数据y_j~i,i=1,2,…,k,j=1,2,…,n_i,n_1+n_2+…+n_i=n,而是Z_j~i=min(y_j~i,τ_j~i是预先给定的截断时间,记δ_j~i=I(y_j~i≤τ_j~i),则我们只能观察到(Z_j~i,δ_j~i),i=1,2,…,k,j=1,…,n_i。本文给出了用(Z_j~i,δ_j~i)来估计未知参数的一种方法,证明了这种参数估计的存在性,无偏性及相合性。