%function A=NRZ_L(a,m)
%非歸零電平碼顯示
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
a=[1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0];
m=20;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(1);
%a為原碼，m為圖顯密集度
%A為顯示序列
n=length(a);%獲取原碼碼元數
pA=ones(1,m);%由圖顯密集度m產生行向量pA
for i=1:n%讀取原碼a的數據，產生相應plot顯示行向量A
    A(m*(i-1)+1:i*m)=a(i)*pA;
end
mn=m*n;
t=1:mn;%產生和A相同元素個數的行向量t
H=plot(t,A);%繪制圖形，并將plot句柄賦值給H
set(H,'linewidth',2);%通過H設置圖形線段的寬度
axis([0 mn -1.2 2.2]);%設置x,y軸的范圍
title('原始序列顯示','FontSize',7.5)
%顯示標題，并設置字體的大小
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(2);
subplot(2,1,1)
%function p_fsk=PSK(a,m)
%PSK相位調制
%a為原碼,m為顯示密集度
%A=NRZ_L(a,m);%調用NRZ_L函數,得到行向量A
n=length(a);
mn=m*n;
B=ones(1,mn);%產生常數"1"行向量B
t=1:mn;
C=A+B;%A+B產生新的向量C,供PSK調制計算
p_psk=sin(1/m*2*pi*t+C*pi);%實現PSK調制
%"1"時,相位改變2*pi,即不改變
%"0"時,相位改變1*pi
H=plot(t,p_psk-1,'b',t,A+1,'b');%在同一副圖中顯示
set(H,'Linewidth',2);%設置圖形顯示屬性
axis([0 mn -3 3]);
title('經PSK調制顯示','FontSize',7.5);
subplot(2,1,2)
p_psk1=fft(p_psk);
H=plot(t,p_psk1);%在同一副圖中顯示
set(H,'Linewidth',2);%設置圖形顯示屬性
axis([0 mn -3 3]);
title('經PSK調制頻率顯示','FontSize',7.5)