through_channel.m

来自「信号处理」· M 代码 · 共 29 行

function rx_signal = Through_channel(c_end,tx_signal,signal_len)

% global sim_consts;
global cir_matrix_t 
global flag


cir_matrix = [];
c_end = c_end/norm(c_end);
 C2=c_end;
cir_len = length(c_end);
    cir_matrix = [cir_matrix C2];
    for k=1:signal_len+cir_len-1
        cir_matrix = [cir_matrix C2+(c_end-C2)/(2*(signal_len+cir_len-1))];
    end 
cir_matrix_t = cir_matrix;

[ cir_len , conv_len] = size(cir_matrix) ;
tx_temp=zeros(1,cir_len);
 for k=1:conv_len
        if k<signal_len+1
            tx_temp = [tx_signal(1,k) tx_temp(1,1:end-1)];
        else tx_temp = [0 tx_temp(1,1:end-1)];
        end
        cir = cir_matrix(:,k);
        cir = cir/norm(cir)';
        rx_signal(1,k) = tx_temp * cir;
 end