应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » Matlab循环的“诡异”问题
北京石油化工学院2026年研究生招生接收调剂公告
31/1返回列表
查看: 528  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

harp

木虫 (正式写手)

[求助] Matlab循环的“诡异”问题

问题：可以运行一次两排“%%%%%%%%%%” 之间的代码，也就是被循环语句，而且结果合理，速度非常快，按下enter就能显示结果，这应该能说明这段代码没有问题吧？

可加上for语句后运行起来却不行，运行的非常慢，一直busy，等不到结果；另外，就算把两排“%%%%%%%%%%” 之间的代码复制重新运行也不行；使用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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-07-11 15:36:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lds588

铁杆木虫 (著名写手)

愚钝

【答案】应助回帖

★ ★ ★
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:
问题：可以运行一次两排“%%%%%%%%%%” 之间的代码，也就是被循环语句，而且结果合理，速度非常快，按下enter就能显示结果，这应该能说明这段代码没有问题吧？

可加上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循环处加一个变量k（52和54行为我所加）
把改后的代码多执行几次你就能发现问题（while语句执行次数的问题）

目的是看看while循环什么时候跳出，我发现当运算时间较长时，始终跳不出while语句
只有生成的随机数合理，才会很容易跳出while循环，所以会执行快一些，这是一个概率问题，与你是否重启matlab没有关系
赞一下(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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

harp

木虫 (正式写手)
引用回帖:
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 harp 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602 找调剂 +4 逆时针快乐 2026-04-02 4/200 2026-04-04 19:32 by 蓝云思雨
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +9 万事宜臻 2026-04-04 9/450 2026-04-04 11:49 by 啵啵啵0119
[考研] 一志愿北京科技大学材料工程085601，求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 本科211，293分请求调剂 +6 莲菜就是藕吧 2026-04-03 7/350 2026-04-04 11:12 by 莲菜就是藕吧
[考研] 一志愿双非085502，267分，过四级求调剂 +3 再忙也要吃饭啊 2026-04-03 3/150 2026-04-04 05:03 by gswylq
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[考研] 化工调剂303分，过四级 +28 栖梧待风 2026-04-02 28/1400 2026-04-03 21:40 by qlm5820
[考研] 303求调剂 +10 DLkz1314. 2026-03-30 10/500 2026-04-03 18:03 by Jimmyandyou
[考研] 366求调剂 +7 sbdnd 2026-04-03 7/350 2026-04-03 12:40 by cymywx
[考研] 一志愿华东理工大学，080500学硕，317分，求调剂 +13 s1145 2026-03-31 15/750 2026-04-03 11:44 by msi123
[考研] 366求调剂一志愿东北大学 +8 运气来得若有似� 2026-04-02 8/400 2026-04-02 21:39 by dongzh2009
[考研] 0832食品科学与工程学硕282调剂 +4 鱼在水中游a 2026-04-02 7/350 2026-04-02 14:12 by baoball
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +9 双马尾痞老板2 2026-04-01 9/450 2026-04-02 12:14 by oooqiao
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
[考研] 一志愿安徽大学计算机科学与技术学硕，331分求调剂 +5 蒋昌鹏qtj 2026-04-01 5/250 2026-04-02 08:10 by fxue1114
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 262求调剂 +9 励志一定发文章 2026-03-31 10/500 2026-04-01 12:22 by sunshine0013
[考研] 【调剂】一志愿厦大生物与医药调剂 +3 Echo虾米 2026-03-31 3/150 2026-04-01 08:40 by JourneyLucky
[考研] 求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 262求调剂 +7 ZZ..000 2026-03-30 8/400 2026-03-31 10:05 by cal0306
信息提示
关闭
请填处理意见
关闭