通信原理 单 双极性码的归零和非归零01等概下的波形和功率谱密度的matlab实现

在实际的基带传输系统中,对在信道中传输的电波形的要求：

（1） 不含直流成分，且低频分量应尽量少。

（2） 含有丰富的定时信息，以便从接收码流提取定时信息。

（3） 功率谱主瓣宽度窄，以节省传输频带。

（4） 不受信源统计特性的影响，能适应信源的变化。

（5） 具有内在的检错能力，码型具有一定规律性，以便进行宏观监测。

（6） 编译码简单，以降低通信延时和成本。

1、单极性归零和非归零的01等概出现时的波形和功率谱密度

Ts = 1; % 码元周期N_sample = 1024; % 单个码元抽样点数dt = Ts / N_sample; % 抽样时间间隔N = 100; % 码元数t = 0 : dt : (N * N_sample - 1) * dt; % 序列传输时间gt1 = ones(1, N_sample); % NRZRanNum = randperm(N , N/2);%随机抽取四分之一的数抬到高电平se = []for i = 1 : Nflag = 0;for j = 1 :N / 2if RanNum(j)==iflag = 1;endendif flag == 1se = [se,gt1];elsese = [se,zeros(1, N_sample)];endendfft_se = fftshift(fft(se)); PE = 10 * log10(abs(fft_se) .^ 2 / (N * Ts));PEL = (-length(fft_se) / 2 : length(fft_se) / 2 - 1) / 100*(1/Ts);subplot(2, 1, 1);plot(t, se);grid on;axis([0 20 -1.5 1.5]);title('单极性NRZ01等概');subplot(2, 1, 2);plot(PEL, PE); grid on; axis([-20 20 -100 100]); title('单极性NRZ01等概');

Ts = 1; % 码元周期N_sample = 1024; % 单个码元抽样点数dt = Ts / N_sample; % 抽样时间间隔N = 100; % 码元数t = 0 : dt : (N * N_sample - 1) * dt; % 序列传输时间gt1 = [ones(1, N_sample / 2), zeros(1, N_sample / 2)]; % NRZRanNum = randperm(N , N/2);%随机抽取二分之一的数抬到高电平se = []for i = 1 : Nflag = 0;for j = 1 :N / 2if RanNum(j)==iflag = 1;endendif flag == 1se = [se,gt1];elsese = [se,zeros(1, N_sample)];endendfft_se = fftshift(fft(se)); PE = 10 * log10(abs(fft_se) .^ 2 / (N * Ts));PEL = (-length(fft_se) / 2 : length(fft_se) / 2 - 1) / 100*(1/Ts);subplot(2, 1, 1);plot(t, se);grid on;axis([0 20 -1.5 1.5]);title('单极性RZ01等概');subplot(2, 1, 2);plot(PEL, PE); grid on; axis([-20 20 -100 100]); title('单极性RZ01等概');

2、双极性归零和非归零的01等概出现时的波形和功率谱密度

Ts = 1; % 码元周期N_sample = 1024; % 单个码元抽样点数dt = Ts / N_sample; % 抽样时间间隔N = 100; % 码元数t = 0 : dt : (N * N_sample - 1) * dt; % 序列传输时间gt1 = ones(1, N_sample); % NRZRanNum = randperm(N , N/2);%随机抽取二分之一的数抬到高电平se = []for i = 1 : Nflag = 0;for j = 1 :N / 2if RanNum(j)==iflag = 1;endendif flag == 1se = [se,gt1];elsese = [se,(-1)*ones(1, N_sample)];endendfft_se = fftshift(fft(se)); PE = 10 * log10(abs(fft_se) .^ 2 / (N * Ts));PEL = (-length(fft_se) / 2 : length(fft_se) / 2 - 1) / 100*(1/Ts);subplot(2, 1, 1);plot(t, se);grid on;axis([0 20 -1.5 1.5]);title('双极性NRZ01等概');subplot(2, 1, 2);plot(PEL, PE); grid on; axis([-20 20 -100 100]); title('双极性NRZ01等概');

Ts = 1; % 码元周期N_sample = 1024; % 单个码元抽样点数dt = Ts / N_sample; % 抽样时间间隔N = 100; % 码元数t = 0 : dt : (N * N_sample - 1) * dt; % 序列传输时间gt1 = [ones(1, N_sample / 2), zeros(1, N_sample / 2)]; % RZRanNum = randperm(N , N/2);%随机抽取二分之一的数抬到高电平se = []for i = 1 : Nflag = 0;for j = 1 :N / 2if RanNum(j)==iflag = 1;endendif flag == 1se = [se,gt1];elsese = [se,(-1)*[ones(1, N_sample / 2), zeros(1, N_sample / 2)]];endendfft_se = fftshift(fft(se)); PE = 10 * log10(abs(fft_se) .^ 2 / (N * Ts));PEL = (-length(fft_se) / 2 : length(fft_se) / 2 - 1) / 100*(1/Ts);subplot(2, 1, 1);plot(t, se);grid on;axis([0 20 -1.5 1.5]);title('双极性RZ01等概');subplot(2, 1, 2);plot(PEL, PE); grid on; axis([-20 20 -100 100]); title('双极性RZ01等概');