A instabilidade torna um sistema LTI não linear (ou variante no tempo)?

7

Estou tirando essa pergunta da pergunta de johnny . Matt L. e eu tivemos conclusões diretamente opostas à pergunta de johnny.

Quero dissociar a questão de questões de causalidade e outras coisas bobas.

Portanto, temos um sistema recursivo simples de primeira ordem descrito com a equação de E / S no domínio do tempo:

y[n]=py[n1] + x[n]nZ

Obviamente, a transformação Z disso é

Y(z)=pz1Y(z) + X(z)

e função de transferência

H(z)Y(z)X(z)=zzp

Normalmente, identificaríamos isso como um sistema de LTI simples e realizável, com zero em e pólo em . Mas na outra questão, existe uma questão relativa à linearidade e invariância no tempo para o caso em que .0pp=1 

Para quais valores p esse sistema é linear? Para quais valores p esse sistema é invariante no tempo?

Creio que este é o núcleo do desacordo que tenho com o Dr. Matt L.

Robert Bristow-Johnson
fonte
Se o expandirmos para sua forma inicial. Teoricamente, o sistema é sempre linear.
aluno
ok @learner, então se , então MattL sugere que você pode ter para a entrada na saída incluindo diferente de zero . a entrada é identicamente zero, a saída é diferente de zero. esse sistema pode ser linear? p=1x[n]=0 n
y[n]=A(1)nAR
A
Robert Bristow-johnson
11
Esclareça: para quais números inteiros o relacionamento válido? Para todos os números inteiros ? Para ? Para ? E nos dois últimos casos, especifique qual é o valor de (respectivamente ), pois concordamos que ou . n
y[n]=py[n1] + x[n]
nn0n>0y[1]y[0]y[0]=py[1]+x[0]y[1]=py[0]+x[1]
precisa saber é o seguinte
11
Eu nunca tinha colocado uma restrição sobre isso. na outra questão, eu estava explícito que . Eu deveria ter dito isso aqui, mas eu negligenciei. então eu vou modificar a pergunta. na minha opinião, se não for especificado, a suposição é que não há restrição. nZ
Robert bristow-johnson
e como você pode derivar ou qualquer outro específico disso. se você assume causalidade, recebe um lado direito . se você assume anti-causalidade (que não é a maneira como afirmei a equação da diferença), você obtém um lado esquerdo . nos dois casos, isso afeta um limite ou outro para o somatório da convolução.
y[n]=i=+h[i]x[ni]
y[0]y[n]h[n]h[n]
Robert Bristow-johnson

Respostas:

3

Não sei como a discussão chegou a esse ponto, mas é bastante complicado seguir esse período de tempo. É uma equação de diferença comum com coeficientes constantes, portanto, define um sistema linear, invariante no tempo. Não há necessidade de avançar mais se houver uma solução única . Aqui temos um problema com a parte enfatizada.

Vamos primeiro escrever as equações do espaço de estado do descritor: o sistema é descrito por

[1100][s1[k+1]s2[k+1]]=[p000][s1[k]s2[k]]+[01]x[k]y[k]=[10][s1[k]s2[k]]+[1]x[k]

Agora é aqui que acho que é a parte problemática. Esse sistema, apesar de ser LTI, não é regular (as palavras-chave relevantes são sistemas descritores regulares , isentos de impulso e de índice 1). Em outras palavras, não existe para o qual a expressão é diferente de zero e, portanto, um dos modos é e, de fato, é para este exemplo. Isso significa que nosso sistema possui problemas de exclusividade da solução. Ao contrário do sistema causal de LTI, não há garantia da existência de uma solução exclusiva. Não há garantia de uma solução admissível para esse assunto (sistemas impulsivos de chavão). Portanto, o raciocínio LTI das outras respostas não será suficiente.λdet(λEA)000

E é isso que causa o problema, até onde posso dizer pelo argumento de Matt L, que ele encontrou duas soluções não triviais distintas para o mesmo sistema e concluiu que esse não pode ser um sistema linear. Mas isso também pressupõe a singularidade e a existência de uma solução e condições iniciais.

Ele difere apenas dos sistemas regulares na maneira em que a exclusividade e a existência de garantias dos sistemas LTI padrão não podem ser assumidas. Os modelos não podem mais ser considerados como tendo trajetórias admissíveis para todos os sinais possíveis.

percusse
fonte
1

