Processamento digital de sinais

...

C^(..) + 4c + 3c = 3R

X1 = c

X2 = c^.

dx1/dt = dc/dt

X^(..) 1 = c^.

X^(..) = x²

dx2/dt = (d c^. )/dt

X^. 2 = c^(..)

X^.2 = -4C – 3C = 3 R

X2 = -4x2 – 3x1 + 3R

X^. = AX+ Bu

X^. 1 0 1 X1 0 × R

X^. 2 -4 -3 X2 3

X1

X2

Y = [1 0]

Filtros FIR

Exemplos:

Filtro de interpolação digital

Projetar um filtro de interpolação digital para upsample um sinal por sete, usando o método bandlimited. Especifique um fator "bandlimitedness" de 0,5 e use amostras na interpolação..

upfac = 7;

alpha = 0.5;

h1 = intfilt(upfac,2,alpha);

O filtro funciona melhor quando o sinal original é bandlimited para alfa vezes a freqüência de Nyquist. Criar um sinal de ruído bandlimited gerando 200 números aleatórios gaussianos e filtrando a seqüência com um filtro lowpass FIR de ordem 40. Redefinir o gerador de números aleatórios para obter resultados reprodutíveis.

Lowp = fir1 (40, alfa);

Rng (’default’)

X = filtro (lowp, 1, randn (200, 1));

Aumentar a taxa de amostragem do sinal inserindo zeros entre cada par de amostras de x.

Xr = upsample (x, upfac);

Use a função de filtro para produzir um sinal interpolado.

Y = filtro (h1,1, xr);

Compensar o atraso introduzido pelo filtro. Trace os sinais originais e interpolados.

Atraso = média (grpdelay (h1));

Y (1: atraso) = [];

Haste (1: upfac: upfac * comprimento (x), x)

aguente

Parcela (y)

Xlim ([400 700])

Intfilt também executa a interpolação polinomial de Lagrange.

Documentation

cfirpm

Complex and nonlinear-phase equiripple FIR filter design

Projeto de filtro FIR com e sem filtro de fase complexa e não-linear

Descrição

O cfirpm permite que restrições arbitrárias de domínio de frequência sejam especificadas para o projeto de um filtro FIR possivelmente complexo. O erro de filtro Chebyshev (ou minimax) é otimizado, produzindo projetos de filtro FIR equiripple.

As funções de resposta de frequência de fresp predefinidas estão incluídas para um número de projetos de filtro comuns, conforme descrito abaixo. Para obter todas as funções predefinidas de resposta em frequência, a opção de simetria ’sym’ assume como padrão ’even’ se nenhuma freqüência negativa estiver contida em f e D = 0; Caso contrário, ’sym’ assume como padrão ’none’. Para todas as funções de resposta de frequência predefinidas, dspecifies um deslocamento de atraso de grupo de tal forma que a resposta do filtro tenha um atraso de grupo de n / 2 + d em unidades do intervalo de amostra. Os valores negativos criam menos atraso; Valores positivos criam mais atraso. Por padrão d = 0:

Nota: Para filtros @bandpass, o primeiro elemento no vetor de freqüência deve ser menor ou igual a zero eo último elemento deve ser maior ou igual a zero.

Exemplos

recolher todos

Filtro Passa-vazio Equiripple

Projetar um filtro lowpass de fase linear de 31 derivações. Mostre suas respostas de magnitude e fase.

B = cfirpm (30, [- 1 -0,5-0,4 0,7 0,8 1], @ passagem baixa);

FIR Aproximação para Allpass resposta

Desenhe um filtro FIR allpass de fase não linear de ordem 22 com resposta em frequência dada aproximadamente por, onde.

N = 22; % Ordem do filtro

F = [-1 1]; % Limites da banda de frequência

W = [1 1]; % Pesos para otimização

Gf = linspace (-1,1,256); Grade de pontos de freqüência

D = exp (-1i * pi * gf * n / 2 + 1i * pi * pi * sinal (gf). * Gf. * Gf * (4 / pi));

Resposta de freqüência desejada

Use cfirpm para calcular o filtro FIR. Trace as respostas de magnitude real e aproximada em dB e as respostas de fase em graus.

B = cfirpm (n, f, ’allpass’, w, ’real’); % Aproximação

Freqz (b, 1,256, ’inteiro’)

Subplot (2,1,1)% Resposta de sobreposição

aguente

Plot (pi * (gf + 1), 20 * log10 (abs (fftshift (d))), ’r--’)

Subtrama (2,1,2)

aguente

Plot (pi * (gf + 1), unwrap (ângulo (fftshift (d))) * 180 / pi, ’r--’)

Lenda (’Aproximação’, ’Desejado’, ’Localização’, ’Sudoeste’)

Introduced before R2006a

Processamento Digital de Sinais

segunda-feira,