matlab 生成 尖峰信号,MATLAB的FFT波形图只含有一个尖峰

我有一些初级的,基本的物理网络数据,我试图应用FFT,但我有点困惑的结果,我认为他们是正确的。我主要是用这个代码

example

fftshift documentation page

有一些调整,但我不太确定我的代码或结果有什么问题。我已经附上了我的代码和我得到的结果的快照。

load aami3am.mat

Fs = 720; % Sampling frequency

T = 1/Fs; % Sample time

L = 60000; % Length of signal

t = (0:L-1)*T; % Time vector

plot(t(1:43081),val(1:43081))

title('aami4b_h Signal')

xlabel('Seconds')

ylabel('ECG Amplitude')

NFFT = 2^nextpow2(L);

val = detrend(val);

Y = fft(val,NFFT)/L;

f = Fs/2*linspace(0,1,NFFT/2+1);

plot(f,2*abs(Y(1:NFFT/2+1)))

title('FFT')

yY = fftshift(Y);

fshift = (-L/2:L/2-1)*(Fs/L);

powershift = abs(yY).^2/L;

plot(fshift(1:L),powershift(1:L))

title('FFT Shifted')

FbE7M.png

lIQnW.png

UWjOy.png

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值