could anyone help me to fix the issue

I am using the command randi(4,4,1) it gives the output

 ans =

     1
     1
     1
     3

Could anyone tell me how to get the output in the following manner

ans=
    1
    2
    3
    4

[Walter Roberson]

Are you asking to find a random seed such that randi(4,4,1) will happen to emit [1; 2; 3; 4] ?

[jaah navi]

yes.

[Stephan]

Why dont you use:

x = (1:4)'

[jaah navi]

Is there any other possibility to have the output as

1
2
3
4.

Answer 1

[Walter Roberson]

s = struct('Type', 'twister', 'Seed', 1, 'State', ...
  uint32([2629073562; 2983301384; 681580311; 4033622241; 792772838; 3306961981; 92883131; 1785085746; 3364128315; 2402025379; 3224868746; 1145213362; 3784365245; 1948434636; 2667646161; 2598854474; 921967201; 1345782310; 4019597455; 2906395199; 1349669984; 2676817993; 4201769589; 2002781766; 3540177092; 4224925813; 3661313599; 1709930435; 1812273278; 2973452884; 1592291796; 3452239013; 3588672187; 3228651068; 3191454495; 3286135343; 2640545275; 3096953148; 3746505897; 2292163827; 2164382601; 1581410039; 2413832827; 2536571847; 179684232; 1638923698; 3155158821; 1454330362; 4050607484; 607322300; 2216566078; 597866774; 1036426282; 732815996; 3131865469; 2440339870; 2814550949; 1479383443; 2449469876; 3810238677; 2923086221; 437801529; 2891199990; 1886893516; 3898673786; 376204646; 1392372379; 4123661669; 1140754642; 3539167101; 2386309702; 3740957436; 4033654965; 1720449988; 3434980330; 4213508374; 3576835843; 2818865106; 1653162115; 2935588114; 3870616539; 102614847; 1476834675; 1220770796; 1233652508; 2385138085; 1300608482; 2753953039; 262993567; 4009374062; 2143978386; 1613109469; 3072671496; 1223816410; 4088822114; 3382188205; 2250056281; 1926821318; 3806317775; 2882166470; 94227745; 2877123406; 1030225246; 1555072155; 95009460; 1855512191; 2840453856; 2478087736; 836340051; 1566383306; 2625414289; 3519123538; 840509954; 2373484829; 676528503; 2783662465; 3557034492; 2566048980; 2347785709; 3566819907; 1311855742; 198269976; 2693520819; 2127070362; 961491174; 3932714317; 837664826; 4277891831; 3535515583; 2831416447; 3505045078; 3763313683; 367436315; 3614057572; 3780746374; 2693039652; 2297021184; 2224934154; 698822522; 2718629137; 1175446314; 2603507610; 2067589016; 2280810156; 2037033584; 3956938481; 1112874779; 3264939860; 2054107185; 2354026721; 1958640221; 2844284824; 3775753525; 2462549847; 3562644229; 3683686884; 3714884555; 2266356233; 2808583945; 980888698; 1137581788; 2771236582; 1975939317; 1605707990; 614167064; 767063856; 4227905160; 3590303986; 2932373212; 2230415839; 127157074; 2328724316; 3356372094; 3215726425; 282321962; 4226412442; 106823192; 3925701436; 610765913; 281952627; 1832011890; 2670621135; 1012800992; 1489632964; 1371755819; 2529629289; 607643288; 3941535311; 3202770816; 284461833; 3696778854; 325625733; 2671656400; 1391137252; 1240723705; 3132941411; 314202987; 1301784708; 2575857120; 313287791; 3569720512; 981744121; 2986286440; 1051168756; 3881027887; 1088809168; 3421075971; 2655923113; 2577977181; 3968201444; 1406585576; 4001025594; 1854233928; 3825832114; 954907921; 3109212504; 1706316388; 1550715292; 2934259476; 1892992132; 4050317728; 4110743408; 359920673; 3542060425; 597068400; 843058885; 3799566677; 2064802063; 860032883; 2732673041; 3457529378; 3660165513; 1125914118; 4212233956; 3757086488; 2632116865; 462233764; 161423763; 141910303; 2265109737; 1983126147; 2231395445; 2147383559; 611750005; 2194081082; 1699618892; 3841025952; 1294478167; 936120505; 4102008957; 4268194620; 1994019922; 1187332833; 3953675561; 1917283619; 896867387; 1634068959; 3950876764; 485558549; 45383166; 2795959904; 1717914918; 3449856475; 3796449494; 3318166191; 1007512487; 480498390; 1730605673; 3972660041; 2287008439; 539639323; 1569683418; 3795864463; 453065119; 2449404120; 265049099; 3514892628; 2665563904; 2534523080; 1969344934; 4294466702; 1478959417; 1858551310; 3029422620; 4121720519; 2933843154; 3209557296; 3111170404; 4264711934; 2448282623; 1388503581; 1999436060; 818418468; 3283991819; 2356924521; 2684567658; 1424193429; 1187340812; 3847120877; 1987859863; 1879502046; 3422594099; 1419478013; 4148487199; 2538837125; 3694851522; 636350498; 1097832595; 3779331880; 351970715; 2534774459; 3389311029; 3762283879; 3742425828; 3882821767; 2683353; 1981273229; 4068324016; 539226467; 2411256222; 1780609115; 2059099269; 2889497980; 1123848930; 540086248; 3467353606; 3362203896; 4058078927; 93044840; 426751932; 824266620; 3444590461; 2122918776; 3339845861; 4233923286; 3733051177; 2957657929; 3110908772; 2551930098; 1294126893; 2213336821; 2124571119; 49780221; 2901327722; 2493306969; 2545470627; 102527300; 3876142393; 4097726412; 3510695954; 660912408; 2033930425; 1601509561; 969180562; 635252598; 177954239; 3054207519; 4122051269; 2787463443; 1664731394; 2907371963; 1484884283; 1560546623; 2902374922; 2395942225; 451352804; 3346805556; 3459298550; 3482428591; 10753957; 4101820340; 3306178891; 1122941824; 355877597; 1683964498; 238724805; 3926649337; 3197046734; 4277634633; 2288745211; 4202067531; 526022968; 4017453944; 1499184106; 2677441952; 2227353703; 1995296581; 1690255681; 2920887680; 605849491; 3072795503; 746910746; 2709796449; 1225135658; 657841564; 4070363626; 144842260; 3718575695; 3159187032; 88291794; 3129049475; 410962657; 1728726693; 2397606939; 4126386549; 751549633; 4226219908; 1549973222; 3060733996; 3741110422; 2530947598; 2627897488; 2317706652; 1170828427; 2671701715; 1153351468; 3762293788; 4093330405; 1641962571; 699324101; 3173743570; 1798831929; 3467616712; 4198420524; 2448981354; 499920867; 968642107; 2140815539; 4193124145; 1639223168; 284638153; 396985542; 3543438633; 237854258; 3938010494; 24441053; 2947436871; 1273496002; 719279415; 3574242559; 1040109604; 3849196601; 3250223302; 3411729501; 3031943234; 2932285520; 2932420675; 2011314805; 2480850074; 3207806491; 462404995; 3279042455; 4270524229; 1064389665; 1894847490; 721365878; 357178131; 2827490451; 2604438657; 242514037; 678802395; 1322770750; 2747624534; 4246466163; 4188936761; 1207204018; 980275996; 841637218; 1468131552; 4102349079; 586888764; 3105466755; 1628818384; 2991889790; 2801191520; 2114916962; 1124291831; 3242113092; 1082871720; 3625937786; 2796251125; 1651820702; 3427511545; 2035120316; 1024058911; 4209506140; 2527167744; 213886228; 2514956543; 2450260579; 327684603; 3444379103; 3884997363; 2844468873; 2261078634; 3926825101; 439487268; 3789435080; 1212963762; 4259079565; 2772611204; 2534236055; 2430244594; 916922266; 903950702; 3381351589; 2268543712; 3616954837; 1273083041; 1682465785; 1342921678; 2593265787; 3033724173; 2988544460; 1824668777; 1214999983; 257453352; 4187931679; 3523379959; 1481153225; 4290295859; 220376185; 4136013972; 511679284; 3510589272; 1404047266; 3712771231; 2125374725; 32037606; 3601685135; 3433623522; 841138647; 1610171318; 3920699442; 1084892922; 3146108732; 3672652561; 1148331655; 2473777375; 3039860130; 1170979324; 575756423; 1389255297; 339011744; 2351938991; 4050094885; 2773634239; 3715040333; 3920910597; 562407139; 29887881; 1623822350; 358193390; 1822261341; 1243290919; 755541153; 153529770; 994467513; 1339524978; 3174283928; 2782204324; 182251010; 1833692038; 3477775846; 2820237500; 2165639585; 4089211432; 4010345846; 1238643345; 3710224584; 2251039304; 1196168985; 3165387311; 3920626153; 1796963839; 2112227260; 3358845388; 593715887; 3046897033; 2968478428; 3846781604; 3565923316; 2452128692; 68338106; 1427007580; 2192917968; 75235680; 2134869635; 1857807303; 3745016485; 2601327385; 35512535; 559792668; 2093376088; 939608650; 2259549051; 2605376692; 2058599240; 2980581379; 3160415220; 2739135905; 254886981; 1652380747; 371107437; 1123937393; 4185309254; 864314942; 2739416220; 2185572068; 1163546293; 3491702910; 996401156; 1198755052; 2898003956; 745796080; 4127642404; 4237523457; 1274635091; 3144139009; 1600421663; 4226154574; 538248802; 373236455; 116925273; 256]) );