A resposta é não. Para rotular um sistema de tempo discreto como LTI, procuramos apenas suas propriedades de linearidade e invariância do tempo e não precisamos nos preocupar se é estável ou não. Essa é outra propriedade independente de um sistema que pode coexistir mutuamente com outras propriedades. E, de fato, muitos sistemas de LTI são instáveis ​​e ainda são sistemas de LTI. Para uma abundância de exemplos, consulte o livro de Alan Oppenheim: Signals & Systems, 2ed, capítulo 2. (ou qualquer outro livro de faculdade sobre sinais e sistemas, ou processamento de sinal digital). Considere, por exemplo, filtros IIR não estáveis ​​que ainda são lineares e tempo invariável. (de fato, seu exemplo é um desses)

Chegando ao seu LCDDE que supostamente define um sistema de tempo discreto recursivo, como você deve saber, o próprio LCDDE não é suficiente para especificar exclusivamente uma solução, pois você também precisa de um conjunto de condições auxiliares (condições iniciais). Sem essas condições iniciais definidas explicitamente, você não pode resolver a equação nem determinar se o sistema que representa será LTI ou causal. Porque, para algumas condições iniciais, pode não ser causal, não linear e variar no tempo, enquanto para outro conjunto (ou seja, as condições iniciais de descanso) será linear, invariante no tempo e causal. Portanto, para que um único LCCDE represente exclusivamente um sistema LTI, suas condições iniciais devem ser definidas corretamente como repouso inicial e não arbitrariamente ...

Fat32
fonte
0

Linearidade e invariância no tempo não dependem do valor de . Os dois possíveis sistemas lineares e invariantes no tempo (LTI) descritos pela equação da diferençap

(1)y[n]=py[n1]+x[n],p0

são dados pela transformação inversa da função de transferência, conforme formulado na pergunta:Z

(2)H(z)=zzp

Observe que (2) não define exclusivamente uma resposta de impulso, a menos que a região de convergência (ROC) de seja fornecida. Para a função de transferência (2), existem dois ROCs possíveis:e. No primeiro caso, a resposta de impulso correspondente é do lado direito e corresponde a um sistema causal de LTI:H(z)|z|>|p||z|<|p|

(3)h1[n]=pnu[n]

onde é a função de etapa da unidade. Se escolhermoscomo o ROC, obtemos uma resposta de impulso do lado esquerdo correspondente a um sistema anti-causal:u[n]|z|<|p|

(4)h2[n]=pnu[n1]

Se , o sistema LTI causal caracterizado por é estável, e o sistema LTI anti-causal caracterizado por é instável e, se o oposto se mantém.|p|<1h1[n]h2[n]|p|1

Até agora, vimos que a equação da diferença (1) define dois sistemas de LTI, um causal e o outro anti-causal. As duas respostas de impulso correspondentes e são soluções da equação da diferença (1) para . No entanto, essas não são as únicas soluções de (1). É bem conhecido da teoria das equações de diferença linear que a solução geral é dada por uma solução específica para alguns , e por uma solução para a equação homogênea definida por . Para a equação da diferença dada, a equação homogênea correspondente éh1[n]h2[n]x[n]=δ[n]x[n]x[n]=0

(5)y[n]=py[n1]

A solução de (5) é

(6)yh[n]=cpn,cR

Agora podemos expressar a solução geral de (1) (com ) combinando uma solução específica ( ou ) com :x[n]=δ[n]h1[n]h2[n]yh[n]

(7)y[n]=h1[n]+yh[n]=h1[n]+cpn

Observe que a solução pode ser obtida em (7) escolhendo . Observe também que (7) é válido para todos os .h2[n]c=1nZ

Se escalássemos o sinal de entrada e com algum , a saída resultante seriax1[n]=aδ[n]aR

(8)y1[n]=ah1[n]+cpn

Esse sinal de saída geralmente não é igual a (com fornecido por (7)), a menos que ou . Consequentemente, o sistema correspondente não é linear. Para o sinal de entrada a saída éay[n]y[n]c=0c=1x2=δ[nn0]

(9)y2[n]=h1[nn0]+cpn

que geralmente não é igual a (novamente, a menos que ou ). Portanto, o sistema correspondente também não é invariante no tempo.y[nn0]c=0c=1

Em suma, a equação da diferença (1) descreve infinitamente muitos sistemas com respostas ao sinal de entrada dado por (7). Apenas dois desses sistemas são LTI, os outros não. Os dois sistemas de LTI são descritos pelas respostas de impulso e fornecidas por (3) e (4), respectivamente.x[n]=δ[n]h1[n]h2[n]

Matt L.
fonte
então, basicamente, o sistema descrito por esta equação não é LTI. ou não necessariamente LTI. direita? agora, para simplificar a discussão, você pode concordar que o motivo de afirmar que não é LTI é porque com uma entrada que é identicamente zero , que o a saída pode ser diferente de zero e ?
y[n]=py[n1] + x[n]nZ
 
