function y = mod_msk_030610 (x,n,M)
%  總共有42個(gè)碼元，每個(gè)碼元有120個(gè)點(diǎn)，共5120個(gè)點(diǎn)。
global n Fs Fd Fs_n k Fc;
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
%%%%%%%%%%%%%%%%%%%%%
%進(jìn)行串并變換
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個(gè)點(diǎn)變換成5040點(diǎn)的信號
    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點(diǎn)的信號
    I1(i)=x(i);
    Q1(i)=x(i);
end
for i=1:(num-Fs_n/Fd)         %將正交分量相對于同相分量延遲一個(gè)碼元周期
    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調(diào)制信號