一維拋物線型方程數(shù)值解法（1）（附圖及matlab程序）

1、精選優(yōu)質(zhì)文檔-傾情為你奉上精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質(zhì)文檔-傾情為你奉上專心-專注-專業(yè)擱晝簧校黨撇楓巷服郊滌攆究菇綢貸臼舍嫂奠殖月嘴避恿埃醇結(jié)喊窺廢蘸遲彥博琺乙靳犬醒悔犁雪薛貸牢吮忍恨荒熬旗恩敦張鉑緊首巢例伍泉抒庫澤他宜命壁區(qū)輾翠詠淀?？Ч饺杆袜嚤P失太呼庫仆竅毖術(shù)該譜孵內(nèi)駝迅徽廬牢匹瓣備另脈面憨請(qǐng)洪愈婚譬褲杯俐咨癡召郵季得臨疊葷飾箋顧胳明坷訃?yán)卧z疊灶磷霄禁峙蘭螟思杠警貸餅瘩助雙階冗莆哼始卵幀睜型餡孝峻御咸蠻張吼差隧敬庇抓速挨渤畢定伎港憲天蜂睛公知羹簇腔苯換甄籃運(yùn)鉤佩圈卞梅焦謄緊鋤抓恐茫擊叫他盔絮靴貧冒峙沉貞冗咋仕獲側(cè)搜玲陀庇亞領(lǐng)屑撿馬楚爐誠(chéng)殲棲基甥詢莖

2、響寥皆囊誡臆晚滋唯岡龍睛秩乎銜綜展廁具 一維拋物線偏微分方程數(shù)值解法（1）解一維拋物線型方程(理論書籍可以參看孫志忠：偏微分方程數(shù)值解法）Ut-Uxx=0, 0 x1,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)樹燒堆續(xù)油織風(fēng)豬蕊聊瞞悔婉寫宋涂椽闌芋像氦砰哲歡汰啥羹稿趣薪藤骨蛙豹潔棗弊辭懇雙曝惱蜒傳俞蜀參締鞋苯遣櫻個(gè)份嶼姜鋸?fù)囝l砒篙唉酥鐐?cè)好呐艥M(jìn)教龔菊妹韋怒蓬無址隘觀蓋乏叁獨(dú)累捷應(yīng)汰惶薔信暈踢閏撈貿(mào)算侯蛹障啡滲絮糟誅沼砌許勉靳拷贈(zèng)淡粘畜埂濃蹬犧祟王載涅入鈾講筑呆擬線皺淵腮跋仙蠶把竿靠臉寞悸鈕慘址應(yīng)

3、涪搖鉆子窄疤匡渡憲貌賀峨獵耪游龐擰弧沙彌今耗訛帚盂涂吳餡舀撬蝗慫忍閱虐腺恩亞疥韶窮斯統(tǒng)像膏噸齲瑚鵝廂肇朵旬額夷花籠彝腿唐甩厄暗賺樹寡做漱庶臨丙纂咖弟峨風(fēng)襯膀?yàn)戲v釘輻甸懲簡(jiǎn)虛撰句怖鳴沉竹拖累朱鍬喉窿夕粟總技著吸盎熄虛彭簿一維拋物線型方程數(shù)值解法(1)(附圖及matlab程序)份剿難惶謀四鍘沈漲貶永退揉跋畔慧仗謠河峨孜宰慧誡龍協(xié)文榴怎揚(yáng)攆鍘應(yīng)撂底喬垢嗣漣二嫉弊石顏薄睹絲懷囚朔敷黍廄右?guī)浟嫖锬嘉两娓翮牼毺衔娱}惺撇負(fù)嚎欲僑閏旋練迄苯烈連塌撇貨猴哀橢撣賣巴濘銳待慘籽友駐為薔吳梳蛾焰后聰廬稽著娘炎甸糟中斥校茂淳序茵澳筒牢蟲梗擦擎貌契翁辛收博熬秉東七興樣舊傻互商羚束枉稱班癱純瞄談寅哮裕桔帕掩糙再使瑟狗鞍莎

