MATLAB: Combinations of a vector with replacement

#combinatronics

Hello everyone. A = [2 3 4]. I am trying to generate these 2 combinations:

22, 23 ,24 ,32 ,33 ,34, 42 ,43 ,44 (order matters so 23 not equal to 32 so nchoosek(A,2) does not work)
222, 233, 232, 234, 243, 333, 322, 323, 344, 343 , 324,342, 422, 424, 433 ,434 , 432, 423, 444. perms(A) does not do the job here because it is without replacement.

how can I do these?

Best Answer

A = [2 3 4]
[A1,A2] = ndgrid(A);
arrayfun(@(k) polyval([A2(k),A1(k)],10), 1:numel(A1))
[A1,A2,A3] = ndgrid(A);
arrayfun(@(k) polyval([A3(k),A2(k),A1(k)],10), 1:numel(A1))

Related Solutions

MATLAB: How to solve the equation quickly

try 'fsolve'

fplot('(47.457/b1)-(76.3792)+(3.0499)-(((2.7272^b1)*(1.0033)+(1.59^b1)*(0.46376)+(2.4506^b1)*(0.89632)+(2.6129^b1)*(0.96046)+(1.367^b1)*(0.31264)+(0.39284^b1)*(-0.93436)+(1.5651^b1)*(0.44795)+(1.5083^b1)*(0.41096)+(0.81499^b1)*(-0.20458)+(1.6526^b1)*(0.50236)+(1.7819^b1)*(0.5777)+(2.4563^b1)*(0.89867)+(0.88014^b1)*(-0.12768)+(2.5081^b1)*(0.91954)+(2.741^b1)*(1.0083)+(1.8194^b1)*(0.59851)+(2.033^b1)*(0.70951)+(1.6526^b1)*(0.50236)+(2.4283^b1)*(0.88718)+(1.2561^b1)*(0.22802)+(0.67065^b1)*(-0.39951)+(1.4447^b1)*(0.36793)+(0.49793^b1)*(-0.69729)+(1.2089^b1)*(0.18969)+(2.3449^b1)*(0.85226)+(0.49793^b1)*(-0.69729)+(0.39284^b1)*(-0.93436)+(0.49793^b1)*(-0.69729)+(1.8122^b1)*(0.59453)+(0.12951^b1)*(-2.044)+(0.58896^b1)*(-0.5294)+(2.6787^b1)*(0.98533)+(0.39284^b1)*(-0.93436)+(2.7946^b1)*(1.0277)+(1.7509^b1)*(0.56012)+(1.9784^b1)*(0.68228)+(2.4229^b1)*(0.88499)+(1.2561^b1)*(0.22802)+(1.742^b1)*(0.55506)+(2.0081^b1)*(0.69717)+(1.2864^b1)*(0.25188)+(2.3003^b1)*(0.83305)+(0.19421^b1)*(-1.6388)+(1.616^b1)*(0.47995)+(2.521^b1)*(0.92467)+(0.88014^b1)*(-0.12768)+(1.5459^b1)*(0.43562)+(2.0526^b1)*(0.7191)+(0.74549^b1)*(-0.29371)+(2.505^b1)*(0.91829)+(1.8026^b1)*(0.58922)+(1.671^b1)*(0.51344)+(2.0037^b1)*(0.69501)+(2.5104^b1)*(0.92042)+(1.6863^b1)*(0.52251)+(0.74549^b1)*(-0.29371)+(1.4695^b1)*(0.38492)+(1.719^b1)*(0.54174)+(1.9258^b1)*(0.65532)+(1.8122^b1)*(0.59453)+(0.39284^b1)*(-0.93436)+(0.49793^b1)*(-0.69729)+(2.2284^b1)*(0.80128)+(2.2822^b1)*(0.82514)+(1.3017^b1)*(0.26365)+(2.0761^b1)*(0.73051)+(2.4563^b1)*(0.89867)+(2.7046^b1)*(0.99495)+(1.6818^b1)*(0.51988)+(0.26185^b1)*(-1.34)+(0.39284^b1)*(-0.93436)+(0.39284^b1)*(-0.93436)+(0.12951^b1)*(-2.044)+(0.74549^b1)*(-0.29371)+(0.39284^b1)*(-0.93436)+(0.58896^b1)*(-0.5294)+(0.19421^b1)*(-1.6388)+(0.88014^b1)*(-0.12768)+(0.58896^b1)*(-0.5294)+(0.26185^b1)*(-1.34)+(0.26185^b1)*(-1.34)+(0.26185^b1)*(-1.34)+(0.39284^b1)*(-0.93436)+(0.26185^b1)*(-1.34)+(2.794^b1)*(1.0275)+(1.8122^b1)*(0.59453)+(0.26185^b1)*(-1.34)+(0.39284^b1)*(-0.93436)+(0.26185^b1)*(-1.34)+(0.67065^b1)*(-0.39951)+(1.1598^b1)*(0.14828)+(0.74549^b1)*(-0.29371)+(1.2561^b1)*(0.22802)+(0.58896^b1)*(-0.5294)+(1^b1)*(0)+(0.26185^b1)*(-1.34)+(2.0761^b1)*(0.73051)+(1^b1)*(0)+(1.3017^b1)*(0.26365)+(1.0556^b1)*(0.054134)+0)/(5^(b1*47.457)))+((((2.7272^b1)+(1.59^b1)+(2.4506^b1)+(2.6129^b1)+(1.367^b1)+(0.39284^b1)+(1.5651^b1)+(1.5083^b1)+(0.81499^b1)+(1.6526^b1)+(1.7819^b1)+(2.4563^b1)+(0.88014^b1)+(2.5081^b1)+(2.741^b1)+(1.8194^b1)+(2.033^b1)+(1.6526^b1)+(2.4283^b1)+(1.2561^b1)+(0.67065^b1)+(1.4447^b1)+(0.49793^b1)+(1.2089^b1)+(2.3449^b1)+(0.49793^b1)+(0.39284^b1)+(0.49793^b1)+(1.8122^b1)+(0.12951^b1)+(0.58896^b1)+(2.6787^b1)+(0.39284^b1)+(2.7946^b1)+(1.7509^b1)+(1.9784^b1)+(2.4229^b1)+(1.2561^b1)+(1.742^b1)+(2.0081^b1)+(1.2864^b1)+(2.3003^b1)+(0.19421^b1)+(1.616^b1)+(2.521^b1)+(0.88014^b1)+(1.5459^b1)+(2.0526^b1)+(0.74549^b1)+(2.505^b1)+(1.8026^b1)+(1.671^b1)+(2.0037^b1)+(2.5104^b1)+(1.6863^b1)+(0.74549^b1)+(1.4695^b1)+(1.719^b1)+(1.9258^b1)+(1.8122^b1)+(0.39284^b1)+(0.49793^b1)+(2.2284^b1)+(2.2822^b1)+(1.3017^b1)+(2.0761^b1)+(2.4563^b1)+(2.7046^b1)+(1.6818^b1)+(0.26185^b1)+(0.39284^b1)+(0.39284^b1)+(0.12951^b1)+(0.74549^b1)+(0.39284^b1)+(0.58896^b1)+(0.19421^b1)+(0.88014^b1)+(0.58896^b1)+(0.26185^b1)+(0.26185^b1)+(0.26185^b1)+(0.39284^b1)+(0.26185^b1)+(2.794^b1)+(1.8122^b1)+(0.26185^b1)+(0.39284^b1)+(0.26185^b1)+(0.67065^b1)+(1.1598^b1)+(0.74549^b1)+(1.2561^b1)+(0.58896^b1)+(1^b1)+(0.26185^b1)+(2.0761^b1)+(1^b1)+(1.3017^b1)+(1.0556^b1)+0)*76.3792)/(5^(b1*47.457)))',[0 100]);