rng(s);
randi(4, 4, 1)

[Stephan]

 o o
  |
 \_/

[jaah navi]

ok.could you please tell me what are the those values mentioned in uint32([......]).If I use all the values i am getting the output or else i am unable to get the output.

[jaah navi]

is there any other way inorder to have the output without mentioning those numbers in uint32([.....])

[Walter Roberson]

The Twister random number generate uses an internal state of 625 uint32 values. The internal state is not as simple as a simple uint32 seed. When you use rng() to seed with a simple number, it uses the number given and a particular algorithm to initialize the 625 internal state variables.

But if you are looking for something simpler, then rng() with any of the following seeds:

238 366 503 518 723 767 846 867 1029 1220 1273 1291 1456 1468 1704 1952 2163 2587 2674 2957 3632 4373 4519 4569 4583 5053 5062 5163 5352 6543 6587 6637 6708 6836 7152 7216 7347 7390 7403 7423 7640 7656 7884 7976 7993 8253 8304 8452 8733 8772 8999 9619 9691 10290 10699 11026 11304 11548 11593 12078 12269 12425 12635 12885 12935 13549 13733 13855 13858 14237 14674 15103 15473 15738 16139 16314 16343 16699 16958 17249 17669 18413 18812 19012 19308 19322 19551 19580 19696 19744 20607 20669 20808 21126 21298 21345 21522 21819 22390 22674 22719 22802 22943 23216 23420 23590 23695 23972 24193 24397 24690 24776 25030 25129 25283 25467 25615 25717 25724 25755 25796 25849 26304 26321 26349 26761 27307 27929 28293 28434 28716 29021 29296 29565 31036 31077 31121 31236 31330 31344 31730 32133 32200 32577 33116 33130 33162 33389 33665 34064 34316 34478 34541 35358 36141 37069 37239 37523 37588 37742 37748 38221 38300 38351 38442 38940 39229 39708 40043 40249 40459 40703 41640 42118 42137 42403 42624 42633 42969 43417 43557 43686 43846 43917 44050 44484 44685 44827 45708 45919 46210 46231 46358 46408 46424 46442 46455 46466 46490 46667 46751 46821 47327 47903 47949 48029 48248 48286 48432 48745 49199 49375 49377 49383 49578 49644 49880 50444 50509 50727 50835 51076 51156 51698 52214 52331 53086 53106 53414 53630 53730 53741 53973 54111 54177 54285 54384 54996 55176 55648 55717 56331 56599 56840 57557 57699 57712 57883 58234 58436 58441 58596 58639 58761 58763 59032 59430 60014 60166 60167 60179 60234 60487 60877 61081 62145 62472 62954 62985 63237 63440 63531 63593 63635 63754 64571 64724 64894 65020 65092 65108 65140 65395

