零误差计算

研究采用有误差的数值计算来获得无误差的准确值具有重要的理论价值和应用价值.这种通过近似的数值方法获得准确结果的计算被称为零误差计算.本文首先指出,只有一致离散集合中的数才能够开展零误差计算,即有非零隔离界的数集,这也是"数"可以进行零误差计算的一个充要条件.以此为基本出发点,本文分析代数数零误差计算的最低理论精度,该精度对应于恢复近似代数数的准确值时必要的误差控制条件,但由于所采用恢复算法的局限性,这一理论精度往往不能保证成功恢复出代数数的准确值.为此,本文给出采用PSLQ (partial-sum-LQ-decomposition)算法进行代数数零误差计算所需的精度控制条件,与基于LLL (Lenstra-Lenstra-Lovász)算法相比,该精度控制条件关于代数数次数的依赖程度由二次降为拟线性,从而可降低相应算法的复杂度.最后探讨零误差计算未来的发展趋势.     

阳树脂再生度对水汽系统氢电导率监测准确性的影响研究

水汽系统氢电导率监测用氢交换柱附加误差超标是造成氢电导率测量不准确的主要原因，而氢交换柱附加误差超标主要源于阳树脂的再生度低。理论分析和研究结果表明，氢交换柱附加误差与阳树脂再生度及水样Cl^-浓度均有关系，再生度越低，Cl^-浓度越高，氢交换柱附加误差就越大。对于超临界水汽系统氢电导率测量用氢交换柱，为了保证附加误差满足不大于5%的标准要求，建议控制树脂再生度在85%以上。阳树脂再生度取决于再生工艺，动态再生和电再生的效果要比静态浸泡再生好，可以保证树脂的再生度和氢交换柱的附加误差满足要求。     

Hamilton系统下基于相位误差的精细辛算法

Hamilton系统是一类重要的动力系统，辛算法（如生成函数法、SRK法、SPRK法、多步法等）是针对Hamilton系统所设计的具有保持相空间辛结构不变或保Hamilton函数不变的算法.但是，时域上，同阶的辛算法与Runge-Kutta法具有相同的数值精度，即辛算法在计算过程中也存在相位误差，导致时域上解的数值精度不高.经过长时间计算后，计算结果在时域上也会变得“面目全非”.为了提高辛算法在时域上解的精度，将精细算法引入到辛差分格式中，提出了基于相位误差的精细辛算法（HPD-symplectic method），这种算法满足辛格式的要求，因此在离散过程中具有保Hamilton系统辛结构的优良特性.同时，由于精细化时间步长，极大地减小了辛算法的相位误差，大幅度提高了时域上解的数值精度，几乎可以达到计算机的精度，误差为O（10-13）.对于高低混频系统和刚性系统，常规的辛算法很难在较大的步长下同时实现对高低频精确仿真，精细辛算法通过精细计算时间步长，在大步长情况下，没有额外增加计算量，实现了高低混频的精确仿真.数值结果验证了此方法的有效性和可靠性.     

急性髓细胞白血病基因筛选模型的贝叶斯分析

基因表达数据蕴含着大量的生物信息,在生物基因信息研究中,筛选表达水平发生显著变化的差异基因是认识疾病形成机理和辅助靶点药物研究的关键问题.根据急性髓细胞白血病(AML)的基因表达数据,构造基因均值差序列,建立贝叶斯分层混合模型,并为模型的参数赋予具有基因生物特征的先验信息.采用马尔可夫链蒙特卡洛(MCMC)算法对模型参数进行估计,并筛选出急性髓细胞白血病差异表达基因.在实际数据分析中,从美国生物信息中心(NCBI)的高通量基因表达数据库中获取急性髓细胞白血病基因数据集,从经过非特异滤波预处理的14688个急性髓细胞白血病基因中筛选出711个差异表达基因,差异表达基因数仅占急性髓细胞白血病基因总数的4.84%,这一结果与基因差异表达的生物学原理相吻合.     

基于模型的小域估计方法及其拓展——由单变量推广到多变量情形

小域估计是当今抽样调查领域的研究热点,通过建立小域估计模型,将所有调查区域通过公共的参数连接起来,利用相关区域数据作为辅助信息,解决多层次抽样调查中某些域样本量不足的问题,提高参数的估计精度。本文首先在目标变量是单变量的情形下,研究了域层次模型和单元层次模型的参数估计方法,进而推广到目标变量是多变量的情形,考虑变量之间的相关性,研究了多变量Fay-Herriot模型参数的经验最优线性无偏预测及其均方误差矩阵的性质,最后通过一个实例,验证了参数估计具有良好效果。     

基于转台误差隔离的捷联惯导系统尺寸效应参数标定

