function output_sig=multipath_chann(input_sig,num,var_pow,delay,fd,t_interval,counter,count_begin)
%input_sig輸入信號(hào)矩陣,加了cp后的信號(hào)，大小為NL×(子載波個(gè)數(shù)+cp長(zhǎng)度lp)；
%num多徑數(shù);
%var_pow各徑相對(duì)主徑的平均功率,單位dB；
%delay各徑延時(shí),單位s；
%fd最大dopple頻率；
%t_interval為離散信道抽樣時(shí)間間隔，等于OFDM符號(hào)長(zhǎng)度/(子載波個(gè)數(shù)+cp長(zhǎng)度lp)；
%output_sig為經(jīng)過多徑信道的輸出信號(hào)矢量
%counter各徑間隔記錄
%count_begin本次產(chǎn)生信道開始記錄的初始位置

t_shift=floor(delay/t_interval);%歸一化各徑延時(shí)
%theta_shift=2*pi*fc*delay;
[nl,l]=size(input_sig);
output_sig=zeros(size(input_sig));

chann_l=nl*l;%信道采樣點(diǎn)數(shù)，若一個(gè)調(diào)制符號(hào)采樣一個(gè)信道點(diǎn)，則信道采樣點(diǎn)數(shù)等于輸入信號(hào)中的調(diào)制符號(hào)個(gè)數(shù)
selec_ray_chan=zeros(num,chann_l);%初始化頻率選擇性信道，徑數(shù)＝num
pow_per_channel=10.^(var_pow/10);%各徑功率線性化，從dB轉(zhuǎn)變成線性
total_pow_allchan=sum(pow_per_channel);%各徑功率之和
%以下for循環(huán)產(chǎn)生相互獨(dú)立的num條rayleigh信道
for k=1:num
    atts=sqrt(pow_per_channel(k));
    selec_ray_chan(k,:)=atts*rayleighnew(chann_l,t_interval,fd,count_begin+k*counter)/sqrt(total_pow_allchan);
end
for k=1:l
    input_sig_serial(((k-1)*nl+1):k*nl)=input_sig(:,k).';%輸入信號(hào)矩陣轉(zhuǎn)變成串行序列
end
delay_sig=zeros(num,chann_l);%初始化延時(shí)后的送入各徑的信號(hào)，每徑所含符號(hào)數(shù)為chann_l
%以下for循環(huán)為各徑的輸入信號(hào)做延遲處理
for f=1:num
    if t_shift(f)~=0
        delay_sig(f,1:t_shift(f))=zeros(1,t_shift(f));
    end
    delay_sig(f,(t_shift(f)+1):chann_l)= input_sig_serial(1:(chann_l-t_shift(f)));
end
output_sig_serial=zeros(1,chann_l);%初始化輸出信號(hào)串行序列
%得到各徑疊加后的輸出信號(hào)序列
for f=1:num
        output_sig_serial= output_sig_serial+selec_ray_chan(f,:).*delay_sig(f,:);
end
for k=1:l
    output_sig(:,k)=output_sig_serial(((k-1)*nl+1):k*nl).';%輸出信號(hào)串行序列轉(zhuǎn)變成與輸入信號(hào)相同的矩陣形式，做為本函數(shù)輸出
end
%注意，在本函數(shù)中沒有為信號(hào)疊加白噪聲