读取数据表格，这个数据是一个连续波，对这个连续波行采用goertzel算法提取二次谐波，在mtlab function中实现

在MATLAB中，Goertzel算法是一种快速频率分量检测的技术，特别适用于实时信号处理，如音频分析或测量连续波形的特定频率成分。如果你有一个包含连续波的数据表格，并想要提取其中的二次谐波，你可以按照以下步骤编写函数：

function [harmonic, residual] = goertzel_continuous_wave(data, frequency_of_interest)
% 初始化参数
window_size = length(data); % 数据长度
sample_rate = 1; % 假设采样率是每秒一次，如果实际不是，需要替换为实际值
alpha = 0.54 + 0.46i; % Goertzel算法的系数

% 计算Goertzel滤波器
k = (0:(window_size - 1)) * sample_rate;

% 应用Goertzel算法
y = zeros(window_size, 1);
for n = 1:window_size
    y(n) = alpha .* conj(y(n - 1)) + data(n);
end

% 获取二次谐波（通常为频率_of_interest*2）
harmonic = abs(y(round(frequency_of_interest * window_size))) / sqrt(window_size);

% 计算剩余信号（去除二次谐波后的残余部分）
residual = data - harmonic * cos(2 * pi * frequency_of_interest * k);

% 返回结果
end

在这个函数中，`data`是输入的连续波信号，`frequency_of_interest`是你想要寻找的二次谐波频率。函数首先计算了Goertzel滤波器系数，然后应用该算法并查找指定频率处的幅度值。最后返回的是提取出的二次谐波值和原始信号的残余部分。