 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 670 | 回復(fù): 4 | ||
小魚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 運行結(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]) 求高手指點 |
回復(fù)此樓» 猜你喜歡
2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有138人回復(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ù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有29人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有54人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助：如何用matlab實現(xiàn)一個沿x軸運動的且半徑隨x坐標(biāo)變化的圓的軌跡圖
已經(jīng)有4人回復(fù)
- 如何在matlab中計算對稱矩陣的高精度本征函數(shù)？
已經(jīng)有0人回復(fù)
- MATLAB程序總是出錯?
已經(jīng)有1人回復(fù)
- MATLAB問題 程序怎么會出錯？Error: Missing variable or function. 怎么修改？
已經(jīng)有1人回復(fù)
- matlab作出的圖都縮到一塊，如何調(diào)整？
已經(jīng)有7人回復(fù)
- MATLAB程序設(shè)計與工程應(yīng)用
已經(jīng)有144人回復(fù)
- matlab提前設(shè)定語句運行時間的問題
已經(jīng)有4人回復(fù)
- 求大神幫忙解題�。。�！
已經(jīng)有4人回復(fù)
- 絕熱毛細(xì)管（制冷）校核計算時壅塞流流量的計算
已經(jīng)有17人回復(fù)
- MATLAB程序調(diào)整
已經(jīng)有4人回復(fù)
- 關(guān)于使用lsqcurvefit進行最小二乘曲線擬合時候出現(xiàn)的誤差巨大問題求助�。�！
已經(jīng)有7人回復(fù)
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對論效應(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)整時間關(guān)系
已經(jīng)有2人回復(fù)
- 【beefly個人文集】畫NBO軌道的最簡單、最省錢方法
已經(jīng)有65人回復(fù)
2樓2012-11-10 09:12:24
小魚118
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 5820.5
- 散金: 51
- 紅花: 4
- 帖子: 206
- 在線: 161.5小時
- 蟲號: 1134404
- 注冊: 2010-10-29
- 性別: MM
- 專業(yè): 光學(xué)
3樓2012-11-10 09:28:57
4樓2012-11-10 10:56:31
【答案】應(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]這兩個參數(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)
我要訂閱樓主 小魚118 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 10/500 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
|
[考研]
 0854，計算機類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|