线性模型中M检验原假设分布的随机加权逼近

在线性模型中M-方法可以用于线性假设检验, 其中M检验、Wald检验和Rao的计分型检验是最常用的检验准则. 但是在计算这些检验的临界值时都涉及到未知参数的估计. 在本文中我们利用随机加权的方法来逼近这些检验的原假设分布. 结果表明在原假设和局部对立假设之下随机加权统计量的渐近分布与原检验统计量在原假设之下的渐近分布相同. 因此我们不需要对冗余参数进行估计,利用随机加权的方法就可以得到这些检验的临界值. 而且在局部对立假设之下可以实现对功效的计算. 当取不同的误差分布和不同的随机权时, 我们对本文的方法进行了蒙特卡洛模拟. 结果表明用随机加权方法来逼近原假设分布是非常精确的.     

Strong and electromagnetic J/ψ and ψ（2S） decays into pion and kaon pairs

We discuss interference effects important for the form factors extraction in the vicinity of J/ψ and ψ（2S） resonances in combination with resonance parameters determination. The implementation to the Monte Carlo event generator PHOKHARA of the J/ψ and ψ（2S） contributions to the muon, pion and kaon pairs production associated with a photon at next-to-leading order is also described.     

A Novel Mirror for Cold Molecules with a Semi-Gaussian Beam

We propose a novel mirror for cold molecules with a blue-detuned semi-Gaussian beam and study the dynamic reflection process of cold molecules by Monte Carlo simulation. Our study shows that this mirror can realize a specular reflection of cold iodine molecular beam with a temperature of 30mK by a reflectivity of 58.2% when the semi-Gaussian laser power is 1.0 kW. When a semi-Gaussian CO2 laser beam with a power of 5.8 kW is used, the reflectivity of this mirror can reach about 100%.     

高空离轨发动机流场红外辐射特性研究

高超声速飞行器大攻角机动时,其离轨发动机产生的喷流与高速稀薄的大气来流产生强烈干扰,流场情况复杂,流场红外辐射也是天基红外系统探测的标志性事件.本文针对高超声速飞行器发动机喷流与稀薄来流的相互干扰情况,采用数值求解Navier-Stokes方程模拟干扰流场,采用逐线积分法得到气体红外辐射特性,结合反向蒙特卡洛方法计算得...     

Monte Carlo simulation of the property of a scintillation bar in the multi-neutron correlation spectrometer

To perform a kinematically complete measurement of the dissociation reaction for neutron-rich nuclei, a multi-neutron correlation spectrometer is proposed at Peking University. A Monte Carlo simulation code based on GEANT4 is developed for a single scintillation bar which processes not only the energy deposition but also the light propagation in the scintillator and the light collection and conversion to signal at the end of the bar in a realistic way. The simulating method is described in detail in this paper, and the timing and position resolutions and detector efficiency are studied based on the simulation and compared with the experimental results. A new method of crosstalk rejection has been demonstrated to be important for the design of the whole spectrometer.     

氢氘物态方程研究进展北大核心CSCD

针对近二十多年的氢氘物态方程理论研究工作进行了综述分析,结合本课题组改进的自由能模型、直接量子蒙特卡洛和量子分子动力学方法的模拟结果,对多个研究小组采用不同方法获得的氢氘宽区物态方程数据进行了定量评估分析。结果表明:在当前理论框架下,仅基于第一原理数值模拟得到的氢氘物态方程能够描述的热力学相空间有限;多模型集成的H-REOS.3数据库在105K以下温度与数值模拟结果的相对差别较大,且数据稀少,二者均不能满足工程应用需求。建议采用基于半经验模型的宽区域物态方程研究方法,即结合高精度的实验研究、数值模拟和解析模型,构建满足工程应用需求的氢氘宽区实用物态方程。     

基于多光谱漫反射的牛肉品质参数检测方法研究

为了满足生鲜肉品质参数无损检测领域,对轻便式、低成本设备的开发需求,提出一种基于多光谱漫反射技术的生鲜肉品质检测方法。首先根据漫反射近似理论,结合牛肉样品散射系数、吸收系数及折射率等参数,在无线细垂直光束的蒙特卡洛仿真的基础上,对具有一定发散角度LED光源进行了初始化的校正,分别从光源照射位置概率分布、不同角度的照射概率分布、仰角、方向角的概率分布、不同角度光线入射样品时反射引起能量损失及对光子权重的影响,得到在LED光源发散角情况下,不同源探距下的漫反射率与检测深度,确定了光源与检测器之间的最佳距离为15 mm,然后根据此距离,搭建了多光谱漫反射检测平台,检测平台由8组中心波长为470,535,575,610,650,720,780和960 nm的LED光源组成,与所要检测的生鲜牛肉品质参数相对应。同时利用LED光源的发散角,确定了光源到样品表面的垂直距离与每个光源的安装位置,保证光源照射到样品的区域是均匀的。样品的漫射光强经由信号采集与放大电路的处理后传至上位机,并在上位机完成建模与分析。最后为验证该检测系统的性能,以生鲜牛肉新鲜度参数中的颜色(L*,a*,b*)与pH值为指标,利用60个样品进行了试验,分别得到8个光源下的原始光强值与校正后的反射率值,然后将牛肉样品按照3∶1比例分为校正集与预测集,针对原始光强值与反射率值,分别利用多元线性回归(multiple linear regression, MLR),偏最小二乘回归(partial least squares regression, PLSR)与偏最小二乘支持向量机回归(partial least-squares support vector machine, LS-SVM)三种方法,建立各个参数在原始光强与反射率数据两种情况下的预测模型,并得到最佳模型结果。结果表明,利用反射率数据建模结果均好于光强数据结果,其中参数L*,a*,b*的MLR建模结果优于PLSR与LS-SVR,其预测集相关系数分别为0.983 2,0.907 2及0.935 9,预测集误差分别为1.00,2.14及0.67。参数pH值的LS-SVR建模结果优于PLSR与MLR,其预测集相关系数为0.942 0,误差为0.19。最后利用未参与试验的20块牛肉样品对模型进行了验证,颜色L*,a*,b*及pH参数的预测值与实测值的相关系数均大于0.85,结果证明,利用多光谱漫反射技术以及所搭建的多光谱漫反射检测系统对生鲜牛肉品质参数检测是可行的,该方法能够为设计便携式或微型化生鲜牛肉品质的无损检测仪器提供参考与依据。     

