matlab中对下main的代码进行注释，N=length(x_1); f=fs*(0:(N/2))/N; X_1=fft(x_1)/N; X_2=fft(x_2)/N; X_3=fft(x_3)/N; Y_1=fft(y_1)/N; Y_2=fft(y_2)/N; Y_3=fft(y_3)/N;

% 获取信号长度 N = length(x_1); % 计算频率向量 f = fs * (0:(N/2)) / N; % 对x_1进行FFT，并归一化（除以信号长度） X_1 = fft(x_1) / N; % 对x_2进行FFT，并归一化（除以信号长度） X_2 = fft(x_2) / N; % 对x_3进行FFT，并归一化（除以信号长度） X_3 = fft(x_3) / N; % 对y_1进行FFT，并归一化（除以信号长度） Y_1 = fft(y_1) / N; % 对y_2进行FFT，并归一化（除以信号长度） Y_2 = fft(y_2) / N; % 对y_3进行FFT，并归一化（除以信号长度） Y_3 = fft(y_3) / N;