GPU光线跟踪算法加速结构研究摘要：基于GPU的光线跟踪算法是当前图形学研究的一个热点，也是将来用于广告、电影、游戏等娱乐产业的关键技术。本文论述了如何对基于GPU的光线跟踪算法进行实现，以及利用各种加速结构，加速算法实现，提高算法执行效率，并对各种加速结构的效果进行了比较研究。关键词：GPGPU 光线跟踪 BVH KD-Tree1．引言近年来，CPU无论在运算能力，还是在可编程性上都得到了大幅的提高，GPU已经在需要大量运算的密集运算领域发挥了举足轻重的作用。各种基于CPU的密集运算被移植到GPU上，以利用GPU巨大的运算能力，加速整个算法的运算过程。光线跟踪算法是生成真实感图形的一种非常重要的方法，在电影、游戏、广告等产业，获得广泛的应用，而光线跟踪算法也是典型的密集运算算法，利用原始的基于CPU的光线跟踪渲染一幅图片是非常耗时的操作。因此，如果能够将CPU上的光线跟踪算法，映射到CPU上，加速光线跟踪算法的执行时间，将会带来巨大的经济效益。因此，基于CPU的光线跟踪算法已成为国内外科研人员的研究热点。2．基于GPU的光线跟踪 相关工作当前，主要由两种方法利用CPU来加速光线跟踪算法。第一种是Carr等人提出来的，将CPU转换为一个蛮力的执行光线一三角形求交的计算器，而将任何的光线生成以及着色过程在CPU上完成。这就需要CPU依然执行绝大部分的渲染工作。C arr等人指出，在ATI Radeon 8500上，每秒最快能够执行1亿2千万次的光线一三角形求交。同时，作者也指出，由于GPU的单精度浮点的限制，图片上依然存在一些不太真实的地方。第二种方法由Purcell等人提出的，改种方法将整个光线跟踪器都移植到CPU上进行实现。从光线的产生，加速结构的遍历，到最后的着色过程都在GPU上执行。此后，有很多相同的项目都是基于Purcell的模型上进行的。 GPU上的光线跟踪算法的映射方式将传统的CPU上执行的光线跟踪算法，映射成为一个GPU协助的，或者基于GPU的光线跟踪器有众多方法。下面重点介绍Purcell提出的映射模型，以及在本文的实现中提出的一个基于CPU的Whitted模型的光线跟踪器。该光线跟踪器的布局如图所示：在Purcell的论文中，它将光线一三角形求交，以及遍历过程分离成两个独立的遍历内核和求交内核。本文的实现中，也按照上述模型图，将光线跟踪算法分解成光线生成，光线一三角形求交，着色这三个步骤。在对光线进行跟踪之前，需要生成从视点指向屏幕的原始光线( primary ray)。在一个GPU上，能够使用光栅器的插值的能力，在一个单一的内核调用中，产生所有的原始光线。给定观察矩形（被采样用于产生图片的投影平面的一部分）的四个角，以及视点，首先计算出这个视锥体的四条边线。如果让光栅器在这4条光线之间，按照512×512规格，在这四条光线之间按照方向进行插值，最终就可以获得能够产生一幅512×512图片（一个像素一个采样点）的所有原始光线的方向。同时能够将这些方向存储在一个纹理里，并把它作为求交内核的输入。所有的原始光线具有相同的起始点，但是仍然把它存储在一个同方向纹理具有相同维度的纹理内。因为当生成阴影光线或者反射光线的时候，光线的原点会发生改变。求交内核把光线的原点，方向，以及场景的描述作为输入数据。在内核被调用数次之后，我们对于每一个像素输出一个击中记录。如果一条光线击中了场景中的某个三角形，返回击中点的3个重心坐标，以及相关的被击中的三角形。此外，还将输出被发现的交点沿光线的距离，以及被击中三角形的材质。这就需要使用5个浮点数值组成一个击中记录。纹理只能够支持4个颜色通道( RCBA)，所以，如果能把击中记录裁减到4个值，那么将是非常有益的。观察发现，只需要3个重心坐标的两个，因为在三角形内部，它们相加的和总是1。这就使得在一个单独的RGBA纹理中存储交点记录是可行的，并且它的维度同其它两个光线纹理的维度相同。Moller和Trumbore提出了一个高效的光线一三角形求交算法，使用这个算法，并利用CPU在向量计算上的优势来进行求交计算。下面列出了求交的代码，这个代码也展示了如何利用向量指令来提高效率。当所有的原始光线都已经计算出了相交的状态的时候，就能够查询着色过程所需要的表面法线和材质的信息。每一个击中记录都存储了一个指向材质纹理的索引，这个材质纹理包含了三角形的法线，材质颜色以及类型。三个顶点的法线根据击中记录的中心坐标进行了插值。最终的颜色能够按(N-L)C进行计算，此处Ⅳ是法线，L是光源的方向，G是三角形的颜色。现在根据击中的三角形所具有的材质的类型（漫反射材质，或者镜面反射材质），需要产生二次光线，以此来计算阴影和反射。1)如果一条光线射出场景之外，像素就被赋予全局的背景颜色。2)如果一条光线击中了一个漫反射材质表面，就发射一条阴影射线( shdow ray)。这些光线的起始点在击中点，方向为从击中点指向光源。3)如果一条光线击中了一个镜面反射材质表面。就发射一条镜面反射光线。镜面发射光线的起始点也在击中点，但是它的方向是在击中点处关于入射光线和插值后的法线对称的方向。一个真正的Whitted类型的光线跟踪器也支持透明材质，从而能够产生折射光线。但由于主要是研究加速结构，所以在本文的实现中，没有考虑折射光线。4)如果阴影光线击中了某个几何体，这就说明在光源和击中点之间，存在某个几何体，所以这个像素就应该是黑色（处于阴影中）。当跟踪阴影光线的时候，不关心最近的那个击中点，更加关心的是是否存在这样的击中点。因此，当有一个交点被发现，就可以停止整个求交过程，从而加速算法的处理过程。在本文的实现中，以相同的方式跟踪阴影光线和反射光线，因此，就没有使用到这个优化策略。已经对每一个像素产生了正确二次光线，如果需要，就能够执行另外一趟遍历/求交过程，对上述的二次光线进行跟踪。每一次调用着色程序就能够对每一个像素返回一个颜色值和一条新的光线。着色内核也可以将前一次着色程序的输出当作本次着色程序的输入。这就使得能够在跟踪连续的光线的时候合并这些连续的镜面反射的颜色。同Carr等人的程序不同，本文所采用的程序不存在浮点精度太低的问题，因为Ceforce 7300在整个管线中支持真正的32位浮点操作。3．加速结构的实现和比较均匀栅格均匀栅格是第一个在GPU上实现的加速结构。Purcell给出了很多选择均匀栅格作为加速结构的理由，但是Purcell没有详细的说明为什么均匀网格对于硬件实现而言比其它的加速结构要更加的简单。当在探讨了均匀栅格的一些主要特性的时候，更加清晰的知道了均匀栅格为什么会成为一个好的GPU机速结构。首先，只用使用简单的算术运算，就能够对于每个体素的遍历在常量时间能被定位和存取。这就消除了对树的遍历的需要，以及重复的纹理查找工作，而纹理查找是相当耗时的。其次，体素的遍历是通过递增算术运算来完成的。这就消除了对堆栈的需要，使得我们能够从光线的起始点开始，以距离递增的顺序访问体素成为可能。再其次，由于对于体素的访问是沿着光线，以距离递增的方式遍历的，所以，一旦在一个被访问的体素中报道发现有一个交点，就可以停止这条光线对体素的遍历过程，从而提高整个遍历过程的速度。最后，用于遍历的代码非常适合用向量编写，而向量形式的编码风格又非常适合GPU的指令集。然而，均匀栅格的缺点就是由于它是空间细分结构的一种特殊情况，多个体素可能包含相同三角形的多个引用。由于无法使用mailbox技术，这就意味着需要对于相同的光线和三角形之间进行不止一次的相交测试。 KD-tree最近，Havran等人对基于CPU的光线跟踪算法的加速结构进行了比较，得出的结论是对于众多不同类型的测试场景，平均而言，KD-tree是最快的。所以，有必要考察一下对于基于KD-tree的GPU光线跟踪算法，是否也会有相似的结论。就像均匀栅格一样，KD-tree也是一种空间细分结构。同均匀网格不同的是，KD-tree利用一个二叉树将场景表示成一个层次结构。在二叉树中，我们将内部节点和叶子节点区分开。叶子节点用来表示体素和与之相关的保存在该体素内的三角形的引用。一个内部节点用来表示空间区域的某个部分。所以，内部节点包含一个分裂面的两个子树的引用，而叶子节点只包含一个三角形列表。KD-tree的创建过程从上而下，根据一个评价函数，通过放置一个分离平面，递归的将场景分离成两个体素。我们能够以递归的方式遍历KD-tree，但是由于GPU没有堆栈结构，所以无法应用递归的策略。取而代之的是，我们能够通过记住我们沿着光线前进了多远来向上或者向下遍历树。这种策略消除了需要堆栈的限制，使得用CPU来完成对KD-tree结构的遍历成为可能。当使用GPU对KD-tree进行遍历的时候，KD-tree像均匀栅格那样被表示成一个纹理的集合。这就意味着有一个保存树数据的纹理，一个保存三角形列表的纹理，和一个保存实际的三角形数据的纹理。GPU的遍历首先调用一个初始化内核，然后按照需要，多次调用合并后的遍历和求交内核。 包围体层次(BVH)给定一些随机的光线，通过计算遍历包围体层次的平均花费，就可以测量出该包围体层次的质量。迄今为止，还没有构建最优的包围体层次的算法，也就是说，如何准确的测量一个包围体层次的平均遍历时间还不是很明显。Goldsmith和Salmon提出了一个评价函数，通常被称为表面积启发式函数。他们通过父节点和孩子节点的表面积之比来形式化的表述这个关系，此评价函数如下所示：此处，hit(n)是光线击中节点n的情况，Sn是节点n的表面积，c和p分别表示父节点和孩子节点。这个评价函数给出了，当用一条随机的光线同层次结构求交的时候，成本上的估计。由于没有最优的方法去有效的构造一个最优的BVH，提出了不同的构造技巧。下面，将列出比较通用的方法。在实践中，对于包围体应用的最广泛的就是轴对齐包围盒(AABB)。AABB易于实现，并且同光线的求交测试非常快。大多数有关BVH的论文在描述BVH的创建的时候，通常分别以Kay和Kajiya，或者Goldsmith和Salmon这两种基本的想法为基础。Kay和Kajiaya建议以自上而下递归的方式进行BVH的创建。Goldsmith和Salmon提出了一个更加复杂的自底向上的构造方式。Goldsmith和Salmon指出，BVH的质量同作为输入传人的三角形的顺序有关。因此，他们建议在构造BVH之前，随机打乱三角形的顺序。下述算法就是利用Kay/Kajiya的思想创建某个场景的包围体层次的方法：4．结束语本文成功的在GPU上实现了用于光线跟踪算法中的各种加速结构，并对这些加速结构在GPU上的加速效果进行了比较。均匀栅格作为第一个在CPU上实现的光线跟踪器的加速结构，也被证明是最慢的，除非是只包含一个单独的物体的场景的情况。均匀栅格不适合几何体的密度非常高的场景。另外，对于均匀栅格的CPU上的遍历表示，也需要大量的数据。Foley和Sugerman认为，对于大多数场景，KD-tree的效率要比均匀栅格高。但是，在KD-tree的遍历过程中，无论是重置阶段还是回退阶段，片元程序都非常的复杂，但这种复杂性也使得其能够在场景的几何体的密度改变的时候做出适当的调整。本文实现的BVH被证明在加速效果上要超过均匀栅格和KD-tree，在现阶段，BVH是在GPU上实现的最快的加速结构。并且在GPU上实现BVH加速结构要比实现其他加速结构更加的简单。参考文献：[1]Randima Femado编，姚勇，王小琴译．GPU精粹一实时图形编程的技术，技巧和技艺[M]．北京：人民邮电出版社，2006.[2] Matt Pharr编著，龚敏敏译.GPU精粹2-高性能图形芯片和通用计算编程技巧[M].北京：清华大学出版社．[3]昊恩华，柳有权．基于图形处理器(GPU)的通用计算叨．计算机辅助设计与图形学学报，2004，16(5): 601-[4] Philip ，David 著，周长发译，计算机图形学几何工具算法详解[M].北京：电子工业出版社，2005.[5] Martin Christen. Implementing ray tracing on GPU. Master´sthesis, University of Applied Sciences Basel,
如果那个问题采用我的的话 这个分是不是也给我呢
1 中国古代数学的发展 在古代世界四大文明中，中国数学持续繁荣时期最为长久。从公元前后至公元14世纪，中国古典数学先后经历了三次发展高潮，即两汉时期、魏晋南北朝时期和宋元时期，并在宋元时期达到顶峰。 与以证明定理为中心的希腊古典数学不同，中国古代数学是以创造算法特别是各种解方程的算法为主线。从线性方程组到高次多项式方程，乃至不定方程，中国古代数学家创造了一系列先进的算法(中国数学家称之为“术”)，他们用这些算法去求解相应类型的代数方程，从而解决导致这些方程的各种各样的科学和实际问题。特别是，几何问题也归结为代数方程，然后用程式化的算法来求解。因此，中国古代数学具有明显的算法化、机械化的特征。以下择要举例说明中国古代数学发展的这种特征。 线性方程组与“方程术” 中国古代最重要的数学经典《九章算术》(约公元前2世纪)卷8的“方程术”，是解线性方程组的算法。以该卷第1题为例，用现代符号表述，该问题相当于解一个三元一次方程组： 3x+2y+z＝39 2x+3y+z＝34 x+2y+3z＝26 《九章》没有表示未知数的符号，而是用算筹将x�y�z的系数和常数项排列成一个(长)方阵： 1 2 3 2 3 2 3 1 1 26 34 39 “方程术”的关键算法叫“遍乘直除”，在本例中演算程序如下：用右行(x)的系数(3)“遍乘”中行和左行各数，然后从所得结果按行分别“直除”右行，即连续减去右行对应各数，就将中行与左行的系数化为0。反复执行这种“遍乘直除”算法，就可以解出方程。很清楚，《九章算术》方程术的“遍乘直除” 算法，实质上就是我们今天所使用的解线性方程组的消元法，以往西方文献中称之为“高斯消去法”，但近年开始改变称谓,如法国科学院院士、原苏黎世大学数学系主任教授在他撰写的教科书[4]中就称解线性方程组的消元法为“张苍法”，张苍相传是《九章算术》的作者之一。 高次多项式方程与“正负开方术” 《九章算术》卷4中有“开方术”和“开立方术”。《九章算术》中的这些算法后来逐步推广到开更高次方的情形，并且在宋元时代发展为一般高次多项式方程的数值求解。秦九韶是这方面的集大成者，他在《数书九章》(1247年)一书中给出了高次多项式方程数值解的完整算法，即他所称的“正负开方术”。 用现代符号表达，秦九韶“正负开方术”的思路如下：对任意给定的方程 f(x)=a0xn+a1xn-1+……+an-2x2+an-1x+an=0 (1) 其中a0≠0，an<0,要求(1)式的一个正根。秦九韶先估计根的最高位数字，连同其位数一起称为“首商”，记作c，则根x＝c+h，代入(1)得 f(c+h)=a0(c+h)n+a1(c+h)n-1+……+an-1(c+h)+an=0 按h的幂次合并同类项即得到关于h的方程: f(h)=a0hn+a1hn-1+……+an-1h+an=0 (2) 于是又可估计满足新方程(2)的根的最高位数字。如此进行下去，若得到某个新方程的常数项为0，则求得的根是有理数；否则上述过程可继续下去，按所需精度求得根的近似值。 如果从原方程(1)的系数a0,a1，…,an及估值c求出新方程(2)的系数a0,a1，…,an的算法是需要反复迭代使用的，秦九韶给出了一个规格化的程序，我们可称之为“秦九韶程序”， 他在《数书九章》中用这一算法去解决各种可以归结为代数方程的实际问题，其中涉及的方程最高次数达到10次，秦九韶解这些问题的算法整齐划一，步骤分明，堪称是中国古代数学算法化、机械化的典范。 多元高次方程组与“四元术” 绝不是所有的问题都可以归结为线性方程组或一个未知量的多项式方程来求解。实际上，可以说更大量的实际问题如果能化为代数方程求解的话，出现的将是含有多个未知量的高次方程组。 多元高次方程组的求解即使在今天也绝非易事。历史上最早对多元高次方程组作出系统处理的是中国元代数学家朱世杰。朱世杰的《四元玉鉴》(1303年)一书中涉及的高次方程达到了4个未知数。朱世杰用“四元术”来解这些方程。“四元术”首先是以“天”、“地”、“人”、“物”来表示不同的未知数，同时建立起方程式，然后用顺序消元的一般方法解出方程。朱世杰在《四元玉鉴》中创造了多种消元程序。 通过《四元玉鉴》中的具体例子可以清晰地了解朱世杰“四元术”的特征。值得注意的是，这些例子中相当一部分是由几何问题导出的。这种将几何问题转化为代数方程并用某种统一的算法求解的例子，在宋元数学著作中比比皆是，充分反映了中国古代几何代数化和机械化的倾向。 一次同余方程组与“中国剩余定理” 中国古代数学家出于历法计算的需要，很早就开始研究形如： X≡Ri (mod ai) i=1,2,...,n (1) (其中ai 是两两互素的整数)的一次同余方程组求解问题。公元4世纪的《孙子算经》中已有相当于求解下列一次同余组的著名的“孙子问题”： X≡2(mod3) ≡3(mod5) ≡2(mod7) 《孙子算经》作者给出的解法，引导了宋代秦九韶求解一次同余组的一般算法——“大衍求一术”。现代文献中通常把这种一般算法称为“中国剩余定理”。 插值法与“招差术” 插值算法在微积分的酝酿过程中扮演了重要角色。在中国，早从东汉时期起，学者们就惯用插值法来推算日月五星的运动。起初是简单的一次内插法，隋唐时期出现二次插值法(如一行《大衍历》，727年)。由于天体运动的加速度也不均匀，二次插值仍不够精密。随着历法的进步，到了宋元时代，便产生了三次内插法(郭守敬《授时历》，1280年)。在此基础上，数学家朱世杰更创造出一般高次内插公式，即他所说的“招差术”。 朱世杰的公式相当于 f(n)=n△+ n(n�1)△2+ n(n�1)(n�2)△3 + n(n�1)(n�2)(n�3)△4+…… 这是一项很突出的成就。 这里不可能一一列举中国古代数学家的所有算法，但仅从以上介绍不难看到，古代与中世纪中国数学家创造的算法，有许多即使按现代标准衡量也达到了很高的水平。这些算法所表达的数学真理，有的在欧洲直到18世纪以后依赖近代数学工具才重新获得(如前面提到的高次代数方程数值求解的秦九韶程序，与1819年英国数学家W. 霍纳重新导出的“霍纳算法”基本一致；多元高次方程组的系统研究在欧洲也要到18世纪末才开始在E. 别朱等人的著作中出现；解一次同余组的剩余定理则由欧拉与高斯分别独立重新获得；至于朱世杰的高次内插公式，实质上已与现在通用的牛顿-格列高里公式相一致)。这些算法的结构，其复杂程度也是惊人的。如对秦九韶“大衍求一术”和“正负开方术”的分析表明，这些算法的计算程序，包含了现代计算机语言中构造非平易算法的基本要素与基本结构。这类复杂的算法，很难再仅仅被看作是简单的经验法则了，而是高度的概括思维能力的产物，这种能力与欧几里得几何的演绎思维风格截然不同，但却在数学的发展中起着完全可与之相媲美的作用。事实上，古代中国算法的繁荣，同时也孕育了一系列极其重要的概念，显示了算法化思维在数学进化中的创造意义和动力功能。以下亦举几例。 负数的引进 《九章算术》“方程术”的消元程序，在方程系数相减时会出现较小数减较大数的情况，正是在这里，《九章算术》的作者们引进了负数，并给出了正、负数的加减运算法则，即“正负术”。 对负数的认识是人类数系扩充的重大步骤。公元7世纪印度数学家也开始使用负数，但负数的认识在欧洲却进展缓慢，甚至到16世纪，韦达的著作还回避负数。 无理数的发现 中国古代数学家在开方运算中接触到了无理数。《九章算术》开方术中指出了存在有开不尽的情形：“若开方不尽者，为不可开”，《九章算术》的作者们给这种不尽根数起了一个专门名词——“面”。“面”，就是无理数。与古希腊毕达哥拉斯学派发现正方形的对角线不是有理数时惊慌失措的表现相比，中国古代数学家却是相对自然地接受了那些“开不尽”的无理数，这也许应归功于他们早就习惯使用的十进位制，这种十进位制使他们能够有效地计算“不尽根数”的近似值。为《九章算术》作注的三国时代数学家刘徽就在“开方术”注中明确提出了用十进制小数任意逼近不尽根数的方法，他称之为“求微数法”，并指出在开方过程中，“其一退以十为步，其再退以百为步，退之弥下，其分弥细，则……虽有所弃之数，不足言之也”。 十进位值记数制是对人类文明不可磨灭的贡献。法国大数学家拉普拉斯曾盛赞十进位值制的发明，认为它“使得我们的算术系统在所有有用的创造中成为第一流的”。中国古代数学家正是在严格遵循十进位制的筹算系统基础上，建立起了富有算法化特色的东方数学大厦。 贾宪三角或杨辉三角 从前面关于高次方程数值求解算法(秦九韶程序)的介绍我们可以看到，中国古代开方术是以�c+hn的二项展开为基础的，这就引导了二项系数表的发现。南宋数学家杨辉著《详解九章算法》(1261年)中，载有一张所谓“开方作法本源图”，实际就是一张二项系数表。这张图摘自公元1050年左右北宋数学家贾宪的一部著作。“开方作法本源图”现在就叫“贾宪三角”或“杨辉三角”。二项系数表在西方则叫“帕斯卡三角”�1654年。 走向符号代数 解方程的数学活动，必然引起人们对方程表达形式的思考。在这方面，以解方程擅长的中国古代数学家们很自然也是走在了前列。在宋元时期的数学著作中，已出现了用特定的汉字作为未知数符号并进而建立方程的系统努力。这就是以李冶为代表的“天元术”和以朱世杰为代表的“四元术”。所谓“天元术”，首先是“立天元一为某某”，这相当于“设为某某”，“天元一”就表示未知数，然后在筹算盘上布列“天元式”，即一元方程式。该方法被推广到多个未知数情形，就是前面提到的朱世杰的“四元术”。因此，用天元术和四元术列方程的方法，与现代代数中的列方程法已相类似。 符号化是近世代数的标志之一。中国宋元数学家在这方面迈出了重要一步，“天元术”和“四元术”，是以创造算法特别是解方程的算法为主线的中国古代数学的一个高峰�。 2 中国古代数学对世界数学发展的贡献 数学的发展包括了两大主要活动：证明定理和创造算法。定理证明是希腊人首倡，后构成数学发展中演绎倾向的脊梁；算法创造昌盛于古代和中世纪的中国、印度，形成了数学发展中强烈的算法倾向。统观数学的历史将会发现，数学的发展并非总是演绎倾向独占鳌头。在数学史上，算法倾向与演绎倾向总是交替地取得主导地位。古代巴比伦和埃及式的原始算法时期，被希腊式的演绎几何所接替，而在中世纪，希腊数学衰落下去，算法倾向在中国、印度等东方国度繁荣起来；东方数学在文艺复兴前夕通过阿拉伯传播到欧洲，对近代数学兴起产生了深刻影响。事实上，作为近代数学诞生标志的解析几何与微积分，从思想方法的渊源看都不能说是演绎倾向而是算法倾向的产物。 从微积分的历史可以知道，微积分的产生是寻找解决一系列实际问题的普遍算法的结果�6�。这些问题包括：决定物体的瞬时速度、求极大值与极小值、求曲线的切线、求物体的重心及引力、面积与体积计算等。从16世纪中开始的100多年间，许多大数学家都致力于获得解决这些问题的特殊算法。牛顿与莱布尼兹的功绩是在于将这些特殊的算法统一成两类基本运算——微分与积分，并进一步指出了它们的互逆关系。无论是牛顿的先驱者还是牛顿本人，他们所使用的算法都是不严格的，都没有完整的演绎推导。牛顿的流数术在逻辑上的瑕疵更是众所周知。对当时的学者来说，首要的是找到行之有效的算法，而不是算法的证明。这种倾向一直延续到18世纪。18世纪的数学家也往往不管微积分基础的困难而大胆前进。如泰勒公式，欧拉、伯努利甚至19世纪初傅里叶所发现的三角展开等，都是在很长时期内缺乏严格的证明。正如冯·诺伊曼指出的那样：没有一个数学家会把这一时期的发展看作是异端邪道；这个时期产生的数学成果被公认为第一流的。并且反过来，如果当时的数学家一定要在有了严密的演绎证明之后才承认新算法的合理性，那就不会有今天的微积分和整个分析大厦了。 现在再来看一看更早的解析几何的诞生。通常认为，笛卡儿发明解析几何的基本思想，是用代数方法来解几何问题。这同欧氏演绎方法已经大相径庭了。而事实上如果我们去阅读笛卡儿的原著，就会发现贯穿于其中的彻底的算法精神。《几何学》开宗明义就宣称：“我将毫不犹豫地在几何学中引进算术的术语，以便使自己变得更加聪明”。众所周知，笛卡儿的《几何学》是他的哲学著作《方法论》的附录。笛卡儿在他另一部生前未正式发表的哲学著作《指导思维的法则》(简称《法则》)中曾强烈批判了传统的主要是希腊的研究方法，认为古希腊人的演绎推理只能用来证明已经知道的事物，“却不能帮助我们发现未知的事情”。因此他提出“需要一种发现真理的方法”，并称之为“通用数学”(mathesis universakis)。笛卡儿在《法则》中描述了这种通用数学的蓝图，他提出的大胆计划，概而言之就是要将一切科学问题转化为求解代数方程的数学问题： 任何问题→数学问题→代数问题→方程求解而笛卡儿的《几何学》，正是他上述方案的一个具体实施和示范，解析几何在整个方案中扮演着重要的工具作用，它将一切几何问题化为代数问题，这些代数问题则可以用一种简单的、几乎自动的或者毋宁说是机械的方法去解决。这与上面介绍的古代中国数学家解决问题的路线可以说是一脉相承。 因此我们完全有理由说，在从文艺复兴到17世纪近代数学兴起的大潮中，回响着东方数学特别是中国数学的韵律。整个17—18世纪应该看成是寻求无穷小算法的英雄年代，尽管这一时期的无穷小算法与中世纪算法相比有质的飞跃。而从19世纪特别是70年代直到20世纪中，演绎倾向又重新在比希腊几何高得多的水准上占据了优势。因此，数学的发展呈现出算法创造与演绎证明两大主流交替繁荣、螺旋式上升过程： 演绎传统——定理证明活动 算法传统——算法创造活动 中国古代数学家对算法传统的形成与发展做出了毋容置疑的巨大贡献。 我们强调中国古代数学的算法传统，并不意味中国古代数学中没有演绎倾向。事实上，在魏晋南北朝时期一些数学家的工作中，已出现具有相当深度的论证思想。如赵爽勾股定理证明、刘徽“阳马”�一种长方锥体体积证明、祖冲之父子对球体积公式的推导等等，均可与古希腊数学家相应的工作媲美。赵爽勾股定理证明示意图“弦图”原型，已被采用作2002年国际数学家大会会标。令人迷惑的是，这种论证倾向随着南北朝的结束，可以说是戛然而止。囿于篇幅和本文重点，对这方面的内容这里不能详述，有兴趣的读者可参阅参考文献�3�。 3 古为今用，创新发展 到了20世纪，至少从中叶开始，电子计算机的出现对数学的发展带来了深远影响，并孕育出孤立子理论、混沌动力学、四色定理证明等一系列令人瞩目的成就。借助计算机及有效的算法猜测发现新事实、归纳证明新定理乃至进行更一般的自动推理……，这一切可以说已揭开了数学史上一个新的算法繁荣时代的伟大序幕。科学界敏锐的有识之士纷纷预见到数学发展的这一趋势。在我国，早在上世纪50年代，华罗庚教授就亲自领导建立了计算机研制组，为我国计算机科学和数学的发展奠定了基础。吴文俊教授更是从70年代中开始，毅然由原先从事的拓扑学领域转向定理机器证明的研究，并开创了现代数学的崭新领域——数学机械化。被国际上誉为“吴方法”的数学机械化方法已使中国在数学机械化领域处于国际领先地位，而正如吴文俊教授本人所说：“几何定理证明的机械化问题，从思维到方法，至少在宋元时代就有蛛丝马迹可寻，”他的工作“主要是受中国古代数学的启发”。“吴方法”，是中国古代数学算法化、机械化精髓的发扬光大。 计算机影响下算法倾向的增长，自然也引起一些外国学者对中国古代数学中算法传统的兴趣。早在上世纪70年代初，著名的计算机科学家就呼吁人们关注古代中国和印度的算法�5�。多年来这方面的研究取得了一定进展，但总的来说还亟待加强。众所周知，中国古代文化包括数学是通过著名的丝绸之路向西方传播的，而阿拉伯地区是这种文化传播的重要中转站。现存有些阿拉伯数学与天文著作中包含有一定的中国数学与天文学知识，如著名的阿尔·卡西《算术之钥》一书中有相当数量的数学问题显示出直接或间接的中国来源，而根据阿尔·卡西本人记述，他所工作的天文台中就有不少来自中国的学者。 然而长期以来由于“西方中心论”特别是“希腊中心论”的影响以及语言文字方面的障碍，有关资料还远远没有得到发掘。正是为了充分揭示东方数学与欧洲数学复兴的关系，吴文俊教授特意从他荣获的国家最高科学奖中拨出专款成立了“吴文俊数学与天文丝路基金”，鼓励支持年轻学者深入开展这方面的研究，这是具有深远意义之举。 研究科学的历史，其重要意义之一就是从历史的发展中获得借鉴和汲取教益，促进现实的科学研究，通俗地说就是“古为今用”。吴文俊对此有精辟的论述，他说：“假如你对数学的历史发展，对一个领域的发生和发展，对一个理论的兴旺和衰落，对一个概念的来龙去脉，对一种重要思想的产生和影响等这许多历史因素都弄清了，我想，对数学就会了解得更多，对数学的现状就会知道得更清楚、更深刻，还可以对数学的未来起一种指导作用，也就是说，可以知道数学究竟应该按怎样的方向发展可以收到最大的效益”。数学机械化理论的创立，正是这种古为今用原则的硕果。我国科学技术的伟大复兴，呼唤着更多这样既有浓郁的中国特色、又有鲜明时代气息的创新。
多年研究软起动器，发现软起动器对电动机的过载保护有些简单化，虽然说是反时限保护，但实际是采用定时分段的办法，有时误动作，有时烧电动机。对于电动机断续过载保护时由于电动机早已过热，那么它的过载能力已经减小，对于冷态的电动机来说，它的过载能力要比热态的电动机过载能力大的多。如果要真正反应电动机的过载能力又能对电动机起到过载保护就必需通过热积分，采用热记忆功能。这样才能保正系统的可靠性和保护的灵敏性。 两种典型的数学模型软起动器对电动机具有控制、保护、监测等功能，对电动机的热过载保护采用的反时限保护特性有多种数学模型，其中典型的有两类： （1）等I2t的时间电流特性 （2）IEC 60255-3[1]推荐的数学模型以上式中： Ir — 电流整定值I — 实际电流值t — 动作时间（s）K — 表征特性的常数α— 函数指数 脱扣器的控制方式脱扣器的控制方式可采用：（1）积分法以两种典型的数学模型为例，分别求积分值： 设定K1或K2的动作值，控制动作时间t。（2）查表法设定I—t对照表，根据当前I控制动作时间t。但是在实际运行中两种方法均存在弊端。如用积分法上述的两类数学模型都可能造成在低于动作值时仍能误动作；如用查表法在通常电流不断变化的情况下，很难合理的控制过载脱扣的延时时间。为了较好的解决低压断路器的智能控制器中长延时脱扣器的延时控制，本文试图按热保护的基本原理进行分析和探讨。 2 热保护的基本要求根据热平衡关系，电气设备的发热应等于散热与蓄热之和，即 （1）式中：P — 发热功率；Kr— 散热系数；S — 散热表面积；τ— 温升；c — 比热；G — 发热体重量；t — 时间。微分方程的解为：过载保护元件应在小于被保护电气设备温升允许值的设置值动作，断开电路。3 按热平衡原理整定过载长延时脱扣4 动作值和热时间常数的计算 动作值按电动机起动器和断路器的要求，k2应分别小于和，为同时满足这两种要求，并留有裕度，可取k2=～。由式（11）可取K=k22T（12）以K作为式（6）或（7）的截止值，当A≥K时控制器动作，实现长延时保护功能。式（9）和（10）可转换为： 热时间常数的计算在已知任意—N值下要求的tr值，即可计算T。 延时时间的计算按式（13）计算在不同过载电流下的延时时间，并考虑电流测量误差的影响，计算结果见表1（计算时取T=642s）。 5 动作值的测量和计算为测量智能脱扣器实态通电时的A值，可以采用数值积分的方法等间隔的测量电流和计算A值并与K值比较。设测量间隔为Δt,并且初始温升为0，由式（6）和（7） 上列各式中N可以为变量。逐次计算，逐次与k比较，直至Ax≥k时控制器动作。则 ……在有辅助电源的情况下，A值逐渐递减，直至软起动器重新起动，A值又开始递增；或辅助电源断开，A值清零。为防止过载脱扣后，软起动器在短时内的再接通并在短时内再分断，可设置一定的恢复时间，以保证在恢复时间内，软起动器不得起动。 6 测量误差分析对式（8）微分：对应表1中的计算值tr，在表2中列出p和f的相应值。表2 与表1中计算值tr对应的p和f值表2的误差传递系数f的估算值与表1的计算结果基本相符。由表1及表2可以看出在较低过载倍数下由电流测量误差所引起的延时时间误差较大。7 保护特性的斜率调节 建立数学模型为了满足不同的配合需要，现在有的制造厂提供了改变长延时保护特性斜率的调节功能[2]或参照IEC 60255标准提供了不同数学模型的保护特性。为了实现保护特性的斜率调节，本文推荐两种数学模型并用的方案。（1）基本数学模型经对比分析我们可以以式（7）作为基本保护特性的基本数学模型。（2）用于斜率调节的数学模型可选用国家标准GB （等同IEC 60255-3）推荐的数学模型用于斜率调节。根据GB ： （16）式中：N=I/Ir 指数α可选K为常数现以三种斜率的保护特性为例：● A型反时限 tr=K/（） （17）● B型反时限 tr=K/（N-1） （18）● C型反时限 tr=K/（N4-1） （19）K值可根据保护要求设定，或参照前述基本保护特性NIr（如N=2或N=6）对应的时间tr设定。 动作值的测量和控制将式（17）、（18）、（19）变换为A=t（） （20）A=t（N-1） （21）A=t（N4-1） （22）在实际运行中可每经过一个等间隔Δt进行一次累加，逐次计算A值，逐次与K值比较，直至达到设定值K值，求出延时时间tr。以式（21）为例，设对应式（20）和（22）可以采用同样方法进行计算和控制。但是应用此方法计算有两个问题需要解决：（1）设定N的阈值通常在K的设定值范围，在N=的条件下，计算值tr很可能小于1h，不能满足软起动器要求。为了防止在及以下的误脱扣，需设定阈值,如设定Nd=，当N≤Nd时可仍按基本数学模型控制和计算。（2）阈值上下数学模型的转换如在N>Nd时，按式（20）～（22）的数学模型进行计算和控制。现举例说明如下● 保护特性取式（21），设定K=根据式（12）计算T值，取k2=在N≤Nd时按前面第4节所述方法进行计算和控制。在N>Nd时按式（21）的数学模型进行计算，如果在尚未达到动作值时电流又下降使N≤Nd，并且当前A值为Ay。则此后需按基本数学模型累加计算A值： （24）…………式中初始值Ay为原数学模型下保留的A值。以下按前面第4节所述方法进行计算和控制。如果此后又回复N>Nd条件，应重新按式（21）的数学模型计算和控制。在反复转换数学模型时不需改变K值和当前的A值。● 保护特性取式（22），设定K=1200根据式（12）计算T值，取k2=在N≤Nd时按前面第4节所述方法进行计算和控制。在N>Nd时按式（22）的数学模型进行计算，如果在尚未达到动作值电流又下降至N≤Nd，并且当前A值为Ay。则需按式（24）计算A值。如果此后又回复N>Nd条件，应重新按式（22）的数学模型计算和控制。在反复转换数学模型时不需改变K值和当前A值。 误差分析对式（16）微分 式（19）、（20）和（21）三种数学模型时间相对误差与电流相对误差之间的传递系数计算值见表3。 表3 三种数学模型时间相对误差与电流相对误差之间的传递系数计算值由表3中可见，当α=和α=1时在Nr≥的情况下，要满足延时时间的误差不超过±10%的要求并不困难；但是在α=4时，因特性曲线斜率值大，要达到同样的指标是有一定难度的，即使电流测量误差为±2%，再考虑K的控制误差和数值化整等因素，延时时间的误差也可能大于±10%。 8 结束语本文提出的一套利用数值积分法解决反时限保护特性的实时测量和控制方法，既可比较合理、方便的提供多种保护特性，又可较好的解决负载不断变化情况下的热记忆问题，还有助于提高长延时控制单元的抗干扰能力。由于在实时控制中，微处理器在很短时间内无法完成一些函数的复杂数学运算，本文中的一些计算公式和参数在工程计算中需要进行了变换和处理，在CMC系列软起动器中得到了应用，通过实际运行达到了理想的效果。
在统计学中，统计模型是指当有些过程无法用理论分析 方法 导出其模型，但可通过试验或直接由工业过程测定数据，经过数理统计法求得各变量之间的函数关系。下文是我为大家整理的关于统计模型论文的 范文 ，欢迎大家阅读参考!
统计套利模型的理论综述与应用分析
【摘要】统计套利模型是基于数量经济学和统计学建立起来的，在对历史数据分析的基础之上，估计相关变量的概率分布，并结合基本面数据对未来收益进行预测，发现套利机会进行交易。统计套利这种分析时间序列的统计学特性，使其具有很大的理论意义和实践意义。在实践方面广泛应用于个对冲基金获取收益，理论方面主要表现在资本有效性检验以及开放式基金评级，本文就统计套利的基本原理、交易策略、应用方向进行介绍。
【关键词】统计套利 成对交易 应用分析
一、统计套利模型的原理简介
统计套利模型是基于两个或两个以上具有较高相关性的股票或者其他证券，通过一定的方法验证股价波动在一段时间内保持这种良好的相关性，那么一旦两者之间出现了背离的走势，而且这种价格的背离在未来预计会得到纠正，从而可以产生套利机会。在统计套利实践中，当两者之间出现背离，那么可以买进表现价格被低估的、卖出价格高估的股票，在未来两者之间的价格背离得到纠正时，进行相反的平仓操作。统计套利原理得以实现的前提是均值回复，即存在均值区间(在实践中一般表现为资产价格的时间序列是平稳的，且其序列图波动在一定的范围之内)，价格的背离是短期的，随着实践的推移，资产价格将会回复到它的均值区间。如果时间序列是平稳的，则可以构造统计套利交易的信号发现机制，该信号机制将会显示是否资产价格已经偏离了长期均值从而存在套利的机会 在某种意义上存在着共同点的两个证券(比如同行业的股票)， 其市场价格之间存在着良好的相关性，价格往往表现为同向变化，从而价格的差值或价格的比值往往围绕着某一固定值进行波动。
二、统计套利模型交易策略与数据的处理
统计套利具 体操 作策略有很多，一般来说主要有成对/一篮子交易，多因素模型等，目前应用比较广泛的策略主要是成对交易策略。成对策略，通常也叫利差交易，即通过对同一行业的或者股价具有长期稳定均衡关系的股票的一个多头头寸和一个空头头寸进行匹配，使交易者维持对市场的中性头寸。这种策略比较适合主动管理的基金。
成对交易策略的实施主要有两个步骤：一是对股票对的选取。海通证券分析师周健在绝对收益策略研究―统计套利一文中指出，应当结合基本面与行业进行选股，这样才能保证策略收益，有效降低风险。比如银行，房地产，煤电行业等。理论上可以通过统计学中的聚类分析方法进行分类，然后在进行协整检验，这样的成功的几率会大一些。第二是对股票价格序列自身及相互之间的相关性进行检验。目前常用的就是协整理论以及随机游走模型。
运用协整理论判定股票价格序列存在的相关性，需要首先对股票价格序列进行平稳性检验，常用的检验方法是图示法和单位根检验法，图示法即对所选各个时间序列变量及一阶差分作时序图，从图中观察变量的时序图出现一定的趋势册可能是非平稳性序列，而经过一阶差分后的时序图表现出随机性，则序列可能是平稳的。但是图示法判断序列是否存在具有很大的主观性。理论上检验序列平稳性及阶输通过单位根检验来确定，单位根检验的方法很多，一般有DF，ADF检验和Phillips的非参数检验(PP检验)一般用的较多的方法是ADF检验。
检验后如果序列本身或者一阶差分后是平稳的，我们就可以对不同的股票序列进行协整检验，协整检验的方法主要有EG两步法，即首先对需要检验的变量进行普通的线性回归，得到一阶残差，再对残差序列进行单位根检验，如果存在单位根，那么变量是不具有协整关系的，如果不存在单位根，则序列是平稳的。EG检验比较适合两个序列之间的协整检验。除EG检验法之外，还有Johansen检验，Gregory hansan法，自回归滞后模型法等。其中johansen检验比较适合三个以上序列之间协整关系的检验。通过协整检验，可以判定股票价格序列之间的相关性，从而进行成对交易。
Christian L. Dunis和Gianluigi Giorgioni(2010)用高频数据代替日交易数据进行套利，并同时比较了具有协整关系的股票对和没有协整关系股票对进行套利的立即收益率，结果显示，股票间价格协整关系越高，进行统计套利的机会越多，潜在收益率也越高。
根据随机游走模型我们可以检验股票价格波动是否具有“记忆性”，也就是说是否存在可预测的成分。一般可以分为两种情况：短期可预测性分析及长期可预测性分析。在短期可预测性分析中，检验标准主要针对的是随机游走过程的第三种情况，即不相关增量的研究，可以采用的检验工具是自相关检验和方差比检验。在序列自相关检验中，常用到的统计量是自相关系数和鲍克斯-皮尔斯 Q统计量，当这两个统计量在一定的置信度下，显著大于其临界水平时，说明该序列自相关，也就是存在一定的可预测性。方差比检验遵循的事实是：随机游走的股价对数收益的方差随着时期线性增长，这些期间内增量是可以度量的。这样，在k期内计算的收益方差应该近似等于k倍的单期收益的方差，如果股价的波动是随机游走的，则方差比接近于1;当存在正的自相关时，方差比大于1;当存在负的自相关是，方差比小于1。进行长期可预测性分析，由于时间跨度较大的时候，采用方差比进行检验的作用不是很明显，所以可以采用R/S分析，用Hurst指数度量其长期可预测性，Hurst指数是通过下列方程的回归系数估计得到的：
Ln[(R/S)N]=C+H*LnN
R/S 是重标极差，N为观察次数，H为Hurst指数，C为常数。当H>时说，说明这些股票可能具有长期记忆性，但是还不能判定这个序列是随机游走或者是具有持续性的分形时间序列，还需要对其进行显著性检验。
无论是采用协整检验还是通过随机游走判断，其目的都是要找到一种短期或者长期内的一种均衡关系，这样我们的统计套利策略才能够得到有效的实施。
进行统计套利的数据一般是采用交易日收盘价数据，但是最近研究发现，采用高频数据(如5分钟，10分钟，15分钟，20分钟收盘价交易数据)市场中存在更多的统计套利机会。日交易数据我们选择前复权收盘价，而且如果两只股票价格价差比较大，需要先进性对数化处理。Christian L. Dunis和Gianluigi Giorgioni(2010)分别使用15分钟收盘价，20分钟收盘价，30分以及一个小时收盘价为样本进行统计套利分析，结果显示，使用高频数据进行统计套利所取得收益更高。而且海通证券金融分析师在绝对收益策略系列研究中，用沪深300指数为样本作为统计套利 配对 交易的标的股票池，使用高频数据计算累计收益率比使用日交易数据高将近5个百分点。
三、统计套利模型的应用的拓展―检验资本市场的有效性
Fama(1969)提出的有效市场假说，其经济含义是:市场能够对信息作出迅速合理的反应,使得市场价格能够充分反映所有可以获得的信息,从而使资产的价格不可用当前的信息进行预测,以至于任何人都无法持续地获得超额利润.通过检验统计套利机会存在与否就可以验证资本市场是有效的的，弱有效的，或者是无效的市场。徐玉莲(2005)通过运用统计套利对中国资本市场效率进行实证研究，首先得出结论：统计套利机会的存在与资本市场效率是不相容的。以此为理论依据，对中国股票市场中的价格惯性、价格反转及价值反转投资策略是否存在统计套利机会进行检验，结果发现我国股票市场尚未达到弱有效性。吴振翔，陈敏(2007)曾经利用这种方法对我国A股市场的弱有效性加以检验，采用惯性和反转两种投资策略发现我国A股若有效性不成立。另外我国学者吴振翔，魏先华等通过对Hogan的统计套利模型进行修正，提出了基于统计套利模型对开放式基金评级的方法。
四、结论
统计套利模型的应用目前主要表现在两个方面：1.作为一种有效的交易策略，进行套利。2.通过检测统计套利机会的存在，验证资本市场或者某个市场的有效性。由于统计套利策略的实施有赖于做空机制的建立，随着我股指期货和融资融券业务的推出和完善，相信在我国会有比较广泛的应用与发展。
参考文献
[1] . Burgess：A computational Methodolology for Modelling the Dynamics of statistical arbitrage, London business school,PhD Thesis,1999.
[2]方昊.统计套利的理论模式及应用分析―基于中国封闭式基金市场的检验.统计与决策,2005,6月(下).
[3]马理,卢烨婷.沪深 300 股指期货期现套利的可行性研究―基于统计套利模型的实证.财贸研究,2011,1.
[4]吴桥林.基于沪深 300 股指期货的套利策略研究[D].中国优秀硕士学位论文.2009.
[5]吴振翔,陈敏.中国股票市场弱有效性的统计套利检验[J].系统工程理论与实践.2007,2月.
关于半参统计模型的估计研究
【摘要】随着数据模型技术的迅速发展，现有的数据模型已经无法满足实践中遇到的一些测量问题，严重的限制了现代科学技术在数据模型上应用和发展，所以基于这种背景之下，学者们针对数据模型测量实验提出了新的理论和方法，并研制出了半参数模型数据应用。半参数模型数据是基于参数模型和非参数模型之上的一种新的测量数据模型，因此它具备参数模型和非参数模型很多共同点。本文将结合数据模型技术，对半参统计模型进行详细的探究与讨论。
【关键词】半参数模型 完善误差 测量值 纵向数据
本文以半参数模型为例，对参数、非参数分量的估计值和观测值等内容进行讨论，并运用三次样条函数插值法得出非参数分量的推估表达式。另外，为了解决纵向数据下半参数模型的参数部分和非参数部分的估计问题，在误差为鞅差序列情形下，对半参数数据模型、渐近正态性、强相合性进行研究和分析。另外，本文初步讨论了平衡参数的选取问题，并充分说明了泛最小二乘估计方法以及相关结论，同时对半参数模型的迭代法进行了相关讨论和研究。
一、概论
在日常生活当中，人们所采用的参数数据模型构造相对简单，所以操作起来比较容易;但在测量数据的实际使用过程中存在着相关大的误差，例如在测量相对微小的物体，或者是对动态物体进行测量时。而建立半参数数据模型可以很好的解决和缓解这一问题：它不但能够消除或是降低测量中出现的误差，同时也不会将无法实现参数化的系统误差进行勾和。系统误差非常影响观测值的各种信息，如果能改善，就能使其实现更快、更及时、更准确的误差识别和提取过程;这样不仅可以提高参数估计的精确度，也对相关科学研究进行了有效补充。
举例来说，在模拟算例及坐标变换GPS定位重力测量等实际应用方面，体现了这种模型具有一定成功性及实用性;这主要是因为半参数数据模型同当前所使用的数据模型存在着一致性，可以很好的满足现在的实际需要。而新建立的半参数模型以及它的参数部分和非参数部分的估计，也可以解决一些污染数据的估计问题。这种半参数模型，不仅研究了纵向数据下其自身的t型估计，同时对一些含光滑项的半参数数据模型进行了详细的阐述。另外，基于对称和不对称这两种情况，可以在一个线性约束条件下对参数估计以及假设进行检验，这主要是因为对观测值产生影响的因素除了包含这个线性关系以外，还受到某种特定因素的干扰，所以不能将其归入误差行列。另外，基于自变量测量存在一定误差，经常会导致在计算过程汇总，丢失很多重要信息。
二、半参数回归模型及其估计方法
这种模型是由西方著名学者Stone在上世纪70年代所提出的，在80年代逐渐发展并成熟起来。目前，这种参数模型已经在医学以及生物学还有经济学等诸多领域中广泛使用开来。
半参数回归模型介于非参数回归模型和参数回归模型之间，其内容不仅囊括了线性部分，同时包含一些非参数部分，应该说这种模型成功的将两者的优点结合在一起。这种模型所涉及到的参数部分，主要是函数关系，也就是我们常说的对变量所呈现出来的大势走向进行有效把握和解释;而非参数部分则主要是值函数关系中不明确的那一部分，换句话就是对变量进行局部调整。因此，该模型能够很好的利用数据中所呈现出来的信息，这一点是参数回归模型还有非参数归回模型所无法比拟的优势，所以说半参数模型往往拥有更强、更准确的解释能力。
从其用途上来说，这种回归模型是当前经常使用的一种统计模型。其形式为：
三、纵向数据、线性函数和光滑性函数的作用
纵向数据其优点就是可以提供许多条件，从而引起人们的高度重视。当前纵向数据例子也非常多。但从其本质上讲，纵向数据其实是指对同一个个体，在不同时间以及不同地点之上，在重复观察之下所得到一种序列数据。但由于个体间都存在着一定的差别，从而导致在对纵向数据进行求方差时会出现一定偏差。在对纵向数据进行观察时，其观察值是相对独立的，因此其特点就是可以能够将截然不同两种数据和时间序列有效的结合在一起。即可以分析出来在个体上随着时间变化而发生的趋势，同时又能看出总体的变化形势。在当前很多纵向数据的研究中，不仅保留了其优点，并在此基础之上进行发展，实现了纵向数据中的局部线性拟合。这主要是人们希望可以建立输出变量和协变量以及时间效应的关系。可由于时间效应相对比较复杂，所以很难进行参数化的建模。
另外，虽然线性模型的估计已经取得大量的成果，但半参数模型估计至今为止还是空白页。线性模型的估计不仅仅是为了解决秩亏或病态的问题，还能在百病态的矩阵时，提供了处理线性、非线性及半参数模型等方法。首先，对观测条件较为接近的两个观测数据作为对照，可以削弱非参数的影响。从而将半参数模型变成线性模型，然后，按线性模型处理，得到参数的估计。而多数的情况下其线性系数将随着另一个变量而变化，但是这种线性系数随着时间的变化而变化，根本求不出在同一个模型中，所有时间段上的样本，亦很难使用一个或几个实函数来进行相关描述。在对测量数据处理时，如果将它看作为随机变量，往往只能达到估计的作用，要想在经典的线性模型中引入另一个变量的非线性函数，即模型中含有本质的非线性部分，就必须使用半参数线性模型。
另外就是指由各个部分组成的形态，研究对象是非线性系统中产生的不光滑和不可微的几何形体，对应的定量参数是维数，分形上统计模型的研究是当前国际非线性研究的重大前沿课题之一。因此，第一种途径是将非参数分量参数化的估计方法，也称之为参数化估计法，是关于半参数模型的早期工作，就是对函数空间附施加一定的限制，主要指光滑性。一些研究者认为半参数模型中的非参数分量也是非线性的，而且在大多数情形下所表现出来的往往是不光滑和不可微的。所以同样的数据，同样的检验方法，也可以使用立方光滑样条函数来研究半参数模型。
四、线性模型的泛最小二乘法与最小二乘法的抗差
(一)最小二乘法出现于18世纪末期
在当时科学研究中常常提出这样的问题：怎样从多个未知参数观测值集合中求出参数的最佳估值。尽管当时对于整体误差的范数，泛最小二乘法不如最小二乘法，但是当时使用最多的还是最小二乘法，其目的也就是为了估计参数。最小二乘法，在经过一段时间的研究和应用之后，逐步发展成为一整套比较完善的理论体系。现阶段不仅可以清楚地知道数据所服从的模型，同时在纵向数据半参数建模中，辅助以迭代加权法。这对补偿最小二乘法对非参数分量估计是非常有效，而且只要观测值很精确，那么该法对非参数分量估计更为可靠。例如在物理大地测量时，很早就使用用最小二乘配置法，并得到重力异常最佳估计值。不过在使用补偿最小二乘法来研究重力异常时，我们还应在兼顾着整体误差比较小的同时，考虑参数估计量的真实性。并在比较了迭代加权偏样条的基础上，研究最小二乘法在当前使用过程中存在的一些不足。应该说，该方法只强调了整体误差要实现最小，而忽略了对参数分量估计时出现的误差。所以在实际操作过程中，需要特别注意。
(二)半参模型在GPS定位中的应用和差分
半参模型在GPS相位观测中，其系统误差是影响高精度定位的主要因素，由于在解算之前模型存在一定误差，所以需及时观测误差中的粗差。GPS使用中，通过广播卫星来计算目标点在实际地理坐标系中具体坐标。这样就可以在操作过程中，发现并恢复整周未知数，由于观测值在卫星和观测站之间，是通过求双差来削弱或者是减少对卫星和接收机等系统误差的影响，因此难于用参数表达。但是在平差计算中，差分法虽然可以将观测方程的数目明显减少，但由于种种原因，依然无法取得令人满意的结果。但是如果选择使用半参数模型中的参数来表达系统误差，则能得到较好的效果。这主要是因为半参数模型是一种广义的线性回归模型，对于有着光滑项的半参数模型，在既定附加的条件之下，能够提供一个线性函数的估计方法，从而将测值中的粗差消除掉。
另外这种方法除了在GPS测量中使用之外，还可应用于光波测距仪以及变形监测等一些参数模型当中。在重力测量中的应用在很多情形下，尤其是数学界的理论研究，我们总是假定S是随机变量实际上，这种假设是合理的，近几年，我们对这种线性模型的研究取得了一些不错的成果，而且因其形式相对简洁，又有较高适用性，所以这种模型在诸多领域中发挥着重要作用。
通过模拟的算例及坐标变换GPS定位重力测量等实际应用，说明了该法的成功性及实用性，从理论上说明了流行的自然样条估计方法，其实质是补偿最小二乘方法的特例，在今后将会有广阔的发展空间。另外 文章 中提到的分形理论的研究对象应是非线性系统中产生的不光滑和不可微的几何形体，而且分形已经在断裂力学、地震学等中有着广泛的应用，因此应被推广使用到研究半参数模型中来，不仅能够更及时，更加准确的进行误差的识别和提取，同时可以提高参数估计的精确度，是对当前半参数模型研究的有力补充。
五、 总结
文章所讲的半参数模型包括了参数、非参数分量的估计值和观测值等内容，并且用了三次样条函数插值法得到了非参数分量的推估表达式。另外，为了解决纵向数据前提下，半参数模型的参数部分和非参数部分的估计问题，在误差为鞅差序列情形下，对半参数数据模型、渐近正态性、强相合性进行研究和分析。同时介绍了最小二乘估计法。另外初步讨论了平衡参数的选取问题，还充分说明了泛最小二乘估计方法以及有关结论。在对半参数模型的迭代法进行了相关讨论和研究的基础之上，为迭代法提供了详细的理论说明，为实际应用提供了理论依据。
参考文献
[1]胡宏昌.误差为AR(1)情形的半参数回归模型拟极大似然估计的存在性[J].湖北师范学院学报(自然科学版)，2009(03).
[2]钱伟民，李静茹.纵向污染数据半参数回归模型中的强相合估计[J].同济大学学报(自然科学版)，2009(08).
[3]樊明智，王芬玲，郭辉.纵向数据半参数回归模型的最小二乘局部线性估计[J].数理统计与管理，2009(02).
[4]崔恒建，王强.变系数结构关系EV模型的参数估计[J].北京师范大学学报(自然科学版).2005(06).
[5]钱伟民，柴根象.纵向数据混合效应模型的统计分析[J].数学年刊A辑(中文版).2009(04)
[6]孙孝前，尤进红.纵向数据半参数建模中的迭代加权偏样条最小二乘估计[J].中国科学(A辑：数学)，2009(05).
[7]张三国，陈希孺.EV多项式模型的估计[J].中国科学(A辑)，2009(10).
[8]任哲，陈明华.污染数据回归分析中参数的最小一乘估计[J].应用概率统计，2009(03).
[9]张三国，陈希孺.有重复观测时EV模型修正极大似然估计的相合性[J].中国科学(A辑).2009(06).
[10]崔恒建，李勇，秦怀振.非线性半参数EV四归模型的估计理论[J].科学通报，2009(23).
[11]罗中明.响应变量随机缺失下变系数模型的统计推断[D].中南大学，2011.
[12]刘超男.两参数指数威布尔分布的参数Bayes估计及可靠性分析[D].中南大学，2008.
[13]郭艳.湖南省税收收入预测模型及其实证检验与经济分析[D].中南大学，2009.
[14]桑红芳.几类分布的参数估计的损失函数和风险函数的Bayes推断[D].中南大学，2009.
[15]朱琳.服从几类可靠性分布的无失效数据的bayes分析[D].中南大学，2009.
[16]黄芙蓉.指数族非线性模型和具有AR(1)误差线性模型的统计分析[D].南京理工大学，2009.
猜你喜欢：
1. 统计学分析论文
2. 统计方面论文优秀范文参考
3. 统计优秀论文范文
4. 统计学的论文参考范例
在人类历史发展和社会生活中，数学发挥着不可替代的作用，同时也是学习和研究现代科学技术必不可少的基本工具。关于数学方面的论文我们可以写哪些呢?下面我给大家带来关于数学方向的优秀论文题目有哪些，希望能帮助到大家!
最全组合数学论文题目
1、并行组合数学模型方式研究及初步应用
2、数学规划在非系统风险投资组合中的应用
3、金融经济学中的组合数学问题
4、竞赛数学中的组合恒等式
5、概率 方法 在组合数学中的应用
6、组合数学中的代数方法
7、组合电器局部放电超高频信号数学模型构建和模式识别研究
8、概率方法在组合数学中的某些应用
9、组合投资数学模型发展的研究
10、高炉炉温组合预报和十字测温数学建模
11、证券组合的风险度量及其数学模型
12、组合数学中的Hopf方法
13、PAR方法在组合数学问题中的应用研究
14、概率方法在组合数学及混合超图染色理论中的应用
15、一些算子在组合数学中的应用
16、陀螺/磁强计组合定姿方法的相关数学问题研究
17、高中数学人教版新旧教材排列组合内容的比较研究
18、生物絮凝吸附-曝气生物滤池组合工艺处理生活污水的数学模拟研究
19、基于数学形态学-小波分析组合算法的牵引网故障判定方法
20、证券组合投资的灰色优化数学模型的研究
21、一些算子在组合数学中的应用
22、概率方法在组合数学中的应用
23、组合数学中的Hopf方法
24、概率方法在组合数学中的某些应用
25、概率方法在组合数学及混合超图染色理论中的应用
26、竞赛数学中的组合恒等式
27、Stern-Lov醩z定理及在组合结构中的应用
28、几类特殊图形的渐近估计及数值解
29、Fine格路和有禁错排
30、基于DFL的Agent自主学习模型及其应用研究
31、基于DFL的多Agent自动推理平台设计
32、预应力混凝土斜拉桥施工监控概率方法研究
33、最大概率方法与最近邻准则下的图像标注
34、亚式期权定价的偏微分方程方法和概率方法
35、编目空间碎片的碰撞概率方法研究及应用
36、基于概率方法的机器人定位
37、民用建筑内部给水设计秒流量的概率方法研究
38、图论中的组合方法和概率方法
39、物理概率方法预估贮存寿命研究
40、静载下结构参数识别的误差分析和概率方法
41、概率方法在组合计数证明中的应用
42、基于非概率方法的结构全寿命总费用评估
43、概率方法在组合数学中的应用
44、概率方法与邻点可区别全染色的色数上界
45、既有钢筋混凝土结构耐久性评定的概率方法
46、概率方法在多任务EEG脑机接口中的应用研究
47、应用概率方法对居住小区给水设计秒流量的推求
48、概率方法与图的染色问题
49、概率方法对居住小区设计秒流量的推求
50、概率方法在组合数学中的某些应用
51、概率方法在组合恒等式证明中的应用
52、遗传算法的研究与应用
53、基于空间算子代数理论的链式多体系统递推动力学研究
54、关于Weidmann猜想及具有转移条件微分算子的研究
55、实数编码遗传算法杂交算子组合研究
56、基于OWA算子理论的混合型多属性群决策研究
57、序列算子与灰色预测模型研究
58、具有转移条件的Sturm-Liouville算子和具有点作用的Schrodinger算子谱分析的研究
59、高精度径向基函数拟插值算子的构造及其应用
60、多线性算子加权Hardy算子与次线性算子的相关研究
数学建模论文题目
1、高中数学核心素养之数学建模能力培养的研究
2、小学数学建模数字化教学的设计与实施策略——以“自行车里的数学问题”为例
3、培养低年段学生数学建模意识的微课教学
4、信息化背景下数学建模教学策略研究
5、数学建模思想融入解析几何的实际应用探讨
6、以数学建模为平台培养大学生创新能力的SWOT分析──以内蒙古农业大学为例
7、基于高等数学建模思维的经济学应用
8、以数学建模促进应用型本科院校数学专业的发展
9、高等代数在数学建模中的应用探讨
10、融入数学建模思想的线性代数案例教学研究
11、以“勾股定理的应用”为例谈初中数学的建模教学
12、经管概率统计中的数学建模思想研究——评《经管与 财税 基础》
13、数学建模实例——河西学院校内充电站最佳选址问题
14、基于数学建模探讨高职数学的改革途径
15、大数据时代大学生数学建模应用能力的提升研究
16、“数学写作之初见建模”教学设计及思考
17、大学数学教学过程中数学建模意识与方法的培养简析
18、基于建模思想的高等数学应用研究
19、小学数学建模教学实践
20、依托对口支援平台培养大学生的数学建模能力
21、跨界研究在数学建模教与学中的应用
22、基于结构参数的机织物等效导热率数学建模
23、数学建模对大学生综合素质影响的调查研究
24、计算机数学建模中改进遗传算法与最小二乘法应用
25、数学建模在高中数学课堂的教学策略分析
26、发动机特性数字化处理与数学建模
27、数学建模中的数据处理——以大型百货商场会员画像描绘为例
28、数学建模竞赛对医学生 学习态度 和自学能力的影响
29、数学建模思想与高等数学教学的融会贯通
30、试论数学建模思想在小学数学教学中的应用
31、浅析飞机地面空调车风量测控系统数学建模及工程实施
32、高中数学教学中数学建模能力的培养——基于核心素养的视角
33、注重数学建模 提炼解题思路——对中考最值问题的探究
34、在数学建模教学中培养思维的洞察力
35、刍议数学建模思想如何渗透于大学数学教学中
36、数学建模竞赛背景下对高校数学教学的思考
37、数学建模课程对高职学生创新能力的培养探究
38、高等数学教学中数学建模思想方法探究
39、初中数学教学中数学建模思想的渗透
40、无线激光通信网络海量信息快速调度数学建模
41、基于多元线性回归模型的空气质量数据校准——2019年大学生数学建模竞赛D题解析
42、中学数学建模教学行为探究
43、数学建模竞赛成果诊断倒逼教学资源库优化的机制研究
44、基于数学建模活动的高校数学教学改革
45、数学建模与应用数学的结合研究
46、谈初中数学建模能力的培养
47、数学建模在初中数学应用题解答中的运用
48、基于数学建模思想的高等数学 教学方法 研究
49、数学建模融入高等数学翻转课堂模式研究
50、数学软件融入数学建模课程教学的探讨
最新小学数学教学论文题目
小学数学教材问题探析
小学数学生活化教学研究
小学数学___教学方法有效性分析
小学数学多媒体课件设计研究
小学生数学思维培养探究
小学数学中创新意识的培养
数学作业批改中巧用评语
新课标下小学数学教学改革研究
数学游戏在小学数学教学中的应用
《9和几的进位加法》教学设计
小学数学教学中素质 教育 研究
小学数学学困生的转化策略
小学数学教学中的情感教育
《六的乘法口诀》教学 反思
浅谈数学课堂中学生问题意识的培养
问答式学习课堂教学怎样转向小组合作学习
浅谈农村课堂的有效交流
浅谈在实践活动中提高学生解决实际问题的能力
浅谈小学应用题教学
浅谈学生合作意识的培养
“层次性体验”在数学课堂中的应用
数学课堂教学中学生探索能力的培养
小学数学低段学生阅读能力培养点滴
“观察、 品味、 顿悟” 我谈小学数学空间与图形教学
浅谈小学数学课堂教学中的“留白”
润物细无声--小班化数学作业面批有效策略的尝试
“我的妈妈体重 50 千克” 对培养良好数感的思考
“圆的面积” 教学一得
利用图解法解决逆推题
我教《24 时计时法》
《解简易方程》 教学反思
“可能性” 的反思
折线统计图折射出的“光芒”
《平均数》 教学反思
数学课堂上的“失误“也是一种资源
幽默语言在教学中的应用
“圆的认识” 教学片断与反思
计算机多媒体与小学数学教学的整
充分发挥学生的主体作用
“圆柱的体积” 教学反思
“平行四边形的面积” 听课反思
听“逆向求和应用题” 有感
小学低年级教学策略的实践与反思
“相遇问题” 建立“数学模型”
如何提高课堂语言评价的有效性
“20 以内退位减法” 教学反思
关于数学方向的优秀论文题目相关 文章 ：
★ 关于数学专业毕业论文题目
★ 数学方面毕业论文题目参考大全
★ 关于数学专业毕业论文题目参考
★ 数学的优秀论文
★ 数学优秀论文范文
★ 数学学术论文的题目
★ 数学教育论文题目
★ 数学教育方向的论文范文
★ 数学教育方向相关论文(2)
整理了两篇医学论文开题报告范文，欢迎参考其格式。
医学硕士开题报告范文-《治疗四肢长骨骨折不愈合》
关键词：钢丝环扎 四肢长骨骨折 济南论文 开题报告
一、选题依据、目的和意义：
骨折不愈合是骨科临床常见病症，其中以四肢长骨多发，例如胫骨，股骨，肱骨等，针对四肢长骨骨折不愈合二次手术我院多才用植骨术配合LCP重新内固定。自体髂骨作为植骨材料具有较多的优点：如取材简单、组织相容性好、无移植排斥反应、骨诱导作用强等，这些优点使得髂骨成为一种最佳的植骨供材，这在临床上已形成共识。植骨是治疗骨折不愈合的重要方法,其机制是爬行替代所引起的支架作用与供给矿物质的作用,爬行替代顺利进行的条件要求准确的复位、充分的植骨和坚强的固定。为达到充分的植骨，及早促进骨折愈合，我们采用髓内外360°植骨的方法，外用钢丝环扎，配合LCP坚强内固定，术后3~12个月内进行随访，根据愈合情况和功能恢复情况分析手术的临床疗效。选题目地在于探讨治疗四肢长骨骨折不愈合的手术改进方法和疗效，为临床治疗提供参考。
本课题以导师多年的临床资料为依据，通过对骨折不愈合手术治疗的国内外文献进行系统整理，结合山东中医药大学附属医院骨科病房对四肢长骨骨折不愈合患者的随访调查及回顾性分析，根据骨科特殊生物力学特点和导师治疗骨折不愈合的多年临床体会，分析治疗效果，并对手术中的细节问题做初步探讨与论述。同时也希望可以通过对导师的临床实践的研究、总结，能为今后的临床工作提供一些帮助和指导。
二、本课题目前国内外研究的动态、水平
治疗骨折不愈合，可分为手术治疗和非手术治疗，其中手术治疗最重要的就是植骨术加更改断端内固定。骨折不愈合应用自体骨移植治疗效果显著，已经形成共识。 植骨是治疗骨不连的重要方法,植骨方式临床多采用髓内外联合植骨。沿肌间隙进入, 骨膜下小心剥离显露骨折部位, 取出内固定器械, 清除骨断端间瘢痕, 咬除硬化骨, 打通髓腔, 修整骨折端, 手法复位, 按照骨缺损情况取骨。髓内植骨以比髓腔稍粗的骨棒,贴紧髓腔骨质;髓外上盖植骨宜用螺丝钉固定植骨块;骨碎屑充分填充残余的空隙,这样才能确实达到植骨的目的和要求。自体皮- 松质骨植骨的爬行替代缩短了骨折愈合过程,新鲜的自体骨具有生物活性,不存在免疫排异,无传染疾病的风险,同时存在骨传导和骨诱导能力。
内固定物更换得坚持以下原则，原钢板内固定者,可更换成交锁髓内针或更长的钢板置于张力侧;原交锁髓内针内固定者,可选用更大号髓内针或钢板内固定;原先短钢板内固定者，可改成较长的钢板。所有病例均需植骨。更换内固定物后,，术后石膏外固定者,应及早进行肌肉收缩锻炼活动,骨痂生长良好后,去石膏开始关节屈伸功能锻炼。但是临床上医师应该具体问题具体对待，可以根据骨痂生长情况酌情处理,出院时务必详细医嘱病人注意事项,配合医生,直到骨折完全愈合。LCP钢板内固定适用于四肢长骨骨折不愈合,可用拉力螺钉固定碎骨块及移植骨块, 并对断端行轴向加压锁定。手术关键是将骨折端的瘢痕结缔组织全部切除, 骨端硬化骨全部咬除, 露出正常骨质, 钻通髓腔, 植入的骨块必须牢固的嵌入缺损区, 间隙用松质骨填满,。应积极正确指导术后功能锻炼, 严格定期随访及指导。避免过早的不正确的负重。综上所述，对于骨折不愈合的治疗，自体骨移植疗效确切，安全稳妥，技术成熟，应用广泛，值得提倡。
三、课题研究的主要内容
1.临床资料
病例来源
本研究病例均采集于山东中医药大学附属医院骨科病房
(二)采集时间
2009年5月～2010年12月
(三)病例选择
1.诊断标准[2]
(1)病史：明确的外伤史，骨折后6个月没有愈合，并且没有进一步愈合倾向已有3个月。
(2)症状：患者骨折端成角、旋转、侧移位、短缩畸形或者节段性骨缺损、持重疼痛或不能持重、局部在应力下疼痛等。
(3)体征：局部窦道形成、流脓、假关节形成或伴有局部软组织瘢痕、缺损等
(4)辅助检查：X线表现：骨端硬化，髓腔封闭;骨端萎缩疏松，中间存在较大的间隙;或骨端硬化，相互成为杵臼假关节等这三种形式中的任何一种就可以定为骨折不愈合。
2.纳入病例标准：
(1)符合本病诊断标准;
(2)骨折平均愈合时间超过半年以上，有假关节形成;
(3)骨折平均愈合时间超过半年以上，多次复查X线拍片显示，骨折线
清晰可见，未见内外骨痂或内外骨痂极少;
(4)拍片显示骨折线增宽，骨折端骨面致密性硬化，骨髓腔封闭，骨质疏松，骨痂间无骨小梁形成，或伴有明显的骨缺损;
(5)临床表现有骨的感染、缺损、畸形、肢体不等长、局部窦道形成、流脓等。
3.排除病例标准：
(1)不符合上述诊断标准者
(2)患者有严重的内科疾病，不能够耐受手术者
(3)精神疾病患者
(4)资料不全影响判断者
2.疗效观察方法
对骨不连愈合的评价应包括骨愈合和功能恢复双重评价：
(1)骨愈合评价标准：本评价结果决定于四项指标:骨愈合、感染、畸形和肢体长度，其中骨愈合标准为X线示骨折线模糊，有连续骨痴通过骨折线，拆除或试行松动外固定物后骨折无异常活动，下肢可无痛行走，上肢持物骨折处有稳定感。 评价标准:
优:骨折愈合，无感染，断端畸形<7°，双侧肢体不等长<2 CM。
良:骨折愈合及其他三标准中两项。
可:骨折愈合及其他三标准中一项。
差:骨折未愈合或再骨折或虽愈合但不具备其他三标准中任何一个。
(2)功能评价标准
功能的评价分上肢与下肢的不同，上肢主要考虑其灵活性，而下肢主要功能为负重行走。
将下肢评价指标定为以下五项:①明显跛行;②踝或膝任何一关节僵硬(完全伸膝或踝完全背伸时，活动范围较正常或对侧丧失15°以上):③软组织情况不良;④有限制活动或影响睡眠的疼痛存在:⑤丧失工作能力或生活不能自理。
优:存在工作能力且无其他四项指标。
良:存在工作能力且具以上四指标中一至二项。
可:存在工作能力并具以上指标中三至四项。
差:丧失工作能力或生活不能自理，不考虑是否具备其他指标。
对上肢功能评价参照“Steuart和Hdlly对上肢功能评价标准”[3]
观察指标为三项:疼痛、关节活动范围、日常活动能力。
l:上肢功能评价标准
分数 痛疼 任一关节活动受限 日常活动
优 无 <20° 完全不受限
良 用力或疲劳后 20~40° 轻微受限
差 持续性 >40° 严重受限
5.课题进度及安排：
2009-05——2010-12 收集病例及随访
2010-10——2010-12 资料汇总及数据分析
2011-01——2011-03 撰写论文、定稿
四、本课题特色、预期取得的结果
骨折不愈合应用自体骨移植治疗效果显著已经形成共识，治疗过程中的经验总结需要不断的进行，更要求开展回顾性工作及进行系统的整理。因此，骨折不愈合的临床资料分析就显得尤为重要。
本课题通过搜集整理山东中医药大学附属医院骨科2009至2010年期间的患者临床资料，对于自体骨移植治疗骨折不愈合的相关性问题进行临床研究与总结。应用统计分析评分进行术前、术后及相关方面比较，对自体骨移植治疗骨折不愈合的临床疗效获得客观、真实、准确的评价，并进一步指导临床工作。
五、可行性分析
山东中医药大学附属医院骨科是山东省中医管理局评定的重点学科、重点科室,在省内知名度较高，病人来源广泛。导师王明喜主任医师从事临床工作30余年，具有丰富的临床经验，对治疗骨折不愈合做过大量研究、临床工作，并取得了良好的效果。本课题搜集整理山东中医药大学附属医院骨科近几年的临床资料，并在导师指导下对这些一手资料进行研究与总结。
四肢长骨骨折不愈合由于并发症较多，治愈比较困难，手术后功能恢复过程漫长，因此在治疗过程中，经验的总结是非常必需的，也是可行的。本课题主要研究山东省中医院近年应用钢丝环扎360°植骨配合LCP内固定治疗四肢长骨骨折不愈合的治疗效果分析情况，因此在选题上可行性较强。课题的研究也得到了学校、附院等各部门、科室的大力支持。相信可以圆满地完成课题。
主要参考文献
[1] 胥少汀，葛宝丰，徐印坎，等.实用骨科学[M]. 北京人民军医出版社,2007
[2] 王亦璁，等.骨与关节损伤[M].人民卫生出版社，2007
[3] 夏和桃.组合式外固定器简介[EB/OL].北京骨外固定技术研究所，2005
[4] 蒋协运.骨科临床疗效评价标准[M].人民卫生出版社, 2005
[5] Boyd HB, Lipinski SW, Wiley JH J Bone Joint Surg Am, 1961; 43(2):159—168
[6] Audige L,Griffin D,Bhandari M,et al. Path analysis of factors for delayed healing and nonunion in 416 operatively treated tibial shaft fractures [J].Orthop Relat Res, 2005; 438:221～232.
[7] James J McCarthy, John Nonunions [J/OL]. eMedicine, Apr 9, 2004
[8] 周来喜，林本丹，钟志刚，等.胫骨骨折三种固定器械的生物力学比较和临床研究[J].骨与关节损伤杂志，2000;15 ( 5 ): 428 ～430
[9] 李峰，欧阳跃平.骨不连临床研究进展[J]. 国际骨科学杂志, 2007;28(2)：117～119
[10] Harvey EJ,Henley MB,Swiontkowsid MF,et al. Iryury,2003; 34(2):111～116
[11] 任可，张春才，赵建宁，等.持续动态压应力下骨折愈合时软骨内骨化的特点及其机制[J]. 解剖学杂志，2008;31(4)：570～574
[12] 李兴华.交锁髓内钉治疗胫骨骨折不愈合[J].中医正骨,2007;19(2)：43～44
[13] 吴国华.4种固定方法治疗胫腓骨双骨折的疗效对比[J].现代中西医结合杂志，2004;13( 1): 47～48
[14] 权毅，潘显明，王元山，等.交锁髓内钉断钉与骨不连的力学研究及临床意义[J]. 中国矫形外科杂志，2003;11(3,4)：207～209
[15] Streeker W, Suger G, Kinzl L. Local complications of in-tramedullary nailing[J].Othopade,1996;25:274～291.
[16] Farmanullah, Muhammad Shoaib Khan ，Syed Muhammad OF MANAGEMENT OF TIBIAL NON-UNION DEFECT WITH ILIZAROV FIXATOR [J]. Ayub Med Coll Abbottabad， 2007; 19(3)
[17] Davies R,Holt N,Nayagam care of pin sites with external fixation[J].http:// Journalof bone and Joint Surgery，2005;87:716～719
[18] 李起鸿.骨外固定技术临床应用中的几个问题[J]. 中华骨科杂志，1996;16: 604.
[19] 杨立民.当今骨科感染的特点与对策[J].骨与关节损伤杂志,1999;14: 139
[20] 陈文红，史振满，陈建常，等.感染性胫骨骨折不愈合的外固定架治疗[J].中国骨与关节损伤杂志，2007;22(8)：691～692
[21] Sluzalek M,Gazdzik T S,M rozek S et al. External fixation in thetreatment of severe tibial fractures complicated by soft tissue injury [J]. Ortop traumato1 Rehabil,2004;6( 1): 103 ～112
[22] 胡蕴玉，陆裕朴，刘伟.异种骨移植修复骨缺损实验研究[J].中华骨科杂志，1990;10: 33～36
[23] Lu WJ，L B,B NR, et al. Chin ，2006; 9(5):272～275
医学论文开题报告范文：细胞信号转导与靶向抗肿瘤药物的研发
一、 选题的目的和意义
定量结构活性关系(Quantitative Structure-Activity Relationships，简称 QSAR)是20世纪60年代发展起来的一门新兴学科，是由结构活性关系(Structure-Activity Relationship，简称 SAR )发展而来的。QSAR 是通过对已知结构且有生物活性系列化合物(如一系列有相同药理作用的结构相似的化合物)进行化学信息学的计算, 选用适当的数学模型建立活性与化合物结构之间定量关系，解释由于分子结构的变化影响化合物生物活性的改变,推测其可能的作用机理。然后建立有效的QSAR模型，如果有新化合物的出现，且其结构数据已知，可以预测其生物活性，也可以优化结构改变现有化合物的结构以提高其生物活性。这种方法广泛应用于药物、农药、化学毒剂等生物活性分子的合理设计。在经历40多年的发展过程中，定量构效活性关系在国际上已成为一个相当活跃的研究领域。
尽管肿瘤的化学治疗已取得重大进展，新的抗肿瘤药物不断出现，但肿瘤的化学治疗仍存在着许多问题，这主要是因为实体肿瘤占恶性肿瘤的90%但多数实体瘤如肺癌、肝癌、结肠癌及胰腺癌等还缺乏有效的药物;现有的抗肿瘤药物毒副反应太大，缺乏选择性;肿瘤细胞对抗肿瘤药物产生抗药性[1]。
QSAR主要侧重于药物早期的研究和发展，为新药物分子的筛选的和设计开拓了新的途径[2]，在受体结构已知的情况下，对抗肿瘤药物进行定量构效活性关系研究，用生成与受体结构互补的配体的方法来发现可以针对特定肿瘤、特定靶点的非细胞毒类药物，使之更具有选择性和针对性。随着新QSAR模型的建立，极大地缩短了新药合成的时间，降低了开发成本，并能在某种程度上预测药物对特定肿瘤人群的有效性。为肿瘤治疗起到了积极地推动作用。
二、国内外研究现状
肿瘤的化学治疗药物发展很快,每年都有大量抗肿瘤药物研究文献发表，各国对抗肿瘤药物的开发也予以高度重视和大量投资，美国就此专门成立了美国国立癌症研究(National Cancer Institute，简称NCI)，成为了世界抗肿瘤的权威机构。
国内研发方向主要以含中草药及其活性成分的抗肿瘤药物，可以归纳为以下几个方面：(1)对现有药物进行结构改造以改善其药理学特性，如增加水溶性、降低毒副作用等;(2)以新的作用机理或作用靶点为指导寻找新的活性物质作为先导化合物;(3)发现新的作用靶点。在当前生物学的后基因时代，科学家们要面对数千个潜在的药物靶点，探讨它们与小分了化合物的相互作用;(4)加强定量构效活性构关系研究.
近年来随着分子生物学和计算机技术的迅速发展，使得开发新药的技术路线发生了重大变革。国际上越来越多的研究机构在新抗肿瘤药物的开发中使用计算机辅助分子设计,它大大加快了新药设计的速度，节省了创制新药工作的人力和物力，使药物学家能够以理论作指导，有目的地开发新药。计算机辅助分子设计主要分两种情况:一种是在受体结构已知的情况下,采用生成与受体结构互补的配体的方法来寻找新药物;另一种是在受体结构未知的情况下,采用对一组具有类似活性的化合物建立定量结构活性关系,在此模型基础上进行结构修饰来预测生成新的化合物。
QSAR作为抗肿瘤药物设计研究中的一个重要计算方法和常用手段,在新药的开发和研制过程中占据了重要位置。近半个世纪以来,QSAR研究对有机合成化学、药物化学及药物设计的发展起了巨大的推动作用,已经成为研究物质理化性质与生物活性以寻求分子解释的一个强有力工具。下面就定量活性结构活性关系研究的一些常见方法作简要地介绍如下。
1、二维定量结构活性关系方法(2D-QSAR)传统的二维定量结构活性关系方法很多,有Hansch法、基团贡献法和分子连接性指数法等[3] 。
其中最为著名、应用最为广泛的是Hansch 法。 它假设同系列化合某些生物活性的变化是和它们某些可测量的物理化学性质(疏水性、电性质和空间立体性质等)的变化相联系的,并假定这些因子是彼此孤立的,采用多重自由能相关法,借助多重线性回归等统计方法就可以得到定量结构活性关系模型。
基团贡献法是Free-Wilson 在对有机物亚结构信息和生物毒性的相关研究基础上建立的一种方法。这种模式认为有机物与受体间的毒性效应是该有机物特定位置上不同取代基团毒性贡献的加和。Free-Wilson 法仅适用于具有相同母体结构的有机物,常被用来对有机物进行毒性初评。
分子连接性指数法(Molecular connective index ,MCI) 是由Kier 和Hall 提出的。它是根据分子中各个骨架原子排列或相连接的方式来描述分子的结构性质。MCI 是一种拓扑学参数,有零阶项(0Xv )、易阶项(1Xv )、二阶项(2Xv ) 等等，可以根据分子的结构式和原子的点价(δ) 计算得到,与有机物的毒性数据有较好的相关性。MCI 能较强地反映分子的立体结构,但反映子电子结构的能力较弱,因此缺乏明确的物理意义,但由于其具有方便、简单且不依赖于实验等优点,近年来得到广泛应用和发展[4～8]。
2、三维定量结构活性关系方法(3D-QSAR)随着结构活性关系理论和统计方法的进一步发展,20 世纪80 年代,三维结构信息被陆续引入到定量结构活性关系研究中, 即3D-QSAR。与2D-QSAR 比较,3D-QSAR 方法在物理化学上的意义更为明确,能间接反映药物分子和靶点之间的非键相互作用特征。因此,近十多年来3D - QSAR 方法得到了迅速的发展和广泛的应用，研究方法也很多[9] ,比如分子形状分(molecular shape analysis ,MSA) ，距离几何方法( distance geometry , DG ,比较分子力场分析(comparative molecular field analysis ,CoMFA) ,比较分子相似因子分析( comparative molecular similarityindices analysi CoMSIA) 以及虚拟受体( phesudo receptor) 等方法。其中应用最为广泛的CoMFA 方法。
3、随着技术的发展和生产技术的进步，又出现了一些先进的方法来构建QSAR模型，都具有很好的预测能力。其中又以启发发(heuristic method,简称HM)，支持向量机(Support Vector Machine，简称SVM)，基因表达式编程(Gene Expression Programming,简称GEP)比较常见。支持向量机(Support Vector Machine)是Vapnik[10]等人根据统计学理论提出的一种新的通用学习方法，它是建立在统计学理论的VC维理论和结构风险最小原理基础上的，能较好地解决小样本、非线性、高维数等实际问题[11-12]，已成功地应用于分类、函数逼近和时间序列预测等方面[13-15];基因表达式编程(GEP)是基于生物学遗传思想，保持了生物学的特性，具有良好的结果重现性，同时也能够进行“遗传变异”控制，最终能获得可靠的实验效果。
三、主要研究内容
1、查阅中外文文献选取数据来源。
2、理化参数与结构参数的计算。
3、具体的结构参数的分析。
4、SVM与GEP的方法研究。
5、定量结构关系式的建立。
6、定量结构关系式的验证。
7、得出结论和总结。
四、论文工作计划
3月中旬—4月初：选题。
4月初—4月底：查阅资料，熟悉实验原理及方法，准备开题报告。
5月10日： 开题。
5月初日—5月底日：进行毕业设计实验，记录数据，撰写论文。
6月初日—6月中旬日：进行毕业论文答辩。
五、参考文献
[1] 任华益. 中华综合临床医学杂志(山东) , 2005, 7(2): 28 -33.
[2] 徐娟，王林编译. 计算机辅助药物设计中的QSAR和QSMR研究. 国外医学•药学分册, 2003, 30(3): 135-138.
[3] 郭宗儒. 药物化学总论. 北京:中国医药科技出版社, 1994. 108.
[4] Bakulh H Rao, Shyam R, Asolekar. QSAR models to predict effect of ionic strength on sorption of chlorinated benzenes and phenols at sediment-water interface. Water Research, 200l, 35(14): 3391-3401.
[5] 冯长君, 堵锡华, 唐自强. 取代芳烃对发光菌、大型蚤、呆鲦鱼急性毒性的QSAR研究. 应用化学, 2002, 19(11): 1037 -1042.
[6]秦正龙, 冯长君. 取代苯酚的定量结构-活性P性质相关性研究. 有机化学, 2003, 23(7): 654-658.
[7] 堵锡华. 取代芳香族化合物生物活性的拓扑学 . 南昌大学学报(理学版), 2005, 29(2): 155-160.
[8] Aleksandar Sablji C. QSAR models for estimating properties of persistent organic pollutants required in evaluation of their environmental fate and risk. Chemosphere, 2001, 43(3): 363 -375.
[9] 徐筱杰, 侯廷军，乔学斌，章威. 计算机辅助药物分子设计. 北京: 化学工业出版社, 2004.
[10] Vapnik Nature of Statistical Learning Theory.
NY: Springer-Verlag,1995.
[11] 阎辉，张学工，李衍达. 应用SVM方法进行沉淀微相识别.物探化探计算技术, 2000, (2): 158 -164.
[12] 张学工. 关于统计学习理论与支持向量机. 自动化学报, 2000, (1): 32 -42.
[13] Vapnik V, Golowich S, Smola A. Supportvector method for function approximation, regression estimation, and signal processing. In: Mozer M, Jordan M, Petsche Teds. Neural Information Processing System, MIT Press, 1997-09.
[14]马云潜，张学工. 支持向量机函数拟合在分形插值中的应用.清华大学学报(自然科学版) , 2000, (3): 76- 78.
[15] Muller K-R, Smola A J, Ratsch G . Predicting time series with support vector machines. In:Proc of ICANN 97, Springer Lecture Notes In Computer Science, 1997: 999-1005.
在统计学中，统计模型是指当有些过程无法用理论分析 方法 导出其模型，但可通过试验或直接由工业过程测定数据，经过数理统计法求得各变量之间的函数关系。下文是我为大家整理的关于统计模型论文的 范文 ，欢迎大家阅读参考!
统计套利模型的理论综述与应用分析
【摘要】统计套利模型是基于数量经济学和统计学建立起来的，在对历史数据分析的基础之上，估计相关变量的概率分布，并结合基本面数据对未来收益进行预测，发现套利机会进行交易。统计套利这种分析时间序列的统计学特性，使其具有很大的理论意义和实践意义。在实践方面广泛应用于个对冲基金获取收益，理论方面主要表现在资本有效性检验以及开放式基金评级，本文就统计套利的基本原理、交易策略、应用方向进行介绍。
【关键词】统计套利 成对交易 应用分析
一、统计套利模型的原理简介
统计套利模型是基于两个或两个以上具有较高相关性的股票或者其他证券，通过一定的方法验证股价波动在一段时间内保持这种良好的相关性，那么一旦两者之间出现了背离的走势，而且这种价格的背离在未来预计会得到纠正，从而可以产生套利机会。在统计套利实践中，当两者之间出现背离，那么可以买进表现价格被低估的、卖出价格高估的股票，在未来两者之间的价格背离得到纠正时，进行相反的平仓操作。统计套利原理得以实现的前提是均值回复，即存在均值区间(在实践中一般表现为资产价格的时间序列是平稳的，且其序列图波动在一定的范围之内)，价格的背离是短期的，随着实践的推移，资产价格将会回复到它的均值区间。如果时间序列是平稳的，则可以构造统计套利交易的信号发现机制，该信号机制将会显示是否资产价格已经偏离了长期均值从而存在套利的机会 在某种意义上存在着共同点的两个证券(比如同行业的股票)， 其市场价格之间存在着良好的相关性，价格往往表现为同向变化，从而价格的差值或价格的比值往往围绕着某一固定值进行波动。
二、统计套利模型交易策略与数据的处理
统计套利具 体操 作策略有很多，一般来说主要有成对/一篮子交易，多因素模型等，目前应用比较广泛的策略主要是成对交易策略。成对策略，通常也叫利差交易，即通过对同一行业的或者股价具有长期稳定均衡关系的股票的一个多头头寸和一个空头头寸进行匹配，使交易者维持对市场的中性头寸。这种策略比较适合主动管理的基金。
成对交易策略的实施主要有两个步骤：一是对股票对的选取。海通证券分析师周健在绝对收益策略研究―统计套利一文中指出，应当结合基本面与行业进行选股，这样才能保证策略收益，有效降低风险。比如银行，房地产，煤电行业等。理论上可以通过统计学中的聚类分析方法进行分类，然后在进行协整检验，这样的成功的几率会大一些。第二是对股票价格序列自身及相互之间的相关性进行检验。目前常用的就是协整理论以及随机游走模型。
运用协整理论判定股票价格序列存在的相关性，需要首先对股票价格序列进行平稳性检验，常用的检验方法是图示法和单位根检验法，图示法即对所选各个时间序列变量及一阶差分作时序图，从图中观察变量的时序图出现一定的趋势册可能是非平稳性序列，而经过一阶差分后的时序图表现出随机性，则序列可能是平稳的。但是图示法判断序列是否存在具有很大的主观性。理论上检验序列平稳性及阶输通过单位根检验来确定，单位根检验的方法很多，一般有DF，ADF检验和Phillips的非参数检验(PP检验)一般用的较多的方法是ADF检验。
检验后如果序列本身或者一阶差分后是平稳的，我们就可以对不同的股票序列进行协整检验，协整检验的方法主要有EG两步法，即首先对需要检验的变量进行普通的线性回归，得到一阶残差，再对残差序列进行单位根检验，如果存在单位根，那么变量是不具有协整关系的，如果不存在单位根，则序列是平稳的。EG检验比较适合两个序列之间的协整检验。除EG检验法之外，还有Johansen检验，Gregory hansan法，自回归滞后模型法等。其中johansen检验比较适合三个以上序列之间协整关系的检验。通过协整检验，可以判定股票价格序列之间的相关性，从而进行成对交易。
Christian L. Dunis和Gianluigi Giorgioni(2010)用高频数据代替日交易数据进行套利，并同时比较了具有协整关系的股票对和没有协整关系股票对进行套利的立即收益率，结果显示，股票间价格协整关系越高，进行统计套利的机会越多，潜在收益率也越高。
根据随机游走模型我们可以检验股票价格波动是否具有“记忆性”，也就是说是否存在可预测的成分。一般可以分为两种情况：短期可预测性分析及长期可预测性分析。在短期可预测性分析中，检验标准主要针对的是随机游走过程的第三种情况，即不相关增量的研究，可以采用的检验工具是自相关检验和方差比检验。在序列自相关检验中，常用到的统计量是自相关系数和鲍克斯-皮尔斯 Q统计量，当这两个统计量在一定的置信度下，显著大于其临界水平时，说明该序列自相关，也就是存在一定的可预测性。方差比检验遵循的事实是：随机游走的股价对数收益的方差随着时期线性增长，这些期间内增量是可以度量的。这样，在k期内计算的收益方差应该近似等于k倍的单期收益的方差，如果股价的波动是随机游走的，则方差比接近于1;当存在正的自相关时，方差比大于1;当存在负的自相关是，方差比小于1。进行长期可预测性分析，由于时间跨度较大的时候，采用方差比进行检验的作用不是很明显，所以可以采用R/S分析，用Hurst指数度量其长期可预测性，Hurst指数是通过下列方程的回归系数估计得到的：
Ln[(R/S)N]=C+H*LnN
R/S 是重标极差，N为观察次数，H为Hurst指数，C为常数。当H>时说，说明这些股票可能具有长期记忆性，但是还不能判定这个序列是随机游走或者是具有持续性的分形时间序列，还需要对其进行显著性检验。
无论是采用协整检验还是通过随机游走判断，其目的都是要找到一种短期或者长期内的一种均衡关系，这样我们的统计套利策略才能够得到有效的实施。
进行统计套利的数据一般是采用交易日收盘价数据，但是最近研究发现，采用高频数据(如5分钟，10分钟，15分钟，20分钟收盘价交易数据)市场中存在更多的统计套利机会。日交易数据我们选择前复权收盘价，而且如果两只股票价格价差比较大，需要先进性对数化处理。Christian L. Dunis和Gianluigi Giorgioni(2010)分别使用15分钟收盘价，20分钟收盘价，30分以及一个小时收盘价为样本进行统计套利分析，结果显示，使用高频数据进行统计套利所取得收益更高。而且海通证券金融分析师在绝对收益策略系列研究中，用沪深300指数为样本作为统计套利 配对 交易的标的股票池，使用高频数据计算累计收益率比使用日交易数据高将近5个百分点。
三、统计套利模型的应用的拓展―检验资本市场的有效性
Fama(1969)提出的有效市场假说，其经济含义是:市场能够对信息作出迅速合理的反应,使得市场价格能够充分反映所有可以获得的信息,从而使资产的价格不可用当前的信息进行预测,以至于任何人都无法持续地获得超额利润.通过检验统计套利机会存在与否就可以验证资本市场是有效的的，弱有效的，或者是无效的市场。徐玉莲(2005)通过运用统计套利对中国资本市场效率进行实证研究，首先得出结论：统计套利机会的存在与资本市场效率是不相容的。以此为理论依据，对中国股票市场中的价格惯性、价格反转及价值反转投资策略是否存在统计套利机会进行检验，结果发现我国股票市场尚未达到弱有效性。吴振翔，陈敏(2007)曾经利用这种方法对我国A股市场的弱有效性加以检验，采用惯性和反转两种投资策略发现我国A股若有效性不成立。另外我国学者吴振翔，魏先华等通过对Hogan的统计套利模型进行修正，提出了基于统计套利模型对开放式基金评级的方法。
四、结论
统计套利模型的应用目前主要表现在两个方面：1.作为一种有效的交易策略，进行套利。2.通过检测统计套利机会的存在，验证资本市场或者某个市场的有效性。由于统计套利策略的实施有赖于做空机制的建立，随着我股指期货和融资融券业务的推出和完善，相信在我国会有比较广泛的应用与发展。
参考文献
[1] . Burgess：A computational Methodolology for Modelling the Dynamics of statistical arbitrage, London business school,PhD Thesis,1999.
[2]方昊.统计套利的理论模式及应用分析―基于中国封闭式基金市场的检验.统计与决策,2005,6月(下).
[3]马理,卢烨婷.沪深 300 股指期货期现套利的可行性研究―基于统计套利模型的实证.财贸研究,2011,1.
[4]吴桥林.基于沪深 300 股指期货的套利策略研究[D].中国优秀硕士学位论文.2009.
[5]吴振翔,陈敏.中国股票市场弱有效性的统计套利检验[J].系统工程理论与实践.2007,2月.
关于半参统计模型的估计研究
【摘要】随着数据模型技术的迅速发展，现有的数据模型已经无法满足实践中遇到的一些测量问题，严重的限制了现代科学技术在数据模型上应用和发展，所以基于这种背景之下，学者们针对数据模型测量实验提出了新的理论和方法，并研制出了半参数模型数据应用。半参数模型数据是基于参数模型和非参数模型之上的一种新的测量数据模型，因此它具备参数模型和非参数模型很多共同点。本文将结合数据模型技术，对半参统计模型进行详细的探究与讨论。
【关键词】半参数模型 完善误差 测量值 纵向数据
本文以半参数模型为例，对参数、非参数分量的估计值和观测值等内容进行讨论，并运用三次样条函数插值法得出非参数分量的推估表达式。另外，为了解决纵向数据下半参数模型的参数部分和非参数部分的估计问题，在误差为鞅差序列情形下，对半参数数据模型、渐近正态性、强相合性进行研究和分析。另外，本文初步讨论了平衡参数的选取问题，并充分说明了泛最小二乘估计方法以及相关结论，同时对半参数模型的迭代法进行了相关讨论和研究。
一、概论
在日常生活当中，人们所采用的参数数据模型构造相对简单，所以操作起来比较容易;但在测量数据的实际使用过程中存在着相关大的误差，例如在测量相对微小的物体，或者是对动态物体进行测量时。而建立半参数数据模型可以很好的解决和缓解这一问题：它不但能够消除或是降低测量中出现的误差，同时也不会将无法实现参数化的系统误差进行勾和。系统误差非常影响观测值的各种信息，如果能改善，就能使其实现更快、更及时、更准确的误差识别和提取过程;这样不仅可以提高参数估计的精确度，也对相关科学研究进行了有效补充。
举例来说，在模拟算例及坐标变换GPS定位重力测量等实际应用方面，体现了这种模型具有一定成功性及实用性;这主要是因为半参数数据模型同当前所使用的数据模型存在着一致性，可以很好的满足现在的实际需要。而新建立的半参数模型以及它的参数部分和非参数部分的估计，也可以解决一些污染数据的估计问题。这种半参数模型，不仅研究了纵向数据下其自身的t型估计，同时对一些含光滑项的半参数数据模型进行了详细的阐述。另外，基于对称和不对称这两种情况，可以在一个线性约束条件下对参数估计以及假设进行检验，这主要是因为对观测值产生影响的因素除了包含这个线性关系以外，还受到某种特定因素的干扰，所以不能将其归入误差行列。另外，基于自变量测量存在一定误差，经常会导致在计算过程汇总，丢失很多重要信息。
二、半参数回归模型及其估计方法
这种模型是由西方著名学者Stone在上世纪70年代所提出的，在80年代逐渐发展并成熟起来。目前，这种参数模型已经在医学以及生物学还有经济学等诸多领域中广泛使用开来。
半参数回归模型介于非参数回归模型和参数回归模型之间，其内容不仅囊括了线性部分，同时包含一些非参数部分，应该说这种模型成功的将两者的优点结合在一起。这种模型所涉及到的参数部分，主要是函数关系，也就是我们常说的对变量所呈现出来的大势走向进行有效把握和解释;而非参数部分则主要是值函数关系中不明确的那一部分，换句话就是对变量进行局部调整。因此，该模型能够很好的利用数据中所呈现出来的信息，这一点是参数回归模型还有非参数归回模型所无法比拟的优势，所以说半参数模型往往拥有更强、更准确的解释能力。
从其用途上来说，这种回归模型是当前经常使用的一种统计模型。其形式为：
三、纵向数据、线性函数和光滑性函数的作用
纵向数据其优点就是可以提供许多条件，从而引起人们的高度重视。当前纵向数据例子也非常多。但从其本质上讲，纵向数据其实是指对同一个个体，在不同时间以及不同地点之上，在重复观察之下所得到一种序列数据。但由于个体间都存在着一定的差别，从而导致在对纵向数据进行求方差时会出现一定偏差。在对纵向数据进行观察时，其观察值是相对独立的，因此其特点就是可以能够将截然不同两种数据和时间序列有效的结合在一起。即可以分析出来在个体上随着时间变化而发生的趋势，同时又能看出总体的变化形势。在当前很多纵向数据的研究中，不仅保留了其优点，并在此基础之上进行发展，实现了纵向数据中的局部线性拟合。这主要是人们希望可以建立输出变量和协变量以及时间效应的关系。可由于时间效应相对比较复杂，所以很难进行参数化的建模。
另外，虽然线性模型的估计已经取得大量的成果，但半参数模型估计至今为止还是空白页。线性模型的估计不仅仅是为了解决秩亏或病态的问题，还能在百病态的矩阵时，提供了处理线性、非线性及半参数模型等方法。首先，对观测条件较为接近的两个观测数据作为对照，可以削弱非参数的影响。从而将半参数模型变成线性模型，然后，按线性模型处理，得到参数的估计。而多数的情况下其线性系数将随着另一个变量而变化，但是这种线性系数随着时间的变化而变化，根本求不出在同一个模型中，所有时间段上的样本，亦很难使用一个或几个实函数来进行相关描述。在对测量数据处理时，如果将它看作为随机变量，往往只能达到估计的作用，要想在经典的线性模型中引入另一个变量的非线性函数，即模型中含有本质的非线性部分，就必须使用半参数线性模型。
另外就是指由各个部分组成的形态，研究对象是非线性系统中产生的不光滑和不可微的几何形体，对应的定量参数是维数，分形上统计模型的研究是当前国际非线性研究的重大前沿课题之一。因此，第一种途径是将非参数分量参数化的估计方法，也称之为参数化估计法，是关于半参数模型的早期工作，就是对函数空间附施加一定的限制，主要指光滑性。一些研究者认为半参数模型中的非参数分量也是非线性的，而且在大多数情形下所表现出来的往往是不光滑和不可微的。所以同样的数据，同样的检验方法，也可以使用立方光滑样条函数来研究半参数模型。
四、线性模型的泛最小二乘法与最小二乘法的抗差
(一)最小二乘法出现于18世纪末期
在当时科学研究中常常提出这样的问题：怎样从多个未知参数观测值集合中求出参数的最佳估值。尽管当时对于整体误差的范数，泛最小二乘法不如最小二乘法，但是当时使用最多的还是最小二乘法，其目的也就是为了估计参数。最小二乘法，在经过一段时间的研究和应用之后，逐步发展成为一整套比较完善的理论体系。现阶段不仅可以清楚地知道数据所服从的模型，同时在纵向数据半参数建模中，辅助以迭代加权法。这对补偿最小二乘法对非参数分量估计是非常有效，而且只要观测值很精确，那么该法对非参数分量估计更为可靠。例如在物理大地测量时，很早就使用用最小二乘配置法，并得到重力异常最佳估计值。不过在使用补偿最小二乘法来研究重力异常时，我们还应在兼顾着整体误差比较小的同时，考虑参数估计量的真实性。并在比较了迭代加权偏样条的基础上，研究最小二乘法在当前使用过程中存在的一些不足。应该说，该方法只强调了整体误差要实现最小，而忽略了对参数分量估计时出现的误差。所以在实际操作过程中，需要特别注意。
(二)半参模型在GPS定位中的应用和差分
半参模型在GPS相位观测中，其系统误差是影响高精度定位的主要因素，由于在解算之前模型存在一定误差，所以需及时观测误差中的粗差。GPS使用中，通过广播卫星来计算目标点在实际地理坐标系中具体坐标。这样就可以在操作过程中，发现并恢复整周未知数，由于观测值在卫星和观测站之间，是通过求双差来削弱或者是减少对卫星和接收机等系统误差的影响，因此难于用参数表达。但是在平差计算中，差分法虽然可以将观测方程的数目明显减少，但由于种种原因，依然无法取得令人满意的结果。但是如果选择使用半参数模型中的参数来表达系统误差，则能得到较好的效果。这主要是因为半参数模型是一种广义的线性回归模型，对于有着光滑项的半参数模型，在既定附加的条件之下，能够提供一个线性函数的估计方法，从而将测值中的粗差消除掉。
另外这种方法除了在GPS测量中使用之外，还可应用于光波测距仪以及变形监测等一些参数模型当中。在重力测量中的应用在很多情形下，尤其是数学界的理论研究，我们总是假定S是随机变量实际上，这种假设是合理的，近几年，我们对这种线性模型的研究取得了一些不错的成果，而且因其形式相对简洁，又有较高适用性，所以这种模型在诸多领域中发挥着重要作用。
通过模拟的算例及坐标变换GPS定位重力测量等实际应用，说明了该法的成功性及实用性，从理论上说明了流行的自然样条估计方法，其实质是补偿最小二乘方法的特例，在今后将会有广阔的发展空间。另外 文章 中提到的分形理论的研究对象应是非线性系统中产生的不光滑和不可微的几何形体，而且分形已经在断裂力学、地震学等中有着广泛的应用，因此应被推广使用到研究半参数模型中来，不仅能够更及时，更加准确的进行误差的识别和提取，同时可以提高参数估计的精确度，是对当前半参数模型研究的有力补充。
五、 总结
文章所讲的半参数模型包括了参数、非参数分量的估计值和观测值等内容，并且用了三次样条函数插值法得到了非参数分量的推估表达式。另外，为了解决纵向数据前提下，半参数模型的参数部分和非参数部分的估计问题，在误差为鞅差序列情形下，对半参数数据模型、渐近正态性、强相合性进行研究和分析。同时介绍了最小二乘估计法。另外初步讨论了平衡参数的选取问题，还充分说明了泛最小二乘估计方法以及有关结论。在对半参数模型的迭代法进行了相关讨论和研究的基础之上，为迭代法提供了详细的理论说明，为实际应用提供了理论依据。
参考文献
[1]胡宏昌.误差为AR(1)情形的半参数回归模型拟极大似然估计的存在性[J].湖北师范学院学报(自然科学版)，2009(03).
[2]钱伟民，李静茹.纵向污染数据半参数回归模型中的强相合估计[J].同济大学学报(自然科学版)，2009(08).
[3]樊明智，王芬玲，郭辉.纵向数据半参数回归模型的最小二乘局部线性估计[J].数理统计与管理，2009(02).
[4]崔恒建，王强.变系数结构关系EV模型的参数估计[J].北京师范大学学报(自然科学版).2005(06).
[5]钱伟民，柴根象.纵向数据混合效应模型的统计分析[J].数学年刊A辑(中文版).2009(04)
[6]孙孝前，尤进红.纵向数据半参数建模中的迭代加权偏样条最小二乘估计[J].中国科学(A辑：数学)，2009(05).
[7]张三国，陈希孺.EV多项式模型的估计[J].中国科学(A辑)，2009(10).
[8]任哲，陈明华.污染数据回归分析中参数的最小一乘估计[J].应用概率统计，2009(03).
[9]张三国，陈希孺.有重复观测时EV模型修正极大似然估计的相合性[J].中国科学(A辑).2009(06).
[10]崔恒建，李勇，秦怀振.非线性半参数EV四归模型的估计理论[J].科学通报，2009(23).
[11]罗中明.响应变量随机缺失下变系数模型的统计推断[D].中南大学，2011.
[12]刘超男.两参数指数威布尔分布的参数Bayes估计及可靠性分析[D].中南大学，2008.
[13]郭艳.湖南省税收收入预测模型及其实证检验与经济分析[D].中南大学，2009.
[14]桑红芳.几类分布的参数估计的损失函数和风险函数的Bayes推断[D].中南大学，2009.
[15]朱琳.服从几类可靠性分布的无失效数据的bayes分析[D].中南大学，2009.
[16]黄芙蓉.指数族非线性模型和具有AR(1)误差线性模型的统计分析[D].南京理工大学，2009.
猜你喜欢：
1. 统计学分析论文
2. 统计方面论文优秀范文参考
3. 统计优秀论文范文
4. 统计学的论文参考范例
牛顿插值法是插值法利用函数f(x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f(x)的近似值。
牛顿插值法相对于拉格朗日插值法具有承袭性的优势，即在增加额外的插值点时，可以利用之前的运算结果以降低运算量。牛顿插值法的特点在于：每增加一个点，不会导致之前的重新计算，只需要算和新增点有关的就可以了。
如果这特定函数是多项式，就称它为插值多项式。利用插值基函数很容易得到拉格朗日插值多项式，公式结构紧凑，在理论分析中甚为方便，但当插值节点增减时全部插值基函数均要随之变化，整个公式也将发生变化，这在实际计算中是很不方便的，为了克服这一缺点，提出了牛顿插值。
插值有点像拟合，通过拟合后的公式来计算缺失的点，但是拟合可能不会要求拟合的曲线一定要通过样本点，满足本身的指定的条件即可。插值在满足曲线穿过样本点的基础上可能还有其他的技术指标，单纯的穿过，分段线性插值即可。
将你假设的数字代入，得到方程（▲Z）/（250-291）=（▲Z-69）/(291-300)等式变换，化简，得到（▲Z-69）*41=9*（▲Z）所以解得▲Z=
一、含义不同：
两者都是通过给定n+1个互异的插值节点，求一条n次代数曲线近似地表示待插值的函曲线，这就叫做代数插值；Lagrange插值代数和Newton法插值都属于代数插值的范畴。
Lagrange插值和Newton法插值的结果和余项都是一致的，因为都是利用n次多项式插值，所以一致。
二、计算不同：
Lagrange插值法是通过构造n+1个n次基本多项式，线性组合而得到的。而Newton法插值是通过求各阶差商，递推得到的一个f（x）=f（x0）+（x-x0）f［x0，x1］+（x-x0）（x-x1）f［x0，x1，x2］+（x-x0）（x-x（n-1））f［x0，x1，xn］这样的公式，代进去就可以得到。
牛顿插值法的特点在于：
每增加一个点，不会导致之前的重新计算，只需要算和新增点有关的就可以。
假设已知n+1n+1个点相对多项式函数ff的值为：（x0，f（x0）），（x1，f（x1）），（x2，f（x2）），（xn，f（xn）），求此多项式函数f。
先从求满足两个点（x0，f（x0）），（x1，f（x1））的函数f1（x）说起：
假设f1（x）=f（x0）+b1（x−x0）f1（x）=f（x0）+b1（x−x0），增加一个点，（x0，f（x0）），（x1，f（x1）），（x2，f（x2）），求满足这三个点的函数f2（x）：
假设f2（x）=f1（x）+b2（x−x0）（x−x1）
以上内容参考：百度百科-牛顿插值法
你写的式子谈不上“二阶”，二阶均差是在一阶均差的基础上再作一次差商，类似于二阶导数是一阶导数的导数。而你的式子还是在“一阶”的级别。另外，均差的定义方式直接和插值多项式的形式相关，牛顿插值方式的均差是建立在相邻节点之间；如果你改变了均差的定义方式，相应也得改变插值多项式的形式，使得插值节点上的函数值等于给定值。换句话，定义了你自己的插值方法，当然能不能成功还不一定。
大家知道吗?致谢的言语行为，通过多样化的组合与变换承担了交际中的礼貌功能，那么 毕业 论文的致谢词怎么写呢?下面是我给大家带来的毕业论文致谢词怎么写_毕业论文致谢词模板5篇，以供大家参考，我们一起来看看吧!
毕业论文致谢词模板(一)
三年的博士学习生涯即将结束,博士学位论文是三年学习成果的一个 总结 ,而致谢又是对三年经历的一个总结,今将其发于此,以飨众友.
“雄关漫道真如铁，而今迈步从头越”。在论文即将梓刻之际，回想三年博士生涯，给我提供帮助的人颇多。最应感谢我的导师孔昭君研究员。首先，感谢他将我纳入门墙。正因于他的接纳，才使我完成了一个科学的门外汉向一个科学的探索者的蜕变。其次，感谢他在我入学之际尽力帮我申请到公费名额。正因为这个公费名额，让我减轻了许多经济上的压力，使我在整整的三年攻博期间有心情坐下来真正做点学问。最后感谢他在我找工作之际，为帮我找一份合适的工作费神费力。师恩如海，衔草难报。孔老师深厚的学术造诣、严谨的治学风格、严肃的科学态度、乐观开朗的性格、幽默机智的谈吐深深地令我折服。尤其是他渊博的学识、敏锐的洞察力，为我论文的选题、撰写提供了关键启发和帮助。
感谢我的家人常年对我的支持和理解!
感谢在百忙之中评审我博士学位论文的各位专家和学者!
特别感谢所有支持过我、帮助过我、批评过我、鼓励过我和理解过我的人们!
最后，感谢岁月与困难对我的磨砺
“长风破浪会有时，直挂云帆济沧海”，这是我最喜欢的诗句。以此作为结尾，与所有要感谢的人共勉，相信自己，追逐最初的梦想，永不放弃!
毕业论文致谢词模板(二)
本研究及学位论文是在我的导师___老师的亲切关怀和悉心指导下完成的。他严肃的科学态度，严谨的治学精神，精益求精的工作作风，深深地感染和激励着我。_老师不仅在学业上给我以精心指导，同时还在思想、生活上给我以无微不至的关怀，在此谨向_老师致以诚挚的谢意和崇高的敬意。我还要感谢在一起愉快的度过毕业论文小组的同学们，正是由于你们的帮助和支持，我才能克服一个一个的困难和疑惑，直至本文的顺利完成。
在论文即将完成之际，我的心情无法平静，从开始进入课题到论文的顺利完成，有多少可敬的师长、同学、朋友给了我无言的帮助，在这里请接受我诚挚的谢意!最后我还要感谢培养我长大含辛茹苦的父母，谢谢你们!
最后，再次对关心、帮助我的老师和同学表示衷心地感谢!
毕业论文致谢词模板(三)
本论文是在导师詹怀宇教授和付时雨研究员的悉心指导下完成的。导师渊博的专业知识，严谨的治学态度，精益求精的工作作风，诲人不倦的高尚师德，严以律己、宽以待人的崇高风范，朴实无华、平易近人的人格魅力对我影响深远。不仅使我树立了远大的学术目标、掌握了基本的研究 方法 ，还使我明白了许多待人接物与为人处世的道理。本论文从选题到完成，每一步都是在导师的指导下完成的，倾注了导师大量的心血。在此，谨向导师表示崇高的敬意和衷心的感谢!
本论文的顺利完成，离不开各位老师、同学和朋友的关心和帮助。在此感谢岳保珍高工、张曾教授、李兵云老师、何婉芬老师的指导和帮助;感谢重点实验室的??等老师的指导和帮助;感谢福建农林大学谢拥群教授、陈礼辉教授、黄六莲高工的关心、支持和帮助，在此表示深深的感谢。没有他们的帮助和支持是没有办法完成我的博士学位论文的，同窗之间的友谊永远长存。
毕业论文致谢词模板(四)
岁月如梭，如歌。转眼间，三年的研究生求学生活即将结束，站在毕业的门槛上，回首往昔，奋斗和辛劳成为丝丝的记忆，甜美与欢笑也都尘埃落定。交通大学以其优良的学习风气、严谨的科研氛围教我求学，以其博大包容的情怀胸襟、浪漫充实的校园生活育我成人。值此毕业论文完成之际，我谨向所有关心、爱护、帮助我的人们表示最诚挚的感谢与最美好的祝愿。
本论文是在导师李霞副教授的悉心指导之下完成的。三年来，导师渊博的专业知识，严谨的治学态度，精益求精的工作作风，诲人不倦的高尚师德，朴实无华、平易近人的人格魅力对我影响深远。导师不仅授我以文，而且教我做人，虽历时三载，却赋予我终生受益无穷之道。本论文从选题到完成，几易其稿，每一步都是在导师的指导下完成的，倾注了导师大量的心血，在此我向我的导师李霞副教授表示深切的谢意与祝福!
本论文的完成也离不开其他各位老师、同学和朋友的关心与帮助。在此也要感谢王保民等各位老师在论文开题、初稿、预答辩期间所提出的宝贵意见，感谢管理学院能源课题组为本论文提供的数据和建议，还要感谢同门的师兄师妹们，在科研过程中给我以许多鼓励和帮助。回想整个论文的写作过程，虽有不易，却让我除却浮躁，经历了思考和启示，也更加深切地体会了法学的精髓和意义，因此倍感珍惜。
毕业论文致谢词模板(五)
本人的本科毕业设计论文一直是在导师__教授的悉心指导下进行的。__教授治学态度严谨，学识渊博，为人和蔼可亲。并且在整个毕业设计过程中，__教授不断对我得到的结论进行总结，并提出新的问题，使得我的毕业设计课题能够深入地进行下去，也使我接触到了许多理论和实际上的新问题，使我做了许多有益的思考。在此表示诚挚的感谢和由衷的敬意。
__教授在__方面具有丰富的实践 经验 ，对我的实验工作给予了很多的指导和帮助，使我能够将理论中的结果与实际相结合。另外，他对待问题的严谨作风也给我留下了深刻的印象。在此表示深深的谢意。
___老师在___方面经验非常丰富。尽管我是第一次__，难免遇到许多比较低级的问题，___老师却都极其耐心地予以解答，在此表示深深的谢意。
此外，我还要感谢许多学长和同学在整个过程中的帮助和配合。他们包括
本文是在__老师精心指导和大力支持下完成的。_老师以其严谨求实的治学态度、高度的敬业精神、兢兢业业、孜孜以求的工作作风和大胆创新的进取精神对我产生重要影响。她渊博的知识、开阔的视野和敏锐的思维给了我深深的启迪。同时，在此次毕业设计过程中我也学到了许多了关于微生物发酵方面的知识，实验技能有了很大的提高。
另外，我还要特别感谢师姐对我实验以及论文写作的指导，她为我完成这篇论文提供了巨大的帮助。还要感谢，__和__同学对我的无私帮助，使我得以顺利完成论文。同时实验室的__老师也时常帮助我，在此我也衷心的感谢他。
最后，再次对关心、帮助我的老师和同学表示衷心地感谢!
毕业论文致谢词模板相关 文章 ：
★ 毕业论文致谢词模版4篇
★ 最新大学生毕业论文致谢词范文
★ 本科毕业论文致谢词范文5篇
★ 毕业论文致谢词致谢老师5篇
★ 毕业论文致谢语五篇
★ 毕业论文致谢信范文精选5篇2021
★ 毕业论文最后的致谢词10篇
★ 毕业论文的致谢5篇
★ 毕业论文致谢词部分范文4篇
★ 大学毕业论文答辩致谢词范文
毕业论文致谢10篇（精选）最后深深的感谢呵护我成长的父母。每当我遇到困难的时候，父母总是第一个给我鼓励的人。回顾20多年来走过的路，每一个脚印都浸满着他们无私的关爱和谆谆教诲，10年的在外求学之路，寄托着父母对我的殷切期望。他们在精神上和物质上的无私支持，坚定了我追求人生理想的信念。父母的爱是天下最无私的最宽厚的爱。大恩无以言报，惟有以永无止境的奋斗，期待将来辉煌的事业让父母为之骄傲。我亦相信自己能达到目标。 最后，向所有关心我的亲人、师长和朋友们表示深深的谢意 【6】弹指一挥间，大学四年已经接近了尾声。当自己怀着忐忑不安的心情完成这篇毕业论文的时候，自己也从当年一个从山里走出的懵懂孩子变成了一个成熟青年，回想自己的十几年的求学生涯，虽然只是一个本科毕业，但也实属不容易。首先，从小学到大学的学费和生活费就不是一个小数目，这当然要感谢我的爸爸妈妈，他们都是农民，没有他们的勤勤恳恳和细心安排，我是无论如何也完成不了我的大学生活。没有他们的支持和鼓励，我也不可能完成我经济学第二专业的学习。当然，一个农民家庭要同时供两个大学生上学，没有别人的帮助和接济是相当困难的。因此我要感谢那些在我求学时对我经济和精神上帮助的亲戚、朋友、老师和同学们，我的生活因你们而精彩和充实。 这里严谨的学风、优美的校园环境使我大学四年过的很充实和愉快。我第一学期是在哲学系度过的。在这短短的半年时间里，我有幸和许多优秀的同学一起学习，听睿智的哲学老师讲授哲学。后来由于一些本可抗拒但一时冲动的原因，我依依不舍的离开了哲学系，离开了亲爱的哲学系的老师和同学们。虽然只有短短的半年时光，但我确实学到了很多有用的知识，尤其是对我思想和方法上的指导。这些有用的东西一直对我大学的学习和生活有很重要的指导作用我相信，这些东西将伴随我走完整个人生的道路。现在回想起在哲学班的日子，还是那么的温馨和惬意，我不能不感谢当时哲学班的每一位同学和老师，跟你们在一起学习、生活，那真是其乐融融，妙不可言! 大一下学期，我转专业来到了一个新集体，开始了新的专业。虽然老师和同学都是新的，专业是新的，但我仍然感受到了那种来自老师和同学们的热情和融洽!这对我融入新环境，对新专业产生兴趣是非常重要的。公管院的老师更是让我难忘，他们严谨的学术态度，幽默风趣的授课方式给我留下了深刻的影响。在这篇论文构思和写作过程，我的论文指导老师王国华教授，对我论文的完成更是功不可没，王老师每次给我的疑问给予细心的解答并给出写作建议，对我的论文进行细心的修改，使得我的论文结构一步一步的完善，内容日趋丰满。没有王老师的细心指导，这篇论文是不可能完成的。 书到用时方恨少，在这篇论文的写作过程中，我深感自己的水平还非常的欠缺。生命不息，学习不止，人生就是一个不断学习和完善的过程，敢问路在何方?路在脚下! 【7】岁月如梭，如歌。转眼间，三年的研究生求学生活即将结束，站在毕业的门槛上，回首往昔，奋斗和辛劳成为丝丝的记忆，甜美与欢笑也都尘埃落定。交通大学以其优良的学习风气、严谨的科研氛围教我求学，以其博大包容的情怀胸襟、浪漫充实的校园生活育我成人。值此毕业论文完成之际，我谨向所有关心、爱护、帮助我的人们表示最诚挚的感谢与最美好的祝愿。 本论文是在导师李霞副教授的悉心指导之下完成的。三年来，导师渊博的专业知识，严谨的治学态度，精益求精的工作作风，诲人不倦的高尚师德，朴实无华、平易近人的人格魅力对我影响深远。导师不仅授我以文，而且教我做人，虽历时三载，却赋予我终生受益无穷之道。本论文从选题到完成，几易其稿，每一步都是在导师的指导下完成的，倾注了导师大量的心血，在此我向我的导师李霞副教授表示深切的谢意与祝福! 本论文的完成也离不开其他各位老师、同学和朋友的关心与帮助。在此也要感谢王保民等各位老师在论文开题、初稿、预 答辩 期间所提出的宝贵意见，感谢管理学院能源课题组为本论文提供的数据和建议，还要感谢同门的师兄师妹们，在科研过程中给我以许多鼓励和帮助。回想整个论文的写作过程，虽有不易，却让我除却浮 躁 ，经历了思考和启示，也更加深切地体会了法学的精髓和意义，因此倍感珍惜。 【8】时光的流逝也许是客观的，然而流逝的快慢却纯是一种主观的感受。当自己终于可以从考研、找工作、毕业论文的压力下解脱出来，长长地吁出一口气时，我忽然间才意识到，原来四年已经过去，到了该告别的时候了。一念至此，竟有些恍惚，所谓白驹过隙、百代过客云云，想来便是这般惆怅了。 可是怅然之后，总要说些什么。大学四年，生活其实很简单，只是一些读书、写字和考试的周而复始。如果把这种单调的生活看作一场场循环的演出，提供那么我只是一个安静的演员。这篇毕业论文也称不上什么精彩的台词，只不过是这种循环演出即将告一段落时的谢幕词。但是无论多么蹩脚的演员，无论台下有多少观众，即使是只说给自己听，在他谢幕时也总要感激一些人，是这些人帮助他走上舞台，成功或者不那么成功地“演出”。 我在这里首先要感谢的是我的学位论文指导老师——张杰/刘淑芹老师。这篇毕业论文从开题、资料查找、修改到最后定稿，如果没有她的心血，尚不知以何等糟糕的面目出现。我很自豪有这样一位老师，她值得我感激和尊敬。 感谢和我共度四年美好大学生活的2004级工商管理与法学双学位班的全体同学。感谢管理学院和法学院的所有授课老师，你们使我终身受益。感谢所有关心、鼓励、支持我的家人、亲戚和朋友。 【9】三年的博士学习生涯即将结束,博士学位论文是三年学习成果的一个总结,而致谢又是对三年经历的一个总结,今将其发于此,以 飨 众友. “雄关漫道真如铁，而今迈步从头越”。在论文即将梓刻之际，回想三年博士生涯，给我提供帮助的人颇多。最应感谢我的导师孔昭君研究员。首先，感谢他将我纳入门墙。正因于他的接纳，才使我完成了一个科学的门外汉向一个科学的探索者的蜕变。其次，感谢他在我入学之际尽力帮我申请到公费名额。正因为这个公费名额，让我减轻了许多经济上的压力，使我在整整的三年攻博期间有心情坐下来真正做点学问。最后感谢他在我找工作之际，为帮我找一份合适的工作费神费力。师恩如海，衔草难报。孔老师深厚的学术造诣、严谨的治学风格、严肃的科学态度、乐观开朗的性格、幽默机智的谈吐深深地令我折服。尤其是他渊博的学识、敏锐的洞察力，为我论文的选题、撰写提供了关键启发和帮助。