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

triangle=TRUE;
rectangle=FALSE;

image=50*ones(50,50);
if triangle==0&&rectangle==1
    image(20:30,20:30)=100;
elseif triangle==1&&rectangle==0
Rs          = 50;
image       = 50*ones(Rs,Rs);
for i = 1:50
    for j = 1:50
        if (i-round(Rs/2))^2+(j-round(Rs/2))^2 <=64
            image(i,j)= 100;
        end
    end
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(8,8);

for j=9:6:51
    figure;
    for k=9:6:51
        subplot(3,3,(k-3)/6);
        tdst_freq_all_over{(j-3)/6,(k-3)/6}=abs(matrix_to_be_printed{j,k});
        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(1,15);
figure;
for l=9:3:51
    subplot(4,4,(l-6)/3);
    tdst_freq_along{1,(l-6)/3}=abs(matrix_to_be_printed{9,l});
    contourf(tdst_freq_along{1,(l-6)/3});
    title(['kx=8 ky=',num2str(l-1)]);
end

联系：highspeedlogic

QQ ：1224848052

微信：HuangL1121

邮箱：1224848052@qq.com

网站：http://www.mat7lab.com/

网站：http://www.hslogic.com/