v1.0 可编辑可修改
2121
v1.0 可编辑可修改
五、实验总结
通过本次实验,可以知道由于高斯滤波器的作用,GMSK信号的功率谱密度非常平滑,而MSK信号的功率谱密度则有较大的旁瓣。结果显示随BT值的减小,GMSK信号的功率谱密度衰减速度明显加快,而当BT增大时,其谱密度逐渐向MSK信号的谱密度靠近。当BT取值为无穷大时,GMSK就退化为MSK。通常在码元宽度一定的情况下,滤波器的等效带宽B越大,基带信号经滤波后损失的能量就越小,已调波就会越接近MSK信号,也就是说,BT值越小,GMSK信号的频谱利用率就越高。
附录:源程序代码(重要部分) %绘制调制波形00101010
Ak=[0 0 1 0 1 0 1 0]; %产生8个基带信号 Ak=2*Ak-1;
gt=ones(1,B_sample); %每码元对应的载波信号
2222
v1.0 可编辑可修改 Akk=sigexpand(Ak,B_sample); %码元扩展 temp=conv(Akk,gt); %码元扩展 Akk=temp(1:length(Akk)); %码元扩展 tt=*Tb:Dt:*Tb-Dt;
%g(t)=Q[2*pi*Bb*(t-Tb/2)/sqrt(log(2))]-Q[2*pi*Bb*(t+Tb/2)/sqrt(log(2))]; %Q(t)=erfc(t/sqrt(2))/2;
gausst=erfc(2*pi*Bb*(tt-Tb/2)/sqrt(log(2))/sqrt(2))/2-erfc(2*pi*Bb*(tt+Tb/2)/sqrt(log(2))/sqrt(2))/2;
J_g=zeros(1,length(gausst)); %使J_g 的长度和Gausst的一样 for i=1:length(gausst) if i==1
J_g(i)=gausst(i)*Dt; else
J_g(i)=J_g(i-1)+gausst(i)*Dt; end; end;
J_g=J_g/2/Tb; %计算相位Alpha
2323
v1.0 可编辑可修改 Alpha=zeros(1,length(Akk)); k=1; L=0;
for j=1:B_sample
J_Alpha=Ak(k+2)*J_g(j);
Alpha((k-1)*B_sample+j)=pi*J_Alpha+L*pi/2; end; k=2; L=0;
for j=1:B_sample
J_Alpha=Ak(k+2)*J_g(j)+Ak(k+1)*J_g(j+B_sample); Alpha((k-1)*B_sample+j)=pi*J_Alpha+L*pi/2; end; k=3; L=0;
for j=1:B_sample
J_Alpha=Ak(k+2)*J_g(j)+Ak(k+1)*J_g(j+B_sample)+Ak(k)*J_g(j+2*B_sample); Alpha((k-1)*B_sample+j)=pi*J_Alpha+L*pi/2; end;
2424
v1.0 可编辑可修改 k=4; L=0;
for j=1:B_sample
J_Alpha=Ak(k+2)*J_g(j)+Ak(k+1)*J_g(j+B_sample)+Ak(k)*J_g(j+2*B_sample)+Ak(k-1)*J_g(j+3*B_sample);
Alpha((k-1)*B_sample+j)=pi*J_Alpha+L*pi/2; end; L=0;
for k=5:B_num-2 if k==5 L=0; else
L=L+Ak(k-3); end;
for j=1:B_sample
J_Alpha=Ak(k+2)*J_g(j)+Ak(k+1)*J_g(j+B_sample)+Ak(k)*J_g(j+2*B_sample)+Ak(k-1)*J_g(j+3*B_sample)+Ak(k-2)*J_g(j+4*B_sample);
Alpha((k-1)*B_sample+j)=pi*J_Alpha+mod(L,4)*pi/2;
2525