标签:FFT fft 小波换 加窗 time dt data out
记录一下matlab中时历曲线分析的相关函数
FFT
function out = myFFT(time,data)
% 功能:FFT
% 输入:time:时间
% data:数据
% 输出:out.Fre:频率
% out.Ampti:幅值
% out.Phase:相位
N=length(data);
fft_Data=fft(data);
Ampti = abs(fft_Data) / (N)*2;
Phase = atan2(imag(fft_Data),real(fft_Data));
dt=time(2)-time(1);
Fre=(0:1:(length(time)-1))/(N*dt)*2*pi;
out.Fre=Fre;
out.Ampti=Ampti;
out.Phase=Phase;
end
加窗FFT
function [omega,res] = getWelch(data,dt,N)
% 功能:FFT加窗处理
% 输入:data 数据段
% dt 时间间隔
% N 每段处理数据长度N
fs=1/dt;
% 每段数据之间的重叠区域宽度为:N*1/2
[p,f] = pwelch(data,hamming(N),N/2,N,fs);
omega = f*2*pi;
res = p/2/pi;
% res = sqrt(res); %转化为幅值
end
小波变换
% data时历
% fs采样频率(Hz)
% wt:为f[Hz]与时间t[s]相关(二维矩阵)
[wt,f]=cwt(data,fs);
测试案例
待添加。。。
标签:FFT,fft,小波换,加窗,time,dt,data,out 来源: https://www.cnblogs.com/chetwin/p/16415974.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。