155GIT1 / 7. cvičení: Porovnání verzí

Z GeoWikiCZ
mBez shrnutí editace
 
(Není zobrazeno 37 mezilehlých verzí od 3 dalších uživatelů.)
Řádek 1: Řádek 1:
{{Geoinformatika}}
{{Geoinformatika}}
{{Cvičení|155GIT1|7|Grafy funkcí, uživatelské funkce}}
{{Cvičení|155GIT1|7|Uživatelské funkce}}


== Náplň cvičení ==
== Náplň ==


# grafické okno <code>figure</code>
# grafy funkcí <code>plot(), subplot(), axis()</code>
# uživatelské funkce
# uživatelské funkce
# chybová hlášení, nápověda k funkci
# drobné úlohy


== Ukázky ==
== Ukázky ==


=== Grafické okno ===
=== Uživatelské funkce ===
 
'''Syntax'''
function [vystup1,vystup2,...,vystupn] = nazev(vstup1,vstup2,...,vstupn)
end  % end neni povinne
 
Každá funkce musí být uložena ve vlastním souboru s příponou <tt>.m</tt>. Název souboru musí odpovídat názvu funkce! Klíčové slovo <tt>function</tt> musí být na prvním řádku souboru.


<source lang=octave>
<source lang=octave>
figure
function [min_x, max_x, mean_x] = stat(x)
 
    min_x = min(min(x));
    max_x = max(max(x));
    mean_x = mean(x(:));    % stejný výsledek jako mean_x = mean(mean(x))
 
end
</source>
</source>


=== Grafy funkcí ===
'''Příklad volání'''