针对捷联惯导系统尺寸效应参数标定问题,提出了一种基于转台误差隔离的改进标定方法。先标定出捷联惯导系统中加速度计间的相对位置关系,再基于尺寸效应误差最小原则,对载体坐标系原点位置进行优化,得出尺寸效应参数。克服了传统分立式标定方法的一些不足,从原理上降低了标定对转台的依赖。对于陀螺零偏稳定性优于0.005°/h,加速度计零偏稳定性优于2×10~(-5) g的捷联惯导系统:尺寸效应参数标定重复性小于0.01 mm(转速60°/s);在摇摆条件下(幅度30°,频率0.15 Hz),4000 s纯惯性导航平均定位误差由补偿前的3100 m以上减小到补偿后的150 m以内。     

核电结构土-结相互作用分析分区混合计算方法

土-结构相互作用分析是核电结构进行抗震设计和安全评估的重要环节.在核电结构的土-结相互作用分析中,阻尼和非线性是影响结构反应的重要因素. 若采用频域分析,可以方便考虑阻尼,但需通过等效线性化来考虑非线性,不适合于强震作用下的土体非线性.若采用时域分析的逐步积分方法,适合于考虑非线性,但材料阻尼一般采用瑞利阻尼模型,除了紧靠指定阻尼比的少数几个振型外,其他振型的反应将受到瑞利阻尼模型所确定的大阻尼所抑制,造成地震反应与真实情形有较大差异.若采用时域分析的模态叠加法,可合理计入阻尼效应,但模态叠加法不能考虑非线性.因此,如何合理考虑阻尼和非线性是核电结构土-结相互作用分析需要关注的问题.基于此,本文提出一种模态叠加和时步积分结合的土-结相互作用分区算法.其中,出于安全性考虑,地震作用下核电主体结构一般不允许进入非线性,因此结构可采用模态叠加方法,以便合理考虑结构阻尼;土体和基础采用显式时步积分法,可考虑土体非线性;通过人工边界条件考虑无限域的影响 (辐射阻尼).通过简单算例对该方法进行了验证,并用于CAP1400核电结构的土-结相互作用分析中,对比分析了采用模态阻尼和瑞利阻尼时核电结构和场地反应的差异,结果表明结构阻尼模型对场地的反应影响不大,但对结构反应影响明显,在实际工程中应合理选取阻尼模型.      

CL-20/TNT共晶炸药的分子动力学研究

运用分子动力学(MD)方法,选择凝聚态分子势能优化力场(COMPASS),对六硝基六氮杂异伍兹烷(ε-CL-20)、2,4,6-三硝基甲苯(TNT)晶体及其等摩尔比的CL-20/TNT混合炸药和共晶炸药进行不同温度下恒定粒子数等压等温(NPT)系综模拟研究.结果表明,CL-20/TNT共晶的内聚能密度(CED)和结合能随温度的升高逐渐减小;共晶的CED比混合炸药的大,结合能是混合炸药的2倍多,预示其稳定性明显增强.对相关函数和局部放大结构显示共晶中组分分子间作用主要来自TNT中H和CL-20中O以及CL-20中H和TNT中O之间形成的氢键.通过波动法求得的弹性力学性能结果表明,CL-20/TNT共晶的拉伸模量(E)、体积模量(K)和剪切模量(G)介于ε-CL-20和TNT晶体之间,且随温度的升高而下降,符合一般预期;但共晶炸药的柯西压(C12-C44,Cij弹性系数)、K/G和泊松比(ν)均比其组分炸药ε-CL-20和TNT高得多,预示该共晶具有异常高的延展性和弹性伸长,主要是二组分呈层状交替排列且之间存在较强相互作用所致.     

足绑式行人导航偏航角误差自观测算法（英文）

基于低成本MEMS惯性传感器的足绑式惯性导航系统(INS)和零速修正(ZUPT)算法广泛应用于行人导航中。由于MEMS惯性传感器零漂误差较大,零速修正时偏航角误差的可观测性差,INS偏航角误差不能被有效约束,成为行人导航的主要误差源。行人徒步行走特别是在室内楼道等结构化道路上行走时,行走的轨迹大多情况下可认为是近似直线,基于这个事实,提出了一种减小偏航角误差的算法,称为偏航角误差自观测(YESO)算法。当判定行人以近似直线徒步行走时,由行走轨迹计算出的航向角可近似认为是一个常值,那么由于各种误差引起该航向角发生变化时,可以将该变化量作为足绑式INS偏航角误差的观测量,进一步可利用卡尔曼滤波器估计出偏航角误差,对足绑式INS的偏航角进行校准。在室内楼道进行了约350 m的现场实验,实验验证了YESO算法的有效性。实验结果表明,当分别采用ZUPT和ZUPT+YSEO算法进行导航解算时,航向角误差从-29°减小到-2°,南北向最大位置误差从-35.5 m减小到-5.2 m。YESO算法的实现仅依靠系统自身的信息,没有增加额外的传感器,算法具有很好的工程实用价值并能方便地推广应用于车辆导航等领域。