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); %瑞利衰减幅值。