我求出的相频图是与理论图像相差太大。  这是怎么回事呀?  请高手指点一下

你的t是个什么呀,要特别注意点乘的用法.

有没有考虑相位解缠(unwrap)啊

angle求的相位在[-pi,pi]之间.

把帮助给你看看，你应该就清楚了：
P = angle(Z) returns the phase angles, in radians, for each element of complex array Z. The angles lie between . For complex Z, the magnitude R and phase angle theta are given by R = abs(Z)
theta = angle(Z)
and the statement Z = R.*exp(i*theta)
converts back to the original complex Z.