function mod_out = pskmod1(mod_in)
%function用于對輸入的基帶數(shù)字信號進行調(diào)制；
%A1，A0，兩個一唯數(shù)組用于存處對應數(shù)字信號1，0的不同相位的模擬調(diào)制段信號；
%本程序為方便直接從圖形中看出結(jié)果，采用一個周期正弦波來表示1BIT的數(shù)字信號，如國需要改變調(diào)制頻率，只需相對改變A1，A0數(shù)組；
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);