Descreva como você obteria a resposta em frequência de maneira experimental de um sistema dinâmico LTI.
Discuta as dificuldades ou impedimentos presentes no procedimento.
Discuta as dificuldades ou impedimentos presentes no procedimento.
Passo 1
O enunciado pede duas coisas: (1) descrever, de forma prática e passo a passo, como obter experimentalmente a resposta em frequência de um sistema dinâmico LTI; (2) discutir as dificuldades e impedimentos que aparecem nesse procedimento.
Informações assumidas como dadas (o enunciado não pede para perguntar nada): você tem acesso ao sistema físico que é razoavelmente LTI durante o teste, pode aplicar sinais de entrada (gerador de sinais, actuador) e medir a entrada e a saída com instrumentos de aquisição (ADC, sensores). Também existe um equipamento para gerar sinais digitais/analógicos e software para registrar/processar os sinais.
Informações assumidas como dadas (o enunciado não pede para perguntar nada): você tem acesso ao sistema físico que é razoavelmente LTI durante o teste, pode aplicar sinais de entrada (gerador de sinais, actuador) e medir a entrada e a saída com instrumentos de aquisição (ADC, sensores). Também existe um equipamento para gerar sinais digitais/analógicos e software para registrar/processar os sinais.
Passo 2
O que é exatamente a “resposta em frequência” que vamos medir?
Definição matemática direta que usamos para comparar com o que medimos:
Definição matemática direta que usamos para comparar com o que medimos:
onde
é a transformada de Fourier do sinal de entrada
e
é a transformada de Fourier do sinal de saída
. Aqui
é uma função complexa de frequência: sua magnitude
diz quanto o sistema amplifica/atenua naquela frequência e seu ângulo
diz o deslocamento de fase. No domínio discreto, com DFT/FFT, identificamos índices de frequência
e usamos
.
Passo 3
Panorama de métodos práticos (escolha a que fizer mais sentido em cada laboratório):
- Varredura senoidal (swept sine): aplicar senos individuais em várias frequências, medir amplitude e fase em regime permanente.
- Seno varrente (chirp / swept-sine contínuo): aplicar um sinal cuja frequência varia com o tempo e extrair a resposta por análise tempo-frequência ou por sincronização.
- Excitação broadband + FFT: aplicar ruído branco, PRBS ou multisine e estimar a FRF por razões de espectros (cross- e autospectra).
- Resposta ao impulso/step + FFT: medir resposta ao impulso (ou step e derivar) e computar FFT.
Cada método tem vantagens: swept sine é muito preciso em SNR por frequência; multisine/ruído permite medir todo o espectro de uma vez; impulso/step é direto mas ruidoso.
Passo 4
Procedimento prático detalhado — checklist de preparação (coisas que você sempre faz antes de medir):
- Calibre e verifique sensores/atuadores (ganhos, offset).
- Garanta sincronismo: entrada e saída devem ser medidas com a mesma referência de tempo/clock. Se usar ADCs diferentes, sincronize os clocks.
- Defina a banda de interesse
a (Hz). Relacione com pulsação . - Escolha taxa de amostragem
satisfeita: (Nyquist). Recomenda-se algumas vezes maior (ex.: ) para reduzir erro de fase e facilitar processamento. - Determine resolução em frequência desejada
. A relação fundamental é
onde
é o tempo de aquisição (janelas longas dão melhor resolução). Planeje
de acordo com
.
Passo 5
Como fazer usando varredura senoidal (método clássico, intuitivo, muito usado em prova e em laboratório):
- Escolha uma grade de frequências
cobrindo a banda de interesse (pontos logarítmicos são comuns). - Para cada
: gere com amplitude pequena o suficiente para manter linearidade, grande o suficiente para boa SNR. - Aplique o seno e aguarde o transiente morrer (tempo de acomodação ~ várias constantes de tempo do sistema). Registre apenas o regime permanente.
- Meça a amplitude e fase relativos entre entrada e saída. Formas práticas de medir a fase/amplitude: (a) lock-in / synchronous detection, ou (b) coletar
amostras que contenham um número inteiro de ciclos em e usar DFT para obter o coeficiente complexo. Por exemplo, se usar DFT:
e então
- Repita para todas as
. Construa e .
Observações didáticas: quando você escolhepara cada frequência, tente que contenha um número inteiro de períodos para evitar “leakage” (vazamento espectral). Se não for possível, aplique janela (Hann) e compense ganho da janela.
Passo 6
Como fazer usando excitação broadband e análise por espectros (método eficiente para medir tudo de uma vez):
- Gere um sinal de excitação
com conteúdo nas frequências de interesse: ruído branco filtrado, PRBS, multisine ou chirp. Multisine permite concentrar energia em frequências específicas. PRBS é prático quando só tem um gerador digital. - Registre simultaneamente
e por tempo total (ou em blocos para averaging). - Calcule DFT/FFT para cada bloco e estime espectros: auto-espectro de entrada
e cross-espectro . Uma estimativa robusta usa média entre blocos (Welch): dividir o sinal em blocos, aplicar janela em cada bloco, calcular FFTs, e média dos produtos. - Estimador não-viesado comum da FRF:
onde, por exemplo, para cada bloco
,
é a média de
sobre
. Explique símbolos:
é FFT do bloco
da entrada,
é o complexo conjugado.
5. Calcule a coerência como medida de confiabilidade/linearidade:
5. Calcule a coerência como medida de confiabilidade/linearidade:
Valores de
próximos de
indicam boa relação linear sinal-ruído naquela frequência; valores baixos sugerem ruído, não-linearidade ou má excitação.
Passo 7
Detalhes práticos de processamento (o que realmente salva seu resultado na prova/lab):
- Janela: usar janela (Hann, Hamming) por bloco reduz leakage; lembre-se de compensar o ganho da janela quando comparar amplitude.
- Averaging (Welch): dividir em
blocos e média reduz variância do estimador aproximadamente por . - Resolução vs. tempo: lembre-se que
. Se precisar distinguir picos muito próximos (ex.: modos com baixa separação), T precisa ser grande. - Evite aliasing: coloque filtro anti-alias analógico antes do ADC e escolha
adequado. - Unwrap de fase: a fase obtida do FFT pode estar dentro de
; aplique “phase unwrapping” para obter uma curva contínua antes de interpretar retas de atraso. - Converter magnitude para dB:
. Converter fase para graus: .
Passo 8
Como lidar com atraso e deslocamentos de tempo (prática comum que confunde em prova):
- Um atraso puro
no caminho (cabo, processamento) dá fase linear negativa: contém um termo . Você pode estimar a partir da inclinação da fase e compensar: medir a fase em várias frequências e ajustar por regressão linear. - Se existir um deslocamento de tempo entre os canais de medição (ex.: dois ADCs dessíncronizados), suas estimativas de fase ficam erradas. Solução: sincronizar clocks ou medir e corrigir o atraso entre canais.
Passo 9
Como identificar e tratar não linearidade e ruído (crucial para interpretar se o que você mediu é o
verdadeiro):
- Varie a amplitude do sinal de excitação: se
muda com amplitude, há não-linearidade. Faça pelo menos duas amplitudes (baixa e média) e compare. - Use coerência
: se é baixa numa banda específica, os dados dessa banda não são confiáveis. A causa pode ser ruído alto, excitação fraca naquela banda, ou comportamento não-linear. - Para sistemas com não-linearidade forte, medidas por FRF linear não representam bem o sistema; pode ser necessário identificar um modelo não-linear.
Passo 10
Escolha de método segundo restrições práticas (quando usar o quê):
- Use varredura senoidal quando precisar de maior SNR por frequência e puder gastar tempo (muito precisa).
- Use multisine ou ruído se precisa rapidez e medição de todo espectro de uma vez; ideal para sistemas em que você tem pouco tempo experimental.
- Use chirp/swept-sine em estruturas mecânicas onde a resposta é bem comportada e você quer boa cobertura sem ruído branco.
- Use PRBS para sistemas discretos/controle quando só há gerador digital disponível e você quer um sinal com “branco” aproximado.
Passo 11
Resumo passo-a-passo concreto que eu usaria no laboratório (receita pronta):
- Defino banda
– , e baseado em . - Calibro sensores e sincronizo hardware.
- Escolho tipo de excitação (multisine se quero rápido; sine sweep se quero precisão).
- Gero sinal com amplitude que mantenha linearidade.
- Adiciono tempo de acomodação para transientes; gravo vários blocos (para averaging).
- Computo FFTs por bloco, estimo
e então . - Verifico coerência
; se baixa, investigo causa. - Compenso atraso, aplico janela/averaging, converto magnitude para dB e fase para graus, faço unwrapping.
- Se quiser um modelo, faço ajuste paramétrico nos dados
(ex.: identificação de polos/zeros).
Passo 12
Principais dificuldades e impedimentos (lista comentada — o que tende a ferrar seu experimento na prática):
- Ruído de medição (SNR baixo): sinal fraco ou sensores ruins tornam
ruidoso; solução: aumentar amplitude (sem saturar) e average. - Não-linearidade do sistema: FRF linear não descreve comportamento não-linear; detectada por variação com amplitude ou baixa coerência.
- Leakage/espectro finito: se não usar número inteiro de ciclos nem janela apropriada, o FFT “vaza” energia entre bins. Use janelas e blocos com número inteiro de ciclos quando usar senos.
- Aliasing: sem anti-alias analógico e sem Fs adequado, frequências altas se misturam na banda de interesse.
- Sincronização: canais medidos com clocks diferentes dão erro de fase.
- Tempo de aquisição vs. resolução: para boa resolução
precisa de registros longos (tempo alto). Em prova ou tempo limitado, há trade-off. - Resposta do sensor/atuador: se o sensor ou atuador tem sua própria dinâmica, você mede a cascata e deve compensar ou calibrar.
- Limitação dinâmica do gerador/atuador: não consegue excitar amplitudes/frequências necessárias (ex.: alto f em transmissões mecânicas).
- Sistemas instáveis ou muito amortecidos: resposta pode divergir ou ser muito fraca para medir.
- Variação temporal do sistema: se o sistema muda ao longo do tempo (deriva térmica etc.), médias e pressupostos de LTI falham.
- Restrições de segurança/segurança física: algumas excitações podem danificar o sistema ou serem perigosas.
Passo 13
Dicas práticas de prova (como escrever/argumentar rápido e certo na hora da avaliação):
- Sempre escreva a definição
e explique como obteve e (seno em regime, FFT, cross-spectra). - Mencione unidades e conversões: magnitude em dB
, fase em graus. - Cite coerência
para justificar confiabilidade dos pontos medidos. - Explique trade-offs: tempo longo -> melhor
, mas maior custo experimental; swept-sine -> alto SNR, mas demorado. - Se pedir recomendações, proponha método e parâmetros (ex.: “uso multisine,
, s, média em 8 blocos”).
Resposta
Procedimento (resumo): sincronize hardware, escolha banda e
, selecione método (varredura senoidal para precisão ou excitação broadband + FFT/Welch para rapidez), gere excitação apropriada, registre entrada/saída simultaneamente por tempo
, calcule FFTs, estime
, compute
e
, e verifique coerência
para validar os pontos.
Principais dificuldades/impedimentos: ruído/SNR, não-linearidade, leakage e resolução limitada ( ), aliasing e falta de anti-alias, sincronização de canais (fase), limitações de sensores/atuadores e variação temporal do sistema.
Principais dificuldades/impedimentos: ruído/SNR, não-linearidade, leakage e resolução limitada (