Contents

H3teht4

%{
Olkoon
f(x,y)=sin(3y - x^2+1)+cos(2y^2-2x). Piirrä pintakuva ja korkeuskäyräpiirros,
jälkimmäinen sekä contour että ezcontour-funktioilla. Tässä on mahdollisuus
kokeilla yllä esiteltyjä korkeuskäyrien valitsemistapoja, myös clabel.
Ota alueeksi vaikka  [-2 2 -1 1]
%}

1) ezcontour

close all
ezcontour('sin(3*y - x^2+1)+cos(2*y^2-2*x)',[-2 2 -1 1])
% Onpas helppoa!

2) contour % Lisää joustavuutta ja mahdollisuuksia.

figure
x=-2:.01:2; y=-1:.01:1;
[X,Y]=meshgrid(x,y);
Z=sin(3*Y - X.^2+1)+cos(2*Y.^2-2*X);
contour(x,y,Z,20)    % 20 korkeuskäyrää.

% Antaa seuraavan pyyhkiä edellinen.

Z_vali=[min(min(Z)) max(Z(:))]  % Tahallaan kahdella tavalla
korkeudet=linspace(Z_vali(1),Z_vali(2),30);
contour(x,y,Z,korkeudet)   % Korkeuskäyrät näille korkeuksille.

[C,h]=contour(x,y,Z,korkeudet);
clabel(C,h,korkeudet(1:6:30))   % Osalle käyristä numeroarvot

clabel(C,h,'manual')    % Lisätään arvoja manuaalisesti

Z_vali =

   -1.9999    2.0000

 
    Please wait a moment...
 
   Carefully select contours for labeling.
   When done, press RETURN while the Graph window is the active window.

Published with MATLAB® 7.9