cwiczenia6

Wszelkie prawa do zada

, plików i instrukcji s

zastrze

one, nikt poza prowadz

cymi zaj

cia i studentami IA na AGH nie ma prawa ich

wykorzystywa

. A tych co maj

takie prawa, prosz

, aby nie umieszczali ich w innym miejscu w sieci.

Grafika 3D

Prosz

, aby ka

de zadanie było rozwi

zane w innym skrypcie.

Wszystkie skrypty i funkcje maj

zapisane w dowolnym folderze nale

żą

cym

do Pa

stwa.

ęść

I – funkcja plot3

plot3(x,y,z) rysuje funkcje np. w postaci parametrycznej, przy czym x, y, z – s

to wektory o takim

samym rozmiarze (wszystkie wektory maj

mie

tyle samo elementów).

Zad1. Prosz

przekopiowa

tre

ść

poni

szego skryptu na swój komputer.

figure(1)
t=0:.1:12; x=t; y=sin(t); z=cos(t);
plot3(x,y,z)
grid

;

W drugim oknie prosz

narysowa

ten sam wykres, ale dla zrównanych jednostek

na osiach, za pomoc

powy

szych polece

oraz polecenia

axis equal
Wykorzystuj

c edytor wła

ciwo

ci oraz przegl

dark

obiektów

prosz

opisa

osie „x”, ”y” i „z”.

Zad 2. Oto współrz

dne punktów: x = [1 1 -1 1]; y = [4 4 4 4]; z = [-1 1 1 -1];

Prosz

narysowa

trójk

t, który powstanie z powy

szych punktów, czerwon

lini

o grubo

ci 6

(…'LineWidth',6).

ęść

II – wykresy 3D

siatki powierzchni
i powierzchnie

2.1

Funkcja meshgrid i mesh

Funkcja meshgrid

[X , Y] = meshgrid (x, y) przekształca wektory x i y w macierze X i Y, po to, aby powstała
siatka punktów, które u

ywamy do znajdowania warto

ci funkcji dwóch zmiennych

i generowania trójwymiarowych wykresów. Wiersze wynikowej macierzy X s

kopiami wektora

x, kolumny wynikowej macierzy Y s

kopiami wektora y.

Funkcja mesh

mesh (X, Y, Z) rysuje szkielet siatki

Zad 3. Prosz

stworzy

skrypt, który zawiera poni

sze polecenia:

clear

all

x=-5:5;
y=-10:10;
[X, Y]=meshgrid(x,y)

Prosz

sprawdzi

jak została wygenerowana siatka X i Y, a potem wstawi

; po poleceniu

[X, Y]=meshgrid(x,y);

I wykorzystuj

c powy

siatk

punktów prosz

narysowa

wykres funkcji:

݂ሺݔ, ݕሻ = ܺ

ଷ

+ 3ܻܺ

ଶ

+ 12ܻܺ

W drugim oknie dla tych samych danych prosz

narysowa

powierzchni

za pomoc

polecenia surf(X,Y,Z). Prosz

doło

wska

nik colorbar do obu okien.

Funkcja surf(X,Y,Z) – rysuje powierzchni

kolorow

opisan

przez macierze X,Y,Z.

2.2. Rysowanie siatki powierzchni za pomoc

: mesh(x,y,z), meshc(x,y,z) - meshz(x,y,z)

i waterfall(x,y,z)

Zad 4. Dla nast

puj

cych danych: x=0:0.5:10;y=0:.5:5; z=sin(y)'*cos(x);

w jednym oknie i na 4 wykresach prosz

narysowa

funkcje

mesh(x,y,z), meshc(x,y,z) -

meshz(x,y,z) i waterfall(x,y,z). Korzystaj

c z poni

szej tabelki prosz

opisa

wykresy poprzez

dodanie tytułów.

mesh

Wykres siatkowy

meshc Wykres siatkowy + poziomice

meshz Wykres siatkowy + zasłony na ko

cach

ribbon Wykres wst

ąż

kowy

surf

Wykres powierzchniowy

surfc

Wykres powierzchniowy + poziomice

surfl

Wykres powierzchniowy + cieniowanie

waterfall

funkcja waterfall(x, y, z, c) działa podobnie jak meshz z ta ró

nic

e nie

rysuje linii odpowiadaj

cych kolumnom macierzy

2.2. Rysowanie powierzchni za pomoc

: surf, surfc, surfl

Zad 5. Dla nast

puj

cych danych: x=0:0.5:10;y=0:.5:5; z=sin(y)'*cos(x);

w pierwszym oknie na 3 wykresach prosz

narysowa

funkcje surf, surfc, surfl

w dugim oknie na trzech wykresach dla surf prosz

sprawdzi

jak działa przejrzysto

ść

alpha

ustawiaj

c w kolejnych oknach alpha(0) alpha(0.5) i alpha(1). Prosz

sprawdzi

ukrywanie kraw

dzi:

hidden on

Wyświetlanie ukrytych krawędzi

