【通信】基于matlab多径衰落信道仿真【含Matlab源码 338期】

原创

Matlab科研 2021-08-07 09:33:04 博主文章分类：matlab ©著作权

©著作权归作者所有：来自51CTO博客作者Matlab科研的原创作品，请联系作者获取转载授权，否则将追究法律责任

一、简介

衰落信道参数包括多径扩展和多普勒扩展。时不变的多径扩展相当于一个延时抽头滤波器，而多普勒扩展要注意多普勒功率谱密度，通常使用Jakes功率谱、高斯、均匀功率谱。

多径衰落信道由单径信道叠加而成，而单径信道中最重要的就是瑞利（Rayleigh）平坦衰落信道。

二、源代码

clear all;

a=[0.6154 0.7919 0.9218 0.7382 0.1763 0.4057 0.9355];%每路传输的幅度

t=[0.9169 0.4103 0.8936 0.0579 0.3529 0.8132 0.0099];%每路传输的到达时间
i=0; 
for omega=0:0.05:100; 
multipath_arrival=a.*exp(j*omega*t);
i=i+1;
abs_H(i)=abs(sum(multipath_arrival));%第i路信道的传输函数
end

clear all
ht=10;
hr=2;
c=3*10^8;
f=1*10^9;
lambda=c/f;
R=-1;
d=1:0.5:10000;
d1=sqrt(d.^2+(ht-hr)^2);
d2=sqrt(d.^2+(ht+hr)^2);
a1=exp(j*2*pi.*d1/lambda)./d1;
a2=R*exp(j*2*pi.*d2/lambda)./d2;
a=abs(a1+a2);
ld=log10(d);
la=log10(a);

clear all
ht=10;hr=2;
c=3e8;R=-1; f0=1e8; fi= [1:1:1000];fd=5000000;f= f0+fd*fi; lambda=c./f;
da=[50,300,800,2000];
for i=1:length(da)
d=da(i);
d1=sqrt(d.^2+(ht-hr)^2);
d2=sqrt(d.^2+(ht+hr)^2);
Td=(d2-d1)/c
a1=exp(j*2*pi*d1./lambda)/d1;
a2=R*exp(j*2*pi*d2./lambda)/d2;
a(i,:)=abs(a1+a2);
end
clear all;
an=[1,0.3,-0.8,0.5,-0.4,0.2];tn=[0,1,2,3,4,5;0,0.1,0.2,0.3,0.4,0.5];
signal=[0, zeros(1,0),ones(1,501),zeros(1,1000)]; % transmitted signal
for k=1:2; %for two cases
for i=1:6;
ray(i,:)=an(i)*[0, zeros(1,(100*tn(k,i))),ones(1,501),zeros(1,(1000-100*tn(k,i)))];
end
y(k,:)=sum(ray(:,1:end));
end
t=((1:1:length(y(1,:)))-1)*10^(-2);
subplot(2,2,1);plot(t,signal);
ylabel('transmitted signal s(t)');
title('case 1 & case 2')
axis([ 0 20 -0.5 1.5])
subplot(2,2,2);plot(t,y(1,:));
ylabel('received signal y(t)');