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

本文对一类非均匀二维插值样条的研究,得出了矩形域R上双三次插值样条函数x(u,w),当满足条件时,在矩形域R边界上节点处的三阶和四阶混合偏导数的估计式,推广了均匀分划情形的结果,获得了精确的误差估计。     

端点条件对三次样条函数一阶、二阶导数误差估计的影响

本文的主要结果是对文[2,3]的改进,并获得了精确的误差估计。内容包括:设S(x)是三次样条函数,它适合则S(x)在节点处的一阶、二阶导数有估计:(1)(2)其中不等式都是精确的。     

三元二次样条函数及其计算（Ⅰ）——高维数据拟合

半节点二次样条是Ｃ１类连续的，其二阶导函数是阶梯函数，在半节点处产生跳跃。鉴于此，本文利用最小二乘法，获得了一种半节点二次插值样条边界条件确定方法，该方法可以保证二次插值样条在半节点处二阶导数的变化最小，这相当于保证了曲率的变化最小。为了适应四维数据插值的需要，给出了三元二次插值样条的定义及基样条表示，提出了一种边界条件，证明了其存在唯一性。最后把二次插值样条边界条件确定方法推广到三元二次插值样条上去，获得了三元二次插值样条边界条件的确定方法。该方法的意义是，可以直接从插值条件获得边界条件，从而克服了多元样条边界条件难以确定的困难     

三元二次样条函数及其计算（I）：高维数据拟合

半节点二次样条是Ｃ＾１类连续的，其二阶导函数是阶梯函数，在半节点处产生跳跃。鉴于此，本文利用最小二乘法，获者了一种半节点二次插值样条边界条件确定方法，该方法可以保证二次插值样条在半节点处二阶导数的变化最小，这相当于保证了曲率的变化最小。为了适应四维数据插值的需要，给出了三元二次插值样条的定义及其样条表示，提出了一种边界条件，证明了其存在唯一性。最后把二次插值样条边界条件确定方法推广到三元二次插值样     

矩形域上一种非张量积型样条曲面插值

为降低造型所用样条曲面的次数,提出了一种非张量积形式的组合箱样条曲面。其基本思想为利用由4片三角箱样条曲面片组合构成的矩形片作为模块构造曲面。相较于双三次B样条曲面,其总次数为4次且支撑域更小,而在边界处具有和双三次B样条类似的导矢性质。进一步利用其显式表达式求出角点和控制顶点的关系,得到构造插值曲面的方法。最后,文末给出了计算实例,表明了方法的可行性。     

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

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

尖前缘任意翼型的超音速和高超音速俯仰安定性导数

本文推广了文[1]的欧拉方程摄动法,讨论了尖前缘任意翼型俯仰安定性导数的计算方法,并应用该方法计算了双圆弧翼型的俯仰安定性导数,其数值结果与文[2]的超音速线化位流方程有限差分法的数值结果比较,更接近于实验值。在激波附体的条件下,本文的方法可用于计算大迎角尖前缘任意翼型的超音速和高超音速俯仰安定性导数。     

一类可控制的保凸三次H-样条曲线和曲面的构造法

本文发展的用于曲线与曲面设计的三次H-样条曲线的方法,是文献[3]的推广。它简单易行,计算有效。 本文内容包括: (1)构造了一种可控制的保凸三次H三样条曲线和曲面。其方程由显格式表示: H_i(t)=1/(1 λ)[t~3,t~2,t,1](?) (2)讨论了它的几何特性。 (3)证明了保凸性的充要条件。 (4)得出了最佳的误差估计。     

参数三次H-样条在曲线拟合中的应用

本文系文[1]的续篇,主要方面内容包括:严格推导了根据曲线型位点求H-样条的特征多边形顶点的公式;解决了H-样条曲线拟合端点条件的处理方法;对H-样条曲线拟合的保凸条件进行了探讨;通过实践,归纳了H-样条曲线拟合的工程应用特点。最后指出:H-样条曲线由于方法简单,计算量小,适用于微机CAD/CAM系统中的曲线设计和拟合。     

一类有理参数三次曲线的形状分析(英文)

一条有理参数三次 H-样条曲线是由一组控制顶点和两顶点连线上的百分比参数所确定。移动一个顶点仅影响三段曲线。有理 H-样条具有许多类似于 B-样条曲线的性质 ,也有 B-样条不具有的性质。本文是在文 [1 ]基础上的继续和发展 ,主要对有理 H-样条曲线的形状进行分析 ,讨论其诸如拐点和奇点的几何特征 ,给出有理参数平面三次 H-样条曲线在非退化情况下有拐点的充要条件 ,并证明在区间 ( 0 ,1 )内曲线段无奇点的结论。为了便于对参数曲线段的形状控制和几何特征的进一步认识 ,在许多的实际应用中 ,需要分析参数曲线段上有无多余拐点和奇点 ,如果有就要消除它。故本文的研究结果无论对理论或实际应用都非常重要。     

