哈工大移动通信设计题瑞利衰落信道代码

main

clc;

LengthOfSignal=10240; %设置信号长度(由于最好大于两倍fc奈奎斯特采样)

fm=512; %最大多普勒频移

fc=5120; %载波频率

t=1:LengthOfSignal; % SignalInput=sin(t/100);

SignalInput=sin(t/100)+cos(t/65); %时域信号输入

delay=[0 21 62 100 150 250]; %设置不同路径的时延

power=[0 -1 -5 -11 -16 -20]; %功率衰减系数dB

y_in=[zeros(1,delay(6)) SignalInput]; %为时移补零

y_out=zeros(1,LengthOfSignal); %时域输出信号

fori=1:6

Ray;

y_out=y_out+r.*y_in(delay(6)+1-delay(i):delay(6)+LengthOfSignal-delay(i))*10^(power(i)/20); end; %进行输出信号叠加

figure(1);

subplot(2,1,1);

plot(SignalInput(delay(6)+1:LengthOfSignal),'r'); %画出时域信号输入波形

title('时域信号输入');

subplot(2,1,2);

plot(y_out(delay(6)+1:LengthOfSignal),'r'); %画出时域信号输出波形

title('时域信号输出');

figure(2);

plot(Sf1,'r');

title('多普勒滤波器的频率响应特性'); %画出多普勒滤波器的频率响应特性

Ray

f=1:2*fm-1; %设置通频带宽度

y=0.5./((1-((f-fm)/fm).^2).^(1/2))/pi; %多普勒功率谱函数(基带)

Sf=zeros(1,LengthOfSignal);

Sf1=y;%多普勒滤波器的频率响应特性

Sf(fc-fm+1:fc+fm-1)=y; %(把基带映射到载波频率)

x1=randn(1,LengthOfSignal);

x2=randn(1,LengthOfSignal);

nc=ifft(fft(x1+1i*x2).*sqrt(Sf)); %同相分量nc函数表达式

x3=randn(1,LengthOfSignal);

x4=randn(1,LengthOfSignal);

ns=ifft(fft(x3+1i*x4).*sqrt(Sf)); %正交分量ns函数表达式

r0=(real(nc)+1i*real(ns)); %瑞利衰减因子

r=abs(r0); %瑞利衰减幅值