首页->【MATLAB Usage 】

2 1

function [odinmag,odoutmag,fn]=uiFFT(fs,N,dia,dib)
%UIFFT 此处显示有关此函数的摘要
% 此处显示详细说明
%FFT  module
% fs  Sample freq
% N  FFT分析数据样本数

%加入输入的信号是非平稳信号,对非平稳信号进行分段处理

% dia,dib : Input data 1xn Matrix


%去除直流分量 
dia=dia-mean(dia);
dib=dib-mean(dib);


% add the windows function
wn=hann(N);   % N点汉宁窗
di=di(1:N).*wn';   % 加汉宁窗
do=do(1:N).*wn';

din=fft(di,N);
dout=fft(do,N);

dinmag=20*log(abs(din)/max(abs(din)));
doutmag=20*log(abs(dout)/max(abs(dout)));

fn=((0:N/2)*fs/N);
odinmag=dinmag(1:N/2+1);
odoutmag=doutmag(1:N/2+1);
% 将输出的频谱值使用plot进行显示即可。
end

楼主可见

  1. shwnyoo 1#

    仅供参考,