| 查看: 2585 | 回復(fù): 13 | ||
value_1987新蟲 (初入文壇)
|
[求助]
Matlab用龍格庫塔法求二階非線性方程
|
|
|
求助解上式方程,用龍格庫塔方法,用[t,y]=ode45(‘fname’,tspan,y0),如果u(t)是常數(shù)或正弦變化等,我會(huì)編程序,但是u(t)是一個(gè)隨機(jī)的數(shù),與t有關(guān)的,我編了個(gè)小程序來表示隨機(jī)數(shù),但是不知道怎么放到M文件里面。 |
新蟲 (初入文壇)
金蟲 (著名寫手)
新蟲 (初入文壇)
|
function [y,z]=Runge_kutta(a,b,y0,z0,h) x=a:h:b; y(1)=y0; z(1)=z0; n=(b-a)/h+1; for i=2:n K(1,1)=f1(x(i-1),y(i-1),z(i-1)); K(2,1)=f2(x(i-1),y(i-1),z(i-1)); K(1,2)=f1(x(i-1)+h/2,y(i-1)+K(1,1)*h/2,z(i-1)+K(2,1)*h/2); K(2,2)=f2(x(i-1)+h/2,y(i-1)+K(1,1)*h/2,z(i-1)+K(2,1)*h/2); K(1,3)=f1(x(i-1)+h/2,y(i-1)+K(1,2)*h/2,z(i-1)+K(2,2)*h/2); K(2,3)=f2(x(i-1)+h/2,y(i-1)+K(1,2)*h/2,z(i-1)+K(2,2)*h/2); K(1,4)=f1(x(i-1)+h,y(i-1)+K(1,3)*h,z(i-1)+K(2,3)*h); K(2,4)=f2(x(i-1)+h,y(i-1)+K(1,3)*h,z(i-1)+K(2,3)*h); y(i)=y(i-1)+h/6*(K(1,1)+2*K(1,2)+2*K(1,3)+K(1,4)); z(i)=z(i-1)+h/6*(K(2,1)+2*K(2,2)+2*K(2,3)+K(2,4)); end y(2); plot(y,'r') 自己編寫了隨機(jī)數(shù)的M文件,但是每個(gè)K都是調(diào)用一次M文件,就會(huì)產(chǎn)生不同隨機(jī)數(shù),應(yīng)該要產(chǎn)生相同隨機(jī)數(shù)才行啊 |
木蟲 (小有名氣)
新蟲 (初入文壇)
新蟲 (初入文壇)
鐵桿木蟲 (著名寫手)
方丈大師
新蟲 (初入文壇)
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 309求調(diào)劑 | 誰不是少年 2026-03-29 | 17/850 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +13 | minht 2026-03-29 | 13/650 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +10 | 26要上岸 2026-03-28 | 10/500 |
|
|
[考研] 085601 329分調(diào)劑 +6 | yzsa12 2026-03-31 | 6/300 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +12 | 材料學(xué)257求調(diào)劑 2026-03-28 | 13/650 |
|
|
[考研] 一志愿中國科學(xué)院大學(xué)265求調(diào)劑 +6 | 恬淡ye 2026-03-31 | 7/350 |
|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
[考研] 吉大生物學(xué)326分求調(diào)劑 +3 | sunnyupup 2026-03-31 | 3/150 |
|
|
[考研] 一志愿西交288化工專碩求調(diào)劑 +4 | 好運(yùn)好運(yùn)接接 2026-03-24 | 4/200 |
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +5 | pipiver 2026-03-30 | 5/250 |
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|