基于钙钛矿量子点荧光太阳集光器的蒙特卡洛光子追踪模拟

传统光聚集器热效应明显、结构复杂、成本昂贵。作为替代,荧光太阳集光器具有许多显著优势并能够有效降低太阳能电池的发电成本,因此受到广泛关注。本文通过传统热注入法合成了全无机钙钛矿CsPbBr_3量子点,并在此基础上设计了基于CsPbBr_3量子点的荧光太阳集光原型器件。通过TEM测试和必要的光学表征,证实本文合成的CsPbBr_3量子点具有典型的立方体结构、76.8%的荧光量子产率、512 nm的发光中心波长和22 nm的中心波长半高宽。此外,结合蒙特卡洛智能优化算法,建立了基于CsPbBr_3量子点的荧光太阳集光器的理论计算模型,确定了全无机钙钛矿量子点最优掺杂浓度和最佳平均波长集光效率。仿真结果表明,在量子点掺杂浓度为2.1×10~(-5) mol/L时,最优的集光效率达到5.4%。本文提出的蒙特卡洛光子追踪模拟过程将为未来荧光太阳集光器参数设计提供科学的计算方法。     

基于蒙特卡洛模型的脊髓损伤后神经细胞定位方法

交通事故、高空坠落、重物砸伤等原因,可能导致人体脊髓损伤（spinal cord injury,SCI）。SCI中断了人体神经信号的传输信道,患者大多高位截瘫,丧失肢体功能,大小便失禁,导致终生残疾,生活异常艰难。2016年世界卫生组织统计数据显示,全世界SCI患病率为（258～785人）/百万人口,每年新发病率为（13.8～86人）/百万人口。目前全球约有SCI患者600万人,并且每年新增约60万人。脊髓损伤一直是医学界的重难点研究课题,医生和科研人员经过长期的研究与探索,目前仍未找到有效的治疗方案来修复脊髓损伤后的微环境以及促使损伤神经的再生。该方法主要针对以下两方面问题：(1) 常规的检测仪器如X射线、核磁共振MRI检查、CT等传统仪器只是影像和形态学方面的观察,不能对脊髓损伤部位的神经细胞进行定位和实时的活性分析,致使医生和科研人员对患者的病情不能精确把握,极有可能延误患者病情,甚至危及患者生命。(2) 2013年12月,西安邮电大学量子通信团队提出采用量子中继“搭桥”的方式,即利用量子纠缠传态在脊髓损伤上断点和下断点建立“连接”,来实现人体神经信号中继的方法;并联合西安交通大学医学院进行大鼠和白兔的实验,取得了突破性进展。2015年1月16日,中国科学院再生医学研究团队联合中国武警脑科医院,结合间充质干细胞,使用神经细胞再生支架完成世界首例脊髓损伤“搭桥”手术,并取得了良好的治疗效果。如何保证“搭桥”是建立在活性神经细胞上,而不是已经坏死的神经细胞上？ 如何找到合适的“桥墩”位置,使“搭桥”的距离尽可能小以保证神经信号传输的保真度能够尽可能的高？ 该文提出的方法可以实现对脊髓损伤后的神经细胞进行精确定位。光谱分析越来越广泛的被应用在生物医学领域,对某些疾病的早期发现和无创或微创的高效检测起着至关重要的作用。神经特异性核蛋白可以用来特异性识别神经细胞,由于存活的神经细胞才会产生神经递质,因此,在动物实验的基础上,通过蒙特卡洛模型模拟光子在生物组织中的传输,对平面Oxy建立极坐标转换,得出光在生物组织中的衰减系数矩阵的计算方法,对脊髓损伤前后的相同位置进行近红外光检测,用聚类算法对神经元特异性核蛋白和神经递质的检测数据进行处理,通过matlab仿真,得出脊髓损伤前后衰减系数的二维分布图,确定Oxy平面中脊髓损伤部位体素的X和Y坐标,在Oxy平面的异常部位选取异常点W与Z轴建立Ozw平面,最终确定脊髓损伤部位神经细胞的位置坐标。该方法可以为医生和科研人员在脊髓损伤患者肢体功能重建的研究中提供理论依据。     

三维静电场实验仪

为克服传统方法不能真实反映静电场的三维本质,综合应用了机械、电子、MCU、传感器、步进电机和计算机接口等技术研制了静电场三维实验仪,实现了静电场的三维空间坐标量和电量的测控,并用所研制的静电场三维实验仪对部分特殊电极所形成的电场进行了模拟实验研究.