matlab 生成两个随机的LoRa信号并求第一个信号的自相关和两个信号的互相关，输出两个信号的图像和相关性的图像

以下是生成两个随机的LoRa信号并求第一个信号的自相关和两个信号的互相关的Matlab代码：

% 生成两个随机的LoRa信号
Fs = 10e3; % 采样率
t = 0:1/Fs:1; % 时间向量

% 第一个信号
f1 = 100; % 载波频率
BW1 = 50; % 带宽
SF1 = 10; % 扩频因子
signal1 = lora_signal(Fs, f1, BW1, SF1, length(t));

% 第二个信号
f2 = 200; % 载波频率
BW2 = 30; % 带宽
SF2 = 8; % 扩频因子
signal2 = lora_signal(Fs, f2, BW2, SF2, length(t));

% 求第一个信号的自相关
Rxx = xcorr(signal1);
t_Rxx = linspace(-1,1,length(Rxx));

% 求两个信号的互相关
Rxy = xcorr(signal1,signal2);
t_Rxy = linspace(-1,1,length(Rxy));

% 绘制信号和相关性图像
figure;
subplot(2,2,1);
plot(t, signal1);
title('Signal 1');
xlabel('Time (s)');
ylabel('Amplitude');
subplot(2,2,2);
plot(t, signal2);
title('Signal 2');
xlabel('Time (s)');
ylabel('Amplitude');
subplot(2,2,3);
plot(t_Rxx, Rxx);
title('Auto-Correlation of Signal 1');
xlabel('Time Lag (s)');
ylabel('Correlation Coefficient');
subplot(2,2,4);
plot(t_Rxy, Rxy);
title('Cross-Correlation of Signal 1 and Signal 2');
xlabel('Time Lag (s)');
ylabel('Correlation Coefficient');

其中，`lora_signal`是一个自定义函数，用于生成LoRa信号的代码如下：

function signal = lora_signal(Fs, f, BW, SF, N)
% Fs: 采样率
% f: 载波频率
% BW: 带宽
% SF: 扩频因子
% N: 信号长度

t = (0:N-1)/Fs; % 时间向量
M = 2^SF; % 星座点数
symbols = randi([0,M-1],1,N/SF); % 随机生成星座点
data = de2bi(symbols,'left-msb'); % 将星座点转换为二进制数据
data = reshape(data',[],1); % 将数据串转换为列向量

% 生成Chirp信号
t_chirp = linspace(0,1/SF,N/SF);
f_dev = BW/t(end);
chirp_signal = chirp(t_chirp,0,1/SF,f_dev,'quadratic');
chirp_signal = chirp_signal.*exp(1j*2*pi*f*t_chirp);

% 扩频
PN = round(rand(1,N)*2-1); % 生成随机的PN序列
PN = repmat(PN,SF,1); % 将PN序列复制SF次
PN_signal = PN(:)';
PN_signal = PN_signal*sqrt(SF); % 归一化
PN_signal = PN_signal.*exp(1j*angle(chirp_signal)); % 与Chirp信号相乘

% 调制
signal = real(PN_signal.*exp(1j*2*pi*f*t));
signal = signal.*sqrt(SF); % 归一化
signal = signal.*data; % 与二进制数据相乘

运行以上代码，将会得到以下结果：