I need to make a general code

Question

Zeina Abdullah 2021-12-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1619700-i-need-to-make-a-general-code

评论： Zeina Abdullah 2022-1-3

I mean that I have this code

%function const_interleaver(select)
%% constructing the orignal address matrix
Address=zeros(24); s
for i=1:24
    for j=1:24
        Address(i,j)=24*(i-1)+j;
    end
end
%% chosing the matrix size
if select==1 % 8-bit
    for i=1:8
        for j=1:8
            Address(i,j)=Address(i,j)-(i-1)*16     
        end
    end
elseif select==2 % 16-bit
    for i=1:16
        for j=1:16
            Address(i,j)=Address(i,j)-(i-1)*8;
        end
    end
else select==3 % 24-bit
end 

and it works in cases, and these cases have a common code between them. I need to make a general code, and by choosing one of the cases, part of it will work, for example, in the case of (select=1) 8 * 8, the code will work to part 8, but if it is (select=2) 16 * 16, it will work to part 16 * 16 And if it is (select=3) 24 * 24, it will work to part 24.

and the general code is :

%clear all
%clc
dim =24; % define square matrix dimention
x = [];
coun = 1;
for i=1:dim:(dim.^2-dim+1)
    row = [i:i+dim-1];
    x = [x;row];
end
Input_Index = x;  % generate square matrix in order ascending index
% ==========================================
%        Generate Chaotic interleaver Scheme
%==========================================
%-------------------(1)-------------------------
part1 = [];
for i =1:2
    part1_i =x(dim:-1:dim/2 + 1,i)';
    part1 = [part1 part1_i];
end
part1;
%-------------------(2)-------------------------
part2 = [];
for i =1:2
    part2_i =x(dim/2:-1:1,i)';
    part2 = [part2 part2_i];
end
part2;
%-------------------(3)-------------------------
part3 = [];
for i =3:4
    part3_i =x(dim:-1:dim/2 + 1,i)';
    part3 = [part3 part3_i];
end
part3;
%-------------------(4)-------------------------
part4 = [];
for i =3:4
    part4_i =x(dim/2:-1:1,i)';
    part4 = [part4 part4_i];
end
part4;
%-------------------(5)-------------------------
part5 = [];
for i =5:6
    part5_i =x(dim:-1:dim/2 + 1,i)';
    part5 = [part5 part5_i];
end
part5;
%-------------------(6)-------------------------
part6 = [];
for i =5:6
    part6_i =x(dim/2:-1:1,i)';
    part6 = [part6 part6_i];
end
part6;
%-------------------(7)-------------------------
part7 = [];
for i =7:8
    part7_i =x(dim:-1:dim/2 + 1,i)';
    part7 = [part7 part7_i];
end
part7;
%-------------------(8)-------------------------
part8 = [];
for i =7:8
    part8_i =x(dim/2:-1:1,i)';
    part8 = [part8 part8_i];
end
part8;
%-------------------(9)-------------------------
part9 = [];
for i =9:12
    part9_i =x(dim:-1:dim/2 + dim/4+1,i)';
    part9= [part9 part9_i];
end
part9;
%-------------------(10)-------------------------
part10 = [];
for i =9:12
    part10_i =x(dim/2+dim/4:-1:dim/2+1 ,i)';
    part10= [part10 part10_i];
end
part10;
%-------------------(11)-------------------------
part11= [];
for i =9:12
    part11_i =x(dim/2:-1:dim/4+1 ,i)';
    part11= [part11 part11_i];
end
part11;
%-------------------(12)-------------------------
part12= [];
for i =9:12
    part12_i =x(dim/4:-1:1 ,i)';
    part12= [part12 part12_i];
end
part12;
%-------------------(13)-------------------------
part13 = [];
for i =13:16
    part13_i =x(dim:-1:dim/2 + dim/4+1,i)';
    part13= [part13 part13_i];
end
part13;
%-------------------(14)-------------------------
part14= [];
for i =13:16
    part14_i =x(dim/2+dim/4:-1:dim/2+1 ,i)';
    part14= [part14 part14_i];
