mod_msk_030610.asv

msk,qpsk,oqpsk,qam仿真程序

function y = mod_msk_030610 (x,n)
%  总共有42个码元，每个码元有120个点，共5120个点。
  Fs Fd=2 Fs_n  Fc=10;
num=1024*Fs_n/Fs;
N=fix(num/Fs_n/Fd);
for i=2:num
    if x(i)==1
        x(i)=mod((x(i-1)+1),2);
    elseif x(i)==0
        x(i)=x(i-1);
    end
end
x=2*x-1;   %0:-1 1:+1
%%%%%%%%%%%%%%%%%%%%%
%进行串并变换

for i=1:2:N
    for j=1:Fs_n/Fd
        I(120*(i+1)/2+j)=x(2*120*(i-1)+j);
        Q(120*(i+1)/2+j)=x(2*120*i+j);
    end
end
for i=1:N*120/2             %将2520个点变换成5040点的信号
    I1(2*i-1)=I(i);
    I1(2*i)=I(i);
    Q1(2*i-1)=Q(i);
    Q1(2*i)=Q(i);
end
for i=(N*120+1):num        %将同相与正交信号都变换成5120点的信号
    I1(i)=x(i);
    Q1(i)=x(i);
end
for i=1:(num-Fs_n/Fd)         %将正交分量相对于同相分量延迟一个码元周期
    Q11(i)=Q1(i+120);
end
for i=1:Fs_n/Fd
    Q11(num-Fs_n/Fd+i)=Q1(i);
end
y_I=I1'.*cos(pi/2*Fd*n).*cos(2*pi*Fc*n);
y_Q=Q11'.*sin(pi/2*Fd*n).*sin(2*pi*Fc*n);
y=y_I-y_Q;          %最后的MSK调制信号