传真(Fax)编程

开发平台：
Matlab

main2.m：源码内容
							%本题目中载频为900M
%多普勒频移计算公式为fm = v*fc/c
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 j = 1:4;                              %绘制不同车速条件下的自相关函数曲线，振荡器数目固定
     RayLeigh = Jakes(fm(j),fs,M,N0_0);         %调用仿真器
    
     %包络时间分布曲线
    rayAmp=abs(RayLeigh);
                                                 %包络概率分布曲线
    rayAmpMin=min(rayAmp);                          %最小幅值
    rayAmpMax=max(rayAmp);                          %最大幅值
    rayAmpSpan=linspace(rayAmpMin,rayAmpMax,100);   %均匀取100个值 
    rayAmpCount_tmp=hist(rayAmp,rayAmpSpan);        %生成直方图
    rayAmpCount=rayAmpCount_tmp/length(rayAmp);
    figure(5);
    plot(rayAmpSpan,rayAmpCount,S(j));
    stem(rayAmpSpan,rayAmpCount,S(j));              %离散序列图
    grid on
    xlabel('幅度');
    ylabel('概率');
    title('不同车速条件下的包络概率分布曲线');
    hold on
                                                    %相位概率分布曲线
    rayPhase_tmp=angle(RayLeigh);                        %取相位向量
    rayPhase=rayPhase_tmp(1:M);                     %取M个样值
    rayPhaseMin=min(rayPhase);                      %最小相位
    rayPhaseMax=max(rayPhase);                      %最大相位
    rayPhaseSpan=linspace(rayPhaseMin,rayPhaseMax,100);%均匀取100个值
    rayPhaseCount_tmp=hist(rayPhase,rayPhaseSpan);     
    rayPhaseCount=rayPhaseCount_tmp/length(rayPhase);
    figure(6);
    stem(rayPhaseSpan,rayPhaseCount,S(j));
    grid on
    xlabel('相位(rad)');
    ylabel('概率');
    title('不同车速条件下的相位概率分布曲线');
    hold on
end

						