end
part14;
%-------------------(15)-------------------------
part15= [];
for i =13:16
    part15_i =x(dim/2:-1:dim/4+1 ,i)';
    part15= [part15 part15_i];
end
part15;
%-------------------(16)-------------------------
part16= [];
for i =13:16
    part16_i =x(dim/4:-1:1 ,i)';
    part16= [part16 part16_i];
end
part16;
%-------------------(17)-------------------------
part17 = [];
for i =17:20
    part17_i =x(dim:-1:dim/2 + dim/4+1,i)';
    part17= [part17 part17_i];
end
part17;
%-------------------(18)-------------------------
part18= [];
for i =17:20
    part18_i =x(dim/2+dim/4:-1:dim/2+1 ,i)';
    part18= [part18 part18_i];
end
part18;
%-------------------(19)-------------------------
part19= [];
for i =17:20
    part19_i =x(dim/2:-1:dim/4+1 ,i)';
    part19= [part19 part19_i];
end
part19;
%-------------------(20)-------------------------
part20= [];
for i =17:20
    part20_i =x(dim/4:-1:1 ,i)';
    part20= [part20 part20_i];
end
part20;
%-------------------(21)-------------------------
part21 = [];
for i =21:24
    part21_i =x(dim:-1:dim/2 + dim/4+1,i)';
    part21= [part21 part21_i];
end
part21;
%-------------------(22)-------------------------
part22= [];
for i =21:24
    part22_i =x(dim/2+dim/4:-1:dim/2+1 ,i)';
    part22= [part22 part22_i];
end
part22;
%-------------------(23)-------------------------
part23= [];
for i =21:24
    part23_i =x(dim/2:-1:dim/4+1 ,i)';
    part23= [part23 part23_i];
end
part23;
%-------------------(24)-------------------------
part24= [];
for i =21:24
    part24_i =x(dim/4:-1:1 ,i)';
    part24= [part24 part24_i];
end
part24;
Chaotic_Interleaver_index = [part24;part23;part22;part21;part20;part19;part18;part17;part16;part15;part14;part13 ...
    ;part12;part11;part10;part9;part8;part7;part6;part5;part4;part3;part2;part1];
Chaotic_Interleaver_index_frame = [part24,part23,part22,part21,part20,part19,part18,part17,part16,part15,part14,part13 ...
    ,part12,part11,part10,part9,part8,part7,part6,part5,part4,part3,part2,part1];

6 个评论
显示 4更早的评论隐藏 4更早的评论

Rik 2022-1-1

If you don't use numbered variables, but instead use arrays, it looks like a for loop will solve your problem.

Zeina Abdullah 2022-1-2

thank you @Rik i will try but iam beginner in matlab

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Turlough Hughes 2022-1-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1619700-i-need-to-make-a-general-code#answer_866140

在 MATLAB Online 中打开

Try this:

function [idx,x] = chaoticInterleaver(N)
assert(mod(N,8)==0,'N must be divisible by 8.')
idx = zeros(N);
x = N * (0:N-1).' + (1:N); % from Jan
% Lower part
idx(N-7:2:end-1, :) = frf(x(1:N/2, 1:8), N);
idx(N-6:2:end,:)    = frf(x(N/2+1:end,1:8), N);
if N >= 16
    for ii = 1:4    % Upper part
        idx(ii:4:N-8,:) = frf( x(N/4*(ii-1)+1:N/4*ii, 9:end), N);
    end
