Matlab仿真多径信道瑞利衰落

%主函数
%两径瑞利衰落信道仿真
%设定默认参数
NN=256; %传输符号个数
tb=0.5; %一个比特时间
fs=10; %每符号采样数
ebn0db=[1:2:15];%设定Eb/N0

%建立QPSK信号
x=random_binary(NN,fs)+i*random_binary(NN,fs); %x为QPSK信号

%输入功率和延迟
p0=0; %视距LOS分量
p1=20; %第一路径分量
p2=1; %第二路径分量
delay=1; %按照每符号采样数决定的延迟
delay0=0;
delay1=0;
delay2=delay;

%设定复高斯（瑞利）衰减
gain1=sqrt(p1)abs(randn(1,NN)+irandn(1,NN));
gain2=sqrt(p2)abs(randn(1,NN)+irandn(1,NN));
for k=1:NN
for kk=1:fs
index=(k-1)*fs+kk;
ggain1(1,index)=gain1(1,k);
ggain2(1,index)=gain2(1,k);
end
end
y1=x;
for k=1:delay2
y2(1,k)=y1(1,k)sqrt(p0);
end
for k=(delay2+1):(NNfs)
y2(1,k)=y1(1,k)*sqrt(p0)+y1(1,k-delay1)*ggain1(1,k)+y1(1,k-delay2)*ggain2(1,k);
end

%匹配滤波器
b=-ones(1,fs);
b=b/fs;
a=1;
y=filter(b,a,y2);

[cor lags]=vxcorr(x,y);
cmax=max(max(abs(cor)));
nmax=find(abs(cor)==cmax);
timelag=lag