声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 813|回复: 0

[综合讨论] 关于fourier的疑问

[复制链接]
发表于 2007-5-14 10:49 | 显示全部楼层 |阅读模式

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

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

x
简单介绍哈算法,主要还是编程的问题,望高人解答。
w(x,y)是初始相位函数,Δw(x,y)是相位差函数,算法的目的是想通过相位差函数,恢复相位函数。
波前相位重建的fourier方法,Δw(x,y)=w(x+s/2,y)-w(x-s/2,y)
由fourier变化可得到:         F[w(x,y)](v)=T(v)*F[Δw(x,y)](v)     v是空间域x对应的频域变量。
                                       T(v)=1/2*i*sin(pi*s*v)                    T(v)是剪切转换函数
            
以下是程序:  
clear all;
syms x y z v;
s=0.4;%剪切量
a=(x^2+y^2-2/3)*y+4/3;%原始相位
b=(((x+s/2)^2+y^2-2/3)*y+4/3)-(((x-s/2)^2+y^2-2/3)*y+4/3);%相位差函数
b
B=fourier(b,x,v);%对相位差函数进行fourier变换
C=1/(2*i*sin(pi*s*v));%剪切转换函数
D=B*C;
d=ifourier(D,v,x);%反变换
有算法知道,d应该和a是一个函数,如果通过取值来验证的话,相差也不应该太大,我仿真出来, d是一个只关于y的函数,不知道哪里有问题,望达人解答。
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-14 21:37 , Processed in 0.071636 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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