实训四-数字信号频带传输的仿真设计

1、实训四 数字信号频带传输的仿真设计1、 实验内容1、 基带信号采用不归零矩形脉冲，生成信号的时域波形和功率谱密度。 2、 生成信号的时域波形和功率谱密度。 3、 接收信号的星座图。 4、 仿真系统的误码率。2、 程序与仿真图1.clc;clear all;tm=1;fc=3;Nsample=100;Nnum=1000;dt=tm/Nsample;N=Nsample*Nnum;t=0:dt:Nnum*tm-dt;NFFT=2.16;Nloop=50;st1=0;st2=0;%for i=1:Nloop; d=2*randint(1,Nnum)-1; st_bb=rectpulse(d,Nsamp

2、le); st_2psk=st_bb.*sin(2*pi*fc*t); window=boxcar(length(st_bb); pxx1,f=periodogram(st_bb,window,NFFT,1/dt); pxx2,f=periodogram(st_2psk,window,NFFT,1/dt); st1=st1+pxx1;st2=st2+pxx2;endst1=st1/Nloop;st2=st2/Nloop;%figure(1);subplot(211);plot(t,st_bb);axis(0 10 -1.5 1.5);xlabel('t');ylabel(

3、9;幅度');title('双极性基带信号的时域波形');subplot(212);plot(t,st_2psk);axis(0 10 -1.5 1.5);grid on;xlabel('t');ylabel('幅度');title('2PSK已调信号波形');figure(2);subplot(211);plot(-1*flipud(f);f,0.5*flipud(st1);st1);axis(-15 15 0 1);xlabel('f');ylabel('频谱');title('

4、双极性基带信号的频谱');grid on;subplot(212);plot(-1*flipud(f);f,0.5*flipud(st2);st2);axis(-10 10 0 0.5);xlabel('f');ylabel('频谱');title('2PSK已调信号的频谱');grid on;2.clear all;N=106;d1=randint(1,N).*2-1;d2=randint(1,N).*2-1;d=d1+j*d2;Es=mean(abs(d).2);Eb=Es./2;snr_db=-2:10;for snr=-2:10;

5、 n0=Eb/(10(snr/10); sigma=sqrt(n0/2); n=sigma*randn(size(d)+j*sigma*randn(size(d);rt=d+n;y1=sign(real(rt);y2=sign(imag(rt);ber(1,snr+3)=(length(find(y1=d1)+(length(find(y2=d2)/(2.*N);ber1(1,snr+3)=berawgn(snr,'qam',4);endsemilogy(snr_db,ber,'b- *');hold on;semilogy(snr_db,ber1,'r

6、- ');grid on;xlabel('Eb/N0');ylabel('wumalv');legend('ber','ber-theory');title('QPSK系统性能曲线');3.Qm=2;N_num=105;d1=sign(randn(1,N_num);d2=sign(randn(1,N_num);d=d1+d2*j;Es=mean(abs(d).2);Eb=Es/Qm;EbN0_dB=10;ebn0=10.(EbN0_dB/10);sigma=sqrt(Eb/ebn0/2);n=sigma*r

7、andn(size(d)+j*sigma*randn(size(d);rt=d+n;plot(rt,'.');axis(-2.5 2.5 -2 2);grid on;xlabel('I');ylabel('Q');title('信噪比为10dB的QPSK星座图');4.Qm=2;N_num=105;d1=sign(randn(1,N_num);d2=sign(randn(1,N_num);d=d1+d2*j;Es=mean(abs(d).2);Eb=Es/Qm;EbN0_dB=5;ebn0=10.(EbN0_dB/10);sigma=sqrt(Eb/ebn0/2);n=sigma*randn(size(d)+j*sigma*randn(size(d);rt=d+n;plot