tallenna.m

% Määritellään anonyymi funktio
f = @(x) sin(x)./x .* exp(-(0.05*x).^2);

% Esitetään funktio graafisesti
hold off
x = linspace(-40,40,200)'; % Huomaa transpoosi: x on pystyvektori
y = f(x);
plot(x,y)

% Etsitään kolme ensimmäistä lokaalia minimiä. Alkuarvaukset 5, 10 ja
% 15 on valittu graafisesti
xmin = [];
xmin(1) = fminunc( f, 5 );
xmin(2) = fminunc( f, 10 );
xmin(3) = fminunc( f, 15 );

xmin              % Tulosta löydetyt x-koordinaatit
ymin = f(xmin)    % Tallenna myös vastaavat y-koordinaatit

% Esitä minimipisteet samassa grafiikassa
% Zoomilla voidaan huomata, että minimipisteet on ratkaistu tarkemmin
% kuin millä funktio on piirretty
hold on
plot( xmin, ymin, 'ro' )

% Tallenna muuttuja (vektori) y nimettyyn tiedostoon
save sinc_pulssi1.txt y -ascii 

% Tallenna kaksi muuttuja x ja y nimettyyn tiedostoon peräkkäin (ei
% kovin hyödyllistä)
save sinc_pulssi2.txt x y -ascii

% Luodaan tiedosto jossa on kaksi saraketta. Ensimmäisessä sarakkeessa
% on vektorin x elementit ja toisessa sarakkeessa vektorin y
% elementit. Tämä vaatii että x ja y ovat pystyvektoreita.
z = [x y];
save sinc_pulssi3.txt z -ascii