61 lines
970 B
Mathematica
61 lines
970 B
Mathematica
|
%% Inicialização do ambiente
|
||
|
|
clear ; close all; clc
|
||
|
|
|
||
|
|
%% Exercício 10
|
||
|
|
n = 0:250; % janela
|
||
|
|
|
||
|
|
w1 = 0.1*pi;
|
||
|
|
w2 = 0.2*pi;
|
||
|
|
w3 = 0.3*pi;
|
||
|
|
|
||
|
|
x1 = 2*cos(w1*n);
|
||
|
|
x2 = 2*cos(w2*n);
|
||
|
|
x3 = 2*cos(w3*n);
|
||
|
|
|
||
|
|
% Figura com os resultados de x1
|
||
|
|
figure
|
||
|
|
hold on
|
||
|
|
% Plot de x1
|
||
|
|
subplot(2,1,1)
|
||
|
|
stem(n, x1)
|
||
|
|
title('x1')
|
||
|
|
% Plot do output para x1
|
||
|
|
subplot(2,1,2)
|
||
|
|
stem(n, fun_y(x1))
|
||
|
|
title('output para x1')
|
||
|
|
hold off
|
||
|
|
|
||
|
|
% Figura com os resultados de x2
|
||
|
|
figure
|
||
|
|
hold on
|
||
|
|
% Plot de x2
|
||
|
|
subplot(2,1,1)
|
||
|
|
stem(n, x2)
|
||
|
|
title('x2')
|
||
|
|
% Plot do output para x1
|
||
|
|
subplot(2,1,2)
|
||
|
|
stem(n, fun_y(x2))
|
||
|
|
title('output para x2')
|
||
|
|
hold off
|
||
|
|
|
||
|
|
% Figura com os resultados de x3
|
||
|
|
figure
|
||
|
|
hold on
|
||
|
|
% Plot de x3
|
||
|
|
subplot(2,1,1)
|
||
|
|
stem(n, x3)
|
||
|
|
title('x3')
|
||
|
|
% Plot do output para x3
|
||
|
|
subplot(2,1,2)
|
||
|
|
stem(n, fun_y(x3))
|
||
|
|
title('output para x1')
|
||
|
|
hold off
|
||
|
|
|
||
|
|
% Definição do sistema y
|
||
|
|
function y = fun_y(N)
|
||
|
|
x0 = N; %% For x[n]
|
||
|
|
x1 = [N(1, 2:size(N, 2)) 0]; %% For x[n-1]
|
||
|
|
xn1 = [0 N(1, 2:size(N, 2))]; %% For x[n+1]
|
||
|
|
y = x0.^2 - x1 .* xn1;
|
||
|
|
end
|