一阶偏微分方程教程.ppt

表1后方垂直截面*164442132041418848070221441151583522091230272835313138681789212602092111731727213130337651715453表2前方垂直截面1633244322431667124055609810482322137155311974247853304441143114960318230129420611036258188柱在底面相互接触，侧面互不接触。注：一个计量单位表示4.753×10?13mol/l的多巴胺，表中的数字如28353表示中间后部圆柱含有28353个单位的药物。试估计药物在它影响区域中的分布。图1药物含量分布图*假设*忽略样本组织中多巴胺的原始含量；01假设样本组织的大小与其余脑组织的大小相比可以忽略，且样本组织不靠近脑边界；02假设大脑是均匀的，扩散和衰减决定了多巴胺在大脑中迁移过程，忽略对流过程的影响；03假设仅进行一次多巴胺注射，注射位于原点；04假设注射和取样之间有较长时间间隔，可以忽略注射过程和各个柱体取样时间的差别。05方程模型的建立*在假设中，我们认为分子扩散和成分衰减是主要的迁移方式。成分的衰减显然可看作是与多巴胺的含量密度(浓度)C(x,y,z,t)(剂量单位/mm3)成正比的。设该比例系数(即成分衰减系数)为k。下面来考虑分子的扩散。先考虑物质仅沿x轴方向的扩散。如图2，垂直于x轴任作柱体，截面为A。图2一方面，该柱体中物质扩散时位于区间段[x,x+?x]的物质在时间段[t,t+dt]内的增量为01另一方面，扩散理论中的涅恩斯特实验定律告诉我们，在时间?t内，物质沿x轴正向流过x处截面(面积为A)的质量为(其中Ex0称为x方向的扩散系数)：02同理，在时间?t内，物质沿x轴正向流过x+?x处截面的质量为：01于是在时间?t内，流入微元体[x,x+?x]内的物质质量为：0204030102显然，即由于大脑是均匀的，显然沿各方向的扩散是一致的，且扩散系数Ex(Ey,Ez)均为常数，再考虑到成分的衰减，应有数学建模培训一阶偏微分方程模型偏微分方程的相关概念*偏微分方程：一个包含有多元未知函数及其偏导数的等式。方程中所含未知函数偏导数的最高阶数称为该方程的阶。如：等。如果方程关于未知函数及其各阶偏导数是线性的，则称它是线性的；如果它关于所有最高阶偏导数是线性的，则称它是拟线性的。定解问题：定解条件通常包括边界条件和初始条件两种。含有定解条件的方程求解问题称为定解问题，包括初值问题（Cauchy问题）、边值问题和混合问题。方程的解：若函数u连续并具有方程所涉及的连续的各阶偏导数，且该函数代入方程使得方程在某区域内成为恒等式，则称该函数为方程在该区域内的解（古典解）。满足某些特定条件的解称为特解，这些条件称为定解条件。一般情况下，一个具有n个自变量的m阶方程的解可以含有m个n-1元任意函数，这样的解称为通解。一阶线性偏微分方程*01040203一阶齐次线性偏微分方程显然方程有平凡解u=常数。一般求其非平凡解。以下以含有3个自变量的方程为例，一般形式为常微分方程组*0102称为方程(2)的特征方程组，每一条积分曲线03称为方程(2)的特征线。若由特征方程组(3)推出函数恒为常数，则称该函数为方程组(3)的一个首次积分。若特征方程组(3)的3个独立的首次积分为则特征方程组(3)的通解为例1.求解方程组*，因此得到另一个首次积分为于是原方程的隐式通解为积分为1解：由再由2得得3，因此得到一个首次由(3)可得*0102若(4)的一个首次积分为03的一个首次积分。04于是得到方程组(3)的一个等价形式：05，则它也称为(3)对于一阶齐次线性偏微分方程(2)与它的特征方程组(3)或(4)，我们有以下结论：证明从略。定理1：连续可微函数是(2)的解的充分必要条件是是(4)的首次积分。定理2：如果是(4)的两个独立的首次积分，则它们的任意连续可微函数是(2)的通解。例2.求解方程*解：特征方程组为01010203040506或首次积分