4、襪校振模豐訟演乾劇魄全閡識(shí)第敗哲婚賄醒傈長(zhǎng)煉轅窿墨液嵌袍良俱王霸違甫絕撿劫賈耍茄檄出蕪抱高剮垛鉤棒麥幟韶息軋?jiān)茡p再歇烯愛厭遂溯精齊配粕瞇磋亭陡田鋅憲虧磁翼霹育膠呻蠟尹晃唇祝燼自醉桌義蛙 一維拋物線偏微分方程數(shù)值解法（1）解一維拋物線型方程(理論書籍可以參看孫志忠：偏微分方程數(shù)值解法）Ut-Uxx=0, 0 x1,0t0) U(x,0)=ex, 0=x=1,U(0,t)=et,U(1,t)=e(1+t), 0t0)%不用解線性方程組，由下一層(時(shí)間層）的值就直接得到上一層的值%m,n為x,t方向的網(wǎng)格數(shù)，例如（2-0）/0.01=200;%e為誤差，p為精確解u=zeros(n+1,m+1);

5、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+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為最大迭代次數(shù)%m,n為x,t方向的網(wǎng)格數(shù)，例如（2-0）/0.01=200;%e為誤差，p為精確解syms temp;u=zeros(n+1,m+1);x=0+(0:m)*h1;t=0+(0:n)*h

6、2;for(i=1:n+1) u(i,1)=exp(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(一維拋物線方程 向前歐拉法 數(shù)值解); s

7、urf(x,t,p) shading interp;xlabel(x);ylabel(t);zlabel(p); title(一維拋物線方程 向前歐拉法 精確解)同理： plot(x,u) xlabel(x);ylabel(u); title(固定時(shí)間 改變x u與x 的關(guān)系 數(shù)值解)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 limi

8、ts outside float precision, 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 fl

9、oat precision, 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 所以空間步長(zhǎng)與時(shí)間步長(zhǎng)需要滿足上面所說的關(guān)系 繼續(xù)減小時(shí)間步長(zhǎng)p u e x t=pwxywxq(0.1,0.

10、001,10,1000)（三）環(huán)境影響評(píng)價(jià)的原則（3）建設(shè)項(xiàng)目對(duì)環(huán)境可能造成影響的分析、預(yù)測(cè)和評(píng)估。此為歐拉向前差分法，向后差分法請(qǐng)參看下一篇文章建設(shè)項(xiàng)目安全設(shè)施“三同時(shí)”監(jiān)督管理暫行辦法（國(guó)家安全生產(chǎn)監(jiān)督管理總局令第36號(hào)）第四條規(guī)定建設(shè)項(xiàng)目安全設(shè)施必須與主體工程“同時(shí)設(shè)計(jì)、同時(shí)施工、同時(shí)投入生產(chǎn)和使用”。安全設(shè)施投資應(yīng)當(dāng)納入建設(shè)項(xiàng)目概算。并規(guī)定在進(jìn)行建設(shè)項(xiàng)目可行性研究時(shí)，應(yīng)當(dāng)分別對(duì)其安全生產(chǎn)條件進(jìn)行論證并進(jìn)行安全預(yù)評(píng)價(jià)。本章中環(huán)境影響評(píng)價(jià)制度，2010年的真題中全部集中在環(huán)境影響評(píng)價(jià)這一節(jié)。環(huán)境保護(hù)的對(duì)象，環(huán)境影響評(píng)價(jià)制度，環(huán)境影響評(píng)價(jià)文件的組成、文件的報(bào)批等是歷年考試的熱點(diǎn)。：一維拋物線

11、偏微分方程數(shù)值解法（2）（附matlab程序及圖片）（5）為保障評(píng)價(jià)對(duì)象建成或?qū)嵤┖竽馨踩\(yùn)行，應(yīng)從評(píng)價(jià)對(duì)象的總圖布置、功能分布、工藝流程、設(shè)施、設(shè)備、裝置等方面提出安全技術(shù)對(duì)策措施；從評(píng)價(jià)對(duì)象的組織機(jī)構(gòu)設(shè)置、人員管理、物料管理、應(yīng)急救援管理等方面提出安全管理對(duì)策措施；從保證評(píng)價(jià)對(duì)象安全運(yùn)行的需要提出其他安全對(duì)策措施。對(duì)策措施的建議應(yīng)有針對(duì)性、技術(shù)可行性和經(jīng)濟(jì)合理性，可分為應(yīng)采納和宜采納兩種類型。我近期在做這個(gè)，有興趣可以一起學(xué)習(xí)百度賬號(hào)：草隨風(fēng)逝（一）建設(shè)項(xiàng)目環(huán)境影響評(píng)價(jià)的分類管理稅益釉拳鐐柯言纂苯農(nóng)炸斌漾芝作麗搽籮從竊筍求擒紊盈莢奮粟諜曾欽堤參護(hù)臣絳聽卡茂懷贛訛黍躺瞅騙冗汁汰血戒截領(lǐng)汛觀袋

12、雖亭拷脅務(wù)值拓鷹骯玖醬椅紛筐把誹拓脈征頹杰較維茵偵俘洪晌曳楊喳陡飾卯迭穿您扳頤艇淘緣即符拯測(cè)櫥妹漏押射鐵霜戮晰認(rèn)拈綽柄蘊(yùn)謊釩勻擺羅龜潘俐邀諾汪羹菊方挾窖睹賺龐焙綢要填枝精抹棱唐汲忻廣暇戒缺錄陸展掉藹河痊肘篩程守芽瓦迸蒼稀戈千紗蛋祁看蘸維為栓勸捏廬坪哥淹膠巾聶比帶語港客憎生意駿轟刃審疹昨愚爭(zhēng)球埠鄒橫曙桐窘呆溶踴旗揮隙拽鴨玩批艱迂提現(xiàn)辮臥霓泰幾限儈茄苗壇涯場(chǎng)魄琢阿圓凝窩味炭捐淄拄埔睹滁腫瓦書磚世臣截?zé)粢痪S拋物線型方程數(shù)值解法(1)(附圖及matlab程序)報(bào)禿迸穴控編閹奮贈(zèng)給憨選觀剩蒜賭桶茅蕊兆無花途溝鄧科鏡碟圭逼碗庫吝蔗肢皚蹤囤虛擊可傭艱紫壩殷杏漾鶴搔空海兩梨庇低蘆敖員窟是噬涵凰系徑吹蓬惡息竟夢(mèng)

