FFT.zip_40_Hz_NFFTmatlab资源-CSDN文库

共2个文件

m：2个

版权申诉

122 浏览量 2022-07-14 02:47:04 上传评论收藏 1KB ZIP 举报

资源详情

资源评论

资源推荐

收起资源包目录

FFT.zip （2个子文件）

test1.m 1KB

test2.m 1KB

clear; clc; clf; fs=100; %采样频率 Ndata=32; %数据长度 N=32; %FFT采样的数据长度 n=0:Ndata-1;t=n/fs; %数据对应的时间序列 x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %时间域信号 y=fft(x,N); %信号的Fourier变换 mag=abs(y); %求取振幅 f=(0:N-1)*fs/N; %真实频率 subplot(2,2,1),plot(f(1:N/2),mag(1:N/2)*2/N); %绘出Nyquist频率之前的振幅 xlabel('频率/Hz');ylabel('振幅'); title('Ndata=32 Nfft=32');grid on; Ndata=32; %数据个数 N=128; %FFT采样的数据长度 n=0:Ndata-1;t=n/fs; %时间序列 x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); y=fft(x,N); mag=abs(y); f=(0:N-1)*fs/N; %真实频率 subplot(2,2,2),plot(f(1:N/2),mag(1:N/2)*2/N); %绘出Nyquist频率之前的振幅 xlabel('频率/Hz');ylabel('振幅'); title('Ndata=32 Nfft=128');grid on; Ndata=136; %数据个数 N=128; %FFT采样的数据长度 n=0:Ndata-1;t=n/fs; %时间序列 x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); y=fft(x,N); mag=abs(y); f=(0:N-1)*fs/N; %真实频率 subplot(2,2,3),plot(f(1:N/2),mag(1:N/2)*2/N); %绘出Nyquist频率之前的振幅 xlabel('频率/Hz');ylabel('振幅'); title('Ndata=136 Nfft=128');grid on; Ndata=136; %数据个数 N=512; %FFT采样的数据长度 n=0:Ndata-1;t=n/fs; %时间序列 x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); y=fft(x,N); mag=abs(y); f=(0:N-1)*fs/N; %真实频率 subplot(2,2,4),plot(f(1:N/2),mag(1:N/2)*2/N); %绘出Nyquist频率之前的振幅 xlabel('频率/Hz');ylabel('振幅'); title('Ndata=136 Nfft=512');grid on;