end
end
function out = frf(partX,N)
% flipud, reshape with N rows, and then flip again (hence the name frf).
out = flipud(reshape(flipud(partX),N,[]).');
end

It may work for higher N too if the "upper part" rules are the same. The interleaver index frame is then:

interleaverIndexFrame = idx.';
interleaverIndexFrame = interleaverFrame(:).';

13 个评论
显示 11更早的评论隐藏 11更早的评论

Turlough Hughes 2022-1-2

编辑：Turlough Hughes 2022-1-2

在 MATLAB Online 中打开

chaoticInterleaver.m

It seems that you changed the function definition line to

function [idx,x] = chaoticInterleaverN

...

but it should be as follows:

function [idx,x] = chaoticInterleaver(N)

Simply save the first piece of code I provided in my answer as a function (called chaoticInterleaver.m), and then call it as follows:

[Chaotic_Interleaver_index, Address] = chaoticInterleaver(8) % you can use inputs 8, 16, 24, etc as required
Chaotic_Interleaver_index = 8×8
  23    15     7    32    24    16     8
  55    47    39    64    56    48    40
  21    13     5    30    22    14     6
  53    45    37    62    54    46    38
  19    11     3    28    20    12     4
  51    43    35    60    52    44    36
  17     9     1    26    18    10     2
  49    41    33    58    50    42    34
Address = 8×8
   2     3     4     5     6     7     8
  10    11    12    13    14    15    16
  18    19    20    21    22    23    24
  26    27    28    29    30    31    32
  34    35    36    37    38    39    40
  42    43    44    45    46    47    48
  50    51    52    53    54    55    56
  58    59    60    61    62    63    64

Turlough Hughes 2022-1-2

在 MATLAB Online 中打开

Perhaps it's best to explain my thought process with some figures.

First consider the case of N == 8. In this scenario you start with the matrix x as follows:

N = 8;
idx = zeros(N); % Space for Chaotic Interleaver Index
x = N * (0:N-1).' + (1:N);

The frf() function is applied to the top half and then the lower half of x. You can see how the frf() function rearranges the top half of x for the chaotic interleaver index, idx, in the following figure (blue parts show how two columns are reordered into a given row):

However, those output rows don’t go directly into idx, instead, they're organised into idx's odd numbered rows in the last 8 rows of idx:

>> idx(N-7:2:end-1, :) = frf(x(1:N/2, 1:8), N)
idx =
  23    15     7    32    24    16     8
   0     0     0     0     0     0     0
  21    13     5    30    22    14     6
   0     0     0     0     0     0     0
  19    11     3    28    20    12     4
   0     0     0     0     0     0     0
  17     9     1    26    18    10     2
   0     0     0     0     0     0     0

In the case of N = 8, the first 8 rows are the last 8 rows (this matters for N>=16). Following the same procedure, the lower half of x will fill the odd rows in idx. Ok... Thus far we have covered frf() and the following two lines of code from my original answer:

% Lower part
idx(N-7:2:end-1, :) = frf(x(1:N/2, 1:8), N);
idx(N-6:2:end,:)    = frf(x(N/2+1:end,1:8), N);

For the case of N==16, the process is the same on the left hand side of x:

One thing to remember is that all elements in the first 8 columns of x will end up on the bottom 8 rows of idx. The RHS is treated slightly differently. First, it's separated into 4 parts... let's consider the upper part (highlighted in red above). The procedure is similar except four columns in the input make up a single row in the output:

Each part on the RHS of x produces 2 rows, giving 8 rows in total, making the upper half of idx. Of course, they don’t go directly into idx, the two rows go into row 1 and 5 of idx, subsequently it's rows 2 and 6, rows 3 and 7, and finally rows 4 and 8. For the case of N == 24 its the same pattern, except the four parts on the right are 6 by 16 arrays and frf gives 4 rows that go into idx on rows 1 5 9 and 13, and so on.

Zeina Abdullah 2022-1-3

I mean your code is working very well and I am depend on it . But I have another step after this output. The 8×8 and 16×16 have group of numbers shared or existing . I need to do something like a relationship between these shared outputs . Hint the shared outputs not in the same address

Zeina Abdullah 2022-1-3

And thank you for your suggestion I will read this dismember function and I will try to understand it

请先登录，再进行评论。

Answer 2

yanqi liu 2022-1-1

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1619700-i-need-to-make-a-general-code#answer_866180

在 MATLAB Online 中打开

clc; clear all; close all;
%function const_interleaver(select)
%% constructing the orignal address matrix
Address=zeros(24); 
for i=1:24
    for j=1:24
        Address(i,j)=24*(i-1)+j;
    end
end
%% chosing the matrix size
bns = zeros(1,24);
if select==1 % 8-bit
    for i=1:8
        for j=1:8
            Address(i,j)=Address(i,j)-(i-1)*16;     
        end
    end
    % part 1,2,3,4,5,6,7 and part 8
    bns([1,2,3,4,5,6,7,8]) = 1;
elseif select==2 % 16-bit
    for i=1:16
        for j=1:16
            Address(i,j)=Address(i,j)-(i-1)*8;
        end
    end
    % part 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,and part 16
    bns([1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16]) = 1;
else select==3 % 24-bit
    % part1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23 and part 24 
    bns([1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]) = 1;
end 
% select==1 the part 1,2,3,4,5,6,7 and part 8 .
% select==2 the part 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,and part 16 
% select ==3 the part1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23 and part 24 
%clear all
%clc
dim =24; % define square matrix dimention
x = [];
coun = 1;
for i=1:dim:(dim.^2-dim+1)
    row = [i:i+dim-1];
    x = [x;row];
end
Input_Index = x;  % generate square matrix in order ascending index
% ==========================================
%        Generate Chaotic interleaver Scheme
%==========================================
if bns(1)
%-------------------(1)-------------------------
part1 = [];
for i =1:2
    part1_i =x(dim:-1:dim/2 + 1,i)';
    part1 = [part1 part1_i];
end
part1;
end
if bns(2)
%-------------------(2)-------------------------
part2 = [];
for i =1:2
    part2_i =x(dim/2:-1:1,i)';
    part2 = [part2 part2_i];
end
part2;
end
if bns(3)
%-------------------(3)-------------------------
part3 = [];
for i =3:4
    part3_i =x(dim:-1:dim/2 + 1,i)';
    part3 = [part3 part3_i];
end
part3;
end
if bns(4)
%-------------------(4)-------------------------
part4 = [];
for i =3:4
    part4_i =x(dim/2:-1:1,i)';
    part4 = [part4 part4_i];
end
part4;
end
if bns(5)
%-------------------(5)-------------------------
part5 = [];
for i =5:6
    part5_i =x(dim:-1:dim/2 + 1,i)';
    part5 = [part5 part5_i];
end
part5;
end
if bns(6)
%-------------------(6)-------------------------
part6 = [];
for i =5:6
    part6_i =x(dim/2:-1:1,i)';
    part6 = [part6 part6_i];
end
part6;
end
if bns(7)
%-------------------(7)-------------------------
part7 = [];
for i =7:8
    part7_i =x(dim:-1:dim/2 + 1,i)';
    part7 = [part7 part7_i];
end
part7;
end
if bns(8)
%-------------------(8)-------------------------
part8 = [];
for i =7:8
    part8_i =x(dim/2:-1:1,i)';
    part8 = [part8 part8_i];
end
part8;
end
if bns(9)
%-------------------(9)-------------------------
part9 = [];
for i =9:12
    part9_i =x(dim:-1:dim/2 + dim/4+1,i)';
    part9= [part9 part9_i];
end
part9;
end
if bns(10)
%-------------------(10)-------------------------
part10 = [];
for i =9:12
    part10_i =x(dim/2+dim/4:-1:dim/2+1 ,i)';
    part10= [part10 part10_i];
end
part10;
end
if bns(11)
%-------------------(11)-------------------------
part11= [];
for i =9:12
    part11_i =x(dim/2:-1:dim/4+1 ,i)';
    part11= [part11 part11_i];
end
part11;
end
if bns(12)
%-------------------(12)-------------------------
part12= [];
for i =9:12
    part12_i =x(dim/4:-1:1 ,i)';
    part12= [part12 part12_i];
end
part12;
end
if bns(13)
%-------------------(13)-------------------------
part13 = [];
for i =13:16
    part13_i =x(dim:-1:dim/2 + dim/4+1,i)';
    part13= [part13 part13_i];
end
part13;
end
if bns(14)
%-------------------(14)-------------------------
part14= [];
for i =13:16
    part14_i =x(dim/2+dim/4:-1:dim/2+1 ,i)';
    part14= [part14 part14_i];
end
part14;
end
if bns(15)
%-------------------(15)-------------------------
part15= [];
for i =13:16
    part15_i =x(dim/2:-1:dim/4+1 ,i)';
    part15= [part15 part15_i];
end
part15;
end
if bns(16)
%-------------------(16)-------------------------
part16= [];
for i =13:16
    part16_i =x(dim/4:-1:1 ,i)';
    part16= [part16 part16_i];
end
part16;
end
if bns(17)
%-------------------(17)-------------------------
part17 = [];
for i =17:20
    part17_i =x(dim:-1:dim/2 + dim/4+1,i)';
    part17= [part17 part17_i];
end
part17;
end
if bns(18)
%-------------------(18)-------------------------
part18= [];
for i =17:20
    part18_i =x(dim/2+dim/4:-1:dim/2+1 ,i)';
    part18= [part18 part18_i];
end
part18;
end
if bns(19)
%-------------------(19)-------------------------
part19= [];
for i =17:20
    part19_i =x(dim/2:-1:dim/4+1 ,i)';
    part19= [part19 part19_i];
end
part19;
end
if bns(20)
%-------------------(20)-------------------------
part20= [];
for i =17:20
    part20_i =x(dim/4:-1:1 ,i)';
    part20= [part20 part20_i];
end
part20;
end
if bns(21)
%-------------------(21)-------------------------
part21 = [];
for i =21:24
    part21_i =x(dim:-1:dim/2 + dim/4+1,i)';
    part21= [part21 part21_i];
end
part21;
end
if bns(22)
%-------------------(22)-------------------------
part22= [];
for i =21:24
    part22_i =x(dim/2+dim/4:-1:dim/2+1 ,i)';
    part22= [part22 part22_i];
end
part22;
end
if bns(23)
%-------------------(23)-------------------------
part23= [];
for i =21:24
    part23_i =x(dim/2:-1:dim/4+1 ,i)';
    part23= [part23 part23_i];
end
part23;
end
if bns(24)
%-------------------(24)-------------------------
part24= [];
for i =21:24
    part24_i =x(dim/4:-1:1 ,i)';
    part24= [part24 part24_i];
end
part24;
end
Chaotic_Interleaver_index = [part24;part23;part22;part21;part20;part19;part18;part17;part16;part15;part14;part13 ...
    ;part12;part11;part10;part9;part8;part7;part6;part5;part4;part3;part2;part1];
Chaotic_Interleaver_index_frame = [part24,part23,part22,part21,part20,part19,part18,part17,part16,part15,part14,part13 ...
    ,part12,part11,part10,part9,part8,part7,part6,part5,part4,part3,part2,part1];

4 个评论
显示 2更早的评论隐藏 2更早的评论

yanqi liu 2022-1-2

yes，sir，it is one file

thanks

Zeina Abdullah 2022-1-2

ok thank you @yanqi liu

请先登录，再进行评论。

I need to make a general code

6 个评论
显示 4更早的评论隐藏 4更早的评论

采纳的回答

13 个评论
显示 11更早的评论隐藏 11更早的评论

更多回答（1 个）

4 个评论
显示 2更早的评论隐藏 2更早的评论

另请参阅

类别

标签

Community Treasure Hunt

I need to make a general code

6 个评论 显示 4更早的评论隐藏 4更早的评论

采纳的回答

13 个评论 显示 11更早的评论隐藏 11更早的评论

更多回答（1 个）

4 个评论 显示 2更早的评论隐藏 2更早的评论

另请参阅

类别

标签

Community Treasure Hunt

6 个评论
显示 4更早的评论隐藏 4更早的评论

13 个评论
显示 11更早的评论隐藏 11更早的评论

4 个评论
显示 2更早的评论隐藏 2更早的评论