I only searched from 0 to 2^16-1

[Stephan]

for [1; 3; 5; 7] use rng (767)

;-)

[jaah navi]

if i use 767 i am getting the following result in the command window

N_rng =

     767

Below results are for rng = 767

rows =

       2

rows =

       3     4     1

but what i actually need is

rows=1
rows=2
rows=3
rows=4

where as the order can be 2;3;4;1 or 3;4;1;2 or 3;2;1;4...

but i want the four numbers to be

displayed separately as how i shown.

I also tried with all the seeds displayed above but i unable to get the result what i actually need.

[jaah navi]

I would like to check on what basis of 0 to 2^16-1 the seeds are displayed.

[Walter Roberson]

matches = [];
for K = 0 : 2^16-1
    rng(K);
    if randi(4,4,1) == [1;2;3;4]
        matches(end+1) = K;
        nm = length(matches);
        if mod(nm,50) == 0
            fprintf('+\n');
        elseif mod(nm,10) == 0
            fprintf('+');   %no newline
        else
            fprintf('.');   %no newline
        end
    end
end
disp(matches)

[jaah navi]

ok. could you please tell me why K = 0 : 2^16-1 why it cant be 0 : 2^14-1,0 : 2^15-1,0 : 2^13-1,0 : 2^12-1,0 : 2^11-1 or some other.

