信号频谱质心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