求二维傅立叶变换频谱时原函数的构造问题

timedcy

比如现在要生成 一个  双个圆孔的、三角形孔的 频谱  怎么弄

还有跟不会的网格的 频谱  等   怎样生成这些复杂的原函数呢？？

help里有几个二维的频谱图但是没有code

[ 本帖最后由 eight 于 2007-6-9 15:46 编辑 ]

happy

不是很明白楼主的意思，楼主想知道fft2,fftn等函数的使用办法？

[ 本帖最后由 eight 于 2007-6-9 16:07 编辑 ]

timedcy

原帖由 happy 于 2007-6-9 16:00 发表

                               
登录/注册后可看大图

不是很明白楼主的意思，楼主想知道fft2,fftn等函数的使用办法？

我想画出一个规则图形的傅立叶变换频谱

wavelet能行吗 找了半天不会做

如果用fft2做 比如一个三角形的孔的频谱怎么弄出来呢？

谢谢！:@o

timedcy

难一点的是 画出一个“五”字的频谱图等的

happy

不知道是否是这个意思

1. t = (0:1/3:1)'*2*pi;
   
2. x = sin(t);
   
3. y = cos(t);
   
4. figure(1);
   
5. fill(x,y,'k');
   
6. axis([-10,10,-10,10]);
   
7. axis off;
   
8. box off;
   
9. print -djpeg -r85 sanjiao;

复制代码

1. img = imread('sanjiao.jpg');
   
2. figure(2);
   
3. imshow(img);

复制代码

1. img1=double(img)/255;
   
2. img_mean=mean2(img1);
   
3. centered_img=img1-img_mean;
   
4. img1_fft=fft2(centered_img);
   
5. img1_fft=fftshift(img1_fft);
   
6. figure(3);
   
7. imshow(img1_fft);

复制代码

timedcy

怎么出来是这样的谱啊    连六角形都看不出啊