zoukankan      html  css  js  c++  java
  • S变换画图

    The S -transform with windows of arbitrary
    and varying shape

    想画如上的三张图,可是问题

    (1)滤波窗口的平滑实现

    The filter (applied to the full 200£400 S-transform matrix)
    is constructed in Matlab by constructing a two-dimensional
    “boxcar”, using
    filter = zeros(200,400);
    filter(69:128,180:240)=1;
    The filter edges are then smoothed by setting the value of
    each pixel equal to the average of its four adjacent pixels (this
    procedure is repeated 100 times).

    (2)S变换中的频率与短时傅里叶变换频率怎么对应

    我画出来的频率不太对。。边缘也还没处理。。

    clear all
    % x1=readEIS('35006_6_2121_秒.txt');
    % x=x1(:,2);
    x1=load('qzh_3.txt');
    % x=qzh_3;
    
    points=length(x1);
    figure
    subplot(211),plot(0:0.02:0.02*(points-1),x1);title('原始信号');grid on; 
    
    
     [st_matrix,st_times,st_frequencies] = st(x1) ;
    %  imagesc(st_times,st_frequencies,abs(st_matrix))
    % imcontour(abs(st_matrix))
    ST=abs(st_matrix);
    filter = zeros(1501,3000);
    filter(100:800,900:1400)=1;
    STF=ST.*filter;
    STF(STF==0) = NaN;
    subplot(212),h=imagesc(STF);%其中C为包含有nan的数据
    set(h,'alphadata',~isnan(STF))
    ylim([0 1000])
    % set(gca, 'ydir', 'reverse');
    colorbar
    colormapeditor
    

      

     

  • 相关阅读:
    SpringBoot+SpringCloud
    bootstrap-thymeleaf-分页
    排序-Java
    native2ascii运用
    标准W3C盒子模型和IE盒子模型
    在既定状态下截图
    java.util.zip.ZipException: error in opening zip file
    安装 haproxy
    mysql集群
    最简redis集群配置
  • 原文地址:https://www.cnblogs.com/gisalameda/p/9090004.html
Copyright © 2011-2022 走看看