时频域极化参数———极化椭圆率matlab绘图

以下是一个简单的matlab代码，用于绘制极化椭圆率参数的时频图：

% 构造信号
fs = 100;                       % 采样率
t = 0:1/fs:1-1/fs;              % 时间向量
f1 = 10;                        % 信号频率
f2 = 30;
s1 = sin(2*pi*f1*t);
s2 = sin(2*pi*f2*t);
s = s1 + s2;

% 计算极化椭圆率参数
A = hilbert(s);                 % 解析信号
theta = angle(A);               % 相位信息
E = exp(1i*2*theta);            % 极化矢量
ellip = abs(sum(E))/length(E);  % 极化椭圆率

% 绘制时频图
win_len = 128;                  % 窗口长度
noverlap = 64;                  % 重叠长度
nfft = 256;                     % FFT点数
[S,F,T,P] = spectrogram(s,win_len,noverlap,nfft,fs);
ellip_tf = zeros(size(S));
for i = 1:size(S,2)
    A_tf = hilbert(S(:,i));
    theta_tf = angle(A_tf);
    E_tf = exp(1i*2*theta_tf);
    ellip_tf(i) = abs(sum(E_tf))/length(E_tf);
end
figure;
surf(T,F,20*log10(abs(S)),'EdgeColor','none');
axis tight;
view(0,90);
xlabel('时间(s)');
ylabel('频率(Hz)');
title('信号的时频图');
colorbar;
figure;
surf(T,F,ellip_tf,'EdgeColor','none');
axis tight;
view(0,90);
xlabel('时间(s)');
ylabel('频率(Hz)');
title('信号的极化椭圆率的时频图');
colorbar;

这个代码首先构造了一个包含两个正弦波的合成信号。然后，它计算了该信号的极化椭圆率参数，并将其用于绘制时频图。最后，它绘制了信号的极化椭圆率的时频图。