How to combine 340 sub matrices of 117 by 1 into one matrix of 9945 by 1. I'm working on ofdm(subcarriers=128, cp length=32) and got 117 data carriers. That's why ofdm requires input matrix of 117 by 1

Question

Sajid Sarwar 2018-11-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/430567-how-to-combine-340-sub-matrices-of-117-by-1-into-one-matrix-of-9945-by-1-i-m-working-on-ofdm-subcar

编辑： Stephen23 2018-11-19

kindly tell me how to combine these subsets in one matrix. for example I'm working on following lines.

st=rand(9945,4)
for n=1:85
 Baserow = (n-1) * 117;
for n1 = 1:4
Subset = st(Baserow+(1:117),n1, :)
txSig = ofdmMod(Subset);                      
>> y = awgn(txSig, SNR(k));
>> sigRx = ofdmDemod(y);
end
end

I'm getting 340 sub matrices of "Sigtx" of 117 by 1. How to combine them in a single matrix?

thanku.

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

Sajid Sarwar 2018-11-18

在 MATLAB Online 中打开

madhan ravi sir,

M = 4;                 % Modulation alphabet
k = log2(M);           % Bits/symbol
numSC = 128;           % Number of OFDM subcarriers
cpLen = 32;            % OFDM cyclic prefix length
ofdmMod = comm.OFDMModulator('FFTLength',numSC,'CyclicPrefixLength',cpLen);
ofdmDemod =        comm.OFDMDemodulator('FFTLength',numSC,'CyclicPrefixLength',cpLen);
channel = comm.AWGNChannel('NoiseMethod','Variance', ...
    'VarianceSource','Input port');
ofdmDims = info(ofdmMod)
numDC = ofdmDims.DataInputSize(1)

Stephen23 2018-11-19

@Sajid Sarwar: please do NOT post the same question multiple times. It does not make it faster or easier for you to get help on this forum.

请先登录，再进行评论。

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

Answer 1

Stephen23 2018-11-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/430567-how-to-combine-340-sub-matrices-of-117-by-1-into-one-matrix-of-9945-by-1-i-m-working-on-ofdm-subcar#answer_347646

编辑：Stephen23 2018-11-18

Add these three lines to your code:

C = cell(4,85);
...
for ii = 1:85
    ...
    for jj = 1:4
        sigRx = ...
        C{jj,ii} = sigRx;
    end
end
V = vertcat(C{:})

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

Sajid Sarwar 2018-11-19

在 MATLAB Online 中打开

CB(:,:,1) = [ 0                  0                  0                  0;...
             -0.1815-1j*0.1318  -0.6351-1j*0.4615   0.6351+1j*0.4615   0.1815+1j*0.1318;...
              0                  0                  0                  0;...
              0.7851            -0.2243             0.2243            -0.7851 ];
CB(:,:,2) = [ 0.7851            -0.2243             0.2243            -0.7851;...
              0                  0                  0                  0;...
             -0.1815-1j*0.1318  -0.6351-1j*0.4615   0.6351+1j*0.4615   0.1815+1j*0.1318;...
              0                  0                  0                  0 ];
CB(:,:,3) = [-0.6351+1j*0.4615   0.1815-1j*0.1318  -0.1815+1j*0.1318   0.6351-1j*0.4615;...
              0.1392-1j*0.1759   0.4873-1j*0.6156  -0.4873+1j*0.6156  -0.1392+1j*0.1759;...
              0                  0                  0                  0;...
              0                  0                  0                  0 ];
CB(:,:,4) = [ 0                  0                  0                  0;...
              0                  0                  0                  0;...
              0.7851            -0.2243             0.2243            -0.7851;...
             -0.0055-1j*0.2242  -0.0193-1j*0.7848   0.0193+1j*0.7848   0.0055+1j*0.2242 ];
CB(:,:,5) = [-0.0055-1j*0.2242  -0.0193-1j*0.7848   0.0193+1j*0.7848   0.0055+1j*0.2242;...
              0                  0                  0                  0;...
              0                  0                  0                  0;...
             -0.6351+1j*0.4615   0.1815-1j*0.1318  -0.1815+1j*0.1318   0.6351-1j*0.4615 ];
CB(:,:,6) = [ 0                  0                  0                  0;...
              0.7851            -0.2243             0.2243            -0.7851;...
              0.1392-1j*0.1759   0.4873-1j*0.6156  -0.4873+1j*0.6156  -0.1392+1j*0.1759;...
              0                  0                  0                  0 ];
