数字基带信号波形仿真

通信原理上机实验报告年级：姓名：学号：时间：数字基带信号波形仿真一、实验目的1.熟悉MATLAB软件的工作环境二、实验原理数字基带信号的波形经常采用方波，其中最基本的二进制基带信号波形有单极性归零波形、单极性不归零波形、双极性归零波形、双极性不归零波形。三、实验内容与结果（1）数字基带信号波形的MATLAB仿真下面通过MATLAB程序来仿真一串随机消息代码的基带信号波形、首先产生1000个随机信号序列，分别用单极性归零码、单极性不归零码、双极性归零码和双极性不归零码编码，并且求平均功率谱密度。源代码（以双极性为例）如下：closeallclearallk=14;L=32;N=2^k;M=N/L;dt=1/L;T=N*dt;df=1.0/T;Bs=N*df/2;t=linspace(-T/2,T/2,N);f=linspace(-Bs,Bs,N);EP1=zeros(size(f));EP2=zeros(size(f));EP3=zeros(size(f));forx=1:1000k=round(rand(1,M));nrz=zeros(L,M);rz=zeros(L,M);fori=1:Mifk(i)==1nrz(:,i)=1;rz(1:L/2,i)=1;elsenrz(:,i)=-1;rz(1:L/2,i)=-1;endendnrz=reshape(nrz,1,N);rz=reshape(rz,1,N);NRZ=t2f(nrz,dt);P1=NRZ.*conj(NRZ)/T;RZ=t2f(rz,dt);P2=RZ.*conj(RZ)/T;EP1=(EP1*(x-1)+P1)/x;EP2=(EP2*(x-1)+P2)/x;endfigure(1)subplot(2,2,1);plot(t,nrz)axis([-5,5,min(nrz)-0.1,max(nrz)+0.1])title('Ë«¼«ÐÔ²»¹éÁãÂë','fontsize',12)xlabel('t(ms)','fontsize',12)ylabel('nrz(t)','fontsize',12)gridonsubplot(2,2,2);plot(t,rz)axis([-5,5,min(rz)-0.1,max(rz)+0.1])title('Ë«¼«ÐÔ¹éÁãÂë','fontsize',12)xlabel('t(ms)','fontsize',12)ylabel('rz(t)','fontsize',12)gridonsubplot(2,2,3);plot(f,EP1)axis([-5,5,0,1.2])title('Ë«¼«ÐÔ²»¹éÁãÂ빦ÂÊÆ×ÃܶÈͼ','fontsize',12)xlabel('f(kHz)','fontsize',12)ylabel('P1(f)','fontsize',12)gridonsubplot(2,2,4);plot(f,EP2)axis([-5,5,0,0.3])title('Ë«¼«ÐÔ¹éÁãÂ빦ÂÊÆ×ÃܶÈͼ','fontsize',12)xlabel('f(kHz)','fontsize',12)ylabel('P2(f)','fontsize',12)gridon调用傅里叶变换的函数t2f,该函数定义如下：functionX=t2f(x,dt)X=fftshift(fft(x))*dt;end(1)实验结果1.数字基带信号波形的Simulink仿真（1）Simulink仿真模型（2）数字基带信号波形的Simulink仿真参数模块名称参数名称参数取值BernoulliBinaryGeneratorProbabilityofazero0.5Initialseed61Sampletime1sSwitchCriteriaforpassingfirstinputU2=ThresholdThreshold0.5Switch1CriteriaforpassingfirstinputU2=ThresholdThreshold0PulseGeneratorPeriod1sPulseWidth50(%ofperiod)（3）实验结果