zhentou.m

64路OFDM系统MATLAB仿真

function [all]=zhentou(Snr,raylei)



%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% t1~t10  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%    all=zeros(1,zhenshu*80+320);

    for jjj=1:2

yuanshixinxi=[0 0 1+i 0 0 0 -1-i 0 0 0 1+i 0 0 0 -1-i 0 0 0 -1-i 0 0 0 1+i 0 0 0 0 0 0 0 -1-i 0 0 0 -1-j 0 0 0 1+i 0 0 0 1+i 0 0 0 1+i 0 0 0 1+i 0 0 ]';


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% IFFT  %%%%%%%%%%%%%%%%%%%%%%%%%
iffthou=ifft(yuanshixinxi,64);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 抑制峰均比处理  %%%%%%%%%%%%%%%%%%%%%%%%%%
iffthou(17:32)=-iffthou(17:32);
iffthou(49:64)=-iffthou(49:64);


[paprchuli]=zhijiexuebo(iffthou);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [16 0],16列0行,添加CP %%%%%%%%%%%%%%%%%%%%
addcp= padarray(paprchuli,[16 0],'circular','pre');
walsh_code=[1 -1 1 -1 -1 1 -1 1 1 -1 1 -1 -1 1 -1 1]';
%for kkk=1:5
    
%addcp(kkk*16-15:kkk*16)=walsh_code; 
%end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 经过高斯信道 %%%%%%%%%%%%%%%%%%%%%%%%%

passawgnchannel= awgn(addcp ,Snr,'measured');
%passawgnchannel= addcp ;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 经过瑞利多径信道 %%%%%%%%%%%%%%%%%%%%%%%%%                                                   
%passchannel = filter(raylei, passawgnchannel); 
passchannel =passawgnchannel;

all(jjj*80-80+1:jjj*80)=passchannel';


    end
    
    
    
    
%%    
    
    
    
    

    %%%%%%%%%%%%%%%%%%%%%%%%%  T1  T2  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    
    lt1t2=[1 1 -1 -1 1 1 -1 1 -1 1 1 1 1 1 1 -1 -1 1 1 ...
            -1 1 -1 1 1 1 1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 -1 1 1 -1 1 -1 1 -1 ...
            -1 -1 -1 -1 1 1 -1 -1 1 -1 1 -1 1 1 1 1 ];
      

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% IFFT  %%%%%%%%%%%%%%%%%%%%%%%%%
iffthou=ifft(lt1t2,64);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 抑制峰均比处理  %%%%%%%%%%%%%%%%%%%%%%%%%%

[paprchuli]=zhijiexuebo(iffthou);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 经过高斯信道 %%%%%%%%%%%%%%%%%%%%%%%%%

passawgnchannel= awgn(paprchuli ,Snr,'measured');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 经过瑞利多径信道 %%%%%%%%%%%%%%%%%%%%%%%%%                                                   
passchannel = filter(raylei, passawgnchannel); 
passchannel =passawgnchannel;

all(193:256)=passchannel';
%all(193:256)=paprchuli;%%%%%%%%%%%%%%%%%   完全去除噪声   %%%%%%%%%%%%%%%%%%%%%%%%%%%%
            %%%%%%%%%T2%%%%%%%%  
 
iffthou=ifft(yuanshixinxi,64);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 抑制峰均比处理  %%%%%%%%%%%%%%%%%%%%%%%%%%

[paprchuli]=zhijiexuebo(iffthou);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 经过高斯信道 %%%%%%%%%%%%%%%%%%%%%%%%%

passawgnchannel= awgn(paprchuli ,Snr,'measured');
%passawgnchannel= addcp ;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 经过瑞利多径信道 %%%%%%%%%%%%%%%%%%%%%%%%%                                                   
passchannel = filter(raylei, passawgnchannel); 
%passchannel =passawgnchannel;

all(257:320)=passchannel';


all(161:192)=passchannel(33:64);