生成带宽为9kHz到15kHz的63阶gold序列matlab代码

以下是生成带宽为9kHz到15kHz的63阶gold序列matlab代码：

clear;clc;
N = 63; % 序列长度
M1 = 6; % 编码器1阶数
M2 = 4; % 编码器2阶数

%% 初始化
s1 = ones(1, M1); % 编码器1初始状态
s2 = ones(1, M2); % 编码器2初始状态
c = ones(1, N); % Gold序列

%% 编码器1、编码器2和组合
for i = 1:N
f1 = xor(s1(2), s1(6)); % 编码器1输出
f2 = xor(s2(3), s2(4)); % 编码器2输出
c(i) = xor(f1, f2); % 组合输出
s1 = [f1 s1(1: M1 - 1)]; % 更新编码器1状态
s2 = [f2 s2(1: M2 - 1)]; % 更新编码器2状态
end

%% 生成正弦波
fs = 48e3; % 采样频率
f0 = 12e3; % 正弦波频率
t = (0: N - 1) / fs; % 时间轴
s = sin(2 * pi * f0 * t); % 正弦波信号

%% 乘上Gold序列
c = 2 * c - 1; % 将0,1序列变为1,-1序列
c = repmat(c, length(s), 1); % 扩展成与正弦波长度相同
s_new = s' .* c(:); % 信号与序列相乘，注意需要转置

%% 播放信号
soundsc(s_new, fs); % 播放结果音频信号

%% 绘制结果信号
figure;subplot(2, 1, 1);plot(s);title('原始信号');
subplot(2, 1, 2);plot(s_new);title('加扰后信号');