 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (3346)
- >導(dǎo)師招生 (276)
- >基金申請(qǐng) (241)
- >考研 (233)
- >文獻(xiàn)求助 (127)
- >蟲(chóng)友互識(shí) (112)
- >考博 (101)
- >休閑灌水 (45)
- >招聘信息布告欄 (27)
- >論文投稿 (26)
- >論文道賀祈福 (25)
- >博后之家 (19)
- >碩博家園 (17)
- >找工作 (13)
- >教師之家 (11)
- >SciFinder/Reaxys (10)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 671 | 回復(fù): 4 | |||
小魚(yú)118銀蟲(chóng) (小有名氣)
|
[求助]
MATLAB程序調(diào)整
|
|
function gaosi() i=1; for N=1:0.1:6 [T,Y]=ode45(@last00,[0.001,1000],[0.0000001,0.0001]) JI(i)=max(Y(:,1)) JIm(i)=max(Y(:,2)); NN(i)=N; i=i+1; end plot(NN,JI) xlabel('N'); ylabel('Φ_{int}'); hold on; end function dy=last00(t,y) dy=zeros(2,1); dy(1)=dy(2); dy(2)=y(2)*(1/60)*(1+4*(exp(-y(1)*0.05^(1/2))+exp(-y(1)*0.15^(1/2))+exp(-y(1)*0.25^(1/2))+... exp(-y(1)*0.35^(1/2))+exp(-y(1)*0.45^(1/2))+exp(-y(1)*0.55^(1/2))+exp(-y(1)*0.65^(1/2))+... exp(-y(1)*0.75^(1/2))+exp(-y(1)*0.85^(1/2))+exp(-y(1)*0.95^(1/2)))+2*(exp(-y(1)*0.1^(1/2))+... exp(-y(1)*0.2^(1/2))+exp(-y(1)*0.3^(1/2))+exp(-y(1)*0.4^(1/2))+exp(-y(1)*0.5^(1/2))+... exp(-y(1)*0.6^(1/2))+exp(-y(1)*0.7^(1/2))+exp(-y(1)*0.8^(1/2))+exp(-y(1)*0.9^(1/2)))+exp(-y(1)))-... 0.807*log(1/(N*N))/(0.7846+log(1/(N*N)))*y(2)*(1-(1-exp(2.6*y(1)))/(2.6*y(1)))-0.4117*0.1238/(0.7846+log(1/(N*N)))*y(2)-0.5891*0.1238/(0.7846+log(1/(N*N)))*y(2)*((1-exp(-30*y(1)))/(30*y(1)))-... (1-0.1238/(0.7846+log(1/(N*N))))*y(2)-5.8e-16/(8*0.15)*7.444e15*0.1238/(0.7846+log(1/(N*N)))*30*y(2)*y(2) end 運(yùn)行結(jié)果為 ??? Undefined function or variable 'N'. Error in ==> gaosi>last00 at 18 dy(2)=y(2)*(1/60)*(1+4*(exp(-y(1)*0.05^(1/2))+exp(-y(1)*0.15^(1/2))+exp(-y(1)*0.25^(1/2))+... Error in ==> funfun\private\odearguments at 110 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ==> ode45 at 173 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in ==> gaosi at 4 [T,Y]=ode45(@last00,[0.001,1000],[0.0000001,0.0001]) 求高手指點(diǎn) |
回復(fù)此樓» 猜你喜歡
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有179人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有30人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有58人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助：如何用matlab實(shí)現(xiàn)一個(gè)沿x軸運(yùn)動(dòng)的且半徑隨x坐標(biāo)變化的圓的軌跡圖
已經(jīng)有4人回復(fù)
- 如何在matlab中計(jì)算對(duì)稱矩陣的高精度本征函數(shù)？
已經(jīng)有0人回復(fù)
- MATLAB程序總是出錯(cuò)?
已經(jīng)有1人回復(fù)
- MATLAB問(wèn)題 程序怎么會(huì)出錯(cuò)？Error: Missing variable or function. 怎么修改？
已經(jīng)有1人回復(fù)
- matlab作出的圖都縮到一塊，如何調(diào)整？
已經(jīng)有7人回復(fù)
- MATLAB程序設(shè)計(jì)與工程應(yīng)用
已經(jīng)有144人回復(fù)
- matlab提前設(shè)定語(yǔ)句運(yùn)行時(shí)間的問(wèn)題
已經(jīng)有4人回復(fù)
- 求大神幫忙解題啊�。�！
已經(jīng)有4人回復(fù)
- 絕熱毛細(xì)管（制冷）校核計(jì)算時(shí)壅塞流流量的計(jì)算
已經(jīng)有17人回復(fù)
- MATLAB程序調(diào)整
已經(jīng)有4人回復(fù)
- 關(guān)于使用lsqcurvefit進(jìn)行最小二乘曲線擬合時(shí)候出現(xiàn)的誤差巨大問(wèn)題求助！��！
已經(jīng)有7人回復(fù)
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對(duì)論效應(yīng)的影響
已經(jīng)有17人回復(fù)
- 分享matlab的一些代碼（入門貼）
已經(jīng)有19人回復(fù)
- 求高手指教用matlab擬合相變曲線，為什么不管怎么調(diào)整參數(shù)誤差都這么大，有圖
已經(jīng)有2人回復(fù)
- 一種簡(jiǎn)單已用的熱加工圖畫(huà)法
已經(jīng)有22人回復(fù)
- 用MATLAB重新演繹PSCAD圖線的方法
已經(jīng)有6人回復(fù)
- 【求助】怎么使用MATLAB求PID與響應(yīng)的調(diào)整時(shí)間關(guān)系
已經(jīng)有2人回復(fù)
- 【beefly個(gè)人文集】畫(huà)NBO軌道的最簡(jiǎn)單、最省錢方法
已經(jīng)有65人回復(fù)
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
dbb627: 金幣+1, 歡迎交流 2012-11-12 08:56:50
感謝參與，應(yīng)助指數(shù) +1
dbb627: 金幣+1, 歡迎交流 2012-11-12 08:56:50
|
你的函數(shù)表達(dá)式dy=last00(t,y)里面需要把參數(shù)加進(jìn)去，dy=last00(t,y，N)，然后在ode45(@last00,[0.001,1000],[0.0000001,0.0001]) 里面把參數(shù)值傳遞給表達(dá)式，這個(gè)你可以看看ODE45的help文件�，F(xiàn)在的這個(gè)函數(shù)沒(méi)法傳遞函數(shù)的參數(shù)到ode45. |
2樓2012-11-10 09:12:24
小魚(yú)118
銀蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 5821.5
- 散金: 51
- 紅花: 4
- 帖子: 206
- 在線: 161.5小時(shí)
- 蟲(chóng)號(hào): 1134404
- 注冊(cè): 2010-10-29
- 性別: MM
- 專業(yè): 光學(xué)
3樓2012-11-10 09:28:57
4樓2012-11-10 10:56:31
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★
小魚(yú)118: 金幣+5, ★★★★★最佳答案, 非常感謝，我同學(xué)成功了！ 2012-11-10 13:56:25
dbb627: 金幣+1, 感謝應(yīng)助 2012-11-12 08:57:21
小魚(yú)118: 金幣+5, ★★★★★最佳答案, 非常感謝，我同學(xué)成功了！ 2012-11-10 13:56:25
dbb627: 金幣+1, 感謝應(yīng)助 2012-11-12 08:57:21
|
程序如下，,[0.001,1000],[0.0000001,0.0001]這兩個(gè)參數(shù)可能不是很合理，可以自己調(diào)整一下。 function gaosi() NN=1:0.1:6; JI=zeros(size(NN)); JIm=zeros(size(NN)); for i=1:length(NN) N=NN(i); [T,Y]=ode45(@last00,[0.01,100],[0.001,0.01],[],N); JI(i)=max(Y(:,1)); JIm(i)=max(Y(:,2)); end plot(NN,JI) xlabel('NN'); ylabel('Φ_{int}'); hold on; function dy=last00(t,y,N) dy=[y(2);y(2)*(1/60)*(1+4*(exp(-y(1)*0.05^(1/2))+exp(-y(1)*0.15^(1/2))+exp(-y(1)*0.25^(1/2))+... exp(-y(1)*0.35^(1/2))+exp(-y(1)*0.45^(1/2))+exp(-y(1)*0.55^(1/2))+exp(-y(1)*0.65^(1/2))+... exp(-y(1)*0.75^(1/2))+exp(-y(1)*0.85^(1/2))+exp(-y(1)*0.95^(1/2)))+2*(exp(-y(1)*0.1^(1/2))+... exp(-y(1)*0.2^(1/2))+exp(-y(1)*0.3^(1/2))+exp(-y(1)*0.4^(1/2))+exp(-y(1)*0.5^(1/2))+... exp(-y(1)*0.6^(1/2))+exp(-y(1)*0.7^(1/2))+exp(-y(1)*0.8^(1/2))+exp(-y(1)*0.9^(1/2)))+exp(-y(1)))-... 0.807*log(1/(N*N))/(0.7846+log(1/(N*N)))*y(2)*(1-(1-exp(2.6*y(1)))/(2.6*y(1)))-0.4117*0.1238/(0.7846+log(1/(N*N)))*y(2)-0.5891*0.1238/(0.7846+log(1/(N*N)))*y(2)*((1-exp(-30*y(1)))/(30*y(1)))-... (1-0.1238/(0.7846+log(1/(N*N))))*y(2)-5.8e-16/(8*0.15)*7.444e15*0.1238/(0.7846+log(1/(N*N)))*30*y(2)*y(2)]; end end |
5樓2012-11-10 12:41:02
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 小魚(yú)118 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +5 | Ymlll 2026-03-18 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +6 | 小材化本科 2026-03-18 | 6/300 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +9 | rare12345 2026-03-18 | 9/450 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|