clc
close(gcf)
clear all;
disp('Wykreślanie wsp. szeregu Fouriera dowolnej funkci');
funs=input('Podaj wzór funkcji f(t)=','s');
tk=input('Podaj czas końcowy tk[s]=');
Lp=1024;
Tobs=tk-0;
df=1/Tobs;
fmax=df*Lp/2;
t=[0:Tobs/Lp:Tobs-Tobs/Lp];
fun=eval(funs);
subplot(2,2,1)
plot(t,fun);
title('Przebieg funkcji f(t)');
xlabel('Czas t[s]');
grid on
FUN=fft(fun);
FUNp=FUN(1:Lp/2)/Lp*2;
FUNp(1)=FUNp(1)/2;
m=abs(FUNp);
subplot(2,2,2)
f=[0:df:fmax-df];
axis([0 fmax 0 max(m)]);
hold on
grid on
bar(f,m);
title('Widmo amplitudowe funkcji f(t)-współczynniki abs(cn)')
xlabel('Częstotliwość f[Hz]');
{"html5":"htmlmixed","css":"css","javascript":"javascript","php":"php","python":"python","ruby":"ruby","lua":"text\/x-lua","bash":"text\/x-sh","go":"go","c":"text\/x-csrc","cpp":"text\/x-c++src","diff":"diff","latex":"stex","sql":"sql","xml":"xml","apl":"apl","asterisk":"asterisk","c_loadrunner":"text\/x-csrc","c_mac":"text\/x-csrc","coffeescript":"text\/x-coffeescript","csharp":"text\/x-csharp","d":"d","ecmascript":"javascript","erlang":"erlang","groovy":"text\/x-groovy","haskell":"text\/x-haskell","haxe":"text\/x-haxe","html4strict":"htmlmixed","java":"text\/x-java","java5":"text\/x-java","jquery":"javascript","mirc":"mirc","mysql":"sql","ocaml":"text\/x-ocaml","pascal":"text\/x-pascal","perl":"perl","perl6":"perl","plsql":"sql","properties":"text\/x-properties","q":"text\/x-q","scala":"scala","scheme":"text\/x-scheme","tcl":"text\/x-tcl","vb":"text\/x-vb","verilog":"text\/x-verilog","yaml":"text\/x-yaml","z80":"text\/x-z80"}