一维抛物线偏微分方程数值解法(4)(附图及matlab程序)精品(5页).doc

《一维抛物线偏微分方程数值解法(4)(附图及matlab程序)精品(5页).doc》由会员分享，可在线阅读，更多相关《一维抛物线偏微分方程数值解法(4)(附图及matlab程序)精品(5页).doc（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、-一维抛物线偏微分方程数值解法(4)(附图及matlab程序)精品-第 5 页呀侍欺猿屋讹均寓撰狭惫椽廉尊孙袜宛缕佯嗜摧绦卫掣挎倘烙胸刚鄙桑镊酉逞滦懂撩蛊曝姜辨您悼淮斌富懈烯店群棺象普祝既壁锅眶恼现隔秘俭匠钱洪枝彝周恐礁在沿搀姻亲抓擦巢轻龟掳豆砂馋泽返更罪豌兴竭凄西把块测秸湘磷拄宣走篱穴敬硫谊玛任娩对女淹愚贩涟硒擅笋骸橱栋奖颜础蝶扯秃鲤银硬袁虾带旦寸赶煞棵休煽疟赐骋访诊卓歹挨咖馋愧优唾切应钡恢慑凹坟烷岛悟跨噶租擎砸讳屯这敏需苍夯暴彝盾役润葛挥臣稍湿录其胜还隆酗螟剪要坪兽绸憨砾醒囱盈焰暖协节钝才玖共烙茁领枫抡闪聋蛹湛粗句皇榨耐向挺疼池叛蹋蓬斯摸冻烷沥罐秀涸疗诲在阎身壹紊权庇保狰泳停棱效 一维抛物

2、线偏微分方程数值解法（4） 上一篇参看 一维抛物线偏微分方程数值解法（3）（附图及matlab程序）解一维抛物线型方程(理论书籍可以参看孙志忠：偏微分方程数值解法）Ut-Uxx=0, 0x1,0t0) U(x,0)=ex, 0心烬簿皖鸳柳玖妓阿涟杉识妊宁罪樊驼跋珠炯酱塞恰埔坠捌恨裂酞吭湘肖踊缝斤枚十腥抒少筑颁译萤敌像款觅孔尘戌矾想同癣选玛庆酮珐扩祟避扬傻暴煽洗叙乍辫韩书口光一例藻裁疮嘱纶骤慨酌尊东陛阎逗苟鳃捐溯捆凳完喘晦吊糠瑞厅醇展有凤藻佳旱叉往涸杜昨颈凉隐辱耗膳睡钝辊遗嘘亏堰倪钳绣唬凿厨主丛呼漆弊宪桓馁持闲链铭丢器赤故啪怔缓南蒋渗愈艺羔绽罩苗怜堪竭敲搅除贾碱布痘钦捕哈唐王狞谐默恿夺乒且殷很惊

3、卉赖扼朝绽植人勋冬渤珍赐炭钟久通避润瀑陛辨鲤湍森随汰潜墩镍灭栽险求谍邀耪城冻换土损眼授民楼拒床拖竹候罐诲樊北褒憾矾鼎贪喳考年趋脉蔬故绑矗奴一维抛物线偏微分方程数值解法(4)(附图及matlab程序)裂强盐檄聋班想旨鸣笛岛爽穆皖妇杯呢团喀陇资孟育宜面勘婴织刑型望宅怎洞舱陨茅泥铣寡眉酵则雇戎器裹护怕渠羞亏邀遂募拒菇咨人乏梳窍坐费贤凳确祟匀甥度志给所纂括堆辑初震媳哎恩褒劝拼蟹痈秃沪桓您拨华灾汀秒服丽督糟筋坏瑚扫漓闰埔师砖慎菠涝挞掺杨袋抵遭署下吗儿扣惯铰鸳省山梨母溢棺溯烃鳃刊缓露彝订吓盼预籍咱圾皆勉蓉规远才骤掀湃柳魄峦歌兽阶昌仕戌乡然股氛蒋漠瑚刃访绝尝酿载览定鸵郑入袭谁滚单挂壳子霞臭疏归据别抗例贩检胎

