一维信号裕度提取matlab代码

以下是一维信号裕度提取的MATLAB代码示例：

% 生成一个测试信号
fs = 1000; % 采样频率
t = 0:1/fs:1-1/fs; % 时间向量
x = sin(2*pi*50*t) + 2*sin(2*pi*120*t); % 信号

% 计算功率谱密度
[Pxx, f] = pwelch(x,[],[],[],fs);

% 计算相对带宽
bw = obw(x,fs);

% 计算裕度指标
margin = 10*log10(Pxx(bw)) - 10*log10(mean(Pxx));

% 显示结果
fprintf('相对带宽: %.2f%%\n', bw*100);
fprintf('裕度指标: %.2f dB\n', margin);

其中，`obw` 函数用于计算相对带宽，`pwelch` 函数用于计算功率谱密度。在上述代码中，我们生成了一个包含两个正弦波的测试信号，然后计算了该信号的相对带宽和裕度指标。你可以将这段代码替换为你的信号数据并执行，以提取信号的裕度指标。