%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);
%function p_fsk=FSK(a,m)
%FSK頻率調制
%a為原碼,m為顯示密集度
%A=NRZ_L(a,m);%調用NRZ_L函數,得到一維數組A
n=length(a);
mn=m*n;
t=1:mn;
B=1+A;%由A產生新的向量B,供FSK調制計算
p_fsk=sin(1/m*2*pi*(B.*t));% 實現FSK調制
%"1"對應頻率為1/m*2*pi*2
%"0"對應頻率為1/m*2*pi
H=plot(t,p_fsk-1,'b',t,A+1,'b');%在同一副圖中顯示
set(H,'Linewidth',2);%設置圖形顯示屬性
axis([0 mn -3 3]);
title('經FSK調制顯示','FontSize',7.5)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(3);
%function p_ask=ASK(a,m)
%ASK幅度調制
%a為原碼，m 為顯示密集度
%A=NRZ_L(a,m);%調用NRZ_L函數，得到行向量A
n=length(a);
mn=m*n;
t=1:mn;
p_ask=A.*sin((1/m)*2*pi*t);%設置頻率為1/m*2*pi使正弦信號的一個周期等于一個碼元的長度
H2=plot(t,p_ask-1,'b',t,A+1,'b');
%在同一圖中顯示調制結果
set(H2,'Linewidth',2);%設置圖形顯示屬性
axis([0 mn -3 3]);
title('經ASK調制顯示','FontSize',7.5)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(4);
%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);