O que é densidade cruzada espectral - CSD?

16

Fiz uma pergunta mais cedo, mas não obtive resposta. Então agora estou simplificando: o que são Densidade Espectral Cruzada (CSD) e Sensibilidade Espectral de Potência (PSD)? Qual é a aplicação deles? Como posso obtê-los no MATLAB?

Skk(ω)=limT1

Skeu(ω)=limT1TE{Yk(ω)Yeu(ω)}
Skk(ω)=limT1TE{Yk(ω)Yk(ω)}

Skeu(ω) é a função de densidade espectral cruzada (CSD) entre os sinais gerais yk(t) e yeu(t) , Skk(ω) é a densidade espectral de potência (PSD) do sinal yk(t) , Yk(ω) é a transformada de Fourier finita do sinal yk(t) na frequência ω , Yk(ω) é o complexo conjugado de Yk(ω) e E{} é o operador de expectativa.


Minha pergunta anterior foi: O que significa 'espectro de potência wavelet', 'espectro de potência automática', 'espectro de potência cruzada' significa na aplicação de wavelet? Eu estava estudando sobre a identificação da forma de modo com o método wavelet e esses termos me confundiram.

Homem elétrico
fonte
Você pode postar uma referência a alguns dos materiais que você estudou? É mais fácil ajudá-lo se você o fizer.
Phonon
@Phonon Oi phanon. Eu edito minha pergunta e publico o link. Vocês podem acessar o jornal ou querem que eu faça o upload em algum lugar? tnx
Electricman

Respostas:

15

Densidade espectral de potência é a distribuição de energia ao longo do eixo de frequência. Geralmente é usado para sinais de energia não finitos (principalmente não limitados em sinais de tempo), que não são quadráveis. O PSD do sinal é a autocorrelação da Transformada de Fourier do sinal, conforme declarado pelo teorema de Wiener-Khinchin. No Matlab:

N = length(S);
F = fft(S);
F = F(1:N/2+1);
PSD = (1/(2*pi*N)) * abs(F).^2;
PSD(2:end-1) = 2*PSD(2:end-1);
freq = 0:(2*pi)/N:pi;

consulte: https://de.mathworks.com/help/signal/ug/power-spectral-density-estimates-using-fft.html

A densidade espectral cruzada é a mesma, mas usando a correlação cruzada, para que você possa encontrar a energia compartilhada por uma determinada frequência para os dois sinais usando seu módulo quadrado e a mudança de fase entre os dois sinais nessa frequência usando seu argumento.

A densidade espectral cruzada pode ser usada para identificar a resposta em frequência de um sistema LTI barulhento: se o ruído não estiver correlacionado à entrada ou saída do sistema, sua resposta em frequência poderá ser encontrada no CSD da entrada e saída.

Florian Castellane
fonte
Obrigado por sua resposta, você também escreveria código matlab para CSD? E você escreveria um exemplo de CSD para identificar a resposta em frequência de um sistema LTI barulhento?
Eletricista
@Electricman O MATLAB Signal Processing Toolbox já tem funções para fazer isso. Em particular, cpsd()faz o que você precisa.
Phonon
@Phonon, acho que usa FFT. como posso executar um CSD por transformada wavelet? Obrigado Phonon
Electricman
@ Eletricman Você deve fazer isso como uma pergunta separada.
Phonon
@Honon, se alguém escrever o código CSD baseado em FFT no matlab. Eu posso fazer a função CSD myself.cpsd () baseada em wavelet não me ajuda. Obrigado cargas
Electricman
5

Para acrescentar à explicação bem descrita acima, no caso de wavelets, que são finitas no tempo, é mais correto não usar o termo "poder", mas "energia". Para Fourier, que tem como base funções o senoide que se estende infinitamente no tempo, densidade espectral de potência é o termo correto. Para wavelets, que tem funções básicas como desvios finitos no tempo, devemos usar 'energia'.

forsker_for_dsp
fonte