简介:本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。
简介:Wewilldefineandcharacterizeε-pseudoChebyshevandε-quasiChebyshevsubspacesofBanachspaces.WewillprovethataclosedsubspaceWisε-pseudoChebyshevifandonlyifWisε-quasiChebyshev.
简介:使用Chebyshev-Gauss(CG)伪谱法研究带动量轮和推力器的欠驱动航天器姿态最优控制问题.基于欧拉姿态角和动量矩定理导出两类航天器姿态运动模型,采用Clenshaw-Curtis积分近似得到性能指标函数中的积分项,应用重心拉格朗日插值逼近状态变量和控制变量,将连续最优控制问题离散为具有代数约束的非线性规划(NLP)问题,通过序列二次规划(SQP)算法求解.数值仿真结果表明,对两类欠驱动航天器的姿态机动最优控制均能达到设计控制要求,得到的姿态最优曲线与验证得到的曲线几乎完全重叠.
简介:摘要目的利用体模评价能谱CT去金属伪影(MAR)技术减低脊柱不同植入物伪影的效果。方法选取离体的羊椎体骨模型1具,分别构建椎弓根钉、3D椎体及钛网内固定模型,行能谱CT检查。以10 keV为间隔,重建出70~140 keV之间共8组MAR和非MAR图像。测量置入物旁ROI内的CT值和SD值,计算伪影指数(AI)和信噪比(SNR),同时对骨质、椎弓根钉伪影、3D椎体伪影及钛网伪影严重程度进行主观评分。绘制AI和SNR随keV变化的曲线,采用配对样本t检验比较低能级(70 keV)、中等能级(100 keV)和高能级(130 keV)条件下MAR和非MAR图像AI值和SNR的差异,主观评分的差异比较采用Wilcoxon符号秩和检验。结果在70~140 keV范围内,MAR和非MAR图像椎弓根钉的前方、后方、侧方和3D椎体以及钛网旁AI值随keV的增加而减低,SNR值改善。在70、100、130 keV条件下,椎弓根钉前方、后方、侧方及钛网旁的AI值在MAR组非较MAR组均明显减低(P<0.05)。3D椎体后方,仅在70 keV条件下,MAR图像AI值较非MAR组减低,差异有统计学意义(P<0.001)。在椎弓根钉前方、后方及3D椎体侧方SNR值随keV增加而增加,在其他部位SNR值随keV增加而减低。主观评价中,MAR图像评分均高于非MAR图像(P<0.05)。结论能谱CT MAR技术能够减低脊柱植入物伪影,可明显减少椎弓根钉及钛网伪影,而对3D椎体后方的减少伪影的效果有限。
简介:目的探讨能谱CT金属伪影消除(MARs)技术与普通CT骨算法重建在消减脊柱内固定术后伪影的价值.方法回顾性分析2013年2月—2014年2月解放军第一一七医院22例行脊柱金属内固定置入患者的影像学资料.其中男14例、女8例,年龄30~77岁.22例患者术后不同时间分别行能谱CTMARs技术检查(MARs组)、骨算法成像检查(骨算法组).测量并计算各组感兴趣区的伪影指数(AI),采用骨窗图像质量主观评估标准对骨窗图像质量进行评分,比较MARs组和骨算法组患者骨窗图像质量评分情况,以及不同单能量值MARs组和骨算法组间AI的分布差异.结果MARs组中,随着单能量值的减低,AI值逐渐增高.MARs组各单能量值和骨算法组间比较,AI值差异有统计学意义(F=8.696,P〈0.05);进一步组间LSD-t检验显示,在60~140keV条件下,MARs组AI值均小于骨算法组,差异均有统计学意义(P值均〈0.05).骨窗图像质量评分比较,骨算法组(2.81&#177;0.39)分,优于MARs组(2.00&#177;0.53)分,差异有统计学意义(t=4.856,P〈0.01).结论能谱CTMARs技术和骨算法成像对消减脊柱金属内固定伪影各有优点,应当结合临床需要选择合适的CT成像技术.
简介:StabilizedorChebyshevexplicitmethodshavebeenwidelyusedinthepasttosolvestiffordinarydifferentialequations.MakinguseofspecialpropertiesofChebyshev-likepolynomials,thesemethodshavefavorablestabilitypropertiescomparedtostandardexplicitmethodswhileremainingexplicit.Anewclassofsuchmethods,calledROCK,introducedin[Numer.Math.,90,1-18,2001]hasrecentlybeenextendedtostiffstochasticdifferentialequationsunderthenameS-ROCK[C.R.Acad.Sci.Paris,345(10),2007andCommun.Math.Sci,6(4),2008].InthispaperwediscusstheextensionoftheS-ROCKmethodstosystemswithdiscretenoiseandproposeanewclassofmethodsforsuchproblems,theT-ROCKmethods.Onemotivationforsuchmethodsisthesimulationofmulti-scaleorstiffchemicalkineticsystemsandsuchsystemsarethefocusofthispaper,butournewmethodscouldpotentiallybeinterestingforotherstiffsystemswithdiscretenoise.TwoversionsoftheT-ROCKmethodsarediscussedandtheirstabilitybehaviorisanalyzedonatestproblem.ComparedtotheT-leapingmethod,asignificantspeed-upcanbeachievedforsomestiffkineticsystems.Thebehavioroftheproposedmethodsaretestedonseveralnumericalexperiments.
简介:
简介:Weraiseandpartlyanswerthequestion:whetherthereexistsaMarkovsystemwithrespecttowhichthezerosoftheChebyshevpolynomialsaredense,butthemaximumlengthofazerofreeintervalofthenthChebyshevpolynomialdoesnottendstozero.Wealsodrawtheconclu-tionthataMarkovsystem,underanadditionalassumption,isdenseifandonlyifthemaxi-mumlengthofazerofreeintervalofthenthassociatedChebyshevpolynomialtendstozero.