- 积分
- 320
- 回帖
- 0
- 西莫币
-
- 贡献
-
- 威望
-
- 存款
-
- 阅读权限
- 20
- 最后登录
- 1970-1-1
该用户从未签到
|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?立即注册
×
参考help中关于fft的资料,编写m文件对电流谐波进行分析,数据导入为外部txt文本倒入的,但m文件得离散傅里叶变换的结果,画图,横坐标为频率,纵坐标为谐波幅值,请问如何才能将横坐标设为真实频率?代码如下:
data=importdata('G:\current.txt');
[N,t,A,B,C] = textread('G:\current.txt','%f %f %f %f %f');
Fs = 805; % 采样频率
T = 1/Fs; % 采样时间
L = 805; % 信号长度
plot(t,A);
xlabel('time (ms)');
ylabel('current');
NFFT = 2^nextpow2(L); % 返回最小的2的指数倍数值(大于等于数据长度L)
Y = fft(A,NFFT)/L;
f = Fs*linspace(0,1,NFFT/2); %横坐标频率
plot(f,2*abs(Y(1:NFFT/2))); |
|