作者sunev (Veritas)
看板MATLAB
标题Re: [讨论] 请教N个for回圈如何实现呢
时间Tue Dec 6 18:50:07 2022
※ 引述《balancer12 (Do Something)》之铭言：
: 想请教各位大大
: 我有一个变数A里面可能有三个变异量
: 在A变异量里面又各有两个变数B
: 如此总共会有2*2*2=2^3=8组
: 写法如下
: GroupA_numPort=[ 1 2 3 ];
: GroupB_CapVariable=[ 10e-6 100e-6 ];
: iTest=0;
: for i= 1:length(GroupB_CapVariable)
: for j= 1:length(GroupB_CapVariable)
: for k= 1:length(GroupB_CapVariable)
: A_Port=[ GroupB_CapVariable(i) , GroupB_CapVariable(j) , GroupB_CapVaria
: ble(k) ] ;
: iTest=iTest+1;
: end
: end
: end
: 因变数A有三个，必须写三次for loop
: 若变数A变成5、10、99......
: 不太可能写这麽多个for loop@@
: 请问有什麽方式可以实现N个for 回圈吗？
: 谢谢大家指教！
ndgrid 可以做
结果会像是
[x1,x2,x3]=ndgrid(GroupB_CapVariable);
A_Port=[x1(:),x2(:),x3(:)];
[x1,x2,x3] 这个语法可以选择再用eval硬干字串
或是直接去魔改内建的ndgrid.m
(把nout 改掉，并直接对varagrout做再处理)
这边虽然避开了for-loop
但把结果一次存下来的结果就是会占用记忆体
尤其你的结果是n的d次方 (n,d定义请参照下方)
所以要格外小心记忆体会不会爆掉
当然这也是matlab 避开for-loop的通病就是了
另一个可能是你只是烦语法，不是想加速
那麽可以选择
iTest=0;
n = length(GroupB_CapVariable);
d = length(GroupA_numPort);
for i = 0:n^d-1
% 手动把i转成d进位，并把各位数算出来
% 请自行小心跑太久或overflow，毕竟n^d很可怕
remain = i;
A_Port = zeros(1,d);
for digit=d:-1:1
ind = mod(remain,d^(digit-1))+1;
A_Port(digit) = GroupB_CapVariable(ind);
remain = remain - d^(digit-1);
end
iTest=iTest+1;
end
以上code都还没debug过，请小心使用
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.54.137 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1670323809.A.E71.html
※ 编辑: sunev (140.112.54.137 台湾), 12/06/2022 18:52:07
※ 编辑: sunev (140.112.54.137 台湾), 12/06/2022 18:54:00
※ 编辑: sunev (140.112.54.137 台湾), 12/06/2022 19:04:07
1F：推 YoursEver: 好清流的pseudo-code，好久没见到这种讨论了。 12/06 22:02
2F：推 balancer12: 谢谢sunev大的建议！小的会在试试看～ 12/06 22:19
3F：→ sunev: matlab的TA大多不在意这种讨论.... 12/07 10:22