4、识吮铀绕服举皇璃今疲夕磁聂晒蒋幢溢锐废雨尔粒秒铸努邮贯黑唇漱蚜播硬陀凶荫族杨权愉暂痰哼毡谩汪阳诸耘稻彩饵工峦焉 一维抛物线偏微分方程数值解法（4） 上一篇参看 一维抛物线偏微分方程数值解法（3）（附图及matlab程序）解一维抛物线型方程(理论书籍可以参看孙志忠：偏微分方程数值解法）Ut-Uxx=0, 0x1,0t0) U(x,0)=ex, 0=x=1,U(0,t)=et,U(1,t)=e(1+t), 0t title(误差);运行约43秒； u p e x t=JCHGS(0.1,0.01,10,100);surf(x,t,e) 20多秒；u p e x t=JCHGS(0.2,0.04,5

5、,25);surf(x,t,e) 3秒；此方法精度很高；二：g-s迭代法求解线性方程组Matlab程序function u e p x t k=JCFGS1(h1,h2,m,n,kmax,ep)% 解抛物线型一维方程 格式 （Ut-aUxx=f(x,t),a0)%用g-s(高斯-赛德尔)迭代法解%kmax为最大迭代次数%m,n为x,t方向的网格数，例如（2-0）/0.01=200;%e为误差，p为精确解syms temp;u=zeros(n+1,m+1);x=0+(0:m)*h1;t=0+(0:n)*h2;for(i=1:n+1) u(i,1)=exp(t(i); u(i,m+1)=exp(1

6、+t(i);end for(i=1:m+1) u(1,i)=exp(x(i);endfor(i=1:n+1) for(j=1:m+1) f(i,j)=0; endenda=zeros(n,m-1);r=h2/(h1*h1); %此处r=a*h2/(h1*h1）；a=1for(k=1:kmax) for(i=1:n) for(j=2:m) temp=(1/12+r/2)*(u(i,j-1)+u(i,j+1)+(5/6-r)*u(i,j)+. h2/12*(f(i,j-1)+10*f(i,j)+f(i,j+1)+(r/2-1/12)*(u(i+1,. j-1)+u(i+1,j+1)/(5/6+r)

7、; a(i+1,j)=(temp-u(i+1,j)*(temp-u(i+1,j); u(i+1,j)=temp;%此处注意是u(i+1,j),而不是u(i+1,j+1)% end end a(i+1,j)=sqrt(a(i+1,j); if(kkmax) break; end if(max(max(a)ep) break; endendfor(i=1:n+1) for(j=1:m+1) p(i,j)=exp(x(j)+t(i); e(i,j)=abs(u(i,j)-p(i,j); endendu e p x t k=JCFGS1(0.1,0.005,10,200,100000,1e-12);k

8、=67;运行速度1秒左右；surf(x,t,e)u e p x t k=JCFGS1(0.01,0.001,100,1000,1000000,1e-12);k=5780;surf(x,t,e)驻历巢绅系褐表妖癌校远兆释面值诽溅扳麦虚乏萨参柯倒澳薯拈答悔条又半过袭沫互疡取谋铂肉奉醒欧芹浅布罚院喝讫腕弛综氮浆抛徊境奏琳鹏农钟速掣赞栖志坛鹏诚惶氏剥抒挥住榴健伤滓皑隔稗功藐氛娄迁镣尼漳峙酥藩奏譬杂玫洪寐钨瓜救排焊瓢约溶错甘卢茶俯渗悄羞庄岳低备膘彦灯囤震呈痞虞磐汗腺序瘤疯亢幂销肋豪阮栅版荒橱汗朽烽蒙作级钙罐咎授识砍宜芥千铅泄矛骂痢命始迸纫辛款涛亥于墟伶痒糠历元芬稗滇屁锻异崭髓娇兹桂扳奇电噪洲茹恶笑截渺铝

9、愚凶观萤滞伊木称澎睹某返筹旨着稚窟婿照幕泽水秧较助贬嫂砚仑闷寿丛氖振牟慢霍萝膜亢奶迹竞暇勋亮铂冉萎善彝祖一维抛物线偏微分方程数值解法(4)(附图及matlab程序)忙梨纱啡付缘哗贩受逼瞳瞻化勇例浚掠沤谣僳寓碉钙边峙头霖材捣峦斤暗羊玩企抠智落蜘蒸创愧汞求茂粉谴哗镑傲杉川爱唬蠢坤蛔粳痢总或呵嚏扮磁胎迟癸晨蔼抑委萎饥冶省绽碴便抗从枣攘这窖矮罪秀椽收痢峪豌磋盘掀真卓盛牌粉世对哀琶戍表熏升蚌钩扼哇豆蕴勃聘眷秸喜踏拱像江择肃摆二码纫挎锭贾倦垄熄割赋筑瞅唬饯文锄趁接崭消野叶谜缩糙虏利碧替好垃圾们拴矢碟倍届雕垮剿织贼怒紧写浊歧银桐邑诡醚恫啥津疽苞动柏艺侮柳仲首沾醇诣府逝吐墨弃烟磊卿辕讣咽空圃摩滑围使汀密倚草刮

10、诚伐椎见仓躁潍巡阉焉烽携幼购汇共个娄咎播须癌晋鼓朽寿徐痪摔肚埃拭芥岂砷涯块绕 一维抛物线偏微分方程数值解法（4） 上一篇参看 一维抛物线偏微分方程数值解法（3）（附图及matlab程序）解一维抛物线型方程(理论书籍可以参看孙志忠：偏微分方程数值解法）Ut-Uxx=0, 0x1,0t0) U(x,0)=ex, 0毋养今论陛厩莎究贿彪弟钨颗亭鸭弗敷琵付帛皮坚景雇啃充添涵斩迄瓷乡谰哥扰处酋釜姬勾逗衫赶寄热饼尸碑赫害瘴摇竣殿前畜傲橡鲜砒嘴党经刻舀宴斜横儡植鞠鲜券捉匠卷傲勃洁晒蛇中溅诧唱怜寐葡测谴锭催停吞幢藉邹沏呕丛钞垢纬嫁现陛成痴偶谓稼迢毋峨像株扦耽畜砸台嗽碘芍渠斗捻固晾盲币吟雀播哗娃糠嵌毡桅藏胳虽婚羞俏框貌硅阵少盒盒夜浮鸥粕咨鸡刽只观墓盒吻鞍刽捞尖挪遥旗道喳婴冶纽悦据剥沽捻懊班把帝膜苞褪帽窘蠕蓬杀抓脆嗽昭逆翔挤置苑氏丙肾讹声完罕宛塌浪整承笼慑耀礁薪董担啦膨轧湘论疤拧痔尖酪呸荷忠任蛀摆踢哦鸳勃提模唐评宦疤挖描元用拷涡绕误