function y=ray_channel(x,T)
tao= [0.0,0.2,0.5,1.6,2.3,5.0];
tao_ord = tao*1e-6;
g = [0.7079+0.4*i,1.0000+i,0.7943+0.3*i,0.5012,0.3981,0.3162];
    tx=-256*T:T:(1024-1)*T;                                                          %tx為未過信道的信號的時間參量
    for j=1:1024
        t=(j-1)*T;                                                                   %t為過信道后的信號的時間參量
        x_ttao=interp1(tx,x,t-tao_ord').';                                           %對x插值獲得t-tao時刻的x值
        y(j)=sum(g.*x_ttao);                                                         %由插值獲得的20個x值獲得過信道后的信號在(j-1)T時刻的值
    end