%實驗三：畫出單極性歸零碼及其功率譜
global dt t df N
close all
k=input('取樣點數=2^k,k=[14]');
if k==[],k=14;end
N=2^k;			%采樣點數
L=64;			%每碼元的采樣點數
M=N/L;			%碼元數
Rb=2;			%碼速率為2Mb/s
Ts=1/Rb;		%碼元間隔
dt=Ts/L;		%時域采樣間隔
Rt=0.5;		%占空比
df=1/(N*dt);	%頻域采樣間隔
T=N*dt;		%截短時間
Bs=N*df/2;		%系統帶寬
t=linspace(-T/2,T/2,N);	%時域橫坐標
f=linspace(-Bs,Bs,N);	%頻域橫坐標
EP=zeros(1,N);
for jj=1:100
a=round(rand(1,M));		%產生M個取值0，1等概的隨機碼
s=zeros(1,N);				%產生一個N個元素的零序列
for ii=1:Rt*Ts/dt
   s(ii+[0:M-1]*L)=a;	%產生單極性歸零碼
end
Q=t2f(s);					%付氏變換
P=Q.*conj(Q)/T; 			%P為單極性歸零碼的功率  
EP=(EP*(ii-1)+P)/ii;    %累計平均
aa=30+10*log10(EP+eps);	%加eps以避免除以零
end 

figure(1)
set(1,'position',[10,50,300,200]) %設定窗口位置及大小
figure(2)
set(2,'position',[350,50,300,200])%設定窗口位置及大小
figure(1)
plot(f,aa,'g')
xlabel('f(MHZ)')
ylabel('Ps(f)(MHZ)')
axis([-15,+15,-50,50])
grid
figure(2)
plot(t,s,'r')
xlabel('t(ms)')
ylabel('s(t)(V)')
axis([-10,10,-0.5,1.5])
grid