harp

木蟲 (正式寫手)

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

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

問題：可以運行一次兩排“%%%%%%%%%%” 之間的代碼，也就是被循環(huán)語句，而且結(jié)果合理，速度非常快，按下enter就能顯示結(jié)果，這應該能說明這段代碼沒有問題吧？

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

這是怎么回事��？？？

求高人指教

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

回復此樓

» 猜你喜歡

0856求調(diào)劑已經(jīng)有13人回復
0703化學已經(jīng)有11人回復
329求調(diào)劑已經(jīng)有7人回復
321求調(diào)劑已經(jīng)有8人回復
0856材料化工調(diào)劑總分330 已經(jīng)有12人回復
309求調(diào)劑已經(jīng)有3人回復
070300化學354求調(diào)劑已經(jīng)有4人回復
330分求調(diào)劑已經(jīng)有3人回復
299求調(diào)劑已經(jīng)有8人回復
一志愿北京理工大學本科211材料工程294求調(diào)劑已經(jīng)有6人回復

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

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

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

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

harp

木蟲 (正式寫手)

應助: 1 (幼兒園)
金幣: 3239.6
散金: 166
紅花: 2
帖子: 510
在線: 180.5小時
蟲號: 633829
注冊: 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 ...

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

贊一下

回復此樓

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

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

查看全部 3 個回答

lds588

鐵桿木蟲 (著名寫手)

愚鈍

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

【答案】應助回帖

★ ★ ★
xiegangmai(金幣+3): 謝謝應助！ 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:
問題：可以運行一次兩排“%%%%%%%%%%” 之間的代碼，也就是被循環(huán)語句，而且結(jié)果合理，速度非�？�，按下enter就能顯示結(jié)果，這應該能說明這段代碼沒有問題吧？

可加上for語句后運行起來卻不行，運行的非常慢， ...

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)處加一個變量k（52和54行為我所加）
把改后的代碼多執(zhí)行幾次你就能發(fā)現(xiàn)問題（while語句執(zhí)行次數(shù)的問題）

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

贊一下(1人)

回復此樓

[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

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

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

[考研] 300求調(diào)劑，材料科學英一數(shù)二 +8	leaflight 2026-03-24	8/400	2026-03-29 01:31 by fmesaito
[考研] 081200-11408-276學碩求調(diào)劑 +6	崔wj 2026-03-26	6/300	2026-03-29 01:11 by hanserlol
[考研] 2026年華南師范大學歡迎化學，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學子加入！ +3	llss0711 2026-03-28	5/250	2026-03-29 00:23 by MRMFYTL
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 一志愿中南大學化學0703總分337求調(diào)劑 +5	niko- 2026-03-27	5/250	2026-03-28 14:25 by 唐沐兒
[考研] 藥學105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +6	腳滑的守法公民 2026-03-27	6/300	2026-03-28 11:02 by gjlllb
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 330一志愿中國海洋大學化學工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 348求調(diào)劑 +4	小懶蟲不懶了 2026-03-27	5/250	2026-03-27 12:47 by 果果媽咪
[考研] 一志愿陜師大生物學071000，298分，求調(diào)劑 +5	SYA！ 2026-03-23	5/250	2026-03-27 09:29 by 不吃魚的貓
[考研] 0703化學求調(diào)劑 +3	丹青奶蓋 2026-03-26	5/250	2026-03-26 20:11 by macy2011
[考研] 機械學碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3	@357 2026-03-25	3/150	2026-03-26 16:34 by by.MENG
[考研] 材料科學與工程 317求調(diào)劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚的貓
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3	askeladz 2026-03-26	3/150	2026-03-26 09:08 by 給你你注意休息
[考研] 296求調(diào)劑 +4	汪！？！ 2026-03-25	7/350	2026-03-25 16:41 by 汪��？！
[考研] 086003食品工程求調(diào)劑 +6	淼淼111 2026-03-24	6/300	2026-03-25 10:29 by 3Strings
[考研] 求調(diào)劑 +6	研研，接電話 2026-03-24	7/350	2026-03-24 17:01 by barlinike
[考研] 292求調(diào)劑 +4	鵝鵝鵝額額額額?/a> 2026-03-24	4/200	2026-03-24 16:41 by peike