%觀察二維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