 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 687 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
小魚118銀蟲 (小有名氣)
|
[求助]
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ù)此樓» 猜你喜歡
計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有91人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有29人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有26人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
UJN物理學(xué)專業(yè)調(diào)劑
已經(jīng)有10人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名
已經(jīng)有0人回復(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問題 程序怎么會(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ùn)行時(shí)間的問題
已經(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)的誤差巨大問題求助！��！
已經(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ù)
- 一種簡單已用的熱加工圖畫法
已經(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è)人文集】畫NBO軌道的最簡單、最省錢方法
已經(jīng)有65人回復(fù)
4樓2012-11-10 10:56:31
【答案】應(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ù)沒法傳遞函數(shù)的參數(shù)到ode45. |
2樓2012-11-10 09:12:24
小魚118
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 5834.5
- 散金: 51
- 紅花: 4
- 帖子: 206
- 在線: 161.5小時(shí)
- 蟲號(hào): 1134404
- 注冊(cè): 2010-10-29
- 性別: MM
- 專業(yè): 光學(xué)
3樓2012-11-10 09:28:57
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★
小魚118: 金幣+5, ★★★★★最佳答案, 非常感謝，我同學(xué)成功了！ 2012-11-10 13:56:25
dbb627: 金幣+1, 感謝應(yīng)助 2012-11-12 08:57:21
小魚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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑一志愿華中師范大學(xué) +5 | 小江小江江江 2026-04-02 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 化工調(diào)劑303分，過四級(jí) +8 | 棲梧待風(fēng) 2026-04-02 | 8/400 |
|
|
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +6 | 1段時(shí)間 2026-03-31 | 7/350 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +7 | Ong3 2026-03-27 | 7/350 |
|
|
|
[考研] 289求調(diào)劑 +23 | 新時(shí)代材料 2026-03-27 | 26/1300 |
|
|
|
[考研] 求調(diào)劑，一志愿南京師范大學(xué)計(jì)算機(jī)專碩，初試373，六級(jí)通過， +3 | 計(jì)算機(jī)追夢人 2026-04-01 | 3/150 |
|
|
|
[考研] 11408 321分求調(diào)劑 +3 | huchun12138 2026-03-30 | 4/200 |
|
|
|
[考研] 295材料工程專碩求調(diào)劑 +19 | 1428151015 2026-03-27 | 19/950 |
|
|
|
[考研] 286求調(diào)劑 +16 | PolarBear11 2026-03-26 | 16/800 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +21 | 荷包蛋hyj 2026-03-26 | 21/1050 |
|
|
|
[考研] 調(diào)劑申請(qǐng) +8 | 張張張張zy 2026-03-31 | 9/450 |
|
|
|
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +6 | 小小機(jī)靈蟲 2026-03-29 | 12/600 |
|
|
|
[考研] 一志愿大連理工大學(xué)，機(jī)械工程學(xué)碩，341 +3 | 西瓜田的守望者 2026-03-30 | 3/150 |
|
|
|
[考研] 11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
|
[考研]
 304材料求調(diào)劑
+5
|
鐘llll 2026-03-26 | 5/250 |
|
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +6 | niko- 2026-03-27 | 6/300 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|