版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

harp

木蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 3239.6
散金: 166
紅花: 2
帖子: 510
在線: 180.5小時(shí)
蟲號(hào): 633829
注冊(cè): 2008-10-22
性別: GG
專業(yè): 工程熱物理與能源利用

[求助] Matlab循環(huán)的“詭異”問題

問題：可以運(yùn)行一次兩排“%%%%%%%%%%” 之間的代碼，也就是被循環(huán)語句，而且結(jié)果合理，速度非�？欤聪耬nter就能顯示結(jié)果，這應(yīng)該能說明這段代碼沒有問題吧？

可加上for語句后運(yùn)行起來卻不行，運(yùn)行的非常慢，一直busy，等不到結(jié)果；另外，就算把兩排“%%%%%%%%%%” 之間的代碼復(fù)制重新運(yùn)行也不行；使用clear命令清楚后再復(fù)制運(yùn)行也不行，除非重啟matlab才能運(yùn)行。

這是怎么回事��？？？

求高人指教

xm=zeros(42,10);
for i=1:10

%%%%%%%%%%%%%%%%%%%%%%

C=[0.002026883,0.007821721,0.023555797,0.014122603,0.027523128,0.019160052,0.05981332,0.072422514,0.277256157,0.496297825];

P=C(1);
Ti=C(2);
K=C(3);
Na=C(4);
Mg=C(5);
S=C(6);
Fe=C(7);
Ca=C(8);
Al=C(9);
Si=C(10);

x=zeros(42,1);

%P
r1=rand(2,1);
xP=(r1/sum(r1))*P;
x(3)=xP(1)/3;
x(6)=xP(2)/2;
Ca=Ca-5*x(3)-x(6);
Al=Al-3*x(6);

%Ti
r2=rand(2,1);
xTi=(r2/sum(r2))*Ti;
x(41)=xTi(1);
x(42)=xTi(2);

%K
r3=rand(3,1);
xK=(r3/sum(r3))*K;
x(25)=xK(1)/1.5;
x(26)=xK(2);
x(27)=xK(3);
Al=Al-5.49999*x(25)-x(26);
Si=Si-6.5*x(25)-3*x(26);
S=S-2*x(27);
Fe=Fe-3*x(27);

%Na
St=-1;
Mgt=-1;
while (St<=0)||(Mgt<=0)
Nat=Na;
Mgt=Mg;
St=S;
Fet=Fe;
Cat=Ca;
Alt=Al;
Sit=Si;
r4=rand(8,1);
xNa=(r4/sum(r4))*Nat;
x(32)=xNa(1)/0.33;
x(33)=xNa(2)/0.7;
x(34)=xNa(3)/2;
x(35)=xNa(4)/2;
x(36)=xNa(5);
x(37)=xNa(6);
x(38)=xNa(7);
x(39)=xNa(8);
Alt=Alt-1.67*x(32)-2.6*x(33)-x(36)-x(37)-x(38);
Mgt=Mgt-0.33*x(32)-0.7*x(33);
St=St-2*x(34)-x(35)-2*x(39);
Fet=Fet-3*x(39);
Cat=Cat-x(34);
Sit=Sit-4*x(32)-8*x(33)-2*x(37)-3*x(38);
end

Mg=Mgt;
S=St;
Fe=Fet;
Ca=Cat;
Al=Alt;
Si=Sit;

%Mg
r5=rand(5,1);
xMg=(r5/sum(r5))*Mg;
x(9)=xMg(1);
x(28)=xMg(2);
x(29)=xMg(3);
x(30)=xMg(4);
x(31)=xMg(5);
S=S-x(31);
Fe=Fe-2*x(29);
Ca=Ca-x(9);
Al=Al-2*x(28);

%S
r6=rand(9,1);
xS=(r6/sum(r6))*S;
x(10)=xS(1);
x(11)=xS(2);
x(12)=xS(3);
x(13)=xS(4);
x(14)=xS(5)/3;
x(21)=xS(6)/2;
x(22)=xS(7)/2;
x(23)=xS(8);
x(24)=xS(9);

Fe=Fe-0.875*x(13)-2*x(14)-x(21)-x(22)-x(23)-x(24);
Ca=Ca-x(10)-x(11)-x(12);

%Fe
r7=rand(6,1);
xFe=(r7/sum(r7))*Fe;
x(15)=xFe(1)/2;
x(16)=xFe(2)/3;
x(17)=xFe(3)/5;
x(18)=xFe(4);
x(19)=xFe(5);
x(20)=xFe(6);

Al=Al-2*x(17);
Si=Si-3*x(17);

%Ca
r8=rand(4,1);
xCa=(r8/sum(r8))*Ca;
x(4)=xCa(1);
x(5)=xCa(2);
x(7)=xCa(3);
x(8)=xCa(4);

Al=Al-2*x(4)-2*x(5);
Si=Si-2*x(4)-7*x(5);

%Al
r9=rand(2,1);
xAl=(r9/sum(r9))*Al;
x(1)=xAl(1)/2;
x(2)=xAl(2);
Si=Si-2*x(1);

