Implementacija kompleksnih števil za FT

Morda?

https://www.nayuki.io/page/how-to-imple...

marjan_h je 19. okt 2016 ob 22:34 izjavil:

Kako implementirati računanje z kompleksnimi števili, v nekem standardnem jeziku C/C++/Java?

Namreč želim napisati funkcijo, ki bi izračunala FT (Fourierjevo transformaicjo) v diskretni obliki. Vendar kaj naj naredim ko množim z e^(i*omega*t)?

Pa brez uporabe knjižnic.

Hvala.

Isuse, pac nek struct z dvema doubloma nad katerim definiras operacije.

Najbrž ma za vajo na kakem faksu, da pokaže razumevanje delovanja FT in ne sme knjižnice uporabit.
Da zračunaš diskretno FT, sploh ne rabiš nobenih operacij nad imaginarnimi števili.
Za vsako frekvenco v zanki primerjaš vhodni signal z sinusoido in s kosinusoido trenutne frekvence. V matematičnih izpeljavah se potem ta koeficienta povežeta z imaginarnimi števili, zgolj za to da lahko računaš z njimi in da si lahko predstaviš vrednosti kot točke v 2D grafu. V diskretnem niti ne rabiš računat, lahko čist ločeno zračunaš realne in imaginarne koeficiente, tisti i te naj ne zanima. Na koncu samo če te zanima s kako amplitudo in s kakšno fazo je prisotna določena frekvenca v signalu, moraš pa iz obeh koeficientov to zračunat - dejansko gre za pretvorbo iz kartezičnega v polarni koordinatni sistem - razdalja ti predstavlja amplitudo, kot (v radianih seveda) pa fazo.

Lahko delas loceno sin in cos (real in imag), je pa veliko bolj komot delati s kompleksnimi spremenljivkami.
GCC (in C99) jih ima vgrajene

#include [complex.h]
.....
.....
complex float a,b,c;
a=1.2+3.4I;
b=5.6+7.8I;
c=a*b;
printf("\n\nc = %f %fi\n\n",crealf(c),cimagf(c));

itd....

imas tudi cexpf(), cpowf(), csinf().....

pri include morata namesto [] biti < in > (slotech se upre)

https://www.nayuki.io/page/how-to-imple...