fsolve(@(b1)(47.457/b1)-(76.3792)+(3.0499)-(((2.7272^b1)*(1.0033)+(1.59^b1)*(0.46376)+(2.4506^b1)*(0.89632)+(2.6129^b1)*(0.96046)+(1.367^b1)*(0.31264)+(0.39284^b1)*(-0.93436)+(1.5651^b1)*(0.44795)+(1.5083^b1)*(0.41096)+(0.81499^b1)*(-0.20458)+(1.6526^b1)*(0.50236)+(1.7819^b1)*(0.5777)+(2.4563^b1)*(0.89867)+(0.88014^b1)*(-0.12768)+(2.5081^b1)*(0.91954)+(2.741^b1)*(1.0083)+(1.8194^b1)*(0.59851)+(2.033^b1)*(0.70951)+(1.6526^b1)*(0.50236)+(2.4283^b1)*(0.88718)+(1.2561^b1)*(0.22802)+(0.67065^b1)*(-0.39951)+(1.4447^b1)*(0.36793)+(0.49793^b1)*(-0.69729)+(1.2089^b1)*(0.18969)+(2.3449^b1)*(0.85226)+(0.49793^b1)*(-0.69729)+(0.39284^b1)*(-0.93436)+(0.49793^b1)*(-0.69729)+(1.8122^b1)*(0.59453)+(0.12951^b1)*(-2.044)+(0.58896^b1)*(-0.5294)+(2.6787^b1)*(0.98533)+(0.39284^b1)*(-0.93436)+(2.7946^b1)*(1.0277)+(1.7509^b1)*(0.56012)+(1.9784^b1)*(0.68228)+(2.4229^b1)*(0.88499)+(1.2561^b1)*(0.22802)+(1.742^b1)*(0.55506)+(2.0081^b1)*(0.69717)+(1.2864^b1)*(0.25188)+(2.3003^b1)*(0.83305)+(0.19421^b1)*(-1.6388)+(1.616^b1)*(0.47995)+(2.521^b1)*(0.92467)+(0.88014^b1)*(-0.12768)+(1.5459^b1)*(0.43562)+(2.0526^b1)*(0.7191)+(0.74549^b1)*(-0.29371)+(2.505^b1)*(0.91829)+(1.8026^b1)*(0.58922)+(1.671^b1)*(0.51344)+(2.0037^b1)*(0.69501)+(2.5104^b1)*(0.92042)+(1.6863^b1)*(0.52251)+(0.74549^b1)*(-0.29371)+(1.4695^b1)*(0.38492)+(1.719^b1)*(0.54174)+(1.9258^b1)*(0.65532)+(1.8122^b1)*(0.59453)+(0.39284^b1)*(-0.93436)+(0.49793^b1)*(-0.69729)+(2.2284^b1)*(0.80128)+(2.2822^b1)*(0.82514)+(1.3017^b1)*(0.26365)+(2.0761^b1)*(0.73051)+(2.4563^b1)*(0.89867)+(2.7046^b1)*(0.99495)+(1.6818^b1)*(0.51988)+(0.26185^b1)*(-1.34)+(0.39284^b1)*(-0.93436)+(0.39284^b1)*(-0.93436)+(0.12951^b1)*(-2.044)+(0.74549^b1)*(-0.29371)+(0.39284^b1)*(-0.93436)+(0.58896^b1)*(-0.5294)+(0.19421^b1)*(-1.6388)+(0.88014^b1)*(-0.12768)+(0.58896^b1)*(-0.5294)+(0.26185^b1)*(-1.34)+(0.26185^b1)*(-1.34)+(0.26185^b1)*(-1.34)+(0.39284^b1)*(-0.93436)+(0.26185^b1)*(-1.34)+(2.794^b1)*(1.0275)+(1.8122^b1)*(0.59453)+(0.26185^b1)*(-1.34)+(0.39284^b1)*(-0.93436)+(0.26185^b1)*(-1.34)+(0.67065^b1)*(-0.39951)+(1.1598^b1)*(0.14828)+(0.74549^b1)*(-0.29371)+(1.2561^b1)*(0.22802)+(0.58896^b1)*(-0.5294)+(1^b1)*(0)+(0.26185^b1)*(-1.34)+(2.0761^b1)*(0.73051)+(1^b1)*(0)+(1.3017^b1)*(0.26365)+(1.0556^b1)*(0.054134)+0)/(5^(b1*47.457)))+((((2.7272^b1)+(1.59^b1)+(2.4506^b1)+(2.6129^b1)+(1.367^b1)+(0.39284^b1)+(1.5651^b1)+(1.5083^b1)+(0.81499^b1)+(1.6526^b1)+(1.7819^b1)+(2.4563^b1)+(0.88014^b1)+(2.5081^b1)+(2.741^b1)+(1.8194^b1)+(2.033^b1)+(1.6526^b1)+(2.4283^b1)+(1.2561^b1)+(0.67065^b1)+(1.4447^b1)+(0.49793^b1)+(1.2089^b1)+(2.3449^b1)+(0.49793^b1)+(0.39284^b1)+(0.49793^b1)+(1.8122^b1)+(0.12951^b1)+(0.58896^b1)+(2.6787^b1)+(0.39284^b1)+(2.7946^b1)+(1.7509^b1)+(1.9784^b1)+(2.4229^b1)+(1.2561^b1)+(1.742^b1)+(2.0081^b1)+(1.2864^b1)+(2.3003^b1)+(0.19421^b1)+(1.616^b1)+(2.521^b1)+(0.88014^b1)+(1.5459^b1)+(2.0526^b1)+(0.74549^b1)+(2.505^b1)+(1.8026^b1)+(1.671^b1)+(2.0037^b1)+(2.5104^b1)+(1.6863^b1)+(0.74549^b1)+(1.4695^b1)+(1.719^b1)+(1.9258^b1)+(1.8122^b1)+(0.39284^b1)+(0.49793^b1)+(2.2284^b1)+(2.2822^b1)+(1.3017^b1)+(2.0761^b1)+(2.4563^b1)+(2.7046^b1)+(1.6818^b1)+(0.26185^b1)+(0.39284^b1)+(0.39284^b1)+(0.12951^b1)+(0.74549^b1)+(0.39284^b1)+(0.58896^b1)+(0.19421^b1)+(0.88014^b1)+(0.58896^b1)+(0.26185^b1)+(0.26185^b1)+(0.26185^b1)+(0.39284^b1)+(0.26185^b1)+(2.794^b1)+(1.8122^b1)+(0.26185^b1)+(0.39284^b1)+(0.26185^b1)+(0.67065^b1)+(1.1598^b1)+(0.74549^b1)+(1.2561^b1)+(0.58896^b1)+(1^b1)+(0.26185^b1)+(2.0761^b1)+(1^b1)+(1.3017^b1)+(1.0556^b1)+0)*76.3792)/(5^(b1*47.457))),1)

MATLAB: All combinations not being generated

update: In line 3 of the code A3(:) is repeated twice, resulting in the 'error'

Follwing is what it should read :

my_matrix= [A1(:),A2(:),A3(:)]

Best Answer

Related Solutions

MATLAB: How to solve the equation quickly

MATLAB: All combinations not being generated

Related Question