gen_cpfsk_gmsk.m

来自「是一个对mimo连续相位调制系统的均衡程序」· M 代码 · 共 32 行

function q = gen_cpfsk_gmsk(fs, L, BT)
% This script generates g and q pulse 
% Input:  
%       L: memory length
%       fs: Number of samples per symbol
%       BT: -3dB bandwidth of thd Gaussian pulse 
% Output: 
%       q: q pulse
%fs=16;
%L=1;
%BT = 1;
T=1;
Ts = T/fs;
N = L*fs;

beta = 2*pi/sqrt(log(2));

g = zeros(1,N);
g_tmp = zeros(1,N/2);
q = zeros(1,N);

for i = 1:N/2
    g(i+N/2) = (1/2/T)*(qfunc(beta*BT/T*((i-1)*Ts-T/2)) - qfunc(beta*BT/T*((i-1)*Ts+T/2)));
    g_tmp(i) = (1/2/T)*(qfunc(beta*BT/T*(i*Ts-T/2)) - qfunc(beta*BT/T*(i*Ts+T/2)));
end
g(1:N/2) = g_tmp(N/2:-1:1);

for n = 1:N
    for i = 1:n
       q(n) = q(n) + g(i)*Ts;
    end
end