求解抛物型方程反问题的一个稳定的数值算法_第1页
求解抛物型方程反问题的一个稳定的数值算法_第2页
求解抛物型方程反问题的一个稳定的数值算法_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

求解抛物型方程反问题的一个稳定的数值算法摘要在本文中，在热传导反问题（IHCP）的条件下，我们考虑一个未知边界的数值分析方法。给出的热传导方程，边界条件，和初始条件都是一个无限维的形式。为解决反问题的数值算法基于有限差分法和最小二乘法。为了规范所得的奇异线性方程组，我们采用吉洪诺夫正则化方法来获得稳定的数值近似解。关键词：热传导反问题，有限差分法，一致性，稳定性，正则化方法。1。简介到目前为止，为了分析反问题和热传导反问题，通过测量内部温度以及对热通量的估计，我们已经研究出了各种方法本文旨在确定一个未知函数在IHCP中的计算。通过使用一个位于身体内一个点的传感器并测量在这一点的温度，然后对这一问题应用有限差分法计算，从而确定了这一问题的一个稳定的数值解。本文的安排如下：在第2节中，我们制定了一个一维的IHCP。在第3节中，用有限差分法离散来求解IHCP。最小二乘法和吉洪诺夫正则化方法将在第4节中讨论。最后的数值实验将在第5节给出。2。IHCP的一个计算公式在这一节中，让我们考虑以下IHCP：，并且：其他的都已知，只有q(t)未知，也就是待求的。3。数值方法概述：为了解决上述问题，我们假定q（t）已知，用如下公式离散上述问题：，其中：则上述离散公示可写为：用这一迭代公式，可以一次求出所有点的值。定理1：若U是上述问题的精确解，则当N趋于无穷时，误差是趋于0的。定理2：上述离散公式是无条件稳定。4。最小二乘法和吉洪诺夫正则化上文中提到的可由最小二乘法来确定，当和在时是很小的，如0.001。则E可表示为，要将其最小化，转而解决以下方程组：。为了确定，可以解决如下方程：（13），其中A是的矩阵，。在数学上，IHCPs属于不适定问题类，即测量数据的小误差可导致的估计数量的大偏差。不适定性估计问题的物理原因是，热传导的扩散在固体表面引起的阻尼的变化。结果是，表面的大振幅的变化源于测量数据的小振幅的变化。数据中的误差和噪声可以被误认为是引起通计过程中表面状态变化的原因。由于矩阵是奇异的，方程（13）的解可以由噪声数据的一个大变化而产生，而正则化方法正是用于控制噪声的传播。在我们的计算采用吉洪诺夫正则化方法[1]解矩阵我们用吉洪诺夫正则化来解决方程（13）,吉洪诺夫正则化中定义为如下问题的解：，而对于（13）式定义为：

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论