[Walter Roberson]

Arbitrary. I knew that rng() required that the seed be a nonnegative integer, and I knew that it accepted at least 16 bits, but I did not bother to research whether it accepted 32 bit or 64 bit integers. I figured that I would program it to 16 bits and let it run, and if it created too many entries I could trim back, and if it found nothing then I could try to increase it.

[jaah navi]

ok. I tried with all the combination of seeds and i am getting result in the following manner Below results are for rng = 84030

rows =

     2

rows =

     3     4     1

Below results are for rng = 84543

rows =

     2

rows =

     3

rows =

     4     1
but i am unable to get the result in rows=2
rows=1
rows=4
rows=1
why it is so.Could you help me on this.

[Walter Roberson]

matchtable = cell(4,4,4,4);
for K = 0 : 2^16-1
    rng(K);
    thisrnd = num2cell(randi(4,4,1));
    idx = sub2ind([4, 4, 4, 4], thisrnd{:});
    if isempty(matchtable{idx}); matchtable{idx} = K; end
    if ~any(cellfun(@isempty,matchtable(:))); break; end
end
if K == 2^16-1
    fprintf('Did not find some combinations\n');
else
    fprintf('Done at seed %d\n', K);
end

Afterwards, if all goes well, then matchtable{A,B,C,D} is the lowest-number seed for rng() such that rng(matchtable{A,B,C,D}); randi(4,4,1) gives [A;B;C;D] . Thus if the task is to prove that all possible combinations can be generated at some point, then this finds the seeds that do that.

Answer 2

[Stephan]

 rng (84030)
 randi (4,4,1)

Gives your result - like all other numbers Walter told you...

 rng (767)
 randi (7,4,1)

Gives you:

1
3
5
7

You got it?

Best regards

Stephan

[jaah navi]

I am getting the following result N_rng =

         84030

Below results are for rng = 84030

rows =

       2

rows =

       3     4     1

when i am using 84030.

Answer 3

[Stephan]

See attached file...

Did you use:

 rng = 84030

???

try brackets like shown:

 rng (84030)

Best regards

Stephan

[jaah navi]

ok.when using bracket as how you have mentioned the result is right.But in my code i am trying with different values of rng. say for example N_rng=[1 2 3 4 5 ].The reason why i am using this is to check with all combinations of 4 numbers inorder to get the following result.

N_UE =

     4

Below results are for rng = 1

rows =

     2     4

rows =

     1     3

Below results are for rng = 2

rows =

     1     4

rows =

     3     2

Below results are for rng = 3

rows =

     2     3

rows =

     1     4

Below results are for rng = 4

rows =

     1     4

rows =

     2     3

Below results are for rng = 5

rows =

     4

rows =

     1

rows =

     3     2

But i am unable to get the result rows=1 rows=2 rows=3 rows=4 as separately I tried with several values but i can get the other combinations other than this.

I want all the combinations along with this to be displayed.If i use the command randi(4,4,1) it gives only 1;2;3;4 but i am unable to get the other combinations.Could you please help me what value should i need to use in N_rng=[.......] so that i can have all combinations.

[Dennis]

Maybe you should provide clearer informations about your actual goal.

If you want all possible combinations of vectors with elements 1:4 you shouldnt use randi.

If you want to manipulate randi to get a certain result, you can use the code Walter provided. Change the if condition and it will find you all values you can use with rng to get that result.

[jaah navi]

ok. I am trying with the following code:

N_UE=[ 4 ];
N_rng=[1 2 5 6 9 11 12 19 47 74 84 91 153 ];
for t=1:length(N_UE)
for s=1:length(N_rng)
this_seed = N_rng(s);
fprintf('Below results are for rng = %d\n', this_seed)
rng(this_seed);
unused_rows=1:N_UE(t);
while ~isempty(unused_rows)
N_UE_rows=ceil(sqrt(randi([numel(unused_rows)])));
if (N_UE_rows+1)==numel(unused_rows)
N_UE_rows=numel(unused_rows);
end
rows=unused_rows(randsample(length(unused_rows),N_UE_rows))
[~,idx]=find(ismember(unused_rows,rows));
unused_rows(idx)=[];
end
end
end

with respect to N_rng i am getting all combinations of 4, except

rows=1
rows=2
rows=3
rows=4.

could you please help me what value do i need to include in N_rng to get the above result.

[jaah navi]

Could anyone help me what value do i need to include in N_rng to get the above result as i tried with many number of values.