x[n]=0 nZ
y[n]=cpn nZ
c0 
Robert Bristow-johnson
Além disso, se , isso poderia ser uma exceção à sua afirmação inicial de que "Linearidade e invariância no tempo não dependem do valor de " ? p=0p
Robert Bristow-johnson
e, finalmente, você está aparentemente negando que, mesmo para um escolhido corretamente , que não pode ser uma solução geral para ou sistemas descritos de maneira semelhante? h[n]
y[n]=i=+h[i]x[ni]
y[n]=py[n1] + x[n]nZ
Robert Bristow-johnson
MattL, você está interessado em abordar alguma dessas perguntas?
22675 Robert De Bristow-johnson #
@ robertbristow-johnson: Sim, mas você parece fazer muitas perguntas ao mesmo tempo, e eu sinto que a maioria delas já foi respondida explicitamente pela minha resposta à sua pergunta original. Mas vou dar uma chance ...
Matt L.
0

Acho essa discussão intrigante e queria acrescentar outro ponto de vista à mistura:

O sistema em consideração ( y[n]=py[n1]+x[n]) pode ser pensado como um mapeamento de um espaço vetorial (dimensional infinito) para outro. Vamos chamar isso de mapeamentoM, e podemos (inicialmente) defini-lo como:

M:RZRZ

Essa terminologia diz que M é um mapeamento de RZ (o espaço de todas as funções com valor real de uma variável inteira) para RZ.

Se o sistema tiver zeros (e o sistema em consideração aqui tiver zero em z=1), isso significa que nosso mapeamento Mnão é um para um, porque dois sinais de entrada diferentes levarão ao mesmo sinal de saída. Por exemplo, para qualquer sinal de entrada,x[n], Nós podemos dizer que M(x)=M(x+λ) para qualquer real λ.

O conjunto de funções que são "zeros" do nosso sistema pode ser definido como:

Kzeros={f[n]=λ:λR}

Da mesma forma, observamos que, se nosso sistema tiver pólos (e o sistema em consideração aqui tiver zero em z=1), isso significa que o mapeamento inverso, M1não é um para um. Especificamente,M1(x)=M1(x+λ(1)n) para qualquer real λ.

O conjunto de funções que são "pólos" do nosso sistema pode ser definido como:

Kpoles={f[n]=λ(1)n:λR}

Agora, RZ é um espaço vetorial Kzeros é um espaço vetorial e Kpoles é um espaço vetorial.

Agora podemos definir dois espaços quocientes (consulte a Wikipedia para obter mais informações sobre espaços quocientes):

Qinput=RZ/Kzeros

Qoutput=RZ/Kpoles

Você pode pensar em Qoutput como sendo o subconjunto de RZ que não contém nenhum componente de sinal no formulário λ(1)nou, alternativamente, você pode pensar em Qoutput como sendo idêntico a RZcom classes de equivalência que nos dizem "para nossa aplicação atual, consideraremos qualquer funçãoy[n]ser equivalente ay[n]+λ(1)n para qualquer real λ"

Ao fazer isso, agora podemos redefinir um novo mapeamento M como um mapeamento de Qinput para Qoutput. Esse novo mapeamento é realmente igual ao nosso antigo mapeamento,M, exceto que reduzimos os espaços vetoriais nos quais ele opera. Além disso, esse novo mapeamento agora é uma bijeção (é "um para um" e "para"), portanto, é garantido que também seja invertível.

Finalmente, esse mapeamento, Mé linear .

Portanto, o objetivo de toda essa explicação desmedida é que, definindo as classes de equivalência apropriadas (ou alternativamente, restringindo nosso espaço de funções permitidas a um subespaço de RZ), podemos manter a propriedade de que nosso mapeamento deve ser linear (e invariável no tempo).

Por exemplo, as regras de linearidade nos dizem que, sex[n] é um sinal de entrada e α é escalar real, então M(αx)=αM(x). Portanto, isso implica que, ao definirα=0, devemos, portanto, esperar que M(0×x)=y[n]=0 (ou seja, se introduzirmos o sinal zero em nosso filtro, a saída deve ser y[n]=0)

No entanto, sabemos que é possível ter uma situação em que a entrada para o filtro seja zero, mas a saída seja da forma y[n]=(1)n, então podemos ficar tentados a dizer "que prova que nosso sistema não é linear, porque y[n]é zero ". No entanto, você deve se lembrar que a classe de equivalência que aplicamos no espaço vetorial de saída diz que" para nosso aplicativo atual, consideraremos qualquer funçãoy[n]ser equivalente ay[n]+λ(1)n para qualquer real λ", o que significa que y[n]=(1)né equivalente a zero!

dave_mc
fonte