用matlab计算单相瞬时电压的实时频率并画出频率关于时间变化的图

假设已经获取到单相瞬时电压的时间序列数据，可以使用matlab中的fft函数对其进行傅里叶变换，然后计算出其实时频率。

以下是一个示例代码：

% 假设已经获取到单相瞬时电压的时间序列数据 voltage_data，采样率为 fs
N = length(voltage_data); % 采样点数
f = (0:N-1)*fs/N; % 频率向量
Y = fft(voltage_data); % 进行傅里叶变换
P = abs(Y/N).^2; % 计算功率谱密度
P(2:end-1) = 2*P(2:end-1); % 修正幅度
[max_val, max_idx] = max(P); % 找到功率谱密度最大值
freq = f(max_idx); % 计算实时频率

% 画出频率关于时间变化的图
t = (0:N-1)/fs; % 时间向量
freq_vec = zeros(1, N);
for i = 1:N
    Y = fft(voltage_data(1:i));
    P = abs(Y/i).^2;
    P(2:end-1) = 2*P(2:end-1);
    [max_val, max_idx] = max(P);
    freq_vec(i) = f(max_idx);
end
plot(t, freq_vec);
xlabel('时间（秒）');
ylabel('频率（Hz）');

这段代码中，首先计算出采样率为 fs 的时间序列数据的频率向量 f，然后对其进行傅里叶变换，得到频域信号 Y。接着计算出功率谱密度 P，并修正幅度。然后找到功率谱密度最大值所对应的频率值，即为实时频率 freq。

接着，使用for循环遍历时间序列数据，每次只取前i个数据进行傅里叶变换，并计算出实时频率。将所有实时频率值存储在freq_vec中，最后画出频率关于时间变化的图