%PSK的調制程序
function psk(A,Lt,fc,t,T,d_NRZ);
ht=A*cos(2*pi*fc*t);
figure(2)
d_2psk=2*d_NRZ-1;
s_2psk=d_2psk(1:Lt).*ht;
subplot(221)
plot(t,s_2psk);
axis([0 10 -1.2 1.2]);ylabel('2psk');

subplot(222)
[f,s_2pskf]=T2F(t,s_2psk);
plot(f,10*log10(abs(s_2pskf).^2/T));
axis([-fc-4 fc+4 -50 10]);ylabel('2PSK功率譜密度(dB/Hz)');