Pouzivam typ pisma " Ariel, tucne, size 12" Sirku okna pre editor NotePad upravte podla sirky tohoto riadku ||||||| Mate otvorene jednak command window ML a jednak okno tohto graf_tan.txt Od riadku *start* zacina ukazka. Jej prikazy editujte vo svojom ML-command window a pozorujte odozvu ML * start * clear, clc % zacnime s "cistym" prostredim .... close all % prikaz zavrie vsetky graficke okna... a zaciname :-) x = linspace( -pi/2, pi/2, 401 ); y = tan(x); % zamer: chceme graf funkcie tangens plot( x, y ) % plot dopadol zalostne... ziaden tangens? Preco?? Pretoze % mechanizmus automierky nemohol zvladnut situaciu... Preco? % % opytajme sa na tangens v poslednej zlozke vektora x! y(401) ans = 1.6331e+016 % vlastne tu malo by to byt "inf", zrejme pi/2 nie je "ciste" pi-pol % brat do uvahy hodnoty "okolo" nuly a zaroven radu % 10^16 nie je mozne ! % % definujme premennu "okno", aby sme prikazom "axis(okno)" % vypli automierku ! okno= [-pi/2 pi/2 -20 20]; % je to vektor [ xmin, xmax, ymin, ymax ] axis(okno), shg % tento prikaz zariadi vypnutie automierky, % graf.okno sa inovuje... % graf jednej vetvy je OK. % Teraz by sme radi znazornili graf tangensu na [ -pi/2, 3*pi/2 ] % Urobime nieco prirodzene, co sa vsak ukaze ako nevhodne, % ale poucne! u= x + pi; % vektor x-hodnot pre dalsiu vetvu tangensu % mnozinu povodnych uzlov "posunieme" % o hodnotu pi x2 = [x u]; % zlepime hodnoty povodneho a posunuteho % intervalu % nove hodnoty tangensu nie su potrebne, % platia stare, % pretoze tan( ) je periodicka funkcia y2 = [y y]; % vektor y2 bude mat rovnaku length ako x2 plot(x2, y2), shg % ale automierku znovu treba vypnut! okno2 = [ -pi/2 3*pi/2 -20 20]; % nove okno je nutne ... pre "axis" axis(okno2), shg % dobre to nedopadlo(!), lebo "plot" spojil body "ponad" pi/2 a % tomu treba zabranit ... % ved predsa vetvy tangensu nic nespaja! % % graf nakreslime na dva razy, ked pouzijeme prikaz "hold on" a % ku prvemu grafu sa prida druhy (bez toho, aby sa prvy stratil...) close % zavreli sme graficke okno % avsak premenne ostavaju plot( x, y ), axis( okno2 ) % zatial prva vetva hold on, plot( u, y ), shg % teraz je to uz dobre... grid on % ak chceme siet ... %%-------------------------------END-----------------------------------%%