基于matlab设计的数字基带传输系统

代码能仿真出基础的数字基带信号传输过程，采用双极性不归零编码，其中奈奎斯特滤波器是通过加载由fda工具生成的mat文件。在仿真过程中，加入的高斯白噪声越大，误码率越高。滤波器的滚降系数越小，误码率越高。clc;clearall;N=100;%生成的比特个数n=16;%每个比特的抽样点数signal=rand(1,N)0.5;%生成低密度的个数为N的01比特流Rb=16;Ts=1/(Rb*n);%抽样间隔%画原信号t1=0:1/Rb:(N-1)/Rb;figure(1)subplot(6,1,1)plot(t1,signal,'b.')axis([0,20/Rb,-1.2,1.2])gridontitle('原始信号')%双极性编码code_signal=ones(1,1600);forj=1:100ifsignal(j)==0fori=0:15code_signal(16*j-i)=-1;endendend%扩展了原信号的抽样点数，并生成双极性编码%画双极性编码信号t2=0:Ts:(N*n-1)*Ts;%T=N*Rb为生成的信号的长度subplot(6,1,2)plot(t2,code_signal,'b.')axis([0,20/Rb,-1.2,1.2])gridontitle('双极性编码')transmissionsignal=awgn(code_signal,10);%加入高斯白噪声%画出高斯白噪声后的信号subplot(6,1,3)plot(t2,transmissionsignal)axis([0,20/Rb,-1.2,1.2])gridontitle('加入高斯白噪声后的信号')%奈奎斯特滤波器myNyquist=load('myNyquist0.5.mat');%Rolloff=0.5,fs=256,Band=8den=1;filtersignal=filter(myNyquist.Num,den,transmissionsignal);%经升余弦滚降系统之后信号subplot(6,1,4)plot(t2,filtersignal)axis([0,20/Rb,-2,2])gridontitle('经升余弦滚降系统之后信号')%抽样判决samplesignal=[];fori=n/2:n:length(filtersignal)iffiltersignal(i)0forj=1:nsamplesignal=[samplesignal,1];endelseforj=1:nsamplesignal=[samplesignal,-1];endendend%画抽样判决之后的信号subplot(6,1,5)plot(t2,samplesignal,'b.')gridonaxis([0,20/Rb,-1.2,1.2])title('抽样判决之后的信号')%反向恢复原编码ssignal=[];fori=n/2:n:length(samplesignal)ifsamplesignal(i)0ssignal=[ssignal,1];00.20.40.60.811.2-101原始信号00.20.40.60.811.2-101双极性编码00.20.40.60.811.2-101加入高斯白噪声后的信号00.20.40.60.811.2-202经升余弦滚降系统之后信号00.20.40.60.811.2-101抽样判决之后的信号00.20.40.60.811.2-101恢复信号elsessignal=[ssignal,0];endend%恢复信号subplot(6,1,6)plot(t1,ssignal,'b.')gridonaxis([0,20/Rb,-1.2,1.2])title('恢复信号')%误码率的计算symscntcnt=0;forc=1:100ifsignal(c)==ssignal(c)cnt=cnt+0;elsecnt=cnt+1;endenderror=cnt/n仿真图：