写一个频域恒虚警matlab程序

以下是一个简单的频域恒虚警Matlab程序，用于检测输入信号是否超出了预定的阈值：

% 设置参数
fs = 1000; % 采样率
N = 1024; % FFT长度
fc = 50; % 信号频率
A = 1; % 信号幅值
noise_amp = 0.2; % 噪声幅值
threshold = 0.5; % 阈值

% 生成信号和噪声
t = (0:N-1)/fs;
signal = A*sin(2*pi*fc*t);
noise = noise_amp*randn(size(signal));

% 合成信号
x = signal + noise;

% 计算频域参数
X = fft(x,N);
f = (0:N-1)*fs/N;
Pxx = abs(X).^2/N;

% 检测是否超过阈值
if max(Pxx) > threshold
    disp('警告：频域幅值超过阈值');
else
    disp('通过测试：频域幅值未超过阈值');
end

% 绘制频谱图
plot(f,Pxx);
xlabel('频率 (Hz)');
ylabel('功率谱');
title('频域分析');

该程序生成一个正弦信号，并添加一定的噪声。然后，计算信号的FFT，并计算频率域的功率谱密度。最后，程序检查是否有任何频率的幅值超过了预设的阈值，并输出相应的警告或通过测试的消息。程序还绘制了频谱图，以便更直观地观察信号的频域特性。