广义s变换matlab程序一个以及带注释的matlab一个

广义s变换matlab程序一个以及带注释的matlab例子一个，可以直接运行

%观察二维S变换对特定图形的分析能力
%可选图形：正方形，三角形
%分析结果，给出图像在整个频率平面上均匀分散的各点处的变换结果，同时给出在某一kx处许多ky值处的变换结果
clc;
clear;
TRUE=1;
FALSE=0;

triangle=TRUE;
rectangle=FALSE;

image=50*ones（5050）;
if triangle==0&&rectangle==1
    image（20:3020:30）=100;
elseif triangle==1&&rectangle==0
    for i=20:30
        image（i20:i）=1;
    end
else disp（sprintf（‘Please choose the image either triangle or rectangle‘））;
end
[image_tdst image_tdst_freq]=tdst（image）;
matrix_to_be_printed=image_tdst_freq;

tdst_freq_all_over=cell（88）;

for j=9:6:51
    figure;
    for k=9:6:51
        subplot（33（k-3）/6）;
        tdst_freq_all_over{（j-3）/6（k-3）/6}=abs（matrix_to_be_printed{jk}）;
        contourf（tdst_freq_all_over{（j-3）/6（k-3）/6}）;
        title（[‘kx=‘num2str（j-1）‘ ‘‘ky=‘num2str（k-1）]）;
    end
end

tdst_freq_along=cell（115）;
figure;
for l=9:3:51
    subplot（44（l-6）/3）;
    tdst_freq_along{1（l-6）/3}=abs（matrix_to_be_printed{9l}）;
    contourf（tdst_freq_along{1（l-6）/3}）;
    title（[‘kx=8 ky=‘num2str（l-1）]）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        4718  2019-06-20 14:43  广义s变换.m
     文件        1125  2019-06-20 14:34  例子.m