how to make fft2 surf plot with z axis in dB?

7 views (last 30 days)
i tried to normalize this matrix with the maximum
and: p=10*log10(p);
but i get mess
how can i avoid that?

Accepted Answer

Star Strider
Star Strider on 20 Nov 2014
Edited: Star Strider on 20 Nov 2014
Without access to your data, it’s difficult to say. Two possibilities:
First:
p = 10*log10(abs(p));
Second:
p = 10*log10(abs(p)+1E-2);
  1 Comment
Star Strider
Star Strider on 20 Nov 2014
Interesting, to say the least!
This is what I get with your data and this code:
D = matfile('Kobi_pmat.mat');
P = D.P;
LP = 10*log10(abs(P));
figure(1)
mesh(LP)
grid on
axis([600 800 600 800 -30 0])
xlabel('x')
ylabel('y')
zlabel('F(x,y) (dBi)')
title('Fourier Transform of Rectangular Reflector')
producing this plot:

Sign in to comment.

More Answers (2)

Youssef  Khmou
Youssef Khmou on 20 Nov 2014
This problem does not occur with the following test :
H=exp(-randn(10).^2);
F=abs(fftshift(fft2(H,200,200)));
F=F/max(F(:));
surf(20*log10(F));
try abs() as mentioned earlier .

Kobi
Kobi on 20 Nov 2014
Edited: Star Strider on 20 Nov 2014
here is my data file:
you can load it like this:
data=load('pmat.mat');
P=data.P;

Categories

Find more on Elementary Math in Help Center and File Exchange

Tags

Products

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by