通信原理实验报告
实验三二进制数字信号调制仿真实验
实验四模拟信号数字传输仿真实验
姓名：张力
班级：通信工程三班
学号：2011551326
实验三二进制数字信号调制仿真实验
（1）按照如上介绍的方法，分别产生一组长度为500的二进制单极性不归零信号和归零信号，存档名为Q3_1。并求分别求出它们的功率谱密度。请写出相应的MATLAB程序，将不归零信号波形及功率谱和归零信号波形及功率谱分别画在同一图形的四个子图中，将结果图保存，贴在下面的空白处。
编写的MATLAB程序为：
clear all;
close all;
A=1
fc=2; %2Hz;
N_sample=8;
N=500; %码元数
Ts=1; %1 Baud/s
dt=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_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([0 10 0 1.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,10*log10(abs(d_NRZf).^2/T));
axis([-2 2 -50 10]);ylabel('输入信号功率谱密度（dB/Hz)');
%/////////
subplot(223);%输入NRZ信号波形（单极性）
d_RZ=conv(dd,gt2); plot(t,d_RZ(1:length(t)));
axis([0 10 0 1.2]); ylabel ('输入信号'); 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([-2 2 -50 10]);ylabel('输入信号功率谱密度（dB/Hz)');
输入信号
输入信号功率谱密度（d B /H z )
输入信号
输入信号功率谱密度（d B /H z )
（2）对刚才产生的长度为500的不归零波形对载波频率为2Hz ，幅度为1的余弦信号进行OOK 调制，并求出调制信号的功率谱密度。编写程序实现之，存档名为Q3_2。要求将不归零信号波形及功率谱和OOK 调制信号波形及功率谱分别画在同一图形的四个子图中，图形名为图2，将结果图保存，贴在下面的空白处。 编写的程序为： clear all; close all; A=1
fc=2; %2Hz; N_sample=8;
N=500; %码元数 Ts=1; %1 Baud/s dt=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_sample);%NRZ 波形
figure(1)
subplot(221);%输入NRZ信号波形（单极性）
d_NRZ=conv(dd,gt);
plot(t,d_NRZ(1:length(t)));
axis([0 10 0 1.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,10*log10(abs(d_NRZf).^2/T));
axis([-2 2 -50 10]);ylabel('输入信号功率谱密度（dB/Hz)');
ht=A*cos(2*pi*fc*t);
s_2ask=d_NRZ(1:Lt).*ht;
subplot(223)
plot(t,s_2ask);
axis([0 10 -1.2 1.2]);
ylabel('OOK');
dt=t(2)-t(1);
T=t(end);
df=1/T;
N=length(s_2ask);
f=-N/2*df:df:N/2*df-df;
s_2askf=fft(s_2ask);
s_2askf=T/N*fftshift(s_2askf);
subplot(224)
plot(f,10*log10(abs(s_2askf).^2/T));
axis([-fc-4 fc+4 -50 10]);
ylabel('OOK功率谱密度（dB/Hz）');
