% freq interleaver
% 先求得freq index 'k'
%注意，matlab矩陣元素標(biāo)號(hào)是從1開(kāi)始的，所以共為1-2048個(gè)點(diǎn)
temp(1)=0;
for i_f=2:2048
temp(i_f)=mod((13*temp(i_f-1)+511),2048);
i_f=i_f+1;
end

%產(chǎn)生dn
i_f=1;
j_f=1;
for i_f=1:2048
    if (temp(i_f)>=256)&(temp(i_f)<=1023)
    dn(j_f)=temp(i_f);
    j_f=j_f+1;
elseif (temp(i_f)>=1025)&(temp(i_f)<=1792)
    dn(j_f)=temp(i_f);
    j_f=j_f+1;
    end
    i_f=i_f+1;
end

%產(chǎn)生freq index 'k'
for i_k=1:1536
    k(i_k)=dn(i_k)-255;
    if k(i_k)>769
        k(i_k)=k(i_k)-1;
    end
    i_k=i_k+1;
end
save('k.mat','k');

%現(xiàn)在根據(jù)已知的k，找尋一種更簡(jiǎn)單的freq interleaving的方法
%matrix_k是正變換陣，根據(jù)k中的數(shù)據(jù)，對(duì)1536*1536的單位陣施加相應(yīng)列變換得到
%將待處理的數(shù)據(jù)右乘matrix_k，根據(jù)線性代數(shù)的有關(guān)理論，即得插值后的新數(shù)據(jù)
%逆變換時(shí)只需再右剩matrix_k的逆矩陣即可

%下面這一塊，只要運(yùn)行一次就好

%先定義空矩陣matrix_k和1536階單位陣matrix_eye
matrix_k=[];
matrix_eye=eye(1536);
%根據(jù)插值表k，對(duì)matrix_eye施加相應(yīng)列變換，構(gòu)造matrix_k，即為頻率插值正變換陣
for i_matrixk=1:1536
   matrix_k=[matrix_k,matrix_eye(:,k(i_matrixk))]; 
end
%求matrix_k的逆矩陣：matrix_k_inv，用來(lái)做頻率插值的恢復(fù):de_interleaving
matrix_k_inv=inv(matrix_k);