一类抛物型程的系数反演问题.pdfVIP

摘 要 由于应用技术的需要，对微分方程反问题的研究近年来发展的非常迅速．热 传导方程未知系数的反演问题就是一类重要的偏微分方程反问题，所反演的参数 一般为热传导系数，热源系数及其源项等等．所利用的反演输入数据有边界测量 数据，内点测量值，给定时刻的测量值等．一般而言，这是一类不适定的问题． 对下面热传导方程定解问题 I饥一d2蚝z+q(x)u=0(z，t)∈(o，1)×(o，T) {“(o，t)=，o(t)，t‘(1，t)=go(t)t∈(0，T)(o．0．1) l t‘(z，o)=t‘o(功 o∈(0，1) 问题的解u(z，t)依赖于系数口(z)，因此这里的反问题是非线性的不适定问题．另一 方面，解。‰t)关于时间t本质上是以指数衰减的，因此u0，T)包含q(x)的信息非 常弱．这加剧了所考虑的问题的不适定性．在原反问题有饵的情况下，数值反演系 数q(x)的问题包含了问题非线性和不适定性的处理两个方面． 本文提出了一种新的迭代方法来求解此问题．有别于通常的优化迭代方法 (它在每一步同时求解近似的q(≈)(z)和对应的正问题的解Ⅱ(o)(z，￡))，本文通过变换 ”(z，t)=并箸将上述反问题转化为下面的非线性问题 (毛t)∈(0，1)×(o，T) l仇一a2％=知2％(露如出+酱) (0．0．2) {v(o，t)=啪(t)v(i，t)=枷(￡)t∈(0，T) I。(z，0)=口(岛刃+口2％产一822斧z∈(o，1) 的迭代求解问题，然后再利用变换下的关系 如胁2”比榭+z毫拱厂嘶㈡出+』r嘶m出+笔者卜㈦如 取时间tl∈(o，T)足够的接近T，由 删甜笔者叫圳 (咖) 近似得出g(z)．此方法通过直接解一个新的非线性正问题达到求解反问豚的目的， 在求解。(。，t)的过程中无须讨论q(x)的影响．在给出上述变换反演方法可行的条件 以后，对上述转化后的新问题，我们讨论了它与原来反问题的等价性。进而讨论了 迭代解序列的收敛性．最后对每一迭代过程中的线性正问题，我们沿着时间方向 利用Sweep方法得到”(^J涵亡)的离散解．通过上面的工作，最后我们完成了q(x)的 反演．数值试验说明了本文反演方法的有效性．无论是迭代求解非线性正问题的 过程，还是进而求出々(￡)，都需要计算反演输入数据的二阶导致，这实际上是原问 题不适定性的反映．我们利用了正则化的方法来处理此不适定性的过程，这个过 程尤其是对噪音输入数据的处理是必要的．本文提出的反演方法本质上是一种将 问题的非线性和不适定性分开处理的方法，并且求解系数的不适定性本质上反映 在对末始时刻给定数据的微分中． 关键词；反问题，抛物型方程。迭代方法，Sweep方法，正则化方法，数值解 Abstract Withthe of in researches promotionmanyappliedproblemsarisingengineering＆reas．the ontheinverse by differential have in problemsgovernedpartial equationsdevelopedrapidly recent mostofthe in aJ-eascanbedescribedPDE8．The years，since phenomenaapplled by of reconstructionofcoefficientsof isonethe