clear all; clc; 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 k=1:1:m; figure(1) subplot(2,1,1) %做m序列图 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序列的频谱') %% 二进制信息序列 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('扩频后的待发送序列码'); %% BPSK调制波形 figure(3) subplot(2,1,2) fs=2000; ts=0:0.00001:3.5-0.00001;%为了使信号看起来更光滑，作图时采样频率为100kHz % ps=cos(2*pi*fs*ts); s_b=rectpulse(s,1000); %将冲激信号补成矩形信号 s_bpsk=(1-2.*s_b).*cos(2*pi*fs*ts);%扩频后信号BPSK调制时域波形，(1-2.*s_b)是1，-1序列 plot(ts,s_bpsk); xlabel('s'); axis([0.055,0.085,-1.2,1.2]) title('扩频后bpsk信号时域波形'); subplot(2,1,1) s_bb=rectpulse(x,7000); s_bpskb=(1-2.*s_bb).*cos(2*pi*fs*ts);%无扩频信号BPSK调制时域波形 plot(ts,s_bpskb); xlabel('s'); axis([0.055,0.085,-1.2,1.2]); title('扩频前bpsk信号时域波形') %% BPSK调制频谱 figure(4) N=400000; ybb=fft(s_bpskb,N); %无扩频信号BPSK调制频谱 magb=abs(ybb); fbb=(1:N/2)*100000/N; subplot(2,1,1) plot(fbb,magb(1:N/2)*2/N); axis([1700,2300,0,0.8]); title('扩频前调制信号频谱图'); xlabel('Hz'); subplot(2,1,2) yb=fft(s_bpsk,N); %扩频信号BPSK调制频谱 mag=abs(yb); fb=(1:N/2)*100000/N; plot(fb,mag(1:N/2)*2/N); axis([1700,2300,0,0.8]); title('扩频后调制信号频谱图'); xlabel('Hz');

好的，上面是MATLAB代码，可以实现您的需求。您可以将代码复制到MATLAB中执行，以完成仿真实验。

代码中主要实现了以下功能：

1. 生成7位双极性m序列，并画出其频谱图。

2. 生成50个随机的0和1，并将其扩频成100/7的码率，画出扩频前和扩频后的二进制信息序列。

3. 对扩频后的序列进行BPSK调制，并画出扩频前和扩频后的BPSK调制信号时域波形。

4. 画出扩频前和扩频后的BPSK调制信号的频谱图。

注意：代码中的参数和变量需要根据具体情况进行调整，例如N、m、fs、ts等。

希望能对您有所帮助。

相关问题

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('扩频后的待发送序列码');

clc,clear,close all x1=[0,0];x2=[18.9875,0]; f = norm(x2 - x1); v1 = (x2 - x1) / f; u = d*v1; x3 = x1 + u; 函数或变量 'd' 无法识别。 出错 Untitled3 (第 5 行) u = d*v1;

很抱歉，我在之前的回答中遗漏了一个变量d的定义。d代表了x3到x1的距离，你需要根据具体情况给d赋值。在你的代码中，你可以在计算u之前先定义d的值。下面是修改后的代码示例：

clc, clear, close all
x1 = [0, 0];
x2 = [18.9875, 0];
d = 10;  % 假设x3到x1的距离为10，根据实际情况修改该值
f = norm(x2 - x1);
v1 = (x2 - x1) / f;
u = d * v1;
x3 = x1 + u;

disp(x3);  % 输出x3的坐标

在这个示例中，我假设x3到x1的距离为10。你可以根据你的实际情况修改这个值。运行代码后，将会在命令窗口中输出x3的坐标。希望这次能帮到你！