[Sidney1]

Koledzy i Koleżanki z #matlab trochę może #matematyka #fizyka i #elektronika

Robię funkcję FFT:

t = 0:.00001:0.1-0.00001;
Fs = 1e3;
ta1=0.00035;
ta2=0.00021;
nb=1.8;
kb=exp((-ta1/ta2)*((-nb*ta2/ta1)^(1/nb)));
u=(1000./kb).*(((t./ta1).^nb)./(1+(t./ta1).^nb)).*exp((-t./ta2));

udft = fft(u);
udft = udft(1:length(u)/2+1);

freq = 0:Fs/length(u):Fs/2;
plot(freq,(abs(udft)),'r');

I to wychodzi bardzo ładnie zgadza się z założeniami moimi.

Ale jak zrobić, żeby z tego fft odzyskać sygnał w dziedzinie czasu?
Jak zrobię ifft(udft) wychodzi mi zupełnie inny sygnał.

Jakieś pomysły?

[Sidney1]

@stefan_pmp: a bo mi wykop gwiazdki pousuwał ( ͡° ʖ̯ ͡°) Chodzi Ci o to, że rzadko próbkowany?

[Sidney1]

@stefan_pmp: ale ona chyba wychodzi symetrycznie? A co z odwracaniem?

[Sidney1]

@stefan_pmp: trochę nie rozumiem co się tutaj odwaliło, definiujesz tutaj oddzielne funkcje czy to ma pójść w stary skrypt?

[Sidney1]

@stefan_pmp: rozumiem, że n- to czas, N- to liczba próbek, a co to jest kn?

[zboinek]

@Sidney1 fx po dx

[Kolczaneiro]

@Sidney1 możesz też idft na piechotę napisać w sumie xd