通原实验软件报告(8)

接收信号的功率谱

解调信号的功率谱幅频特性

实验八 第一类部分响应系统 源程序

%第一类部分响应系统 global dt df t f N close all N=2^10; L=16;

%采样点数

%每码元的采样点数

M=N/L; Rb=2;

%码元数

%码速率是2Mb/s

Ts=1/Rb; dt=Ts/L;

%码元间隔 %时域采样间隔

df=1/(N*dt); %频域采样间隔 T=N*dt; Bs=N*df/2; Tb=T/256

f=[-Bs+df/2:df:Bs]; t=[-T/2+dt/2:dt:T/2]; for jj=1:100 a=sign(randn(1,M)); s=zeros(1,N); for ii=1:L

s(ii+[0:M-1]*L)=a; end

%产生双极性NRZ码 %码元矢量

%截短时间 %系统带宽

S=T2F(s);

P=S.*conj(S)/T; %信号功率谱

h1=Tb*Tb*sin((pi*t)/Tb)/((pi*t)/(Tb-t)); H1=t2f(h1);

PH=H1.*conj(H1)/T; %功率谱 Y=S.*H1

end figure(1)

set(1,'Position',[30,100,450,300]) as=abs(P); plot(f,as,'g'); grid

axis([-2,+2,min(as),max(as)]) xlabel('f (KHz)') ylabel('P(f) (W/KHz)') title('双极性NRZ码的功率谱') figure(2)

set(2,'Position',[530,100,450,300]) ass=abs(PH); plot(t,s,'g'); grid

axis([-10,+10,-2,+2]) xlabel('t(us)');ylabel('s(t)(V)') title('双极性NRZ码时域波形') figure(3)

set(3,'Position',[30,200,450,300]) ar=abs(Y); plot(f,ar,'g'); grid

axis([-Bs,+Bs,min(ar),max(ar)]) xlabel('f (KHz)') ylabel('P(f) (W/KHz)')

title('未经低通后接收信号的功率谱') fir=lpf(1/(2*Tb));%低通滤波器 FIR=t2f(fir);

Yz=Y.*FIR;%最终解调后的信号 PYZ=Yz.*conj(Yz)/T; %功率谱 figure(4)

set(4,'Position',[530,200,450,300]) aZ=abs(PYZ); plot(f,aZ,'g'); grid

axis([-2,+2,min(aZ),max(aZ)]) xlabel('f (KHz)') ylabel('P(f) (W/KHz)')

title('经低通后的功率谱幅频特性') 仿真图

双极性NRZ码时域波形

双极性NRZ码的功率谱

未经低通后接收信号的功率谱

经低通后的功率谱幅频特性