Fs=8192; %采样频率
T=1/Fs;
N=8192; %序列长度
f1=400;
f2=1200;
f3=2400;
t=(0:1:N-1)T;
signal=1+2sin(2pif1t)+3sin(2pif2t+pi/2)+5sin(2pif3*t+pi/4); %生成输入信号
% 假设你的信号是一个列向量,名为signal。
% 假设你的采样率为Fs。
% 确保你已经加载了信号数据,并将其存储在signal变量中。
% 计算信号长度和FFT长度
N = length(signal);
L = 2^nextpow2(N);
% 进行FFT变换
Y = fft(signal, L);
f = Fs*(0:(L/2))/L;
% 提取幅频特性
A = abs(Y/L);
% 获取频谱中的奇数和偶数索引
oddFreqIdx = mod(1:length(f), 2) == 1; % 奇数索引
evenFreqIdx = mod(1:length(f), 2) == 0; % 偶数索引
% 提取奇数和偶数频率成分
oddFreqs = f(oddFreqIdx);
oddAmplitudes = A(oddFreqIdx);
evenFreqs = f(evenFreqIdx);
evenAmplitudes= A(evenFreqIdx);
% 显示奇数和偶数频率成分
disp(‘奇数频率成分:’);
disp(‘频率(Hz) 幅值’);
for i = 1:length(oddFreqs)
disp([num2str(oddFreqs(i)) ’ ’ num2str(oddAmplitudes(i))]);
end
disp(‘偶数频率成分:’);
disp(‘频率(Hz) 幅值’);
for i = 1:length(evenFreqs)
disp([num2str(evenFreqs(i)) ’ ’ num2str(evenAmplitudes(i))]);
end