Box Funkton mit fft < Matlab < Mathe-Software < Mathe < Vorhilfe
|
Wir sollen die Box-Funkton in Matlab mittels der Funktion fft darstellen. Allso mittels der fast fourier transformation. Nur leider check ich den Befehl nicht wriklich und hatte auch keine eiführende VO zu Fourierreihen - steh also ziemlich im Regen.
Box Funkton
[mm] f(x)=\begin{cases} 1, & |x| \mbox{für } \le \bruch{a}{2} \\ 0, & \mbox{für } n \mbox{ sonst} \end{cases} [/mm]
Was ich bis jetzt hab (leider nicht ganz korrekt):
---------------begin code ---------------
clear
hold on
L = 20; % Intervalllänge
a = 10; % Intervall der Box
N = 3000; % Teilintervalle (Abtastrate)
n = 1:N; % Intervallvektor ??
X = zeros(1,N);
for m = 1:N
X(m) = ((m/N)-1/2)*L;
end
[mm] f_j [/mm] = (X>=-a/2) & (X<=a/2);
k = (2.*pi.*n)./L;
% analytsiche Lösung
f = [];
f_analytical = [];
for x=X
f_analytical = [f_analytical sum(exp(i.*k.*x).*(2./(k.*L)).*sin((k.*a)/2))];
end
plot(X,f_analytical,'g')
fk = [mm] fft(f_j)/N [/mm] .* (exp(-i*k.*X));
% Berechnung mittels der Koeffizienten von [mm] fft(f_j)
[/mm]
for x=X
f = [f sum(exp(i.*k.*x).*fk)];
end
[mm] plot(X,f_j,'r')
[/mm]
plot(X,f,'b-')
hold off
---------------end code ---------------
Wäre auch nett wenn jemand eine Seite (aus Wikipedia) kennt wo man
sich die Basics gut ansehen kann - Danke.
Lg andi
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 14:20 Sa 17.10.2009 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
