matlab教程【6】微分方程-資料下載頁
2025-05-10 18:39本頁面
【正文】 腦。平時，他只需要將溫度、濕度、壓力等氣象數(shù)據(jù)輸入，電腦就會依據(jù)三個內(nèi)建的微分方程式，計算出下一刻可能的氣象數(shù)據(jù)，因此模擬出氣象變化圖。 這一天， Lorenz想更進一步了解某段紀錄的後續(xù)變化，他把某時刻的氣象數(shù)據(jù)重新輸入電腦，讓電腦計算出更多的後續(xù)結(jié)果。當時，電腦處理數(shù)據(jù)資料的數(shù)度不快，在結(jié)果出來之前，足夠他喝杯咖啡并和友人閑聊一陣。在一小時後，結(jié)果出來了，不過令他目瞪口呆。結(jié)果和原資訊兩相比較，初期數(shù)據(jù)還差不多，越到後期，數(shù)據(jù)差異就越大了，就像是不同的兩筆資訊。而問題并不出在電腦，問題是他輸入的數(shù)據(jù)差了 ，而這些微的差異卻造成天壤之別。所以長期的準確預測天氣是不可能的。 例：求解 lorlenz方程 例：求解 lorlenz方程 zyxydtdzzydtdyyzxdtdx??????????35101038zyxydtdyzydtdyyzxdtdy??????????35( 3)1010( 2)38( 1)設 y(1)=x, y(2)=y, y(3)=z %%program axis([10 50 50 50 50 50])。 %設置坐標軸范圍 view(3) %設置觀察三維圖形的視角 hold on title(39。Lonrenz Attractor39。) %圖的標題 options =odeset(39。OutputF39。,39。odephas339。) 。 %設置輸出方式為三維空間三變量圖形 [t,y]=ode23(39。lorfun39。,[0,20],[0,0,eps],options)。 %%odefile function ydot=lorfun(t,y) ydot=[8/3*y(1)+y(2)*y(3)。 10*y(2)+10*y(3)。 y(2)*y(1)+35*y(2)y(3)]。 作業(yè)： [1,3]區(qū)間內(nèi)的數(shù)值解，并將結(jié)果與解析解（ y=x2+(1/x2))進行比較。 2( 1)42????yxxydxdy605005222??????)()(39。yyxdxdyxdx[0,2]的解。 1和例 2。
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號-1