资源描述:
《通信原理实验报告二》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、通信原理实验报告实验三二进制数字信号调制仿真实验实验四模拟信号数字传输仿真实验姓名:田湘班级:学号:2012550732实验三二进制数字信号调制仿真实验(1)按照如上介绍的方法,分别产生一组长度为500的二进制单极性不归零信号和归零信号,存档名为Q3_1。并求分别求出它们的功率谱密度。请写出相应的MATLAB程序,将不归零信号波形及功率谱和归零信号波形及功率谱分别画在同一图形的四个子图中,将结果图保存,贴在下面的空白处。编写的MATLAB程序为:clearall;closeall;A=1fc=2;%2Hz;
2、N_sample=8;N=500;%码元数Ts=1;%1Baud/sdt=Ts/fc/N_sample;%波形采样间隔t=0:dt:N*Ts-dt;Lt=length(t);%产生二进制信源d=sign(randn(1,N));%dd=sigexpand((d+1)/2,fc*N_sample);N1=length((d+1)/2);dd=zeros(fc*N_sample,N1);dd(1,:)=(d+1)/2;dd=reshape(dd,1,fc*N_sample*N1);gt=ones(1,fc*N_
3、sample);%NRZ波形gt2=ones(1,fc*N_sample/2)%RZboxing/////figure(1)subplot(221);%输入NRZ信号波形(单极性)d_NRZ=conv(dd,gt);plot(t,d_NRZ(1:length(t)));axis([01001.2]);ylabel('输入信号');subplot(222);%输入NRZ频谱dt=t(2)-t(1);T=t(end);df=1/T;N=length(d_NRZ(1:length(t)));f=-N/2*df:d
4、f:N/2*df-df;d_NRZf=fft(d_NRZ(1:length(t)));d_NRZf=T/N*fftshift(d_NRZf);plot(f,10*log10(abs(d_NRZf).^2/T));axis([-22-5010]);ylabel('输入信号功率谱密度(dB/Hz)');%/////////subplot(223);%输入NRZ信号波形(单极性)d_RZ=conv(dd,gt2);plot(t,d_RZ(1:length(t)));axis([01001.2]);ylabel('
5、输入信号');subplot(224);%输入NRZ频谱dt=t(2)-t(1);T=t(end);df=1/T;N=length(d_RZ(1:length(t)));f=-N/2*df:df:N/2*df-df;d_RZf=fft(d_RZ(1:length(t)));d_RZf=T/N*fftshift(d_RZf);plot(f,10*log10(abs(d_RZf).^2/T));axis([-22-5010]);ylabel('输入信号功率谱密度(dB/Hz)');(2)对刚才产生的长度为500
6、的不归零波形对载波频率为2Hz,幅度为1的余弦信号进行OOK调制,并求出调制信号的功率谱密度。编写程序实现之,存档名为Q3_2。要求将不归零信号波形及功率谱和OOK调制信号波形及功率谱分别画在同一图形的四个子图中,图形名为图2,将结果图保存,贴在下面的空白处。编写的程序为:clearall;closeall;A=1fc=2;%2Hz;N_sample=8;N=500;%码元数Ts=1;%1Baud/sdt=Ts/fc/N_sample;%波形采样间隔t=0:dt:N*Ts-dt;Lt=length(t);%
7、产生二进制信源d=sign(randn(1,N));%dd=sigexpand((d+1)/2,fc*N_sample);N1=length((d+1)/2);dd=zeros(fc*N_sample,N1);dd(1,:)=(d+1)/2;dd=reshape(dd,1,fc*N_sample*N1);gt=ones(1,fc*N_sample);%NRZ波形figure(1)subplot(221);%输入NRZ信号波形(单极性)d_NRZ=conv(dd,gt);plot(t,d_NRZ(1:leng
8、th(t)));axis([01001.2]);ylabel('输入信号');subplot(222);%输入NRZ频谱dt=t(2)-t(1);T=t(end);df=1/T;N=length(d_NRZ(1:length(t)));f=-N/2*df:df:N/2*df-df;d_NRZf=fft(d_NRZ(1:length(t)));d_NRZf=T/N*fftshift(d_NRZf);plot(f,1