最佳基带系统的源程序

1、最佳基带系统的源程序：(一)，用到的函数1,f2t函数   function x=f2t(X)   global dt df t f T N   %x=f2t(X)   %x为时域的取样值矢量   %X为x的傅氏变换   %X与x长度相同并为2的整幂   %本函数需要一个全局变量dt(时域取样间隔)   X=X(N/2+1:N),X(1:N/2); 

2、;  x=ifft(X)/dt;   end    2,t2f函数。   function X=t2f(x)   global dt df N t f T   %X=t2f(x)   %x为时域的取样值矢量   %X为x的傅氏变换   %X与x长度相同,并为2的整幂。   %本函数需要一个全局变量d

3、t(时域取样间隔)   H=fft(x);   X=H(N/2+1:N),H(1:N/2)*dt;   end (二),主程序   globaldt t f df N T   close all   clear Eb_N0 Pe    k=input('取样点数2k, k取13左右');   if isempty(k),

4、k=13; end   z=input('每个信号取样点数2z, z<k');   if isempty(z), z=3; end   aa=input('滚降系数=0.5');   if aa=,aa=0.5;end   N=2k;   L=2z;M=N/L;   Rb=2;%码速率是2Mb/sTs=1/Rb;%码元间隔dt=Ts/L;%时域采样间

5、隔df=1/(N*dt);%频域采样间隔T=N*dt;%截短时间Bs=N*df/2;%系统带宽f=-Bs+df/2:df:Bs;%频域横坐标t=-T/2+dt/2:dt:T/2;%时域横坐标g=sin(pi*t/Ts).*cos(pi*t*aa/Ts)./pi*t/Ts.*(1-4*t.2*aa2/Ts2);%升余弦脉冲波形GG=t2f(g);GG=abs(GG);%升余弦脉冲的傅式变换GT=sqrt(GG);GR=GT;%最佳系统的发送接收滤波器的傅式变换for l1=1:20;Eb_N0(l1)=(l1-1);%Eb/N0 in dBeb_n0(l1)=10(Eb_N0(l1)/10);E

6、b=1;n0=Eb/eb_n0(l1);%信道的噪声谱密度sita=n0*Bs;%信道中噪声功率n_err=0;%误码计数for l2=1:20;b=sign(randn(1,M);s=zeros(1,N);%产生冲激序列s(L/2:L:N)=b/dt;SS=t2f(s);S=SS.*GG;%信道的傅式变化a=f2t(S);a=real(a);%不加噪声的输出n_ch=sqrt(sita)*randn(size(t); %信道噪声N_CH=t2f(n_ch);nr=real(f2t(N_CH.*GR);%输出噪声sr=a+nr;%接收信号y=sr(L/2-L/4:L:N-L/4);%取样bb=sign(y);%判决n_err=n_err+length(find(bb=b);%错误累计endPe(l1)=n_err/(M*l2);semilogy(Eb_N0,Pe,'g');%PeEb/N0曲线画图xlabel('Eb/N0');ylabel('Pe');title('PeEb/N0曲线');eb_n0=10.(Eb_N0/10);hold onsemilogy(Eb_N0,0.5*erfc(sqr