hidden off

Domyślny, ukrywa niewidoczne krawędzie

w trzecim oknie sterujemy kolorem i wy

wietlaniem siatki oraz kolorem powierzchni

za pomoc

instrukcji:

surf(x,y,z,'FaceColor',[0.5 0.5 0.5], 'EdgeColor','red')
surf(x,y,z,'FaceColor','blue', 'EdgeColor','none')

–prosz

sprawdzi

jak działa dodatkowo

camlight left ( mo

e zamiast left lepiej b

dzie:

headlight, left, right

lighting phong (phong, flat, gouraud, none)

li nie wida

zmian, prosz

zmieni

kolor siatki i powierzchni

w czwartym oknie

wiczymy cienie - potrzebujemy do tego trzy wykresy funkcji

z=x*exp(-x

-y

), pierwszy jest na rysunku poni

ej:

drugi wykres powinien by

narysowany za pomoc

funkcji:

surfl(x,y,z), a trzeci mie

dodatkowo doło

do surfl(x,y,z) funkcj

shading interp;

shading flat

Powierzchnia z dyskretnymi kolorami

shading interp

Powierzchnia z wypełnieniem kolorami interpolowanymi

shading faceted

Powierzchnia z dyskretnymi kolorami i siatka

w pi

tym oknie zobaczymy jakie s

liwo

ci wy

wietlania poziomic. Znowu podzielimy

okno na 4 wykresy. Bazowym wykresem niech b

dzie

x=0:0.1:10;y=0:.1:5; z=sin(y)'*cos(x); meshc(x,y,z)

Korzystaj

c z polecenia:

contour(x,y,z,n) – kre

lenie poziomic (wykres konturowy), gdzie n – liczba poziomic

prosz

wykre

3 poziomice

i w kolejnym okienku trzy poziomice opisane:

[c,h]=contour(x,y,z,n)
clabel(c,h);

Prosz

sprawdzi

jak działa polecenie contourf(x,y,z,5), w ostatnim, czwartym wykresie.

ęść

IV – dodatkowe wykresy w grafice 3D

bar3, ribbon

Zad 6. Dowoln

funkcj

3D prosz

przedstawi

za pomoc

polecenia

ribbon(Y,Z)

Zad 7. Dla wylosowanych danych Y=round(rand(10,4)*100) na trzech ró

nych wykresach prosz

sprawdzi

jak działa:

bar3(Y,

'detached'

)

bar3(Y,

'grouped'

)

bar3(Y,

'stacked'

) W ka

dym wykresie prosz

zmieni

domy

lny punkt obserwacji

view

Zmiana domy

lnego punktu obserwacji

view(azymut, elewacja)

Okre

la punkt obserwacyjny za pomoc

azymutu i

elewacji

view([x,y,z])

Okre

la punkt obserwacji w układzie kartezja

skim

view(2)

Obserwacja azymut=0, elewacja=90

view(3)

Domy

lny punkt obserwacji: azymut=-37.5 ,

elewacja= 30

ęść

V – animacje

funkcja comet(x,y); - to animacja ruchu komety. x,y – to współrz

dne poło

enia.

Zad 8. Prosz

napisa

skrypt, w którym za pomoc

animacji mo

emy zobaczy

ruch

komety po torze, który jest elips

Funkcje getframe i movie, to kolejny przykład animacji. Ka

de polecenie plot generuje

nowe okno, zawarto

ść

tego okna mo

na przy pomocy polecenie getframe zapisa

np.: w tablicy.

Utworzone tablice odtwarzamy wykorzystuj

c polecenie movie. Jako przykład wykorzystania tych

funkcji podaje si

poni

szy skrypt :

for k = 1:16

plot(fft(eye(k+16)))
axis equal
M(k) = getframe;
end
movie(M,30)

Zad 9. Prosz

napisa

skrypt, w którym za pomoc

animacji mo

emy zobaczy

kolejne

etapy rysowania funkcji sinus.

Przy braku pomysłu na rozwi

zanie prosz

przej

c do

kolejnego zadania, jest tam ukryta podpowied

drawnow – aktualizacja okna

Zad 10. Oto przykład rysowania funkcji przy pomocy polecenia drawnow

x = -2*pi:pi/20:2*pi;
for i=1:8
plot(x(1:i*10),cos(x(1:i*10)))
axis([-2*pi 2*pi -1 1])
drawnow
pause(0.1)
end;

Prosz

wykorzystuj

c dranow narysowa

dowolny wykres funkcji oraz dowoln

powierzchni

ęść

VI – walec, kula, wielok

t, elipsoida

cylinder

Generacja walca

elipsoid

Generacja elipsoidy

(nie działa w pracowni)

fill3

Generacja wielok

( nie działa w pracowni)

sphere

Generacja kuli

Zad 11. Prosz

wygenerowa

figury podane w powy

szej tabeli i za pomoc

polecenia slice

przeci

ąć

dwie z tych figur dowoln

płaszczyzn