二维广义S变换

包含二维广义S变换程序和测试程序，用于对图像进行时频分析。matlab版本R2011b,操作系统为32位win7。在2G内存下最大只能对84*84的图像进行处理，更大的图像会出现out of memory的问题。目前网上还没有这个程序，包括本站，pudn，ilovematlab都没有。把程序上传到这三个地方，望大家一起来改进它。

%观察二维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]=tdst4（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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        4974  2013-03-19 09:11  tdst.m
     文件        1126  2013-01-03 15:37  analysis_at_certain_freq.m