请用MATLAB写出 正弦型、余弦型瞬时值采样比相判据算法代码，包括fcn模块，认真思考后回答

瞬时值采样比相判据算法是一种基于采样的电力系统故障检测方法。下面是正弦型、余弦型瞬时值采样比相判据算法的 MATLAB 代码实现，包括 fcn 模块。

%% 瞬时值采样比相判据算法

% 定义采样频率、采样时间、阈值等参数
fs = 1000;
Ts = 1/fs;
K = 1.5; % 阈值

% 生成测试信号
t = 0:Ts:1;
f1 = 50; % 正弦信号频率
f2 = 60; % 余弦信号频率
x1 = sin(2*pi*f1*t);
x2 = cos(2*pi*f2*t);
x = x1 + x2;

% 瞬时值采样比相判据算法
for n = 1:length(x)
    if n < 2 % 跳过前两个采样点
        continue;
    end
    % 计算瞬时值
    v1 = x(n) - x(n-1);
    v2 = abs(x(n)) + abs(x(n-1));
    % 计算采样比
    r = v1/v2;
    % 判断是否超过阈值
    if abs(r) > K
        disp(['第', num2str(n), '个采样点发生故障']);
    end
end

% fcn 模块
function [r] = sampling_ratio(x, n)
% x: 采样信号向量
% n: 当前采样点下标
% r: 采样比
    if n < 2 % 跳过前两个采样点
        r = 0;
    else
        % 计算瞬时值
        v1 = x(n) - x(n-1);
        v2 = abs(x(n)) + abs(x(n-1));
        % 计算采样比
        r = v1/v2;
    end
end

在上述代码中，我们首先定义了采样频率、采样时间和阈值等参数。接着生成了一个包含正弦信号和余弦信号的测试信号，用于模拟电力系统中的实际信号。然后使用 for 循环遍历采样信号的每个采样点，计算该采样点的瞬时值和采样比，并判断是否超过了阈值。最后，我们还定义了一个 fcn 模块，用于计算每个采样点的采样比。

需要注意的是，上述代码只是一个简单的示例，实际应用中还需要考虑信号的噪声、采样频率的选择等因素，以确保算法的准确性和可靠性。