L=8; M=N/L Rb=2;
%每码元的采样点数 %码元数 %码速率是2Mb/s %码元间隔 %时域采样间隔
Ts=1/Rb; dt=Ts/L;
df=1/(N*dt) T=N*dt
%频域采样间隔 %截短时间 %系统带宽 %滚降系数=0.5
Bs=N*df/2 alpha=0.5
t=linspace(-T/2,T/2,N); %时域横坐标
f=linspace(-Bs,Bs,N)+eps; %频域横坐标
figure(1)
set(1,'Position',[30,100,450,300])%设定窗口位置及大小
hr1=sin(pi*t/Ts)./(pi*t/Ts);
hr2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2); hr=hr1.*hr2; HR=abs(T2F(hr)); GT=sqrt(HR); GR=GT;
for loop1=1:20
Eb_N0(loop1)=(loop1-1) %分贝值变为真值
eb_n0(loop1)=10^(Eb_N0(loop1)/10); Eb=1;
n0=Eb/eb_n0(loop1);
%信道噪声谱密度 sita=n0*Bs; %噪声功率 n_err=0; %误码计数 for loop2=1:5
a=sign(randn(1,M));%发送码元 imp=zeros(1,N); imp(L/2:L:N)=a/dt; IMP=T2F(imp);
n_ch=sqrt(sita)*randn(size(t)); %信道噪声 nr=real(F2T(T2F(n_ch).*GR)); sr=real(f2t(IMP.*HR))+nr; y=sr(L/4:L:N);
%以L/4为起点，步长L,取样点N
aa=sign(y); %接收码元 n_err=n_err+length(find(aa~=a)) ; end
Pe(loop1)=n_err/(M*loop2);%误码率Pe'=n_err/(M*loop2) figure(1)
semilogy(Eb_N0,Pe,'g'); eb_n0=10.^(Eb_N0/10); hold on
semilogy(Eb_N0,0.5*erfc(sqrt(eb_n0))); axis([0,9,1e-4,1]) xlabel('Eb/N0') ylabel('Pe') title('误码率曲线') end 仿真图 误码率曲线
%还原为真值
实验结论
由图可以看出，由于判决时刻有Ts/4的误差，滚降系数为0.5,所以实际的信噪比（绿线）比理想信噪比（蓝线）差。
实验五 Pe~Eb\\No，升余弦滚降系数a=0.5，取样时间无偏差，但信道是多径信道，C(f)=abs(1-0.5*exp(-j*2*pi*f*dt))，dt=Ts/2。 源程序
global dt t f df N T close all clear Eb_N0 Pe
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 %系统带宽 alpha=0.5 %滚降系数
t=[-T/2+dt/2:dt:T/2]; %时域横坐标 f=[-Bs+df/2:df:Bs]; %频域横坐标
figure(1)
set(1,'Position',[30,100,450,300])%设定窗口位置及大小 %升余弦
hr1=sin(pi*t/Ts)./(pi*t/Ts);
hr2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
hr=hr1.*hr2;
HR=abs(t2f(hr)); %取模是为了忽略时延 GT=sqrt(abs(HR));
GR=GT; %发送和接收滤波器模型 tao=Ts/2;
C=1-0.5*exp(-j*(2*pi*f*tao)) ; %多径信道模型
G=C.*HR; %总体特性
for loop1=1:20
Eb_N0(loop1)=(loop1-1) ;?/N0 in dB eb_n0(loop1)=10^(Eb_N0(loop1)/10); Eb=1;
n0=Eb/eb_n0(loop1); %信道的噪声谱密度
sita=n0*Bs; %信道中噪声功率 n_err=0; %误码计数 for loop2=1:5
a=sign(randn(1,M));
imp=zeros(1,N); %产生冲激序列 imp(L/2:L:N)=a/dt; IMP=t2f(imp);
n_ch=sqrt(sita)*randn(size(t)); %信道噪声
nr=real(f2t(t2f(n_ch).*GR.*C)); %输出噪声
sr=real(f2t(IMP.*G))+nr; %接收信号 y=sr(L/2:L:N);
%取样，取样时刻无偏差
aa=sign(y); %判决
n_err=n_err+length(find(aa~=a)); %误码个数 end
Pe(loop1)=n_err/(M*loop2); %实际误码率 figure(1)
semilogy(Eb_N0,Pe,'g'); %实际误码率曲线 eb_n0=10.^(Eb_N0/10); hold on
semilogy(Eb_N0,0.5*erfc(sqrt(eb_n0))); % axis([0,9,1e-4,1]) xlabel('Eb/N0 in db')
理论误码率曲线 ylabel('Pe') title('误码率曲线') end 仿真图
误码率曲线：绿色线为实际曲线，蓝色为理论误码率曲线
实验结论
由上图可以看出,当取样时间偏差为Ts/4的时候,误码率明显高于无偏差的理论误码率。原因可以从信号的眼图看出,没有偏差时的取样点在眼睛睁开最大处,判决效果最佳,而偏差Ts/4之后取样信号的幅度下降,比较容易受噪声干扰,所以误码率上升。
实验六 仿真数字基带传输系统，包括输入、输出信号波形及其功率谱，眼图（升余弦滚降系数a=0.5），Pe~Eb\\No曲线，取样时间无偏差。（单极性归零码的占空比可调） 源程序
%单极性归零码通过基带系统后输出及其功率谱和眼图 global dt df t f N%全局变量 close all clear Eb_N0 Pe N=2^14; L=16;
%采样点数
%每码元的采样点数
M=N/L; Rb=2;
%码元数
%码速率是2Mb/s %码元间隔
Ts=1/Rb;
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库通原实验软件报告(4)在线全文阅读。
相关推荐：