《通信原理》多径衰落信道仿真1

一、实验内容

利用matlab实现不同条件下多径衰落信道的信号。

（1）根据已知条件移动台距离基站初始距离r0=1计算直射径信号和反射径信号并画出其图形和总接收信号；

（2）其他条件不变时，当移动台距离基站初始距离r0=3计算直射径信号和反射径信号并画出其图形和总接收信号；

（3）其他条件不变时，当移动台距离基站初始距离r0=6计算直射径信号和反射径信号并画出其图形和总接收信号；

（4）其他条件不变时，当移动台距离基站初始距离r0=9计算直射径信号和反射径信号并画出其图形和总接收信号；

（5）其他条件不变时，当f=10^8,r0=3计算直射径信号和反射径信号并画出其图形和总接收信号；

（6）其他条件不变时，当v=1,f=2,r0=3,d=15,t1=0.1:0.001:12计算直射径信号和反射径信号并画出其图形和总接收信号。

二、实验结果

 1.实验代码

%多径衰落信道仿真1
 
%其他条件不变，r0=1时的情况
f=1;    %发射信号频谱。
v=0;    %移动台速度，静止情况为0。
c=3d8;  %电磁波速度，光速。
r0=1;   %移动台距离基站初始距离。
d=10;   %基站距离反射墙的距离。
t1=0.1:0.0001:10;   %时间。
E1=cos(2*pi*f*((1-v/c).*t1-r0/c))./(r0+v.*t1);  %直射径信号。
E2=cos(2*pi*f*((1+v/c).*t1+(r0-2*d)/c))./(2*d-r0-v*t1);  %反射径信号。
figure('name','条件变化时的接收信号图','position',[0,40,1000,500]);
subplot(2,4,1);plot(t1,E1,'g-',t1,E2,'b-',t1,E1-E2,'r-');     %画出直射径、反射径和总的接收信号。
xlabel('移动台距离基站的初始距离r0','fontname','黑体');
ylabel('直射径/反射径/总接收信号','fontname','黑体');
title('r0=1的信号图','fontname','黑体');
legend('直射径信号','反射径信号','移动台接受的合成信号','location','b');
axis([0,10,-1.5,2.5]);
 
%其他条件不变，r0=3时的情况
f=1;    %发射信号频谱。
v=0;    %移动台速度，静止情况为0。
c=3d8;  %电磁波速度，光速。
r0=3;   %移动台距离基站初始距离。
d=10;   %基站距离反射墙的距离。
t1=0.1:0.0001:10;   %时间。
E1=cos(2*pi*f*((1-v/c).*t1-r0/c))./(r0+v.*t1);  %直射径信号。
E2=cos(2*pi*f*((1+v/c).*t1+(r0-2*d)/c))./(2*d-r0-v*t1);  %反射径信号。
subplot(2,4,2);plot(t1,E1,'g-',t1,E2,'b-',t1,E1-E2,'r-');     %画出直射径、反射径和总的接收信号。
xlabel('移动台距离基站的初始距离r0','fontname','黑体');
ylabel('直射径/反射径/总接收信号','fontname','黑体');
title('r0=3的信号图','fontname','黑体');
axis([0,10,-2,2]);
 
%其他条件不变，r0=6时的情况
f=1;    %发射信号频谱。
v=0;    %移动台速度，静止情况为0。
c=3d8;  %电磁波速度，光速。
r0=6;   %移动台距离基站初始距离。
d=10;   %基站距离反射墙的距离。
t1=0.1:0.0001:10;   %时间。
E1=cos(2*pi*f*((1-v/c).*t1-r0/c))./(r0+v.*t1);  %直射径信号。
E2=cos(2*pi*f*((1+v/c).*t1+(r0-2*d)/c))./(2*d-r0-v*t1);  %反射径信号。
subplot(2,4,3);plot(t1,E1,'g-',t1,E2,'b-',t1,E1-E2,'r-');     %画出直射径、反射径和总的接收信号。
xlabel('移动台距离基站的初始距离r0','fontname','黑体');
ylabel('直射径/反射径/总接收信号','fontname','黑体');
title('r0=6的信号图','fontname','黑体');
axis([0,10,-2,2]);
 
