haodawei123的博客

代码实现了快速傅里叶变换及其逆变换的计算，幅值、初始相位的计算import numpy as npimport matplotlib.pyplot as pltFs=500#采样频率ts=1/Fs#采样时间间隔L=500#信号长度t=np.linspace(0,L-1,L)*tsx=2+0.7*np.cos(2*np.pi*50*t+np.pi/2)+2*np.cos(2*np.pi*200*t+np.pi/6)y=x#+2*np.random.rand()fft_y=np.fft.ff

clcFs=1000;%采样频率T=1/Fs;%一个信号的采样时间,1/1000=0.001sL=1000;%信号长度t=(0:L-1)T;%时间向量：生成0-999之间的数乘以0.001就是0-1之间生成0-0.999x=0.7sin(2pi50t)+sin(2pi120t);%产生一个数据，作为我们的采样数据y=x+2randn(size(t));%将上面的数据加上噪声,size(t)返回向量t的行列数%y=x;%去掉上面的噪声，加上噪声会幅值会变化plot(Fst(1:50),y(1: