xiangwei.m

来自「GMSK 调制方式 通信专业学生通用 功能很强大」· M 代码 · 共 34 行

%%相位路径程序%%
a=sign(randn(1,2^23-1));
for j=1:50
    suml=0;
    hh=1;
    a((2+10*(j-1)):(2+10*(j-1)+4))=-1;
for i=(2+10*(j-1)):(2+10*(j-1)+15)
    
    suml=suml+a(i-1);
    sum0=0;
        for jj=-2:1:2
            temp=faifai((jj-2)*T,0,jj,0.0000036923,0.3);
            sum0=sum0+a(i+2+jj)*temp;
        end
        fai0=sum0*pi+suml*pi/2;
        for ii=1:8
            t=T/8*ii;
            sum0=0;
            for jj=-2:1:2
                sum0=sum0+a(i+2+jj)*faifai(0,t,jj,0.0000036923,0.3);
            end
            deta_fai=pi*sum0;
            dian(hh)=fai0+deta_fai;
            hh=hh+1;
        end
end
plot(dian,'-','MarkerSize',4.6);
hold on;
end