%其他条件不变，r0=9时的情况
f=1;    %发射信号频谱。
v=0;    %移动台速度，静止情况为0。
c=3d8;  %电磁波速度，光速。
r0=9;   %移动台距离基站初始距离。
d=10;   %基站距离反射墙的距离。
t1=0.1:0.0001:10;   %时间。
E1=cos(2*pi*f*((1-v/c).*t1-r0/c))./(r0+v.*t1);  %直射径信号。
E2=cos(2*pi*f*((1+v/c).*t1+(r0-2*d)/c))./(2*d-r0-v*t1);  %反射径信号。
subplot(2,4,4);plot(t1,E1,'g-',t1,E2,'b-',t1,E1-E2,'r-');     %画出直射径、反射径和总的接收信号。
xlabel('移动台距离基站的初始距离r0','fontname','黑体');
ylabel('直射径/反射径/总接收信号','fontname','黑体');
title('r0=9的信号图','fontname','黑体');
axis([0,10,-2,2]);
 
%其他条件不变，f=10^8,r0=3时的情况
f=10.^8;    %发射信号频谱。
v=0;    %移动台速度，静止情况为0。
c=3d8;  %电磁波速度，光速。
r0=3;   %移动台距离基站初始距离。
d=10;   %基站距离反射墙的距离。
t1=0.1:0.0001:10;   %时间。
E1=cos(2*pi*f*((1-v/c).*t1-r0/c))./(r0+v.*t1);  %直射径信号。
E2=cos(2*pi*f*((1+v/c).*t1+(r0-2*d)/c))./(2*d-r0-v*t1);  %反射径信号。
subplot(2,4,5);plot(t1,E1,'g-',t1,E2,'b-',t1,E1-E2,'r-');     %画出直射径、反射径和总的接收信号。
xlabel('移动台距离基站的初始距离r0','fontname','黑体');
ylabel('直射径/反射径/总接收信号','fontname','黑体');
title('r0=3,f=10^8的信号图','fontname','黑体');
axis([0,10,-2,2]);
 
%v=1,f=2,r0=3,d=15,t1=0.1:0.001:12时的情况
f=2;    %发射信号频谱。
v=1;    %移动台速度，静止情况为0。
c=3d8;  %电磁波速度，光速。
r0=3;   %移动台距离基站初始距离。
d=15;   %基站距离反射墙的距离。
t1=0.1:0.001:12;   %时间。
E1=cos(2*pi*f*((1-v/c).*t1-r0/c))./(r0+v.*t1);  %直射径信号。
E2=cos(2*pi*f*((1+v/c).*t1+(r0-2*d)/c))./(2*d-r0-v*t1);  %反射径信号。
subplot(2,4,6);plot(t1,E1,'g-',t1,E2,'b-',t1,E1-E2,'r-');     %画出直射径、反射径和总的接收信号。
xlabel('移动台距离基站的初始距离r0','fontname','黑体');
ylabel('直射径/反射径/总接收信号','fontname','黑体');
title('v=1,f=2,r0=3,d=15,t1=0.1:0.001:12的信号图','fontname','黑体');
axis([0,12,-0.5,0.5]);
 
%v=1,f=4,r0=3,d=15,t1=0.1:0.001:12时的情况
f=4;    %发射信号频谱。
v=1;    %移动台速度，静止情况为0。
c=3d8;  %电磁波速度，光速。
r0=3;   %移动台距离基站初始距离。
d=15;   %基站距离反射墙的距离。
t1=0.1:0.001:12;   %时间。
E1=cos(2*pi*f*((1-v/c).*t1-r0/c))./(r0+v.*t1);  %直射径信号。
E2=cos(2*pi*f*((1+v/c).*t1+(r0-2*d)/c))./(2*d-r0-v*t1);  %反射径信号。
subplot(2,4,8);plot(t1,E1,'g-',t1,E2,'b-',t1,E1-E2,'r-');     %画出直射径、反射径和总的接收信号。
xlabel('移动台距离基站的初始距离r0','fontname','黑体');
ylabel('直射径/反射径/总接收信号','fontname','黑体');
title('v=1,f=4,r0=3,d=15,t1=0.1:0.001:12的信号图','fontname','黑体');
axis([0,12,-0.5,0.5]);

2.实验结果