bfskmain.m

来自「一个关于仿真的程序」· M 代码 · 共 42 行

clc
clear
echo on
%x表示信噪比
x=0:15;
% y表示信号的误比特率，它的长度与x相同
y=x;
% BFSK调制的频率间隔等于24kHz
FrequencySeparation=24000;
% 信源产生信号的bit率等于10kbit/s
BitRate=10000;
% 仿真时间设置为10秒
SimulatonTime=10;
% BFSK调制信号每个符号的抽样数等于2
SamplesPerSymbol=2;
%发送端和接收端的相对运动速度
Velocity=40;
% 光速
LightSpeed=3*10^8;
% 载波频率
Frequency=825*10^6;
%计算载波波长
WaveLength=LightSpeed/Frequency;
% 根据运动速度和波长计算多普勒频移
% 注意要把运动速度的单位转换成米/秒
Fd=Velocity*10^3/3600/WaveLength;
for i=1:length(x)
%     信道的信噪比依次取x中的元素；
    SNR=x(i);
%     运行仿真程序，得到的误比特率保存在工作区变量BitErrorRate中
    sim('BFSK_rayleigh');
%     计算BitErrorRate的均值作为本次仿真的误比特率
    y(i)=mean(BitErrorRate);
end

% 准备一个空白的空间
% hold off;
figure

% 绘制x和y的关系曲线图，纵坐标采用对数表示
semilogy(x,y,'-*');
grid on;