基于PIV技术的压力场重构算法实现与研究

介绍了有限容积法、直接积分法和Poisson方程法3种基于PIV瞬时速度场重构压力场的基本原理以及相应的计算方法，选取管流突扩流场和偏置方块绕流流场两个不可压缩流场的瞬时速度场数据，采用上述3种压力场重构算法，分别研究了图像噪声、速度场精度、插值算法以及边界条件的类型与精度对重构压力场的影响。最后针对管流突扩过程第20ms的流场，给出了3种重构算法下的压力场云图以及对应的CFD模拟结果。研究表明，有限容积法和直接积分法容易受到噪声的影响而产生剧烈震荡，但是可以在较大的速度场误差范围内保持较高的精度，通过采用双线性插值可以获得更高精度的重构压力场；Poisson方程法不易受到噪声的影响而产生震荡，同时在高精度PIV速度场下的优势较为突出，通过采用双三次差值可以获得更高精度的重构压力场；混合边界条件仅仅测定边界上有限个点的压力值，就获得了接近狄利克雷边界条件下重构压力场的精度，远高于诺依曼边界条件；边界条件的误差严重降低重构压力场的精度，其影响程度比速度场误差还要大。     

基于超分辨率重构方法的湍流流场重构

从低分辨率流场数据中获取精细流场信息具有重要的研究意义。基于卷积神经网络的超分辨率重构方法是近年来发展的一种较为有效的精细流场重构方法。本文采用高效亚像素卷积神经网络（Efficient Sub-Pixel Convolutional Neural Network,ESPCN）,对Rayleigh–Bénard（RB）对流的数值模拟数据和湍流边界层（Turbulent Boundary Layer,TBL）的实验测量数据进行了超分辨率重构,并与双三次插值方法（Bicubic Interpolation）的重构结果进行对比。对比结果表明:在较小的下采样比下,ESPCN方法和Bicubic方法的重构精度相当;在较大的下采样比下,ESPCN方法的重构精度明显优于Bicubic方法。此外,ESPCN方法对数据梯度较大区域的超分辨率重构效果优于Bicubic方法。     

数字地形数据的二维三次卷积插值

由于数字地形数据库是以网格的形式给出网格点上的地形高度信息，在低空突防轨迹优化过程中，往往需要非网格蹼地形的高度信息和地形曲面的导数，这些都是离散的地形数据库所无法提供的。对二维的地形数据进行插值是解决问题的有效手段。本文提出一种二维离散数据的插值方法－－二维三次卷积插值。该插值函数是地形曲面的三次逼近，具有连续的三阶导数，可以非常好地满足低空突防轨迹优化的需要。仿真结果表明，它具有比双线性插值和     

顺拉法的最大拉直力

按照相对地球的惯性坐标系，根据动能定理，首先推导得出了回收物为“有限质量”情况下的拉直力计算公式，而后又根据回收物为“无限质量”的假设，得出了与文（１），（２）相同的拉直力计算公式。在推导中，文（３）方程中一项∫＾ＬｍａｘＬｓＰｄＬ修正为∫＾ｘｍａｘⅡ０Ｐｄｘ。     

亚、超音速细长翼身组合体大迎角气动特性的计算研究

本文提出了一种适于初步设计使用、具有良好精度的亚、超音速细长翼身组合体大迎角气动特性的综合性计算方法。对大迎角情况下的涡升力,采用吸力比拟原理计算;位流升力的计算,采用基本解的数值计算方法。关于机翼翼剖面头部圆度和涡破碎对涡升力的影响,进行经验性修正。翼身干扰的贡献,通过翼身干扰系数进行计算。并按文[4]原理,将亚音速计算方法推广到亚音速前缘的超音速情况。对几种机翼与翼身组合体的计算结果表明,本文方法具有方法简便、计算快速和符合设计精度要求的优点。     

关于二元风洞侧壁抽气的讨论

本文研究了二元风洞侧壁抽气问题。文[1]曾考虑了抽气对风洞内无粘主流和有粘边界层的双重作用。本文研究则进一步包含了侧壁开孔率和开孔方向的影响,更为全面地描述了有侧壁抽气时风洞内的流动,并从理论上探讨了消除侧壁干扰的抽气控制问题。     

三维弹塑性边界元1/r~3积分奇异性的处理

本文将文[2]的基本思想与方法,推广应用于处理三维弹塑性边界元法塑性单元中的1/r~3积分奇异性问题,推导了有关的具体公式,并从理论上证明了该方法处理1/r~3奇异性的有效性。针对奇异性单元。本文提出旋转与坐标转换的分割技术,可大大减少工作量。并使计算机编程更简洁。数值算例研究表明,文中采用的方法是十分有效的。