查看: 88|回复: 0|关注: 0

[未答复] 归一化能量谱密度怎么求

[复制链接]

新手

9 麦片

财富积分


050


9

主题

30

帖子

0

最佳答案
发表于 5 天前 | 显示全部楼层 |阅读模式
w1=0.23;
w2=1;
f0=3.4*10^11;
fs=100;
c=3*10^8;
r=0.0025;
r1=0.005;
r2=0.0008;
t=0:0.01:30;
a=0.25-r*w2;
t3=w2.*t;
t4=floor(t3);
aaa=1/(1-2*a).*(abs(t3-0.5-t4)-a+abs(abs(t3-0.5-t4)-a));
s=r1.*sin(2.*pi.*w1.*t)+r2.*aaa;
sr=exp(j*2*pi*f0.*(t-2*s/c));
L=length(sr);
[tfr,t,f] = tfrspwv(sr.',1:L,L,hamming(63));  
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),:)));
view([0,90]); %时间和瞬时频率
xlim([0,30])
centroid=zeros(1,L);
w=f(1:length(f))*fs;
w=w.';
spe_tfr=abs(tfr);
for k=1:L
    centroid(k)=w*spe_tfr(:,k)/sum(spe_tfr(:,k));
end
figure(2)
plot(t/fs,centroid,'r-.')
xlabel('Time/s')
ylabel('Fre/Hz')
legend('Centroid')
grid on
xlim([0,30])
cen=resample(centroid,5,fs)
LL=length(cen);
[tfr1,t1,f1] = tfrspwv(cen.',1:LL,LL,hamming(63));   
figure(3);
%contour(t/fs,f(1:length(f))*fs,abs(tfr(1:length(f),:)));
mesh(t1/5,f1(1:length(f1))*5,abs(tfr1(1:length(f1),:)));
xlabel('time/s')
ylabel('frequency/Hz')
view([0,90]); % 时间和瞬时频率
xlim([0,30])
请问怎么理解归一化能量谱密度,如何求归一化能量谱密度,以及对应的频率和归一化的高能量密度均值?谢谢

1.png
2.png
3.png
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /2 下一条

快速回复 返回顶部 返回列表