Ω上三重积分优化复化Simpson数值算法

设空间区域 Ω={(x,y,z)|α≤x≤b,φ_1(x)≤y≤φ_2(x),φ_1(x,y)≤z≤φ_2(x,y)}。(1)f(x,y,z)在Ω及其邻域内具有四阶连续偏导数,φ_1(x)与φ_2(x)在[α,a]内可导,φ_1(x,y)与φ_2(x,y)在Ω的投影(xoy面)区域上具有连续偏导数。下面介绍三重积分 I=∫∫∫f(x,y,z)dxdydz (2)的优化复化Simpson数值积分算法。首先将Ω进行划分,把[α,b]分为2m等分,步长与分点为 h_1=(b-α) /2m,x_i=α+ih_1(i=0,1,2,…,2m)。 (3)在x_(2i+1)处把[φ_1(x_(2i+1)),φ_2(x_(2i+1))分为2n等分,步长与分点为 g_(1,2i+1) =((φ_2(x_(2i+1)))-(φ_1(x_(2i+1))))/2n (i=o,1,2,…,m-1), (4) y_(2i+1,j)=φ_1(x_(2i+1))+jg_(1,2i+1) (j=0,1,2,…,2n)。     

一类最佳二次样条插值

若已知区间[a,b]的一个分划△:a=x_0相似文献   

布尔函数线性等价的分析与应用

对于g(x) =f(xA +b) +l·x +c ,给定f(x) ,g(x) ,如何求取等价关系A ,b ,l,c是一个有用的问题 .该文利用Walsh谱和自相关函数谱作为工具 ,给出的算法 1可以求取g(x) =f(xA)型的等价关系 .针对g(x) =f(xA +b)+l·x +c类型的等价关系 ,当b已知时 ,基于Fuller Millan算法给出的算法 2比Fuller Millan算法至少要快k - 1倍 ,其中k为函数绝对自相关函数谱含有的谱类个数 .应用于AES的S 盒的 8个布尔函数间等价关系的求取 ,算法 2比Fuller Millan算法提高速度近 2 0倍 .应用于IP (IsomorphismofPolynomials)问题的分析 ,指出Patarin所给参数的IP问题是可解的 ,因此基于IP问题的密码体制是不安全的 .     

拟合实验数据的新方程--非整数幂多项式方程

将实验数据拟合成方程在科学研究和工程计算上具有十分重要的作用。目前，拟合实验数据最常用的是整数幂多项式f(x)=c0 c1x c2x^2… cnx^n但是用此多项式拟合各类数据时有时误差较大。本文提出一个新方程，即双系列非整数幂多项式g(x)=c0 c1x^a…… ckx^ka ck 1x^(k 1)b ……cnx^nb式中a，b为参数ci(i=0,1,2,3,……n)为待定系数。在拟合各类实验数据时，新方程总是优于整数幂多项式。     

奇妙的自生成对称图案

在下面的程序中, 通过很简单的迭代运算就生成了许多奇妙的对称图案,想信你看见了这些“杰作”,一定会惊叹不已！ 程序中的迭代公式为: x=(x+k＊sin(y))cos(2＊Pi/q)+y＊sin(2＊Pi/q);y=-(x+k＊sin(y))sin(2＊Pi/q)+y＊cos(2＊Pi/q); 其中Pi=3.14159;q=3,4,4,…,12。给定一个q值和一个x,y的初始值,就可以得到一个图     

GFT及离散卷积的并行算法及其实现

一、GFT的计算 GFT是离散富里叶变换DFT的一种推广.它在许多方面有实际应用,其定义为: 设a,b为二个实数,x_n(n=0,1,…,N—1)为一实序列,称 X_k=sum from n=0 to N-1 x_nW_N~((n+a)(k+b)),k=0,1…,N-1,为具有时间参数a及频率参数b的广义DFT.简记为GFT(a,b),其中W_N=e~(-i2π/N)。可以证明其逆变换为     

关于满足[I(x,y)=I(x,I(x,y))] D-蕴涵的解

研究了[I(x,y)=I(x,I(x,y))]方程,[I]为由连续三角模[T]、连续三角余模[S]和强否定[n]生成的D-蕴涵,即[I(x,y)=S(T(n(x),n(y)),y)],给出了满足方程[I(x,y)=I(x,I(x,y))]的解。     

魅力无穷的膨胀动画

下面的程序将演示从中心向四周膨胀的一幅美丽的动画过程。程序通过一个迭代公式来实现。在迭代公式中,对坐标的改变通过一个分段函数来实现：当X≥0时,X＝y－;当X＜0时,X＝y－;而改变y则通过y＝a－x实现。其中a,b,C选取不同的参数迭代生成不同图象,其动画过程也有区别。本文通过尝实试,选取了四组参数,可在程序中输入1,2,3,4来来得到这四组参数。程序在运行中,先在屏幕中心生成一个小小的图形,然后向外膨胀,膨胀到一定时候,又开始向内运动,等运动到内部一定程度,又开始向外运动膨胀,并且整个图形膨胀得更大;这样周…     

第二类物质溶液表面张力与浓度关系的数学模型

利用科学计算绘图软件Origin,将"最大气泡法测定乙醇溶液的表面张力"的实验数据,通过多项式:y=A B*x C* X^2和y=A B*x C*x^2 D*x^3;对数式:y=a-b*ln(x c);指数式:y=A_1*exp(-x/t_1) y_0等数学模型进行拟合。对拟合方程分别求导后,代入吉布斯(J.W.Gibbs)吸附等温式。再利用Origin进一步处理,作第二类物质溶液表面吸附等温线(溶液浓度从无限小至无限大)。然后,将所获得的图形参数进行综合对比分析和论证,最终确立第二类物质溶液表面张力σ与浓度c之间的理想方程,建立起相应的数学模型:即σ=b*exp(-c/a) d。     

常微分方程初值问题的并行算法

<正> 本文介绍常微分方程初值问题二个并行算法:一、分段并行积分法.二、块隐合一步法。一、分段拜行积分法 1.算法初值问题: y′=f(x,y) y(x_0)=y_0的分段并行积分法,概要如下:将求解区间[a,b]分成N个等长的于区间[x_i-1,x_i],x_0=α,x_N=b,i=1,2,…,N_0在点x_i给出M_i个初值y_(ij),j=1,2,…M_i,i=1,2,…,N-1。对于每一i,M_i应取多大,M_i个初值应如何选取,这些问题在后面讨论。这样我们得到了1+M_1+M_2+…+M_(N-1)=S个初值问题:     

小型解释型SLBASIC语言开发与应用（六）

八、字符集字符集是ASCII字符集的子集,是构筑SLBASI语言的最基本元素。包括:英文字母(除了字符串常数外,不分大小写):ABCDEFGHIJKLMNOPQRSTUVWXYZ——变量名、字符串常数必须大写a b c d e f g h i j k l m n o p q rs t u v w x y z数字:0 1 2 3 4 5 6 7 8 9(十六进制数字A B C D E F)九、算术运算符SLBASIC语言中有加、减、乘、除四种运算,加号、减号、乘号、除号,分别用" "、"-""*""/"符号表示。一元操作符:-;负号例:-100二元操作符:;例:12*23乘方运算,23可写成2*2*2。     

Periodical Stability for a Class of Uncertain System Under Process Perturbation

Controlled object for a class of m × n flow show with infinite buffers, is a linear state space model in the sense of ordinary max-algebra R , ⊕, ? . S open[1~3]: ( ) ( ) ( )( ) ( ) X k A X k B u kY k C X k??? == ?? ⊕ ? , k = 1,2,L (1) where A ∈ R p × p , B ∈ R p × q ,C ∈ Rq ×p, p = mn, q = m + n. k denotes the processing batch of jobs J 1, J 2, …, J n; X = [ x1 1 , L , x1 n , x2 1 , L , xm 1 , L , xmn]T denotes the state vector; xi j is the time when job J j was processed…     

催化剂反应中一类非线性边值问题的有限元配点法

§1.引 言 Galerkin方法是求解微分方程边值问题应用最广的一类有限元方法.文[1]利用配置点Galerkin方法研究了边值问题Ly=(a(x)y')+c(x)y=f(x),x∈I=(0,1)y(0)=y(1)=0的近似解.本文利用配置点Galerkin方法研究如下催化剂反应中质量转换问题:Ly=xy"(x)+(s-1)y'(x)+xq(x)y=xf(x,y),x∈I, (1)y'(0)=0, -y'(1)=A(y(1)-1) (2)     

On the structure of binary feedforward inverses with delay 2

Let M‘=S(Mα,f)be a semi-input-memory finite automaton with input alphabet Y and output alphabet X.If X=Y={0,1},then M‘ is a feedforware inverse with delay 2 if and only if there exists a cycle c of state diagram of Mαsuch that f(y0,…,yc,λα（t)0 can be expressed in the form of f ^(1)(y0,…,yc-1,λα（t)) yc for any state t in C and y0,y1,…,yc in Y;or of f^(2)(y0,…,yc-2,λα（t)) yc-c for any state t in Cand y0,y1,…,yc in Y;or for any state t in Cand y0,y1,…yc,in Y,y0,y1…yc satisfies the D[t] condition.The socalled y0,y1…yc satisfying the D[t] condition is that:for some i,j,(i,j)∈{(1,2),(1,3),(2,1),(2,2),(3,1),(3,2)},there exists a (c 2-k)-ary function f^(k),k=1,2,3,such that the Equation(1)and Equation (2)hokl simultaneously for all y‘c-2,…,y‘c 1∈Y. Equation (1);f(y0,…,yc-i,y‘c-i 1,…y‘c,λα（t))=f^(j)(y0,…yc-i,λα（t）） y‘c-i 1 Equation (2):f(y1,…,yc-j 1,y‘c-j 2,…,y‘c 1,λα（t））=f^(j)(y1,…,yc-j 1,λα（t）） y‘c-j s where t=δα（t）and if (i,j)=(1,2)then one and only one of the following conditions C1 and C2 holds for all y‘c-1,y‘c,y‘c 1∈Y.Condition C1:there exists a c-ary function g^(1),such that f(y0,…,yc-2,y‘c-1,y‘c,λα（t））=g^(1),(y0,…,yc-2,λα（t）） y‘c-1( )y‘c;Condition C2:there exists a (c-1)-ary functiong g^(2)such that f(y1,…,yc-2,y‘c-1,y‘c,y‘c 1,λα（t））=g^(2)(y1,…,yc-2,λα（t）） y‘c-1 y‘c,where t=δα(t).     

一种改进的用封锁技术消除隐藏线算法的研究

本文探讨了针对Z=f(x,y)型三维曲面图形的投影画法中,一种改进的利用封锁技术消除隐藏线的算法。并对文献[1]中的原始算法和本文的改进算法所产生的图形结果作了比较及讨论。     

非线性系统降价观测器的构造

本文主要基于李雅普诺夫稳定性概念、指数观测器和系统指数可观的概念[1]探讨了如下 式所示系统建立降价观测器的理论及方法.x=Ax+f(x),y=Hx     

IO加速解决系统瓶颈

在目前的企业信息化实际环境中,由于企业网络中数据流量增大,降低了服务器的应用响应能力,其中,服务器的I/O延迟成为制约系统性能的瓶颈。在这样的背景下,英特尔I/O加速技术应运而生,它能够较好地帮助企业解决了这一系统性能问题,不仅加快了数据传输速度,同时也提高了系统效率及可靠性。IOAT解决方案包括4个硬件组件:D e m p s e y或W o o d C r e s t处理器、B l a c k f o r d或G r e e n C r e e k北桥、ESB2E南桥以及GilgalGbE内置网卡芯片或IntelPCI-E网卡。英特尔网络连接产品事业部业务发展经理齐炜表示,在传统的以太网当…     

非解析复映射构造IFS的参数研究

为了采用非解析复映射构造分形或奇怪吸引子,研究了复映射f(z)=e~(iπ/2)z~n+c的广义M集的1周期参数对构造非线性IFS的影响.在该复映射的M集1周期区域随机选取参数;根据M集的对称性,用与所选参数在M集中对称位置的参数构成迭代函数系;在动力平面上构造出迭代函数系中的所有迭代函数的充满Julia集以及它们的公共吸引域;将随机选出参数所构造出的迭代函数的吸引不动点作为初始迭代点,通过在迭代函数系中连续随机选取一个迭代函数,跟踪这个吸引不动点在动力平面上的公共吸引域内的迭代轨道.通过实验,找到了可以生成分形的非线性IFS的参数选取方法.结果表明:当n取不同值时,非解析复映射族f(z)=e~(iπ/2)z~n+c的广义M集的1周期参数可以用于构造非线性IFS,这种IFS可以大量生成分形山以及具有Z_(n+1)和D_(n+1)对称特性的新分形.     

高阶Carotid-Kundalini函数Julia分形集的特性

Julia分形图是研究复动力系统的一种有力工具。文章研究了由高阶Carotid-Kundalini函数f(z)=cos(Nzkarccos(z))+c生成的Julia分形图的性质:当c为实数和N为实数或纯虚数时,分形图具有对称性;当N为实数,c=0时,图形具有2k+1个主瓣,和2k个从主瓣上发出的触角,且触角无界。     

国防科技大学研究生院二○○○年博士生入学考试数值分析试题

一、选择题 (每小题 3分 ,共 1 2分 )在下列各题的备选答案中 ,请把你认为正确的答案的题号填入括号中 ,少选、多选均不给分。1 .下列命题正确的有 (　　　　 )。( 1 )数据 {xi,fi}ni=0 的 m次样条函数 Sm( x)满足条件S( m)m ( xi) =f ( m)i 　　 i =0 ,… ,n　　 ( 2 )利用数表 {xi,fi}ni=0 构造的拉格朗日插值多项式 Ln( x)是一个不超过 n次的多项式。( 3)设 α是 f( x) =0的根。如果 f( x)在区间[α-δ,α+δ]上二次连续可微并且 f′( x)≠ 0 ,则 x0 ∈ [α- δ,α+ δ],牛顿迭代收敛。( 4 )三次样条函数是一个三次多项式。2 .下列求积…