ls_dft.m

程序里面包含各部分算法仿真 可供参考和使用

function [H]=LS_DFT(N_subc,H,L_delay) 


time_domain =  ifft(fftshift(H,1)) * sqrt(N_subc);
        
        % 保留一段时域样点
         
            leakage = 0;       % 能量泄漏的子载波数
            % 保留时域能量集中的L_delay个样点
            time_cut( 1: L_delay,: ,:) = time_domain( 1: L_delay,:,: );
            % 保留能量泄漏的leakage个样点
            time_cut( N_subc - leakage : N_subc,: ,:) = time_domain( N_subc - leakage : N_subc,: ,:);
            
  
        % 变换到频域
        H = fftshift(fft(time_cut) * 1/sqrt(N_subc),1);