如果是这样的话,可以解释为什么正弦激励会产生100Hz的分量了,主要是因为采集的数据中把正弦信号负半轴的部分给截掉了
仿真程序如下:
clear
fs = 5000;
N = 1000;
a1 = ones(N/(50*2),50);
a2 = repmat((1:50)/50, N/(50*2),1);
a3 = repmat(sin(2*pi*50/fs*(1:50)), N/(50*2),1);
z = zeros(N/(50*2),50);
a = [ a1 z ].';
b1 = a(:);
a = [ a2 z ].';
b2 = a(:);
a = [ a3 z ].';
b3 = a(:);
bfft1 = abs(fft(b1-mean(b1)));
bfft2 = abs(fft(b2-mean(b2)));
bfft3 = abs(fft(b3-mean(b3)));
figure(1); hold on;
plot(b1, 'r');
plot(b2, 'g');
plot(b3, 'b');
figure(2); hold on;
plot(0:fs/N:fs/2-fs/N, bfft1(1:N/2), 'r')
plot(0:fs/N:fs/2-fs/N, bfft2(1:N/2), 'g')
plot(0:fs/N:fs/2-fs/N, bfft3(1:N/2), 'b')
To aiyou18:
最好把你测量的空载时候的正弦信号数据也放上来。
[ 本帖最后由 xray 于 2007-10-24 19:56 编辑 ] |