% 理想低通濾波器系統
global dt t df N%定義全局變量
N=2^13; %采樣點數
L=32;   %每碼元的采樣點數
M=N/L   %碼元數
Rb=2;		%碼速率是2 b/s
Ts=1/Rb;        %碼元間隔
dt=Ts/L;        %時域采樣間隔
df=1/(N*dt)     %頻域采樣間隔
T=N*dt          %截短時間
Bs=N*df/2       %系統帶寬
Na=4;           %示波器掃描寬度為4個碼元
t=[-T/2+dt/2:dt:T/2]; %時域橫坐標
f=[-Bs+df/2:df:Bs];   %頻域橫坐標
g=sin(pi*t/Ts)./(pi*t/Ts);%理想低通濾波器系統響應表達式
G=t2f(g);%傅里葉變換（時域->頻域）
EP=zeros(size(f))+eps;
for ii=1:100
 a=sign(randn(1,M));
 imp=zeros(1,N);    %產生沖激序列
 imp(L/2:L:N)=a/dt;
 S=t2f(imp).*G;       
 s=f2t(t2f(imp).*G);  
 s=real(s);
 P=S.*conj(S)/T;                                                 
 EP=(EP*(ii-1)+P+eps)/ii;
 subplot(212)
 plot(f,30+10*log10(EP));%繪制理想低通濾波器特性
 grid on
 axis([-2,+2,0,30]) %設定坐標軸范圍
 xlabel('f (MHz)')
 ylabel('Ps(f)  (dBm/MHz)')
 title('理想低通頻域特性')
end
subplot(211)
plot(t,g);%繪制理想低通濾波器特性
grid on
hold on
line([-4*pi 4*pi],[0 0])
axis([-6 6 -0.4 1.3])
title('理想低通特性的沖擊響應')