% 计算基波频率 fund_freq = 50; % 计算谐波频率 harmonic_freqs = [fund_freq * i for i in 2:6]; % 设计FIR滤波器 nyquist_freq = 0.5 * fs; cutoff_freq = 2 * max(harmonic_freqs); numtaps = 101; taps = fir1(numtaps, cutoff_freq / nyquist_freq); % 滤波数据 filtered_data = filter(taps, 1, y); % 绘制原始数据和滤波后数据 t = 0:1/fs:(length(data)-1)/fs;

这段代码是用来计算基波频率和谐波频率，设计一个FIR滤波器来滤波数据，然后绘制原始数据和滤波后数据的图像。其中，fund_freq变量表示基波频率为50Hz，harmonic_freqs变量表示谐波频率为基波频率的2到6倍，nyquist_freq表示采样频率的一半，cutoff_freq表示截止频率为谐波频率的最大值的两倍，numtaps表示滤波器的阶数，taps表示FIR滤波器的系数，filtered_data表示滤波后的数据，t表示时间序列。