我正在尝试使用MATLAB绘制FFT的相位。我有这个信号实际上是通过设置为半强度的调光器的电流。无论如何,那真的没关系。基本上,在我的代码中,我将信号放在一个矢量中,i。然后我在i上执行和FFT并将其存储在I中。然后我尝试获得I的大小和角度。
幅度光谱似乎是正确的,但是相位/角度没有,我无法弄清楚原因。有什么建议?我确实意识到我的代码在如何低效编写部分方面有点'狡猾'。我不是真正的MATLAB专业人士......
非常感谢帮助。
谢谢。
% ( 40/240 * sin(2*pi*50*t) for a < t < T
% waveform = {
% ( 0 for 0 < t < a
cycles = 2;
a = 90;
clear i t;
aTime = a/360;
dt = 0.0001;
t = 0;
i = 0;
for n = 0 : cycles - 1;
T = 1/50;
t0 = 0 + (n*T) : dt : T*aTime + (n*T) - dt;
t1 = T*aTime + (n*T) : dt : T/2 + (n*T) - dt;
t2 = T/2 + (n*T) : dt : T*(aTime + 1/2) + (n*T) - dt;
t3 = T*(aTime + 1/2) + (n*T) : dt : T + (n*T) - dt;
t = [t [t0, t1, t2, t3]];
i = [i zeros(1, length(t0))];
i = [i 40/240 * sin(2*pi*50*t1)];
i = [i zeros(1, length(t2))];
i = [i 40/240 * sin(2*pi*50*t3)];
end
subplot(3,2,[1 2])
hold on;
plot(t, 40/240 * sin(2*pi*50*t), ':r');
plot(t, i);
xlabel('time (sec)')
ylabel('i(t)')
title('Current through a 40W, 240V dimmed light with alpha = 90^o')
grid on;
hold off;
axis([0, T*(n(end) + 1), -0.2, 0.2]);
fs = 1/dt;
N = length(i);
df = fs/N;
f = (-fs/2) : df : (fs/2)-df;
I = fftshift(fft(i)/N);
subplot(3,2,3)
plot(f, abs(I))
axis([-1000,1000,0,0.055]);
xlabel('frequency (Hz)')
ylabel('|i(t)|')
title('Magnitude Spectrum')
grid on;
subplot(3,2,5)
plot(f, mod(unwrap(angle(I)), 2*pi))
axis([-1000, 1000, -pi, 2.5*pi]);
xlabel('Radians')
ylabel('Arg(i(t))')
title('Frequency')
grid on;
subplot(3,2,4)
hold on;
plot(t, i);
plot(t, real(0.1*exp(1i*(2*pi*50*t + 4.139))), 'm');
plot(t, real(2*0.02653*exp(1i*(2*pi*150*t + 6.268))), 'g');
plot(t, real(2*0.008844*exp(1i*(2*pi*250*t + 3.156))), 'r');
axis([0, T*(n(end) + 1), -0.2, 0.2]);
hold off;
xlabel('Time')
ylabel('Value')
title('Fourier Series Components')
grid on;
subplot(3,2,6)
hold on;
plot(t, i);
plot(t, real(0.1*exp(1i*(2*pi*50*t + 4.139))), 'm');
plot(t, real(2*0.008844*exp(1i*(2*pi*250*t + 3.156)) + 2*0.02653*exp(1i*(2*pi*150*t + 6.268)) + 0.1*exp(1i*(2*pi*50*t + 4.139))), 'm');
plot(t, real(2*0.02653*exp(1i*(2*pi*150*t + 6.268)) + 0.1*exp(1i*(2*pi*50*t + 4.139))), 'g');
axis([0, T*(n(end) + 1), -0.2, 0.2]);
hold off;
xlabel('Time')
ylabel('Value')
title('Fourier Series Sum')
grid on;
编辑: 这样做是为了在角度和幅度上都应用了ffthift。
这就是我得到的:
答案 0 :(得分:2)
必须解开FFT阶段才有意义。否则会出现2pi的不连续性。
答案 1 :(得分:0)
我认为你的问题是你正在使用fftshift作为幅度,而不是将它用于角度。我建议你要么使用它,要么就是这样,混合这样的东西通常是不好的做法......