f6_10.m

来自「《Matlab时频分析技术及其应用》随书源代码」· M 代码 · 共 22 行

%产生单分量信号
sig1=atoms(128,[64,0.25,20,1]);
%计算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig1);
%计算Renyi信息
disp('单分量信号的Renyi信息度量为')
R1=renyi(tfr,t,f)
%产生两分量信号
sig2=atoms(128,[32,0.25,20,1;96,0.25,20,1]);
%计算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig2);
%计算Renyi信息
disp('两分量信号的Renyi信息度量为');
R2=renyi(tfr,t,f)
%产生四分量信号
sig3=atoms(128,[32,0.15,20,1;96,0.15,20,1;32,0.35,20,1;96,0.35,20,1]);
%计算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig3);
%计算Renyi信息
disp('四分量信号的Renyi信息度量为');
R3=renyi(tfr,t,f)