M = 4;                 % Modulation alphabet
k = log2(M);           % Bits/symbol
numSC = 128;           % Number of OFDM subcarriers
cpLen = 32;            % OFDM cyclic prefix length
ofdmMod = comm.OFDMModulator('FFTLength',numSC,'CyclicPrefixLength',cpLen);
ofdmDemod =        comm.OFDMDemodulator('FFTLength',numSC,'CyclicPrefixLength',cpLen);
channel = comm.AWGNChannel('NoiseMethod','Variance', ...
    'VarianceSource','Input port');
errorRate = comm.ErrorRate('ResetInputPort',true);
ofdmDims = info(ofdmMod)
numDC = ofdmDims.DataInputSize(1)
K = size(CB, 1); % number of orthogonal resources
M = size(CB, 2); % number of codewords in each codebook
V = size(CB, 3); % number of users (layers)
N = 9945; % SCMA signals in frame
EbN0 = 0:20;
SNR  = EbN0 + 10*log10(log2(M)*V/K)+ 10*log10(numDC/numSC);
errorStats  = zeros(V, length(SNR));
Nbits = zeros(V, length(SNR));
BER   = zeros(V, length(SNR));
maxNumErrs = 100;
maxNumBits = 1e7;
Niter      = 10;
for k = 1:length(SNR)
    N0 = 1/(10^(SNR(k)/10)); % noise power
    while ((min(errorStats (:,k)) < maxNumErrs) && (Nbits(1,k) < maxNumBits))
x = randi([0 M-1], V, N); 
 K=4;
h = 1/sqrt(2)*(randn(K, V, N)+1j*randn(K, V, N));
s = scmaenc(x, CB, h);
>> st=transpose(s);
for n=1:85
 Baserow = (n-1) * 117;
for n1 = 1:4
Subset = st(Baserow+(1:117),n1, :)
txSig = ofdmMod(Subset)                    
>> y = awgn(txSig, SNR(k));
>> sigRx = ofdmDemod(y)
end
end
LLR = scmadec(sigRx, CB, h, N0, Niter);
        % symbol to bit conversion
        r    = de2bi(x, log2(M), 'left-msb');
        data = zeros(log2(M)*N, V);
        for kk = 1:V
            data(:,kk) = reshape(downsample(r, V, kk-1).',[],1);
        end
        % LLR to bit conversion
        datadec = reshape((LLR <= 0), [log2(M) N*V]).';
        datar   = zeros(log2(M)*N, V);
        for kk = 1:V
            datar(:,kk) = reshape(downsample(datadec, V, kk-1).', [], 1);
        end
    
%BER(k,:) = errorStats;
%errorStats = errorRate(x, sigRx,1);
err        = sum(xor(data, datar));
        errorStats (:,k)  = errorStats (:,k) + err.';
        Nbits(:,k) = Nbits(:,k) + log2(M)*N;
end
    BER(:,k) = errorStats (:,k)./Nbits(:,k);
    k
end
figure
semilogy(EbN0,BER)
% this is compelte code whereas there are separe files for scmadec, scmaenc and log_sum_exp

Sajid Sarwar 2018-11-19

在 MATLAB Online 中打开

Stephen Cobeldick sir, "combine 340 sub matrices of 117 by 1 into one matrix of 9945 by 4" as follows:

st=rand(9945,4)
for n=1:85
 Baserow = (n-1) * 117;
for n1 = 1:4
Subset = st(Baserow+(1:117),n1, :)
txSig = ofdmMod(Subset);                      
>> y = awgn(txSig, SNR(k));
>> sigRx = ofdmDemod(y);
end
end

then output of sigRx(340 matrices of 117 by 1=9945 by 4

(a_1,1 a_1,2 a_1,3 a_1,4

uptill uptill uptill uptill

a_117,1 a_117,2 a_117,3 a_117,4

uptill uptill uptill uptill

a_9945,1 a_9945,2 a_9945,3 a_9945,4)

i have tried to write a matrix of 9945 by 4, which I required

Stephen23 2018-11-19

编辑：Stephen23 2018-11-19

@Sajid Sarwar: why are you just repeating your comments? So far you have totally ignored anything that I have written: your code does not show any attempt to use what I have shown you in my answer or my comments, so I am not quite sure how I can help you, if you just totally ignore what I write.

请先登录，再进行评论。

How to combine 340 sub matrices of 117 by 1 into one matrix of 9945 by 1. I'm working on ofdm(subcarriers=128, cp length=32) and got 117 data carriers. That's why ofdm requires input matrix of 117 by 1

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

回答（1 个）

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

另请参阅

类别

标签

Community Treasure Hunt

How to combine 340 sub matrices of 117 by 1 into one matrix of 9945 by 1. I'm working on ofdm(subcarriers=128, cp length=32) and got 117 data carriers. That's why ofdm requires input matrix of 117 by 1

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

回答（1 个）

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

另请参阅

类别

标签

Community Treasure Hunt

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

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