多径衰落仿真器matlab代码

该代码对一个时间频率双衰落信道进行了仿真，刻画出了其二阶统计特性。

clc;clear;

v=100;% m/s
fc=900;% Mhz
fs=1e5;% Hz
Ns=10000;
N=10;%多经数

t=0:1/fs:Ns/fs-1/fs;
st=sin（2*pi*500*t）;

delay=[0 8 14 23 38 42 55 69 73 86];%/ us
gain=[0 -21 -35 -40 -44 -46 -52 -55 -58 -60];
input=[zeros（1delay（N）） st];
h=zeros（N+1Ns）;

for i=1:N
    temp=[zeros（1delay（N）） rayleigh_channel（fcvfsNs1）];
    
    h（i:）=temp（delay（N）+1-delay（i）:delay（N）+Ns-delay（i））*10^（gain（i）/20）;
    h（N+1:）=h（N+1:）+h（i:）;
end;
rt=abs（h（N+1:））.*st;

figure（1）;

subplot（221）;
plot（tst）;
title（‘原始符号的波形‘）;
xlabel（‘time/s‘）;
ylabel（‘幅度‘）;

subplot（223）;
plot（trt）;
title（‘经过信道后的符号波形‘）;
xlabel（‘time/s‘）;
ylabel（‘幅度‘）;

subplot（222）;
plot（t20*log10（abs（h（1:））））;
title（‘多径信道幅度响应‘）;
xlabel（‘time/s‘）;
ylabel（‘|h（t）|/dB‘）;

subplot（224）;
plot（tangle（h（N+1:）））;
title（‘多径信道的相位响应‘）;
xlabel（‘time/s‘）;
ylabel（‘angle/rad‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        428  2016-10-19 21:27  rayleigh_channel.m

     文件        916  2016-10-19 21:23  main.m

----------- ---------  ---------- -----  ----

                 1344                    2