Re: [讨论] 请教N个for循环如何实现呢

楼主: sunev (Veritas)   2022-12-06 18:50:07
※ 引述《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过,请小心使用
作者: YoursEver (猪是妳)   2022-12-06 22:02:00
好清流的pseudo-code,好久没见到这种讨论了。
作者: balancer12 (Do Something)   2022-12-06 22:19:00
谢谢sunev大的建议!小的会在试试看~
楼主: sunev (Veritas)   2022-12-07 10:22:00
matlab的TA大多不在意这种讨论....

Links booklink

Contact Us: admin [ a t ] ucptt.com