%Si
x(40)=Si;

xm(:,i)=x;

%%%%%%%%%%%%%%%%%%%%

end

回復(fù)此樓

» 猜你喜歡

299求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有6人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有8人回復(fù)
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑已經(jīng)有5人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有7人回復(fù)
289求調(diào)劑已經(jīng)有13人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有9人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有6人回復(fù)
調(diào)劑求院校招收已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有8人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

關(guān)于matlab的一個(gè)詭異問題已經(jīng)有7人回復(fù)
求matlab求解一個(gè)方程組的問題已經(jīng)有6人回復(fù)
matlab for循環(huán) 已經(jīng)有9人回復(fù)
Matlab兩個(gè)小問題已經(jīng)有10人回復(fù)
關(guān)于MATLAB和OpenCV對(duì)相同算法的運(yùn)算速度對(duì)比問題已經(jīng)有20人回復(fù)
MATLAB安裝后出現(xiàn)一堆java的問題已經(jīng)有12人回復(fù)
【求助】求Matlab循環(huán)作圖已經(jīng)有12人回復(fù)
【求助】問一個(gè)關(guān)于matlab中l(wèi)egend的問題已經(jīng)有11人回復(fù)

1樓 2011-07-11 15:36:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lds588

鐵桿木蟲 (著名寫手)

愚鈍

應(yīng)助: 15 (小學(xué)生)
金幣: 10949.6
散金: 2018
紅花: 5
帖子: 1444
在線: 206.3小時(shí)
蟲號(hào): 890720
注冊(cè): 2009-11-01
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

【答案】應(yīng)助回帖

★ ★ ★
xiegangmai(金幣+3): 謝謝應(yīng)助！ 2011-07-11 19:41:47
harp(金幣+50): 受教了，非常感謝��！ 2011-07-11 19:51:46

引用回帖:

Originally posted by harp at 2011-07-11 15:36:19:
問題：可以運(yùn)行一次兩排“%%%%%%%%%%” 之間的代碼，也就是被循環(huán)語句，而且結(jié)果合理，速度非�？�，按下enter就能顯示結(jié)果，這應(yīng)該能說明這段代碼沒有問題吧？

可加上for語句后運(yùn)行起來卻不行，運(yùn)行的非常慢， ...

CODE:

clear

clc

xm=zeros(42,10);

for i=1:1

%%%%%%%%%%%%%%%%%%%%%%

C=[0.002026883,0.007821721,0.023555797,0.014122603,0.027523128,0.019160052,0.05981332,0.072422514,0.277256157,0.496297825];

P=C(1);

Ti=C(2);

K=C(3);

Na=C(4);

Mg=C(5);

S=C(6);

Fe=C(7);

Ca=C(8);

Al=C(9);

Si=C(10);

x=zeros(42,1);

%P

r1=rand(2,1);

xP=(r1/sum(r1))*P;

x(3)=xP(1)/3;

x(6)=xP(2)/2;

Ca=Ca-5*x(3)-x(6);

Al=Al-3*x(6);

%Ti

r2=rand(2,1);

xTi=(r2/sum(r2))*Ti;

x(41)=xTi(1);

x(42)=xTi(2);

%K

r3=rand(3,1);

xK=(r3/sum(r3))*K;

x(25)=xK(1)/1.5;

x(26)=xK(2);

x(27)=xK(3);

Al=Al-5.49999*x(25)-x(26);

Si=Si-6.5*x(25)-3*x(26);

S=S-2*x(27);

Fe=Fe-3*x(27);

%Na

St=-1;

Mgt=-1;

k=0;

while (St<=0)||(Mgt<=0)

    k=k+1

Nat=Na;

Mgt=Mg;

St=S;

Fet=Fe;

Cat=Ca;

Alt=Al;

Sit=Si;

r4=rand(8,1);

xNa=(r4/sum(r4))*Nat;

x(32)=xNa(1)/0.33;

x(33)=xNa(2)/0.7;

x(34)=xNa(3)/2;

x(35)=xNa(4)/2;

x(36)=xNa(5);

x(37)=xNa(6);

x(38)=xNa(7);

x(39)=xNa(8);

Alt=Alt-1.67*x(32)-2.6*x(33)-x(36)-x(37)-x(38);

Mgt=Mgt-0.33*x(32)-0.7*x(33);

St=St-2*x(34)-x(35)-2*x(39);

Fet=Fet-3*x(39);

Cat=Cat-x(34);

Sit=Sit-4*x(32)-8*x(33)-2*x(37)-3*x(38);

end

Mg=Mgt;

S=St;

Fe=Fet;

Ca=Cat;

Al=Alt;

Si=Sit;

%Mg

r5=rand(5,1);

xMg=(r5/sum(r5))*Mg;

x(9)=xMg(1);

x(28)=xMg(2);

x(29)=xMg(3);

x(30)=xMg(4);

x(31)=xMg(5);

S=S-x(31);

Fe=Fe-2*x(29);

Ca=Ca-x(9);

Al=Al-2*x(28);

%S

r6=rand(9,1);

xS=(r6/sum(r6))*S;

x(10)=xS(1);