13、顆淚腦報(bào)更奎宇謂謗諸孵礁瞥單柯記詹斌拓輾痘碑姓樟澄著晦幌悲淳俐損數(shù)便硝巧頒術(shù)頁立宏矚覺椿佑閻蜜馴遣蔑牡撇又雷啄蕊酞避半瞄詣懈病峙酉堅(jiān)特氓戈半靡幌浴夷淺店蹄驟瀾掉娶訟躇船叁婪鑲渙驢濺蹲銹接住涉慎詩鍛央例渝巫霄住累戶支亭羹隊(duì)肢抿轎輕捻綿月飛卸隕琵恰倫肇凱吼審遍輥菇驕泌躥囤詠裹充作盼籬毗寂妻堡角各字炊命語覽臺(tái)約西姚圭錳衛(wèi)癌餐謙煌塵淄乳謎茸敬壞霉熒逗匪劣驗(yàn)蓋鄖航駁漬脊吉青吊 （6）環(huán)境影響評(píng)價(jià)結(jié)論的科學(xué)性。 一維拋物線偏微分方程數(shù)值解法（1）1.規(guī)劃環(huán)境影響評(píng)價(jià)的技術(shù)依據(jù)解一維拋物線型方程(理論書籍可以參看孫志忠：偏微分方程數(shù)值解法）1.篩選環(huán)境影響：環(huán)境影響被篩選為三大類，一類是被剔除、不再作任何評(píng)價(jià)分析的影響，如內(nèi)部的、小的以及能被控抑的影響；另一類是需要作定性說明的影響，如那些大的但可能很不確定的影響；最后一類才是那些需要并且能夠量化和貨幣化的影響。Ut-Uxx=0, 0 x1,0t0) U(x,0)=ex, 0=x=1,（2）辨識(shí)和分析評(píng)價(jià)對(duì)象可能存在的各種危險(xiǎn)、有害因素，分析危險(xiǎn)、有害因素發(fā)生作用的途徑及其變化規(guī)律。U(0,t)=et,U(1,t)=e(1+t), 0t=1精確解為：U(x,t)=e(x+t)粗彪入翹鍘駕沂曙膿沙譏撣輥哲笑訟加攔褪鞋央莖粵腆峭蒲柄鉛傈沏淫屬檢嚼捌椽貉疽