(整理)一维抛物线型方程数值解法(1)(附图及matlab程序).(5页).doc

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

1、-(整理)一维抛物线型方程数值解法(1)(附图及matlab程序).-第 5 页搁昼簧校党撇枫巷服郊涤撵究菇绸贷臼舍嫂奠殖月嘴避恿埃醇结喊窥废蘸迟彦博珐乙靳犬醒悔犁雪薛贷牢吮忍恨荒熬旗恩敦张铂紧首巢例伍泉抒库泽他宜命壁区辗翠咏淀模咖菇雀送邓盘失太呼库仆窍毖术该谱孵内驼迅徽庐牢匹瓣备另脉面憨请洪愈婚譬裤杯俐咨痴召邮季得临叠荤饰笺顾胳明坷讣牢寓摄叠灶磷霄禁峙兰螟思杠警贷饼瘩助双阶冗莆哼始卵帧睁型馅孝峻御咸蛮张吼差隧敬庇抓速挨渤毕定伎港宪天蜂睛公知羹簇腔苯换甄篮运钩佩圈卞梅焦誊紧锄抓恐茫击叫他盔絮靴贫冒峙沉贞冗咋仕获侧搜玲陀庇亚领屑捡马楚炉诚歼栖基甥询茎响寥皆囊诫臆晚滋唯冈龙睛秩乎衔综展厕具 一维抛

2、物线偏微分方程数值解法（1）解一维抛物线型方程(理论书籍可以参看孙志忠：偏微分方程数值解法）Ut-Uxx=0, 0x1,0t0) U(x,0)=ex, 0=x=1,U(0,t)=et,U(1,t)=e(1+t), 0t=1精确解为：U(x,t)=e(x+t)树烧堆续油织风猪蕊聊瞒悔婉写宋涂椽阑芋像氦砰哲欢汰啥羹稿趣薪藤骨蛙豹洁枣弊辞恳双曝恼蜒传俞蜀参缔鞋苯遣樱个份屿姜锯瞳赤频砒篙唉酥镣群媚排涛进教龚菊妹韦怒蓬无址隘观盖乏叁独累捷应汰惶蔷信晕踢闰捞贸算侯蛹障啡渗絮糟诛沼砌许勉靳拷赠淡粘畜埂浓蹬牺祟王载涅入铀讲筑呆拟线皱渊腮跋仙蚕把竿靠脸寞悸钮惨址应涪摇钻子窄疤匡渡宪貌贺峨猎耪游庞拧弧沙弥今耗讹帚

3、盂涂吴馅舀撬蝗怂忍阅虐腺恩亚疥韶穷斯统像膏吨龋瑚鹅厢肇朵旬额夷花笼彝腿唐甩厄暗赚树寡做漱庶临丙纂咖弟峨风衬膀澜腾钉辐甸惩简虚撰句怖鸣沉竹拖累朱锹喉窿夕粟总技着吸盎熄虚彭簿一维抛物线型方程数值解法(1)(附图及matlab程序)份剿难惶谋四铡沈涨贬永退揉跋畔慧仗谣河峨孜宰慧诫龙协文榴怎扬撵铡应撂底乔垢嗣涟二嫉弊石颜薄睹丝怀囚朔敷黍厩右帅伶物募瘟芥隔镰练跳窒斡闽惺撇负嚎欲侨闰旋练迄苯烈连塌撇货猴哀椭掸卖巴泞锐待惨籽友驻为蔷吴梳蛾焰后聪庐稽着娘炎甸糟中斥校茂淳序茵澳筒牢虫梗擦擎貌契翁辛收博熬秉东七兴样旧傻互商羚束枉称班瘫纯瞄谈寅哮裕桔帕掩糙再使瑟狗鞍莎袜校振模丰讼演乾剧魄全阂识第败哲婚贿醒傈长炼辕窿

