clear all; clc; %% 生成双极性7位M序列 X1 = 0; X2 = 0; X3 = 1; m = 350; % 重复50遍的7位单极性m序列 for i = 1:m Y1 = X1; Y2 = X2; Y3 = X3; X3 = Y2; X2 = Y1; X1 = xor(Y3, Y1); L(i) = Y1; end for i = 1:m M(i) = 1 - 2 * L(i); % 将单极性m序列变为双极性m序列 end %% 画出双极性7位M序列频谱 k = 1:1:m; figure(1); subplot(2, 1, 1); stem(k-1, M); axis([0, 7, -1, 1]); xlabel('k'); ylabel('M序列'); title('双极性7位M序列'); subplot(2, 1, 2); ym = fft(M, 4096); magm = abs(ym); % 求双极性m序列频谱 fm = (1:2048) * 200 / 2048; plot(fm, magm(1:2048) * 2 / 4096); title('双极性7位M序列的频谱'); xlabel('Hz'); %% 生成扩频前待发送二进制信息序列和扩频后序列码 N = 50; a = 0; x_rand = rand(1, N); % 产生50个0与1之间随机数 for i = 1:N if x_rand(i) >= 0.5 % 大于等于0.5的取1，小于0.5的取0 x(i) = 1; a = a + 1; else x(i) = 0; end end t = 0:N-1; figure(2); subplot(2, 1, 1); stem(t, x); title('扩频前待发送二进制信息序列'); tt = 0:349; subplot(2, 1, 2); L = 1:7*N; y = rectpulse(x, 7); s(L) = 0; for i = 1:350 % 扩频后，码率变为100/7*7=100Hz s(i) = xor(L(i), y(i)); end tt = 0:7*N-1; stem(tt, s); axis([0, 350, 0, 1]); title('扩频后的待发送序列码');

% 清空变量并清空命令窗口
clear all;
clc;

%% 生成双极性7位M序列
X1 = 0; X2 = 0; X3 = 1; % 初始化三个寄存器
m = 350; % 重复50遍的7位单极性m序列

% 生成双极性7位M序列
for i = 1:m
Y1 = X1; Y2 = X2; Y3 = X3; % 保存当前寄存器值
X3 = Y2; X2 = Y1; % 移位操作
X1 = xor(Y3, Y1); % 异或操作
L(i) = Y1; % 保存当前寄存器最低位的值
end

% 将单极性m序列变为双极性m序列
for i = 1:m
M(i) = 1 - 2 * L(i);
end

%% 画出双极性7位M序列频谱
k = 1:1:m;

figure(1); % 新建图像窗口
subplot(2, 1, 1); % 画图1
stem(k-1, M); axis([0, 7, -1, 1]); % 设置坐标轴范围
xlabel('k'); ylabel('M序列'); % 设置坐标轴标签
title('双极性7位M序列'); % 设置图像标题

subplot(2, 1, 2); % 画图2
ym = fft(M, 4096); % 对双极性m序列进行FFT
magm = abs(ym); % 求双极性m序列频谱
fm = (1:2048) * 200 / 2048; % 计算频率
plot(fm, magm(1:2048) * 2 / 4096); % 画图
title('双极性7位M序列的频谱'); xlabel('Hz');

%% 生成扩频前待发送二进制信息序列和扩频后序列码
N = 50; % 信息序列长度
a = 0;
x_rand = rand(1, N); % 产生50个0与1之间随机数

% 生成待发送二进制信息序列
for i = 1:N
if x_rand(i) >= 0.5 % 大于等于0.5的取1，小于0.5的取0
x(i) = 1;
a = a + 1;
else
x(i) = 0;
end
end

t = 0:N-1; % 生成时间序列

figure(2); % 新建图像窗口
subplot(2, 1, 1); % 画图1
stem(t, x); % 画图
title('扩频前待发送二进制信息序列');

tt = 0:349; % 生成时间序列
subplot(2, 1, 2); % 画图2
L = 1:7*N; y = rectpulse(x, 7); % 7倍扩频
s(L) = 0;
for i = 1:350 % 扩频后，码率变为100/7*7=100Hz
s(i) = xor(L(i), y(i)); % 异或操作
end

tt = 0:7*N-1; % 生成时间序列
stem(tt, s); axis([0, 350, 0, 1]); % 设置坐标轴范围
title('扩频后的待发送序列码');