pskmod1.m

这是PSK在MATLAB上调制的源程序！希望大家可以用的上！

function mod_out = pskmod1(mod_in)
%function用于对输入的基带数字信号进行调制；
%A1，A0，两个一唯数组用于存处对应数字信号1，0的不同相位的模拟调制段信号；
%本程序为方便直接从图形中看出结果，采用一个周期正弦波来表示1BIT的数字信号，如国需要改变调制频率，只需相对改变A1，A0数组；
a1=[      0    0.0628    0.1253    0.1874    0.2487    0.3090    0.3681    0.4258    0.4818    0.5358     0.5878    0.6374    0.6845    0.7290    0.7705    0.8090    0.8443    0.8763    0.9048    0.9298    0.9511    0.9686    0.9823    0.9921    0.9980    1.0000    0.9980    0.9921    0.9823    0.9686    0.9511    0.9298    0.9048    0.8763    0.8443    0.8090    0.7705    0.7290    0.6845    0.6374    0.5878    0.5358    0.4818    0.4258    0.3681    0.3090    0.2487    0.1874    0.1253    0.0628    0.0000   -0.0628   -0.1253   -0.1874   -0.2487    -0.3090   -0.3681   -0.4258   -0.4818   -0.5358   -0.5878   -0.6374   -0.6845   -0.7290   -0.7705   -0.8090   -0.8443   -0.8763   -0.9048   -0.9298   -0.9511   -0.9686   -0.9823   -0.9921   -0.9980   -1.0000   -0.9980   -0.9921   -0.9823   -0.9686    -0.9511   -0.9298   -0.9048   -0.8763   -0.8443   -0.8090   -0.7705   -0.7290   -0.6845   -0.6374   -0.5878   -0.5358   -0.4818   -0.4258   -0.3681   -0.3090   -0.2487   -0.1874   -0.1253   -0.0628  ];
a0=[ 0.0000    -0.0628   -0.1253   -0.1874   -0.2487    -0.3090   -0.3681   -0.4258   -0.4818   -0.5358   -0.5878   -0.6374   -0.6845   -0.7290   -0.7705   -0.8090   -0.8443   -0.8763   -0.9048   -0.9298   -0.9511   -0.9686   -0.9823   -0.9921   -0.9980   -1.0000   -0.9980   -0.9921   -0.9823   -0.9686   -0.9511   -0.9298   -0.9048   -0.8763   -0.8443   -0.8090   -0.7705   -0.7290   -0.6845   -0.6374   -0.5878   -0.5358   -0.4818   -0.4258   -0.3681   -0.3090   -0.2487   -0.1874   -0.1253   -0.0628      0.0000    0.0628    0.1253    0.1874    0.2487    0.3090    0.3681    0.4258    0.4818    0.5358     0.5878    0.6374    0.6845    0.7290    0.7705    0.8090    0.8443    0.8763    0.9048    0.9298     0.9511    0.9686    0.9823    0.9921    0.9980    1.0000    0.9980    0.9921    0.9823    0.9686    0.9511    0.9298    0.9048   0.8763    0.8443    0.8090    0.7705    0.7290    0.6845    0.6374     0.5878    0.5358    0.4818    0.4258    0.3681    0.3090    0.2487    0.1874    0.1253    0.0628  ];
 for i=1:length(mod_in)
       if  mod_in(i)==0
           for j=1:100
               mod_out(i*100+j-100)=a0(j);
           end    
       else
           for j=1:100
               mod_out(i*100+j-100)=a1(j);
           end
       end
 end
 t=0:0.01:length(mod_in)-0.01;
 plot(t,mod_out);