How do I invert nested functions in MATLAB?

1 vue (au cours des 30 derniers jours)
Ak K
Ak K le 14 Août 2019
Modifié(e) : Walter Roberson le 15 Août 2019
If I have a variable (let's call it a) and then I apply the following multiple functions to transform it into a second variable, (b), how would I reverse it?
b = fftshift(fftn(fftshift(a)));
I tried doing the following:
c = ifftshift(ifftn(ifftshift(b)));
If this had worked, c would have been equal to a (the initial, non-transformed data set). Why doesn't it? Is there a better, more universal way of reversing the function initially applied to A? I tried finverse, but it didn't work in the way I used it because fftshift works on doubles, and I was using a sym-type variable when I tried getting the inverse.
Any help is appreciated

Réponse acceptée

KALYAN ACHARJYA
KALYAN ACHARJYA le 14 Août 2019
Modifié(e) : KALYAN ACHARJYA le 14 Août 2019
"If I have a variable (let's call it a) and then I apply the following multiple functions to transform it into a second variable, (b), how would I reverse it?"
a=magic(3)
b = fftshift(fftn(fftshift(a)));
c = ifftshift(ifftn(ifftshift(b)))
Result:
a =
8 1 6
3 5 7
4 9 2
c =
8.0000 1.0000 6.0000
3.0000 5.0000 7.0000
4.0000 9.0000 2.0000
  3 commentaires
Walter Roberson
Walter Roberson le 14 Août 2019
Modifié(e) : Walter Roberson le 15 Août 2019
Floating point roundoff. I suggest that you subtract the two versions and look at the difference.
fft involves the use of transcendental trig functions, so in order for the inverse to be always bit for bit exact, you would need to use indefinite precision.
You can use counting arguments to prove that if you use finite floating point that asin(sin(x)) cannot always exactly equal x even within the primary period.
Ak K
Ak K le 15 Août 2019
Got it, thanks; I tried that, and the difference was indeed extremely small.

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Fourier Analysis and Filtering dans Help Center et File Exchange

Produits


Version

R2018a

Community Treasure Hunt

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

Start Hunting!

Translated by