main3.m

理想低通滤波器特性曲线绘制 频率图分析 升余弦滤波器特性曲线绘制 频谱图分析 部分响应系统滤波器特性曲线

%部分响应系统
global dt t df N
N=2^13; %采样点数
L=32;   %每码元的采样点数
M=N/L   %码元数
Rb=2;		%码速率是2Mb/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];   %频域横坐标
x=input('请输入部分响应系统的类别（1-5）:');
switch(x);
    case 1
        R=[0 0 0 1 1];
    case 2
        R=[0 0 1 2 1];
    case 3
        R=[0 0 -1 1 2];
    case 4
        R=[0 0 -1 0 1];
    case 5
        R=[-1 0 2 0 -1];
end
g=R(1)*sinc(pi*t/Ts)+R(2)*sinc(pi*(t-Ts)/Ts)+R(3)*sinc(pi*(t-2*Ts)/Ts)+R(4)*sinc(pi*(t-3*Ts)/Ts)+R(5)*sinc(pi*(t-4*Ts)/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([-1,1,0,30])
 xlabel('f (MHz)')
 ylabel('Ps(f)  (dBm/MHz)')
 title('第x类部分响应系统频谱图')
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('第x类部分响应系统的冲击响应')