- 积分
 - 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)));  |   
 
 
 
 |