我不认为你的代码可以在矢量化。然而,这应该快一点
为nset = 1:10 0
为nindy = 1:5
如果(nset 2 + (nindy-1)。* 3) = = 1
C17 = C (1 + (nindy-1)。* 2)…
+ ((C (2 + (nindy-1)。* 2) - C (1 + (nindy-1)。* 2)) / 11)。*一个(nset 1 + (nindy-1)。* 3);
为nimp = 1:20 00
如果B (33。* 6 + nindy nimp) < = C17
D (nset nindy nimp) = 1;
结束
结束
结束
结束
结束
警告:
我没有测试的正确性也没有速度。
