- 积分
- 0
- 注册时间
- 2010-2-23
- 仿真币
-
- 最后登录
- 1970-1-1
|
在此恳请讯号分析高手相助 拜托 谢谢请指导小弟一下
像是小弟有段程序
clear all
clc
%%%%%%%时域DATA
length = 1000;
t = 0:0.00001:0.1;
f1 =60;
f2 = 180;
I1 = sin(2*pi*f1*t) ;
I2 = 0.2*sin(2*pi*f2*t) ;
Itotal = I1+I2;
subplot(211)
plot(t,Itotal)
%%%%%%%频域分析
N = 1024;
Fsampling = 0:60:1200
EVud2 = Itotal(1:N)
Vfft = fft(EVud2);
EVudfft = [abs(Vfft(1:N))*1/N]; %波形DATA做FFT后取大小
%%%%%%%看图验证结果
subplot(212)
for ffti = 1:20
frequency(ffti,1) = Fsampling(ffti);
dia_mag_Vfft(ffti,1) = EVudfft(ffti);
end
stem(frequency,dia_mag_Vfft);
%%%%%%% END
做完结果发现小弟频域分析好像有错误
问题与疑点:
(1) 小弟的波形是60Hz与180Hz正旋波做交叠,所以FFT分析结果应该是60与180有值其它很小值,但是结果却很奇特
因为小弟第一次用请帮解答
(2) 小弟基频是60Hz做傅立叶大小后我发现程式结果看Workspace 中的结果(EVudfft)会不知哪个是基频
是要验证画图才知道,但是小弟是要利用数据做计算请问该如何取0Hz 60Hz 60Hz的倍数?
(3) 能否告知該如何使用FFT分析?
(4) 或是有更好的程序能更精确快速解决像是小弟此范例 |
|