信号相关性分析

clear

clc

close all

% 参数设置

fs = 200;

x = 1/fs:1/fs:20;

% 生成非线性信号

y1 = sin(2*pi*1*x) ; % 添加噪声使其非线性

y2 = 2*sin(2*pi*1*x) + 0.5 * randn(size(x)); % 添加噪声使其非线性

% 计算斯皮尔曼相关系数

[rank_corr, p_value_spearman] = corr(y1', y2', 'Type', 'Spearman');

disp(['斯皮尔曼相关系数: ', num2str(rank_corr)]);

disp(['p-value (Spearman): ', num2str(p_value_spearman)]);

% 绘制信号

figure;

subplot(2, 1, 1);

plot(x, y1);

title('Signal y1');

xlabel('Time (s)');

ylabel('Amplitude');

subplot(2, 1, 2);

plot(x, y2);

title('Signal y2');

xlabel('Time (s)');

ylabel('Amplitude');