O PROCESSAMENTO DIGITAL DE SINAIS

UNIVERSIDADE SÃO FRANCISCO                                                               ENGENHARIA ELÉTRICA             PROCESSAMENTO DIGITAL DE SINAIS LAB. - 01 :  SINAIS DISCRETOS

[pic 1]

        

1- Execute no Matlab, analise e comente o programa no MatLab que visualiza a seguinte sequencia. Pode ser feito em linha de comando :

[pic 2] 

       

       (ao invés do comando plot, utilize o comando stem)

[pic 3]

2- Semelhante ao Exe1, porém criar um arquivo .m no MatLab para visualizar as  seqüências x(n) e x(n-2) e x(n+2). Utilizar o comando subplot para visualizar ao mesmo tempo o gráfico original e o gráfico com a sequencia deslocada. (duvida??? >>help subplot)

 [pic 4]

FUNÇÃO IMPULSO

Geração matemática de uma sequencia de amostragem unitária,

[pic 5] 

[pic 6]

3- Analise o programa MatLab que visualiza a seqüência amostra unitária  δ(n) e δ(n-no) para  n1≤ n ≤n2.

Para isso terá que ser criado a função impseq em um arquivo do Matlab chamado também impseq.m

Após criar a função, executar na linha de comando >>impseq(?,?,?)

[pic 7]

function [x,n]=impseq(n0,n1,n2)

%---------------------------------------

% Gera x(n) = delta(n); n1<= n <= n2

% Gera y(n)=x(n-n0)=delta(n-n0); n1<=n<=n2

%---------------------------------------

n=[n1:n2]; x=[n==0];                    % x(n)=1 somente para n=0

m=[n1:n2]; y=[(m-n0)==0];               % y(n)=1 somente pata n=n0

subplot(2,1,1); stem(n,x);

title('Sequencia unitária Delta(n)'); xlabel('amostra [n]'); ylabel('amplitude [x]');

subplot(2,1,2); stem(m,y);

title('Sequencia unitária Delta(n-no)'); xlabel('amostra [n]'); ylabel('amplitude [y]');

 [pic 8]

4-    a) Execute o programa tal que esboce o impulso unitário δ(n) e δ(n-4).

1. Execute o programa tal que esboce o impulso unitário δ(n) e δ(n+4).

[pic 9]

FUNÇÃO DEGRAU

Representação da seqüência degrau unitário

 [pic 10]

e está relacionada a amostra unitária por   [pic 11]

[pic 12]

5-  Analise o programa MatLab a seguir, criadopara permitir a visualização daa seqüência degrau unitária  u(n) e u(n-no) para  n1≤ n ≤n2 e gere uma função degrau exemplo, de forma similar a função impulso.

[pic 13]

function [x,n]=stepseq(n0,n1,n2)

%---------------------------------------

% Gera x(n) = u(n); n1<= n <= n2

% Gera y(n)=x(n-n0)=u(n-n0); n1<=n<=n2

%---------------------------------------

n=[n1:n2]; x=[n>=0];                   % x(n)=1 somente para n>=0

m=[n1:n2]; y=[(m-n0)>=0];               % y(n)=1 somente pata n>=n0

subplot(2,1,1); stem(n,x);

title('Sequencia degrau unitario u(n)'); xlabel('amostra [n]'); ylabel('amplitude [x]');

subplot(2,1,2); stem(m,y);

title('Sequencia degrau unitario u(n-no)'); xlabel('amostra [n]'); ylabel('amplitude [y]');

[pic 14]

6-     a) Execute o programa tal que esboce a sequencia degrau unitário u(n) e δ(n-4).

1. Execute o programa tal que esboce a sequencia degrau unitário u(n) e δ(n+4).

[pic 15]

7-    Repita os exercícios anteriores para outros valores de deslocamento. Aproveite e varie também o número de amostras.

[pic 16]

Seqüência exponencial real

A função [pic 17], pode ser gerada pelo programa MatLab abaixo. Edite e rode verificando os valores seqüenciais gerados.

[pic 18]

function [x,n] = expseq(n0,n1,n2)

% gera x(n) = a^n ; n real

%--------------------------------------

% [x,n] = expseq(n0,n1,n2)

%

n1 =  0;

n2 = 10;

n = [n1:n2]; x = (0.9).^n;

plot(n(:),x(:), '.') %uma alternativa para o uso do stem

title('LABORATÓRIO DE DSP-Sequência Exponencial');

grid;

xlabel ('n');

ylabel('Sinal Gerado');

[pic 19]

8- Para testar:

[pic 20]

[pic 21]

Geração de uma seqüência exponencial complexa

Seja x(n) dado pela seguinte relação

[pic 22] para um n qualquer.

Onde σ e o fator de amosrtecimento e[pic 23] é a freqüência angular em radianos por segundo. Vamos gerar essa função no intervalo [pic 24] com o,seguinte script do MatLab,