锁相放大器 一次谐波及二次谐波 matlab

以下是用MATLAB实现锁相放大器提取一次谐波和二次谐波的示例代码：

% 生成参考信号和输入信号
t = linspace(0, 1, 1000);
f_in = 10;  % 输入信号频率
ref_freq = [f_in * 2, f_in * 4];  % 参考信号频率
ref_signal = [sin(2 * pi * freq * t) for freq = ref_freq];
input_signal = sin(2 * pi * f_in * t);

% 相乘并积分
products = [input_signal .* ref for ref = ref_signal];
integrals = [trapz(product, t) for product = products];

% 计算幅度和相位
amplitudes = [2 * abs(integral) / length(t) for integral = integrals];
phases = [angle(integral) for integral = integrals];

disp('Amplitude of 1st harmonic: '); disp(amplitudes(1));
disp('Phase of 1st harmonic: '); disp(phases(1));
disp('Amplitude of 2nd harmonic: '); disp(amplitudes(2));
disp('Phase of 2nd harmonic: '); disp(phases(2));

以上代码与Python实现的代码类似，生成参考信号和输入信号，并将参考信号设置为输入信号的两倍和四倍频率。将输入信号和参考信号相乘并积分，然后计算出一次谐波和二次谐波的幅度和相位信息。需要注意的是，在MATLAB中，数组下标是从1开始的，因此需要将Python实现代码中的下标从0开始改为从1开始。