传真(Fax)编程

开发平台：
Matlab

main1.m：源码内容
							%本题目中载频为900M
%多普勒频移计算公式为fm = v*fc/c
close all
clear all                                %清零
fm = [16,45,80,120];                %四组不同的多普勒频移，对应车速为19.2km/h,54km/h,96km/h,144km/h
fm_0=50;                                 %固定多普勒频移  对应车速为60km/h
M = 50000;                               %仿真点数
fs = 5000;                              %时间抽样频率
ts = 1/fs;                               %最小时间间隔
N0 = [6,12,120,480];                  %四组不同的振荡器数目
N0_0 = 64;                              %固定振荡器数
tt=[ts:ts*5:M*ts*5];                          %时间长度
S = ['b','k','r','g'];                   %绘图颜色及方式
   
for i = 1:4;                             %绘制不同振荡器数目条件下的自相关函数曲线，多普勒频移固定
    RayLeigh = Jakes(fm_0,fs,M,N0(i));        %调用仿真器
    [RayLeigh_dt,dt] = xcorr(RayLeigh,'unbiased');%生成自相关函数
    hold on
    figure(1);
    plot(dt*ts*5,RayLeigh_dt,S(i),'LineWidth',1.5);
    grid on
    xlabel('tao(s)');
    ylabel('R(tao)');
    title('不同振荡器数目条件下的自相关函数曲线');
    axis([0,0.3,-1.5,2.2]);
end
for j = 1:4;                              %绘制不同车速条件下的自相关函数曲线，振荡器数目固定
    RayLeigh = Jakes(fm(j),fs,M,N0_0);         %调用仿真器
    [RayLeigh_dt,dt] = xcorr(RayLeigh,'unbiased');%生成自相关函数,计算的是无偏估计
    figure(2);
    hold on
    plot(dt*ts*5,RayLeigh_dt,S(j),'LineWidth',1.5);
    grid on
    xlabel('tao(s)');
    ylabel('R(tao)');
    title('不同车速条件下的自相关函数曲线');
    axis([0,0.3,-1,2.2]);
                                           %包络时间分布曲线
    RayLeighAmp=abs(RayLeigh);                       %求包络幅度
    figure(3);    
    hold on
    plot(tt,20*log10(RayLeighAmp),S(j),'LineWidth',1.5); %幅值取对数
    grid on
    xlabel('时间(s)');
    ylabel('幅度(db)');
    title('不同车速条件下的包络时间分布曲线');
    axis([0,0.3,-28,15]);
    
                                                %相位时间分布曲线
    RayLeighPhase=angle(RayLeigh);                        %取相位向量
    figure(4);
    hold on
    plot(tt,RayLeighPhase,S(j),'LineWidth',1.5);
    grid on
    xlabel('时间（s）');
    ylabel('相位(rad)');
    title('不同车速条件下的相位时间分布曲线');
    axis([0,0.3,-4,4]);
end

						