信号频谱质心matlab,时频分析包tfrstft,然后求频谱质心 | 您所在的位置:网站首页 › 曲线的质心公式推导 › 信号频谱质心matlab,时频分析包tfrstft,然后求频谱质心 |
w1=0.23; w2=1; f0=3.4*10^11; fs=100; c=3*10^8; r=0.00025; r1=0.005; r2=0.0008; t=0:0.01:30; a=0.25-r*w2; t3=w2.*t; t4=floor(t3); y=abs(t3-0.5-t4); aaa=1/(1-2*a).*(y-a+abs(y-a)); s=r1.*sin(2.*pi.*w1.*t)+r2.*aaa; sr=exp(1j*2*pi*f0.*(t-2*s/c)); L=length(sr); %[tfr,t,f] = tfrstft(sr.',1:L,L,hamming(63),0); [tfr,t,f] = tfrstft(sr.',1:L,L-2000,hamming(63),0); % 此处修改了频率的数量,通过tfr的矩阵大小可以分清楚哪个是时间,哪个是频率 figure(); %contour(t/fs,f(1:length(f))*fs,abs(tfr(1:length(f),:))); mesh(t/fs,f(1:length(f))*fs,abs(tfr(1:length(f),:))); xlabel('time/s') ylabel('frequency/Hz') view([0,90]); % 时间和瞬时频率 xlim([0,30]) tfr=abs(tf |
CopyRight 2018-2019 实验室设备网 版权所有 |