h_interp2order.m

来自「代码为对OFDM符号信道估计结果进行线性插值的仿真代码」· M 代码 · 共 34 行

function channels=H_interp2order(pilotin)
len_pilot=length(pilotin);
rate=8;
N=rate*len_pilot/30;
od=zeros(length(pilotin),1);
fake_pilot=pilotin(1)-(pilotin(2)-pilotin(1)); % create pilotin(0) according to other pilots by linear for the purpose of 1st interp
for m=1:len_pilot,
    if m==1
        for i=1:(rate-1) 
            a=i/N;
            c1=a*(a-1)/2;
            c0=-(a-1)*(a+1);
            c_1=a*(a+1)/2;
            channels(i)=c0*fake_pilot+c_1*pilotin(1);
        end
    elseif m==2
         for i=1:(rate-1) 
            a=i/N;
            c1=a*(a-1)/2;
            c0=-(a-1)*(a+1);
            c_1=a*(a+1)/2;
            channels((m-1)*(rate-1)+i)=c1*fake_pilot+c0*pilotin(1)+c_1*pilotin(2);
        end
    else
        for i=1:(rate-1),
             a=i/N;
             c1=a*(a-1)/2;
             c0=-(a-1)*(a+1);
             c_1=a*(a+1)/2;
             channels((m-1)*(rate-1)+i)=c1*pilotin(m-2)+c0*pilotin(m-1)+c_1*pilotin(m);
         end
     end
 end
channels=conj(channels);