* <code>plot()</code> {{bullet}} [http://geo.fsv.cvut.cz/user/gin/git1/matlab/matlab.pdf#86 argumenty]
<source lang=octave>
data = rand(5)
[mi, ma, mn] = stat(data);
fprintf('min  = %.3f\nmax  = %.3f\nmean = %.3f\n', mi, ma, mn)
</source>
 
=== Chybová hlášení, nápověda ===


<source lang=octave>
<source lang=octave>
x = [0:pi/50:2*pi];
function [min_x, max_x, mean_x] = stat(x)
y = sin(x);
% vypocet minimalni, maximalni a prumerne hodnoty matice
plot(x, y)
 
%
    if any(any( imag(x)~=0 ))
% nové okno
        error('stat: Vstupni matice obsahuje komplexni cisla.')
h = figure
    end
plot(x, y 'r+')
    min_x = min(min(x));
%
    max_x = max(max(x));
% více grafů najednou
    mean_x = mean(x(:));
z = cos(x);
 
plot(x, y, 'r+', x, z, 'b*')
end
</source>
</source>


{{fig|octave-graf-1|<code>plot(x, y, 'r+', x, z, 'b*')</code>}}
Vyvolání nápovědy k uživatelské funkci:
 
                         
* <code>subplot()</code>
<source lang=octave>
help stat
</source>


Volání funkce (chybná vstupní data):
<source lang=octave>
<source lang=octave>
x = 1:100;
data1 = [data, 5+2i * ones(5,1)]
figure
[mi, ma, mn] = stat(data1)
% první
subplot(2, 2, 1)
plot(x, x)
% druhý
subplot(2, 2, 2)
plot(x, sqrt(x))
% třetí
subplot(2, 2, 3)
plot(x, log(x))
% čtvrtý
subplot(2, 2, 4)
plot(x, x.^2)
</source>
</source>


{{fig|octave-graf-2|<code>subplot()</code>}}
=== Příkaz return ===
* je-li v běhu funkce vykonán příkaz <tt>return</tt>, je okamžitě činnost funkce ukončena (bez vykonání evtl.dalších příkazů) a na pozici jejího volání jsou vráceny aktuálně vypočtené hodnoty výstupních parametrů funkce
 
=== Drobné úlohy ===
 
'''Funkce pro načtení neformátovaných dat po řádcích'''


* <code>axis()</code> - omezení oblasti grafu
např. záměna čárky za tečku:


<source lang=octave>
<source lang=octave>
x = 0:0.1:5;
function teckazacarku(file)
y = exp(x);
% vyhleda vsechny carky v souboru a zameni je za tecku
% celý graf
% IN: file ... jmeno vstupniho souboru (textovy retezec)
subplot(2, 1, 1)
 
plot(x, y)
if ischar(file) == 0
% výsek x <1, 2>; y <0, 10>
    error('Jmeno vstupniho souboru musi byt textovy retezec.');
subplot(2, 1, 2)
end
plot(x, y)
 
axis([1,2,0,10])
vstup = fopen(file,'r');
 
filenew = [file(1:end-4), '_upr', file(end-3:end)];  % jméno nového souboru
vystup = fopen(filenew,'w');
 
konec = 0;  % indentifikátor konce souboru
while konec==0
    radek = fgets(vstup);  % fgets() načte 1 řádek souboru (se znakem \n)
    for i=1:length(radek)
        if radek(i) == ','
            radek(i) = '.';
        end
    end
    fprintf(vystup,'%s',radek);
    konec = feof(vstup);  % feof() vrátí pro konec souboru hodnotu 1, jinak vrací hodnotu 0  
end
 
fclose all;
fprintf('Data ulozena do souboru %s\n',filenew);
</source>
</source>


{{fig|octave-graf-3|<code>axis()</code>}}
datový soubor [http://athena.fsv.cvut.cz:8000/GIT1/cv7/carky.txt carky.txt] ke stažení


* popisky
Volání funkce:


<source lang=octave>
<source lang=octave>
x = [0:pi/100:pi];
teckazacarku('carky.txt')
figure
plot(x, sin(x))
title('Funkce sinus')
xlabel('hodnota v rad')
ylabel('sin(x)')
legend('prubeh funkce sinus')
</source>
</source>


{{fig|octave-graf-4|Graf s popisky}}
== Úlohy ==
 
* [[155GIT1 / 7. cvičení / Příklady|7. cvičení - příklady]]
<!--
 


=== Uživatelské funkce ===


Syntax:
function [vystup] = nazev(vstup)


== Úlohy ==
{{Cvičení|155GIT1|7|Test Matlab/Octave}}
; Zadání testu
* https://moodle-vyuka.cvut.cz/course/view.php?id=5130
; Výsledky testu
* [https://docs.google.com/spreadsheets/d/1We6CEro7T5q766HukoxxTLYYTlvEdsV0i1B1VcHUD6s/edit?usp=sharing souhrnné výsledky]
-->
<!-- -->

Aktuální verze z 13. 4. 2023, 03:59

Uživatelské funkce

Náplň

  1. uživatelské funkce
  2. chybová hlášení, nápověda k funkci
  3. drobné úlohy

Ukázky

Uživatelské funkce

Syntax

function [vystup1,vystup2,...,vystupn] = nazev(vstup1,vstup2,...,vstupn)

end  % end neni povinne

Každá funkce musí být uložena ve vlastním souboru s příponou .m. Název souboru musí odpovídat názvu funkce! Klíčové slovo function musí být na prvním řádku souboru.

 function [min_x, max_x, mean_x] = stat(x)

    min_x = min(min(x));
    max_x = max(max(x));
    mean_x = mean(x(:));    % stejný výsledek jako mean_x = mean(mean(x))

 end

Příklad volání

 data = rand(5)
 [mi, ma, mn] = stat(data);
 fprintf('min  = %.3f\nmax  = %.3f\nmean = %.3f\n', mi, ma, mn)

Chybová hlášení, nápověda

 function [min_x, max_x, mean_x] = stat(x)
 % vypocet minimalni, maximalni a prumerne hodnoty matice

     if any(any( imag(x)~=0 ))
         error('stat: Vstupni matice obsahuje komplexni cisla.')
     end
     min_x = min(min(x));
     max_x = max(max(x));
     mean_x = mean(x(:));

 end

Vyvolání nápovědy k uživatelské funkci:

 help stat

Volání funkce (chybná vstupní data):

 data1 = [data, 5+2i * ones(5,1)]
 [mi, ma, mn] = stat(data1)

Příkaz return

  • je-li v běhu funkce vykonán příkaz return, je okamžitě činnost funkce ukončena (bez vykonání evtl.dalších příkazů) a na pozici jejího volání jsou vráceny aktuálně vypočtené hodnoty výstupních parametrů funkce

Drobné úlohy

Funkce pro načtení neformátovaných dat po řádcích

např. záměna čárky za tečku:

 function teckazacarku(file)
 % vyhleda vsechny carky v souboru a zameni je za tecku
 % IN: file ... jmeno vstupniho souboru (textovy retezec)

 if ischar(file) == 0
     error('Jmeno vstupniho souboru musi byt textovy retezec.');
 end

 vstup = fopen(file,'r');

 filenew = [file(1:end-4), '_upr', file(end-3:end)];  % jméno nového souboru
 vystup = fopen(filenew,'w');

 konec = 0;  % indentifikátor konce souboru
 while konec==0
     radek = fgets(vstup);  % fgets() načte 1 řádek souboru (se znakem \n)
     for i=1:length(radek)
         if radek(i) == ','
             radek(i) = '.';
         end
     end
     fprintf(vystup,'%s',radek);
     konec = feof(vstup);   % feof() vrátí pro konec souboru hodnotu 1, jinak vrací hodnotu 0 
 end

 fclose all;
 fprintf('Data ulozena do souboru %s\n',filenew);

datový soubor carky.txt ke stažení

Volání funkce:

 teckazacarku('carky.txt')

Úlohy