FIC MATLAB 2017-2/Aula-12

%-- 19-10-2017 19:08:59 --%
clc
% Plots em 2D e 3D
edit
simulink
demo simulink
t = linspace(0, 2*pi, 100)
x = cos(t);
t = linspace(0, 2*pi, 20);
x = cos(t);
plot(t, x)
stem(t, x)
bars(t, x)
bar(t, x)
stairs(t, x)
w = linspace(0, 50*pi, 1000);
H = 1 / sqrt(1 + w^2)
H = 1 / sqrt(1 + w.^2)
H = 1 ./ sqrt(1 + w.^2)
format compact
H = 1 ./ sqrt(1 + w.^2);
plot(w, H)
grid on
grid minor
w = linspace(0, 5000*pi, 1000);
w = linspace(0, 5000*pi, 10000);
H = 1 ./ sqrt(1 + w.^2);
plot(w, H)
grid on
grid minor
xlim([0 50*pi])
semilogx(w, H)
grid on
semilogy(w, H)
grid on
loglog(w, H)
grid on
plot(w, H)
loglog(w, H)
grid on
w = linspace(0, 5000*pi, 100000);
H = 1 ./ sqrt(1 + w.^2);
semilogx(w, H)
w = logspace(-5, 5, 1000);
semilogx(w, H)
H = 1 ./ sqrt(1 + w.^2);
semilogx(w, H)
semilogy(w, H)
loglog(w, H)
grid minor
% Gráficos polares
theta = linspace(0, 2*pi)
theta = linspace(0, 2*pi).
theta = linspace(0, 2*pi);
r = cos(theta);
polar(theta, r)
theta = linspace(0, pi);
r = cos(theta);
polar(theta, r)
theta = linspace(0, 2*pi);
r = cos(2*theta);
polar(theta, r)
r = cos(3*theta);
polar(theta, r)
r = cos((1/3)*theta);
polar(theta, r)
theta = linspace(0, 8*pi);]
theta = linspace(0, 8*pi);
r = cos((1/3)*theta);
polar(theta, r)
theta = linspace(0, 8*pi, 10000);
r = cos((5/6)*theta);
polar(theta, r)
r = cos(17*theta);
polar(theta, r)
doc polar
r = 2*cos(17*theta);
polar(theta, r)
z = 3 + 1j*4
z = 3 + 1i*4
z = 3 + 1k*4
k = sqrt(-1)
z = 3 + k*4
abs(z)
angle(z)
atan(4/3)
atan2(3,4)
atan2(4,3)
angled(z)
angle(z)
angle(z) * (180 / pi)
r = theta
theta = linspace(0, 2*pi, 10000);
r = theta;
polar(theta, r)
% plotyy
x = -1:0.001:1;
y1 = x.^2;
y2 = cos(2*pi*x);
plot(x, y1, 'b', x, y2, 'r')
y2 = 5000*cos(2*pi*x);
plot(x, y1, 'b', x, y2, 'r')
doc plotyy
plottt(x, y1, 'b', x, y2, 'r')
plotyy(x, y1, 'b', x, y2, 'r')
plotyy(x, y1, x, y2)
doc plotyy
x = -1:0.001:1;
close all
tau = x.^2;
plot(x, tau)
ylabel('Torque (tau)')
ylabel('Torque (\tau)')
xlabel('X-men (\upsilon)')
xlabel('X-men (\nu)')
xlabel('X-men (\alpha)')
xlabel('X-men (\alpha_0)')
ylabel('Torque (\int_0^1 \tau)')
titlr('\int_{-\infty}^{\upsilon} \pi^2 d\Omega')
title('\int_{-\infty}^{\upsilon} \pi^2 d\Omega')
close all
% Gráficos em 3D
t = linspace(0,10,1000);
x = cos(t);
y = sin(t);
z = t;
plot3(x,y,z)
grid on
t = linspace(0,50,10000);
x = cos(t);
y = sin(t);
z = t;
plot3(x,y,z)
grid on
xlabel('x')
ylabel('y')
zlabel('z')
comet3(x,y,z)
lorenz
grid on
xlabel('x')
ylabel('y')
zlabel('z')
view(45,90)
view(45,0)
view(0,0)
view(90,0)
% Gráficos de superfície
[X,Y] = meshgrid(0:4:20,0:2:10)
surf(X,Y,0)
surf(X,Y,zeros(size(X)))
[X,Y] = meshgrid(-20:0.5:20);
R = sqrt(X.^2 + Y.^2);
surf(X,Y,R)
Z = sin(R)./R;
surf(X,Y,Z)
axis equal
surf(X,Y,Z)
doc axis
shading faceted
shading flat
shading interp
shading flat
shading interp
R = sqrt(X.^2 + Y.^2) + eps;
R
Z = sin(R)./R;
surf(X,Y,Z)
shading interp
doc light
light('Position',[-1 0 0],'Style','local')
shading interp
light('Position',[1 5 0],'Style','local')
surf(X,Y,Z)
shading interp
light('Position',[1 5 0],'Style','local')
light('Position',[1 5 20],'Style','local')
doc membrane
openExample('matlab_featured/MatlabLogoExample')
surf(X,Y,Z)
mesh(X,Y,Z)
% Curva de niveis
contour(X,Y,Z)
contour(X,Y,Z,5)
contour(X,Y,Z,2)
contour(X,Y,Z,3)
contour(X,Y,Z,4)
surfc(X,Y,Z)
shading interp
meshc(X,Y,Z)
surf(X,Y,Z)
hold on
[X1, X2] = meshgrid(-1:0.1:1);
[X1, Y1] = meshgrid(-1:0.1:1);
G1 = sqrt(X1.^2 + Y1^2) + 1.5;
G1 = sqrt(X1.^2 + Y1.^2) + 1.5;
surf(X1, Y1, G1)
G1 = sqrt(1 - X1.^2 - Y1.^2) + 1.5;
surf(X1, Y1, G1)
G1
iscomplex(Z)
G1 = sqrt(4 - X1.^2 - Y1.^2) + 1.5;
surf(X1, Y1, G1)
surf(X, Y, Z); hold on; surf(X1, Y1, G1); surf(X1, Y1, G2);
G2 = -sqrt(4 - X1.^2 - Y1.^2) + 1.5;
surf(X, Y, Z); hold on; surf(X1, Y1, G1); surf(X1, Y1, G2);
surf(X, Y, Z); hold on; surf(X1, Y1, G1); surf(X1, Y1, G2); axis equal
surf(X, Y, 20*Z); hold on; surf(X1, Y1, G1); surf(X1, Y1, G2); axis equal
G1 = sqrt(4 - X1.^2 - Y1.^2) + 22;
G2 = -sqrt(4 - X1.^2 - Y1.^2) + 22;
surf(X, Y, 20*Z); hold on; surf(X1, Y1, G1); surf(X1, Y1, G2); axis equal
doc sphere
[Gx, Gy, Gz] = sphere
surf(X, Y, 20*Z); hold on; surf(Gx, Gy, Gz + 22); axis equal
shading interp
[Gx, Gy, Gz] = sphere
surf(X, Y, 20*Z); hold on; surf(4*Gx, 4*Gy, 4*Gz + 22); axis equal
a = 1:180;
ii = 1; while true; view(a(mod(ii,length(a))+1),15); end
ii = 1; while true; view(a(mod(ii,length(a))+1),15); pause(0.1) end
ii = 1; while true; view(a(mod(ii,length(a))+1),15); pause(0.1); end
ii = 1; while true; view(a(mod(ii,length(a))+1),15); pause(0.5); end
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.5); end
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.05); end
doc light
xlabel('x'); ylabel('y'); zlabel('z')
light('Position',[30 0 20],'Style','local')
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.05); end
shading interp
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.05); end
a = 1:360;
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.05); end
[X,Y] = meshgrid(-20:0.05:20);
R = sqrt(X.^2 + Y.^2) + eps;
Z = sin(R)./R;
[Gx, Gy, Gz] = sphere(100);
surf(X, Y, 20*Z); hold on; surf(4*Gx, 4*Gy, 4*Gz + 22); axis equal
light('Position',[30 0 20],'Style','local')
shading interp
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.05); end
light('Position',[0 0 40],'Style','local')
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.05); end
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.01); end
transpdemo
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.01); end
surf(X, Y, 20*Z); hold on; surf(4*Gx, 4*Gy, 4*Gz + 22); axis equal
surf(X, Y, 20*Z); hold on; surf(4*Gx, 4*Gy, 4*Gz + 22); axis equal; shading interp
ii = 1; while true; view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.01); end
close all
surf(X, Y, 20*Z); hold on; surf(4*Gx, 4*Gy, 4*Gz + 22); axis equal; shading interp
light('Position',[30 0 20],'Style','local')
light('Position',[0 0 40],'Style','local')
figure
transpdemo
ii = 1; while true; figure(1); view(a(mod(ii,length(a))+1),15); figure(2); view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.01); end
ii = 1; while true; figure(1); view(-a(mod(ii,length(a))+1),15); figure(2); view(a(mod(ii,length(a))+1),15); ii = ii+1; pause(0.01); end
im2 = imread('imagem2.jpeg');
c = 0:0.01:1;
b = [c c(end:-1:1)];
ii = 1; while true;figure(1); view(-a(mod(ii,length(a))+1),15); figure(2); view(a(mod(ii,length(a))+1),15); figure(3); imshow(uint8(b(mod(ii,length(b))+1)*(randi([-130 130],600,600,3) + double(im2)))); pause(0.01); ii = ii+1; end
ii = 1; while true;figure(1); view(-a(mod(ii,length(a))+1),15); figure(2); view(a(mod(ii,length(a))+1),15); figure(3); imshow(uint8(b(mod(ii,length(b))+1)*(randi([-130 130],600,600,3) + double(im2)))); ii = ii+1; end
ii = 1; while true;figure(1); view(-a(mod(ii,length(a))+1),15); figure(2); view(a(mod(ii,length(a))+1),15); figure(3); imshow(uint8(b(mod(ii,length(b))+1)*(randi([-130 130],600,600,3) + double(im2)))); ii = ii+1; pause(eps); end
ii = 1; while true;figure(1); view(-a(mod(ii,length(a))+1),15); figure(2); view(a(mod(ii,length(a))+1),15); figure(3); imshow(uint8(b(mod(ii,length(b))+1)*(randi([-130 130],600,600,3) + double(im2)))); view(-a(mod(ii,length(a))+1),15); ii = ii+1; pause(eps); end
ii = 1; while true;figure(1); view(-a(mod(ii,length(a))+1),15); figure(2); view(a(mod(ii,length(a))+1),15); figure(3); imshow(uint8(b(mod(ii,length(b))+1)*(randi([-130 130],600,600,3) + double(im2)))); view(-a(mod(ii,length(a))+1),90); ii = ii+1; pause(eps); end
animacaum

% animacaum.m

% Obs: rodar o script ao som de 'Shooting stars'

clear all; close all; clc

[X,Y] = meshgrid(-20:0.05:20);
R = sqrt(X.^2 + Y.^2) + eps;
Z = sin(R)./R;

figure(1)
surf(X, Y, 20*Z)
hold on;

[Gx, Gy, Gz] = sphere(100);
surf(4*Gx, 4*Gy, 4*Gz + 22)
axis equal

light('Position',[30 0 20],'Style','local')
light('Position',[0 0 40],'Style','local')

shading interp

figure(2)
transpdemo

im2 = imread('imagem2.jpeg');
%figure(3)
%imshow(im2)

a = 1:360;
c = 0:0.01:1;
b = [c c(end:-1:1)];

ii = 1;

while true;
    
    figure(1);
    view(-a(mod(ii,length(a))+1),15);
    
    figure(2);
    view(a(mod(ii,length(a))+1),15);
    
    figure(3);
    imshow(uint8(b(mod(ii,length(b))+1)*(randi([-130 130],600,600,3) + double(im2))));
    view(-a(mod(ii,length(a))+1),90);
    
    ii = ii+1;
    
    pause(eps);

end