大学毕业设计一席谈v2之二十四 GMSK调制解调器(2)

2021年7月23日，东京奥运会开幕了。24日早上，中国代表团夺得奥运首金。厉害。能站在奥运赛场就是一种成功，如果还能得到奖牌，那更是一种荣誉和加肯定！背后流过的汗水只有自己知晓。运动健儿们加油！我们也要加油学习！本文给出另外一种生成GMSK信号方式的仿真程序，程序中还加入了解调的内容！一起来看程序吧！边看边思考！

2026年2月，冬奥会结束了，从米兰到科尔蒂纳、从瓦尔泰利纳到菲耶梅谷，中国运动员拼搏的身影闪耀在各个赛场。中国健儿最后阶段连续发力，连续夺得金牌，中国代表团以5金4银6铜共15枚奖牌数，创造冬奥会境外参赛最好成绩！真棒，向他们学习，于是接连写了好几篇技术专题付费文章。

此时，有位粉丝读完文章发来留言：如果能直观的观察到MSK和GMSK的码间串扰就好了。博主是否可以在MSK和GMSK的代码里加入眼图来观察和对比。

那我就在本文最后给出仿真程序专门用来比较GMSK方式下的码间串扰情况，至于MSK信号呢？在本人的程序上稍加改动即可分析了。GMSK解调有差分解调、准相干解调和最大似然解调等几种方式。学完下面的仿真程序，大家能看出采用了哪种解调方式吗？

% gmsk_sim1.m

% Simulation program to realize GMSK trans mission system

% 软件版本：2021a

% 2021-6-25

%*************       Preparation part      *************

sr = 256000;    % Symbol rate

br = sr;        % Bit rate

nd = 1000;      % Number of symbols that simulates in each loop

ebn0 = 15;      % Eb/N0

IPOINT = 8;     % Number of oversamples

%**************    Filter initialization    ***************

irfn=21;                  % Number of taps

B=0.25*sr;

B2=0.6*sr;

[xh] = function_gaussf(B,irfn,IPOINT,sr);    

% Trans mitter filter coefficients 

[xh2] = function_gaussf(B2,irfn,IPOINT,sr);  

% Receiver filter coefficients 

%**************  START CALCULATION   ***************

nloop = 100;  % Number of simulation loops

noe = 0;    % Number of error data

nod = 0;    % Number of trans mitted data

for iii=1:nloop     

%%%%%%     Trans mitter     %%%%%%  

%********   Data generation   ********* 

        data1 = randi([0,1],1,nd);  

%*******  GMSK Modulation   ********  

    % 双极性变换

    data11 = 2*data1-1;

    data2 = upsample(data11,IPOINT);    

    % 高斯滤波

    data3 = conv(data2,xh);                

    th = zeros(1,length(data3)+1);

    ich2 = zeros(1,length(data3)+1);

    qch2 = zeros(1,length(data3)+1);

  % 相位累加

    for ii = 2:length(data3)+1