单极性码间干扰与信道均衡

function [p]=smldPe2(snr_in_dB) % [p]=smldPe2(snr_in_dB) % SMLDPE仿真求出以dB为单位的特定信噪比的误码率 E=1; SNR=exp(snr_in_dB*log(10)/10); %信噪比 sgma=E/sqrt(2*SNR); %噪声的标准偏差 M=10000; %码元数 L=32; %每个码元采样点数 N=M*L; %采样点数 Rb=2; %码速率 Ts=1/Rb; %码元间隔 dt=Ts/L; %时域采样间隔 df=1/(N*dt); %频域采样间隔 d=16; %半个码元采样点数 T=N*dt; %截短时宽 Bs=N*df/2; %频带宽度 t=[-T/2+dt/2:dt:T/2]; %时域坐标 t1=[-T/2+dt/2:32*dt:T/2]; %时域坐标 f=[-Bs+df/2:df:Bs]; %频域坐标 alpha=1; %滚降系数 g1=sin(pi*t/Ts)./(pi*t/Ts); g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2); g=g1.*g2; G=sqrt(t2f(g,dt)); %发送和接收滤波器频谱 Cf=1-0.5*exp(-2*i*pi*f*d); %信道传输函数 Ge=1./Cf;%信道均衡 Gf=G.^2.*Cf.*Ge; %系统总的频谱 %求系统后的码元样本序列 a=sign(randn(1,M)); imp=zeros(1,N); imp(L/2:L:N)=a/dt; S=t2f(imp,dt).*Gf; s=f2t(S,dt); s1=zeros(1,N); %此段程序将经傅立叶变换得到的数字序列 s1(N/2+1:N)=s(1:N/2);%在中轴反转,如111-1-11,经变换变成-1-11111 s1(1:N/2)=s(N/2+1:N);%因此应把它反过来,这主要是由 s1=real(s1); %f2t和t2f函数引起的. gr=s1+sgma*randn(1,N); %信号加噪声，实现SNR gr1=gr(16:32:N); %取每个码元中间的一个样本 %随后为判决和误码概率的计算 numoferr=0; for ii=1:M, %随后为判决 if (gr1(ii)<0), decis=-1; %判决为“0” else decis=1; %判决为“1” end; if(decis~=a(ii)), numoferr=numoferr+1; %如果出错，计数器加1 end; end; p=numoferr/M; %误差概率估算