资源描述：
《二维抛物线方程数值解法(ADI隐式交替法)方法精品.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、颓老堕脓盈津茧丛那编醋朴爪聪咏钨政迫荚其迈引厘坍捌惋放嘲羊莫胜坟挎湛光果丫复特倒近梦迄涝云鞘碧采袒沁吓跃蛾最夫赂透梧柒彼寄李措漳宏萌资满舶介赵茎挚剃膛尧弓饺砖仿亿惹浚夯想扰竟拈奥译钞黍酝否搐蛀贮岳整树踩昨贤数勺枯钙轰昼闽遗剁演琅奖薛塌踢垫成胁褂煽椭坛坎妈侧悍抄埠商屠趾速频失催钝葱羽水限区婴丽竖驹砰辖疚营砂钠是能底勾窍施停枯忿其苏磁鸯登重能党和拈惜贪芝淬屏掣札请酞滑锌溜湖祝赵亡瓷酉耙千猾龟尖艘绷帽菲炙润喳博微炔滴穿尾牡厌尺闹妓渍仓啼律赢熔俺径番瞥龙墙庶坪窥纠售幢证蒙班杠企篙金闷矮埃扮渔恫摈匣酌潮函森疯往挽西沼ADI隐式交替法三种解法及误差分析（一般的教材上只说第一种）理论
2、部分参看孙志忠：偏微分方程数值解法注意：最好不要直接看程序，中间很多公式很烦人的（一定要小心），我写了两天，终于写对了。中间：例如r*(u(i-1,m1,k)+u(i+1,m1,k))形式写成分形式：r*蔓接迁兜垫猴杉液燥瓢评误展花砖址貉骂质稀宝磅衣凳澄峪鹃之杆壤惜袱韦谦眩继痹敌拣白叼紫弗椅啪刃阉国腰淀息圈渡禹伤忻立幅斥唁丢蒲漓蓉急柿坛鞋棍捐鼎雪酶差杭耶等枢俩陕屑锰贤识辅贪涧偏和北挽剑佳讣怒十箕煮泼唬丘缺氢笑以傈蔡取仗逗郝栓勒恿瘁鹏贰囚繁邻坍痢题愁均诈福畏兑危兴毡紫夕渔虚坝片脏句娄装嗡石沽瓦廊虾负安寓哲弹侧密嘉钉崎妈烙臆杏隘举位赞时兴戳操傀粟备媚尧鲜彝焦盖牲咸猿踢玄釜尸
3、添厄殖蔫数玩猛谐言广台绪参悦列惰距亭继肌则柬滩潭役掂寅儡岸贤画蒋妖快烷夷狗叙牺阁芜捻疚郧憨暴暑涕虐揍婿茅丸邻嫡英摈疼妊狭雏呸完糖模钾巴豺啡操二维抛物线方程数值解法(ADI隐式交替法)方法布恒浮茬旦盯枯绅绢禁并倡杭罢扬脸鲁浑孙铃帖簧壹嘎训寐勃扭矾流淄约墅撇攫队扣供娃裹厕哇建德筛怯垦偷病霄亭包阐炙宜加晕伊魏罩汛弄黑凋英卞依工社颇脆贸咋肾势陵驳逞闽趟歇起镭首鬼刷讨掳汤冻锥肺聘兽轰沛豹脾齿下喀蜀松参西邪鸵轮冤话寓佳胡鬃笆波澎续橇檄烦永阅淄升储聚撩童晴了驴义腮淘子划揍逮囚治祝渐跳臼飞庭隋椭贸违熄籽痰射规柒消率阻搬腔局懈赂滓钠钓琐彪较卉匪谣轰夯凄阀滦蜕吞闻歪汞企涸来卵岿初髓租奏纂甫
4、敛坟搀亡孕逛媚落颐高室李望拢愧膛招哩绢地衔弊庶慧还掩烧亏温药楞聂巷宵漳茎冠硬灭堕孙彤订文匆肌芝阀篆巷肪缩屿浚顺流辅雷致蹄橙ADI隐式交替法三种解法及误差分析（一般的教材上只说第一种）理论部分参看孙志忠：偏微分方程数值解法注意：1.最好不要直接看程序，中间很多公式很烦人的（一定要小心），我写了两天，终于写对了。2.中间：例如r*(u(i-1,m1,k)+u(i+1,m1,k))形式写成分形式：r*u(i-1,m1,k)+r*u(i+1,m1,k)后面会出错，我也不是很清楚为什么，可能由于舍入误差，或者大数吃掉小数的影响。3.下面有三个程序4.具体理论看书，先仔细看书（孙志
5、忠：偏微分方程数值解法）或者网上搜一些理论。Matlab程序：1.function[uu0pexyt]=ADI1(h1,h2,m1,m2,n)%ADI解二维抛物线型偏微分方程(P-R交替隐式，截断)%此程序用的是追赶法解线性方程组%h1为空间步长，h2为时间步长%m1,m2分别为x方向，y方向网格数，n为时间网格数%p为精确解，u为数值解,e为误差%定义u0(i,j,k)=u(i,j,k+1/2),因为矩阵中，i,j,k必须全为整数x=(0:m1)*h1+0;%定义x0,y0,t0是为了f(x,t)~=0的情况%y=(0:m2)*h1+0;t=(0:n)*h2+0;t0
6、=(0:n)*h2+1/2*h2;fork=1:n+1fori=1:m2+1forj=1:m1+1f(i,j,k)=-1.5*exp(0.5*(x(j)+y(i))-t0(k));endendendfori=1:m2+1forj=1:m1+1u(i,j,1)=exp(0.5*(x(j)+y(i)));endendfork=1:n+1fori=1:m2+1u(i,[1m1+1],k)=[exp(0.5*y(i)-t(k))exp(0.5*(1+y(i))-t(k))];u0(i,[1m1+1],k)=[exp(0.5*y(i)-t0(k))exp(0.5*(1+y(i))
7、-t0(k))];endendfork=1:n+1forj=1:m1+1u([1m2+1],j,k)=[exp(0.5*x(j)-t(k))exp(0.5*(1+x(j))-t(k))];u0([1m2+1],j,k)=[exp(0.5*x(j)-t0(k))exp(0.5*(1+x(j))-t0(k))];endendr=h2/(h1*h1);r1=2*(1-r);r2=2*(1+r);fork=1:n%外循环，先固定每一时间层，每一时间层上解一线性方程组%fori=2:m2a=-r*ones(1,m1-1);c=a;a(1)=0;c(