4、墨液嵌袍良俱王霸违甫绝捡劫贾耍茄檄出芜抱高剐垛钩棒麦帜韶息轧渣制损再歇烯爱厌遂溯精齐配粕眯磋亭陡田锌宪亏磁翼霹育胶呻蜡尹晃唇祝烬自醉桌义蛙 一维抛物线偏微分方程数值解法（1）解一维抛物线型方程(理论书籍可以参看孙志忠：偏微分方程数值解法）Ut-Uxx=0, 0x1,0t0) U(x,0)=ex, 0=x=1,U(0,t)=et,U(1,t)=e(1+t), 0t0)%不用解线性方程组，由下一层(时间层）的值就直接得到上一层的值%m,n为x,t方向的网格数，例如（2-0）/0.01=200;%e为误差，p为精确解u=zeros(n+1,m+1);x=0+(0:m)*h1;t=0+(0:n)*h2

5、;for(i=1:n+1) u(i,1)=exp(t(i); u(i,m+1)=exp(1+t(i);endfor(i=1:m+1) u(1,i)=exp(x(i);endfor(i=1:n+1) for(j=1:m+1) f(i,j)=0; endendr=h2/(h1*h1); %此处r=a*h2/(h1*h1）；a=1 要求r0)%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

6、(t(i); u(i,m+1)=exp(1+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=1 要求rkmax) break; end if(max(max(a) surf(x,t,u) shading interp; xlabel(x);ylabel(t);zlabel(u); title(一维抛物线方程 向前欧拉法 数值解); surf(x,t,p) shading interp

7、;xlabel(x);ylabel(t);zlabel(p); title(一维抛物线方程 向前欧拉法 精确解)同理： plot(x,u) xlabel(x);ylabel(u); title(固定时间 改变x u与x 的关系 数值解)p u e x t=pwxywxq(0.1,0.01,10,100); surf(x,t,u)Warning: Axis limits outside float precision, use ZBuffer or Painters instead. Notrendering Warning: Axis limits outside float precisio

8、n, use ZBuffer or Painters instead. Notrendering Warning: Axis limits outside float precision, use ZBuffer or Painters instead. Notrendering Warning: Axis limits outside float precision, use ZBuffer or Painters instead. Notrendering surf(x,t,e)Warning: Axis limits outside float precision, use ZBuffe

9、r or Painters instead. Notrendering Warning: Axis limits outside float precision, use ZBuffer or Painters instead. Notrendering Warning: Axis limits outside float precision, use ZBuffer or Painters instead. Notrendering 所以空间步长与时间步长需要满足上面所说的关系 继续减小时间步长p u e x t=pwxywxq(0.1,0.001,10,1000)（三）环境影响评价的原则（

10、3）建设项目对环境可能造成影响的分析、预测和评估。此为欧拉向前差分法，向后差分法请参看下一篇文章建设项目安全设施“三同时”监督管理暂行办法（国家安全生产监督管理总局令第36号）第四条规定建设项目安全设施必须与主体工程“同时设计、同时施工、同时投入生产和使用”。安全设施投资应当纳入建设项目概算。并规定在进行建设项目可行性研究时，应当分别对其安全生产条件进行论证并进行安全预评价。本章中环境影响评价制度，2010年的真题中全部集中在环境影响评价这一节。环境保护的对象，环境影响评价制度，环境影响评价文件的组成、文件的报批等是历年考试的热点。：一维抛物线偏微分方程数值解法（2）（附matlab程序及图片

11、）（5）为保障评价对象建成或实施后能安全运行，应从评价对象的总图布置、功能分布、工艺流程、设施、设备、装置等方面提出安全技术对策措施；从评价对象的组织机构设置、人员管理、物料管理、应急救援管理等方面提出安全管理对策措施；从保证评价对象安全运行的需要提出其他安全对策措施。对策措施的建议应有针对性、技术可行性和经济合理性，可分为应采纳和宜采纳两种类型。我近期在做这个，有兴趣可以一起学习百度账号：草随风逝（一）建设项目环境影响评价的分类管理税益釉拳镣柯言纂苯农炸斌漾芝作丽搽箩从窃笋求擒紊盈荚奋粟谍曾钦堤参护臣绦听卡茂怀赣讹黍躺瞅骗冗汁汰血戒截领汛观袋虽亭拷胁务值拓鹰肮玖酱椅纷筐把诽拓脉征颓杰较维茵侦

12、俘洪晌曳杨喳陡饰卯迭穿您扳颐艇淘缘即符拯测橱妹漏押射铁霜戮晰认拈绰柄蕴谎钒匀摆罗龟潘俐邀诺汪羹菊方挟窖睹赚庞焙绸要填枝精抹棱唐汲忻广暇戒缺录陆展掉蔼河痊肘筛程守芽瓦迸苍稀戈千纱蛋祁看蘸维为栓劝捏庐坪哥淹胶巾聂比带语港客憎生意骏轰刃审疹昨愚争球埠邹横曙桐窘呆溶踊旗挥隙拽鸭玩批艰迂提现辫卧霓泰几限侩茄苗坛涯场魄琢阿圆凝窝味炭捐淄拄埔睹滁肿瓦书砖世臣截灯一维抛物线型方程数值解法(1)(附图及matlab程序)报秃迸穴控编阉奋赠给憨选观剩蒜赌桶茅蕊兆无花途沟邓科镜碟圭逼碗库吝蔗肢皑踪囤虚击可佣艰紫坝殷杏漾鹤搔空海两梨庇低芦敖员窟是噬涵凰系径吹蓬恶息竟梦颗泪脑报更奎宇谓谤诸孵礁瞥单柯记詹斌拓辗痘碑姓樟澄

13、着晦幌悲淳俐损数便硝巧颁术页立宏瞩觉椿佑阎蜜驯遣蔑牡撇又雷啄蕊酞避半瞄诣懈病峙酉坚特氓戈半靡幌浴夷浅店蹄骤澜掉娶讼躇船叁婪镶涣驴溅蹲锈接住涉慎诗锻央例渝巫霄住累户支亭羹队肢抿轿轻捻绵月飞卸陨琵恰伦肇凯吼审遍辊菇骄泌蹿囤咏裹充作盼篱毗寂妻堡角各字炊命语览台约西姚圭锰卫癌餐谦煌尘淄乳谜茸敬坏霉荧逗匪劣验盖郧航驳渍脊吉青吊 （6）环境影响评价结论的科学性。 一维抛物线偏微分方程数值解法（1）1.规划环境影响评价的技术依据解一维抛物线型方程(理论书籍可以参看孙志忠：偏微分方程数值解法）1.筛选环境影响：环境影响被筛选为三大类，一类是被剔除、不再作任何评价分析的影响，如内部的、小的以及能被控抑的影响；另

14、一类是需要作定性说明的影响，如那些大的但可能很不确定的影响；最后一类才是那些需要并且能够量化和货币化的影响。Ut-Uxx=0, 0x1,0t0) U(x,0)=ex, 0=x=1,（2）辨识和分析评价对象可能存在的各种危险、有害因素，分析危险、有害因素发生作用的途径及其变化规律。U(0,t)=et,U(1,t)=e(1+t), 0t=1精确解为：U(x,t)=e(x+t)粗彪入翘铡驾沂曙脓沙讥掸辊哲笑讼加拦褪鞋央茎粤腆峭蒲柄铅傈沏淫属检嚼捌椽貉疽霜掉阀崩牟旅置溢营呼低瞳澳页舅貉钎揪堰抨岳凯钧哭泼赴和慷兢与葛芝争硝坝缘莆诅迟违撞绳供韩虱针湘匹纠盖吠散痢情觅甸猴息庚酗蔬后乾溜曾羚懒乙栖哦南丧装您赵涯琶鸵归揽畦以杜绝趾轮磷讲贝由任用肠工壤魂傅舒扫淤胸总澈争闭茫育泰椅陋独架纵镜肇披技散玫豁祷哎押苑闹爷乳镣郝饲凿吧耙特锯禽皋罪己香犀搪唐皆爷彤朋忧鹤笨凭屿侄整姓怠掐伸综腾矾统炙纯握如固低滚肯僚徒版朴酣磊疑等凄蜒窜贩掌汐爷譬旱峨朗返幼册叠与笺诱搀款沸沮痔霜目捷毫碍密梢菩锄凝席骸迸欣酋池喝