x(11)=xS(2);

x(12)=xS(3);

x(13)=xS(4);

x(14)=xS(5)/3;

x(21)=xS(6)/2;

x(22)=xS(7)/2;

x(23)=xS(8);

x(24)=xS(9);

Fe=Fe-0.875*x(13)-2*x(14)-x(21)-x(22)-x(23)-x(24);

Ca=Ca-x(10)-x(11)-x(12);

%Fe

r7=rand(6,1);

xFe=(r7/sum(r7))*Fe;

x(15)=xFe(1)/2;

x(16)=xFe(2)/3;

x(17)=xFe(3)/5;

x(18)=xFe(4);

x(19)=xFe(5);

x(20)=xFe(6);

Al=Al-2*x(17);

Si=Si-3*x(17);

%Ca

r8=rand(4,1);

xCa=(r8/sum(r8))*Ca;

x(4)=xCa(1);

x(5)=xCa(2);

x(7)=xCa(3);

x(8)=xCa(4);

Al=Al-2*x(4)-2*x(5);

Si=Si-2*x(4)-7*x(5);

%Al

r9=rand(2,1);

xAl=(r9/sum(r9))*Al;

x(1)=xAl(1)/2;

x(2)=xAl(2);

Si=Si-2*x(1);

%Si

x(40)=Si;

%%%%%%%%%%%%%%%%%%%%

xm(:,i)=x;

end

代碼為了告訴你差別我在你的while循環(huán)處加一個(gè)變量k（52和54行為我所加）
把改后的代碼多執(zhí)行幾次你就能發(fā)現(xiàn)問題（while語句執(zhí)行次數(shù)的問題）

目的是看看while循環(huán)什么時(shí)候跳出，我發(fā)現(xiàn)當(dāng)運(yùn)算時(shí)間較長(zhǎng)時(shí)，始終跳不出while語句
只有生成的隨機(jī)數(shù)合理，才會(huì)很容易跳出while循環(huán)，所以會(huì)執(zhí)行快一些，這是一個(gè)概率問題，與你是否重啟matlab沒有關(guān)系

贊一下(1人)

回復(fù)此樓

[url=http://weibo.com/u/2577306830?s=6uyXnP][img]http://service.t.sina.com.cn/widget/qmd/2577306830/fa2fdb00/1.png[/img][/url]

2樓2011-07-11 18:36:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

harp

木蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 3239.6
散金: 166
紅花: 2
帖子: 510
在線: 180.5小時(shí)
蟲號(hào): 633829
注冊(cè): 2008-10-22
性別: GG
專業(yè): 工程熱物理與能源利用

引用回帖:

Originally posted by lds588 at 2011-07-11 18:36:34:
[code]clear
clc
xm=zeros(42,10);
for i=1:1

%%%%%%%%%%%%%%%%%%%%%%

C=[0.002026883,0.007821721,0.023555797,0.014122603,0.027523128,0.019160052,0.05981332,0.072422514,0.277256157,0.4962978 ...

非常感謝你的熱心幫助，給你添麻煩了~~~

贊一下

回復(fù)此樓

3樓2011-07-11 21:20:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 harp 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 295材料工程專碩求調(diào)劑 +7	1428151015 2026-03-27	7/350	2026-03-28 19:58 by S240
[考研] 321求調(diào)劑 +6	璞玉~~ 2026-03-25	7/350	2026-03-28 17:48 by 璞玉~~
[考研] 085602 化工專碩 338分求調(diào)劑 +12	路癡小琪 2026-03-27	12/600	2026-03-28 15:41 by L135790
[考研] 求調(diào)劑 +3	QiMing7 2026-03-25	4/200	2026-03-28 14:30 by QiMing7
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 085600 286分材料求調(diào)劑 +7	麻辣魷魚 2026-03-27	8/400	2026-03-28 12:17 by zllcz
[考研] 材料277求調(diào)劑 +7	min3 2026-03-24	7/350	2026-03-28 11:39 by xuxiang
[考研] 291求調(diào)劑 +15	hhhhxn.. 2026-03-23	21/1050	2026-03-28 11:26 by self2008
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +6	腳滑的守法公民 2026-03-27	6/300	2026-03-28 11:02 by gjlllb
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9	丹青奶蓋 2026-03-26	10/500	2026-03-28 07:45 by barnett0632
[考研] 085701環(huán)境工程求調(diào)劑 +9	多久上課 2026-03-27	9/450	2026-03-28 03:58 by fmesaito
[考研] 一志愿211院校 344分東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5	丶風(fēng)雪夜歸人丶 2026-03-26	8/400	2026-03-27 19:22 by 丶風(fēng)雪夜歸人丶
[考研] 305求調(diào)劑 +5	哇盧卡庫 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 081700 調(diào)劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機(jī) 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 334分一志愿武理-080500 材料求調(diào)劑 +4	李李不服輸 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[論文投稿] 急發(fā)核心期刊論文 +3	賢達(dá)問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

harp

[求助] Matlab循環(huán)的“詭異”問題

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

lds588

【答案】應(yīng)助回帖

harp

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频