声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1489|回复: 4

[图像处理] 这个程序如何理解

[复制链接]
发表于 2012-11-9 17:07 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
clc
clear

h(1:128)=cos((0:127)*2*pi*6/256);
h(31:71)=h(31:71)+0.5*cos((-10:30)*2*pi*52/256);
h(129:256)=cos((0:127)*2*pi*25/256);
N = 256;
H = [fft(h) fft(h) fft(h)];
al = [0:(N/2)-1, -(N/2):-1];
for n=1:N
W=exp(-2*pi^2*al.^2./((n-(N/2)-1).^2));
S(n,:) = ifft(H(n+(N/2):n+(N/2)+N-1) .* W);
end

S(N/2+1,:) = ifft(H(N+1:2*N) .* [1 zeros(1,N-1)]);
S_TRANS = abs(S);
figure;
pcolor(S_TRANS)

程序是时间序列的S变换,其公式为
S(t,f)=∫_(-∞)^∞▒〖H(a+f)exp⁡((-2π^2 a^2)/f^2 )  exp⁡(+2πiat)da〗
上面程序中
H = [fft(h) fft(h) fft(h)]; 这句应该怎么理解,为什么要延拓信号的傅里叶变换?
al = [0:(N/2)-1, -(N/2):-1];这句为什么要设置成这种形式?
回复
分享到:

使用道具 举报

 楼主| 发表于 2012-11-12 20:45 | 显示全部楼层
真的木有来帮忙看看啊?
发表于 2012-11-12 22:20 | 显示全部楼层
个人水平有限,感觉这不是编程问题!?
或许应移至信号版块!
发表于 2012-11-13 08:33 | 显示全部楼层
确实是属于信号处理问题,去信号版问问应该能有结果
发表于 2012-11-14 10:02 | 显示全部楼层
楼主应该先说明这个程序是干嘛的?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-24 21:52 , Processed in 0.052900 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表