% s1_.m
% Versao de 2018-01-29

% Isto e' um comentario

clear all
close all

% Apenas no Octave
% Descomente a linha seguinte se estiver a usar o Octave
more off

% Mais um cometario
N=20;
n=0:1:N-1 % qual o efeito de nao se escrever ; no final
% Mais um cometario
n=0:1:N-1;
n

n1 = 0:4:N-1
n2 = 0:0.5:N-1

for i=1:1:N,
    x0(i)=sin(2*pi/N*n(i));
end

i=1;
while i<=N,
    x1(i)=sin(2*pi/N*n(i));
    i=i+1;
end

% e' muito mais: compacto, facil de escrever, rapido a executar
x2 = sin(2*pi/N*n);

% verificar que x0, x1 e x2 sao iguais
sum(abs(x1-x0))
sum(abs(x2-x0))
max(abs(x2-x0)) % mas afinal, o que e' que isto faz ?

% numeros complexos
i = sqrt(-1);
% para evitar definir a variavel i (ou j ou outra qualquer)
% o matlab permite representar o numero sqrt(-1) por 1i;
1i
UmNumeroComplexo = 123+456*1i;
OutroNumeroComplexo = 123*exp(1i*pi/2);
a = 4+5*i;
real(a)
imag(a)
abs(a)
angle(a)
angle(a)*180/pi

a0 = [1+2*i, 3+4*i, -2+5*i, 20-30*i]
a1 = transpose(a0)
a2 = conj(a0)
a3 = a0'

a4 = [a0; a2]

size(a0)
size(a1)
size(a2)
size(a3)
size(a4)

length(a0)
length(a4)

a=3;
if a==2,
    b=1;
elseif a~=3,
    b=2;
else
    b=3;
end

A=[ 0  1  2  3  4;
    5  6  7  8  9;
   10 11 12 13 14;
   20 21 22 23 24;
   30 31 32 33 34;
   ]

B=A*4*5

% qual a diferenc,a ?
C1 = A*B
C2 = A*B'
C3 = A.*B
C4 = A+B*sqrt(-1)

A = [ 1 2 3 4 5 6];
B = A+A*sqrt(-1)

D1=A.*B
D2=A./B
D3=A.^B

who

whos

% capacidades graficas
Np = 50;
n = 0:1:Np-1;
x1 = sin(2*pi/Np*n);

figure(1);
clf;

subplot(2,1,1);
hold on;
stem(n, x1, ...
    'MarkerSize',4, ...
    'MarkerFaceColor',[0 0 1], ...
    'MarkerEdgeColor',[0 0 0], ...
    'Color',[0.75 1 0.5], ...
    'LineWidth',1);


subplot(2,1,2);
hold on;
stem(n, x1, ...
    'MarkerSize',4, ...
    'MarkerFaceColor',[1 0 0],...
    'MarkerEdgeColor',[0 0 0],...
    'LineWidth',1);
disp('Isto esta' em pausa, observe a figura 1 e carregue qualquer tecla');
pause

figure(2); clf;
stem(x1, 'r*')
disp('Mais uma pausa, observe a figura 2 e carregue qualquer tecla');
pause

stem(n, x1, 'g+') % verificar o eixo horizontal
disp('verificar o eixo horizontal (1)');
pause

stem(n+10, x1, 'r*') % verificar o eixo horizontal
disp('verificar o eixo horizontal (2)');
pause

disp('stem, plot, stairs');
stem(n, x1, 'r')
hold on
stairs(n, x1, 'b')
plot(n, x1, 'm')
grid on

title('Titulo do grafico')
xlabel('Etiqueta X')
ylabel('Etiqueta Y')

disp('Observe a figura');
pause

% Execucao de funcoes criadas pelo utilizador
clear all

y=[zeros(1,10), ones(1,15), zeros(1,5)+3, sin(2*pi*(0:1:20)/11)];
n=0:1:length(y)-1;
figure(1); clf;
stem(n,y);

[y(1), fsn(y,1)] % visualize o conteudo de fsn.m
[y(10), fsn(y,10)] % em matlab, os indices de matrizes/vectores come am em 1

func1(y) % o que e' que esta funcao faz ?

help func1 % talvez isto ajude

% Apenas no matlab
% Descomente a linha seguinte se estiver a usar o matlab
% mupad