Attachment 'graf_tan.txt'

Download

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-----------------------------------%%

Attached Files

To refer to attachments on a page, use attachment:filename, as shown below in the list of files. Do NOT use the URL of the [get] link, since this is subject to change and can break easily.

You are not allowed to attach a file to this page.