m06_demodce.m

MATLAB仿真技术与实例应用教程

M = 4;          %设置M的数目
Fd = 1;         %信号采样速率
Fs = 32;        %采样频率
SNRperBit = 5;  %信噪比
adjSNR = SNRperBit-10*log10(Fs/Fd)+10*log10(log2(M));
x = randint(5000,1,M); % 源信号

%正交FSK调制
tone = 0.5;
randn('state',1945724); 
w1 = dmodce(x,Fd,Fs,'fsk',M,tone); 
y1  = awgn(w1, adjSNR, 'measured', [], 'dB');
z1 = ddemodce(y1,Fd,Fs,'fsk',M,tone); 
ser1 = symerr(x,z1)  %输出符号误差率

% 非正交FSK调制
tone = 0.25;   
randn('state',1945724); .
w2 = dmodce(x,Fd,Fs,'fsk',M,tone); 
y2  = awgn(w2, adjSNR, 'measured', [], 'dB');
z2 = ddemodce(y2,Fd,Fs,'fsk',M,tone); 
ser2 = symerr(x,z2)   %输出符号误差率