|
1.相同类似帖没必要发新帖
2.减小采样频率并非单纯将Fs降低而已!
3.for i=1:N; t(i)=wave(i,2); end
这种写法不好, 直接t=wave(:,2);即可- clc; clear;
- wave=load('30002.txt');
- t=wave(:,2); t=detrend(t); N=length(t);
- Fs=10000; Ayy=abs(fft(t))*2/N; Ayy(1)=Ayy(1)/2;
- F=([1:N]-1)*Fs/N; %换算成实际的频率值
- figure; subplot(211); plot(t); title('原始信号');
- subplot(212); plot(F(1:N/2),Ayy(1:N/2)); title('幅度-频率曲线图'); %显示换算后的FFT模值结果
- t=t(1:2:end); t=detrend(t); N=length(t);
- Fs=5000; Ayy=abs(fft(t))*2/N; Ayy(1)=Ayy(1)/2;
- F=([1:N]-1)*Fs/N; %换算成实际的频率值
- figure; subplot(211); plot(t); title('原始信号');
- subplot(212); plot(F(1:N/2),Ayy(1:N/2)); title('幅度-频率曲线图'); %显示换算后的FFT模值结果
复制代码 |
|