Oitava Média Móvel

Depois de juntar os bits desse segmento, criei esta função usando a função de filtro Octaves. Começa com a média móvel simples como base. V é o vetor de coluna de números para calcular a média móvel exponencial. A janela é um número inteiro como um número de dias. Eu usei 12. Aqui está uma explicação matemática dessa função. Observe que a página usa 2 (n1) (onde n é janela ou o número de dias) como alfa. Mas eu uso 1n porque esse valor de alfa se adapta às minhas necessidades. Ajuste o alfa conforme necessário. Alternativamente, às vezes eu preciso das dimensões dos vetores de entrada e saída para combinar. Eu preencho valores inválidos com NaN adicionando meanV NaN (window-1,1) meanV como a última linha na função movingEMean. Você também pode preenchê-lo com SimpleAvg se desejar uma estimativa aproximada.31 Processamento do sinal Este capítulo descreve o processamento do sinal e as funções rápidas de transformação de Fourier disponíveis na Octave. As transformações rápidas de Fourier são calculadas com as bibliotecas FFTW ou FFTPACK dependendo de como a Octave é construída. Calcule a transformada discreta de Fourier de A usando um algoritmo de Transformação de Fourier Rápida (FFT). O FFT é calculado ao longo da primeira dimensão não-singleton da matriz. Assim, se x é uma matriz, fft (x) calcula a FFT para cada coluna de x. Se for chamado com dois argumentos, espera-se que n seja um inteiro especificando o número de elementos de x para usar ou uma matriz vazia para especificar que seu valor deve ser ignorado. Se n for maior que a dimensão ao longo da qual o FFT é calculado, então x é redimensionado e preenchido com zeros. Caso contrário, se n for menor do que a dimensão ao longo da qual o FFT é calculado, então x é truncado. Se for chamado com três argumentos, dim é um número inteiro que especifica a dimensão da matriz ao longo do qual o FFT é executado. Calcule a transformada de Fourier inversa discreta de A usando um algoritmo de Transformação de Fourier Rápida (FFT). O FFT inverso é calculado ao longo da primeira dimensão não-singleton da matriz. Assim, se x é uma matriz, fft (x) calcula a FFT inversa para cada coluna de x. Se for chamado com dois argumentos, espera-se que n seja um inteiro especificando o número de elementos de x para usar ou uma matriz vazia para especificar que seu valor deve ser ignorado. Se n for maior que a dimensão ao longo da qual o FFT inverso é calculado, então x é redimensionado e preenchido com zeros. Caso contrário, se n for menor do que a dimensão ao longo da qual a FFT inversa é calculada, então x é truncado. Se for chamado com três argumentos, dim é um número inteiro que especifica a dimensão da matriz ao longo da qual a FFT inversa é realizada. Calcule a transformada de Fourier discreta bidimensional de A usando um algoritmo de Transformação de Fourier Rápida (FFT). Os argumentos opcionais m e n podem ser usados ​​especifica o número de linhas e colunas de A para usar. Se qualquer um destes for maior do que o tamanho de A. A é redimensionado e preenchido com zeros. Se A é uma matriz multidimensional, cada sub-matriz bidimensional de A é tratada separadamente. Calcule a transformada de Fourier discreta bidimensional inversa de A usando um algoritmo de Transformação de Fourier Rápida (FFT). Os argumentos opcionais m e n podem ser usados ​​especifica o número de linhas e colunas de A para usar. Se qualquer um destes for maior do que o tamanho de A. A é redimensionado e preenchido com zeros. Se A é uma matriz multidimensional, cada sub-matriz bidimensional de A é tratada separadamente. Calcule a transformada de Fourier discreta N-dimensional de A usando um algoritmo de Transformação de Fourier Rápida (FFT). O tamanho opcional do argumento do vetor pode ser usado, especifique as dimensões da matriz a ser usada. Se um elemento de tamanho for menor que a dimensão correspondente de A., a dimensão de A é truncada antes da execução da FFT. Caso contrário, se um elemento de tamanho for maior do que a dimensão correspondente, então A é redimensionada e preenchida com zeros. Calcule a transformada de Fourier discreta N-dimensional inversa de A usando um algoritmo de Transformação de Fourier Rápida (FFT). O tamanho opcional do argumento do vetor pode ser usado, especifique as dimensões da matriz a ser usada. Se um elemento de tamanho for menor que a dimensão correspondente de A., a dimensão de A é truncada antes de executar a FFT inversa. Caso contrário, se um elemento de tamanho for maior do que a dimensão correspondente, então A é redimensionada e preenchida com zeros. Octave usa as bibliotecas FFTW para executar cálculos FFT. Quando a Octave inicia e inicializa as bibliotecas FFTW, elas lêem um arquivo do sistema inteiro (em um sistema Unix, geralmente é um erro) que contém informações úteis para acelerar os cálculos FFT. Esta informação é chamada de sabedoria. O arquivo do sistema permite que a sabedoria seja compartilhada entre todos os aplicativos usando as bibliotecas do FFTW. Use a função fftw para gerar e salvar a sabedoria. Usando os utilitários fornecidos juntamente com as bibliotecas FFTW (fftw-sabedoria em sistemas Unix), você pode até adicionar sabedoria gerada pela Octave para o arquivo de sabedoria do sistema. Gerencie os dados de sabedoria do FFTW. Os dados de sabedoria podem ser usados ​​para acelerar significativamente o cálculo das FFT, mas implica um custo inicial em seu cálculo. Quando as bibliotecas do FFTW são inicializadas, elas lêem um arquivo de sabedoria do sistema (tipicamente em etc), permitindo que a sabedoria seja compartilhada entre outras aplicações além da Octava. Alternativamente, a função fftw pode ser usada para importar sabedoria. Por exemplo, salvará o sabedoria existente usado pela Octave para o conhecimento da cadeia. Esta string pode então ser guardada em um arquivo e restaurada usando os comandos save e load, respectivamente. Essa sabedoria existente pode ser reimportada da seguinte forma. Se a sabedoria é uma corda vazia, então a sabedoria usada é desmarcada. Durante o cálculo das transformações de Fourier é gerada mais sabedoria. A moda em que essa sabedoria é gerada também é controlada pela função fftw. Existem cinco modos diferentes em que a sabedoria pode ser tratada: Especifica que nenhuma medição de tempo de execução do meio ótimo de cálculo de um particular é realizada e uma heurística simples é usada para escolher um plano (provavelmente sub-ótimo). A vantagem deste método é que há pouca ou nenhuma sobrecarga na geração do plano, o que é apropriado para uma transformada de Fourier que será calculada uma vez. Neste caso, uma gama de algoritmos para executar a transformação é considerada e a melhor é selecionada com base em seu tempo de execução. Semelhante a quotmeasurequot. Mas uma gama mais ampla de algoritmos é considerada. Como quotmeasurequot. Mas todos os algoritmos possíveis que podem ser usados ​​para tratar a transformação são considerados. Como a medição de tempo de execução do algoritmo pode ser cara, isso é um compromisso em que quotmeasurequot é usado para transformações até o tamanho de 8192 e além do que o método quotestimatequot é usado. O método padrão é quotestimatequot. O método atual pode ser consultado ou configurado usando Nota que a sabedoria calculada será perdida ao reiniciar o Octave. No entanto, os dados de sabedoria podem ser recarregados se ele for salvo em um arquivo como descrito acima. Os arquivos de sabedoria salva não devem ser usados ​​em diferentes plataformas, pois não serão eficientes e o ponto de cálculo da sabedoria é perdido. O número de threads utilizados para calcular os planos e executar as transformações pode ser definido com Observação: a oitava deve ser compilada com o suporte multi-threaded FFTW para esse recurso. O número de processadores disponíveis para o processo atual é usado por padrão. Convocar dois vetores usando o FFT para computação. C fftconv (x. Y) retorna um vetor de comprimento igual ao comprimento (x) comprimento (y) - 1. Se x e y são os vetores de coeficientes de dois polinômios, o valor retornado é o vetor de coeficientes do polinômio do produto. A computação usa a FFT chamando a função fftfilt. Se o argumento opcional n for especificado, um FFT de N-ponto é usado. Filtre x com o filtro FIR b usando a FFT. Se x for uma matriz, filtre cada coluna da matriz. Dado o terceiro argumento opcional, n. Fftfilt usa o método de sobreposição-adicionar para filtrar x com b usando uma FFT de N-ponto. O tamanho FFT deve ser uma potência igual de 2 e deve ser maior ou igual ao comprimento de b. Se o n especificado não atender a esses critérios, ele é ajustado automaticamente para o valor mais próximo que o faz. Aplique um filtro digital 1-D aos dados x. O filtro retorna a solução para a seguinte equação de diferença linear, tempo-invariante: onde Nlength (a) -1 e Mlength (b) -1. O resultado é calculado sobre a primeira dimensão não-singleton de x ou sobre dim se fornecida. Uma forma equivalente da equação é: onde c aa (1) e d ba (1). Se o quarto argumento si for fornecido, ele é tomado como o estado inicial do sistema e o estado final é retornado como sf. O vetor de estado é um vetor de coluna cujo comprimento é igual ao comprimento do vetor de coeficiente mais longo menos um. Se não for fornecido, o vetor de estado inicial é definido para todos os zeros. Em termos da Transformada Z, y é o resultado de passar o sinal de tempo discreto x através de um sistema caracterizado pela seguinte função de sistema racional: Aplique o filtro FIR 2-D b para x. Se a forma do argumento for especificada, devolva uma matriz da forma desejada. Os valores possíveis são: pad x com zeros em todos os lados antes de filtrar. Unpadded x (padrão) trim x após a filtragem para que os efeitos de borda não estejam incluídos. Observe que esta é apenas uma variação na convolução, com os parâmetros invertidos e b girados 180 graus. Retornar a resposta de freqüência complexa h do filtro racional IIR cujos coeficientes de numerador e denominador são b e a. respectivamente. A resposta é avaliada em n freqüências angulares entre 0 e 2pi. O valor de saída w é um vetor das freqüências. Se um é omitido, o denominador é assumido como sendo 1 (isso corresponde a um filtro FIR simples). Se n for omitido, é assumido um valor de 512. Para a computação mais rápida, n deve influenciar em um pequeno número de primos pequenos. Se o quarto argumento, quotwholeot. É omitido, a resposta é avaliada em freqüências entre 0 e pi. Avalie a resposta nas freqüências específicas no vetor w. Os valores para w são medidos em radianos. Hellip freqz (hellip, Fs) Frequências de retorno em Hz em vez de radianos assumindo uma taxa de amostragem Fs. Se você estiver avaliando a resposta em freqüências específicas w. Essas frequências devem ser solicitadas em Hz, em vez de radianos. Trace a magnitude e a resposta de fase de h em vez de devolvê-las. Trace a magnitude e a resposta de fase de h. Se o argumento freqnorm opcional for verdadeiro, o vetor de freqüência w está em unidades de radianos normalizados. Se Freqnorm for falso, ou não dado, então w é medido em Hertz. Calcule a função sinc. Retornar sin (pix) (pix). Desligue as fases do radian adicionando múltiplos de 2pi conforme apropriado para remover saltos maiores do que tol. Tol padrão para pi. Unwrap funcionará ao longo da dimensão dim. Se dim é não especificado, ele padrão é a primeira dimensão não-singleton. Ajustar um modelo de regressão ARCH às séries temporais e usando o algoritmo de pontuação no documento ARCH original da Englersquos. Em que e (t) é N (0, h (t)). Dado um vetor de séries temporais e ao tempo t-1 e uma matriz de regressores (comuns) x até t. A ordem da regressão da variância residual é especificada por p. Se invocado como archfit (y. K. P) com um inteiro positivo k. Ajuste um processo ARCH (k. P), ou seja, faça o acima com a t-ésima linha de x dada opcionalmente, pode-se especificar o número de iterações iter. O fator de atualização gamma. E os valores iniciais a0 e b0 para o algoritmo de pontuação. Simule uma sequência ARCH do comprimento t com coeficientes AR b e coeficientes CH a. O resultado y (t) segue o modelo onde e (t). Dado o tempo t-1. É N (0, h (t)). Com Para um modelo de regressão linear, execute um teste de Lagrange Multiplicador (LM) da hipótese nula de nenhuma heteroscedesa autônoma condicional contra a alternativa de CH (p). Isto é, O modelo é dado até t-1 e x até t. E (t) é N (0, h (t)) com e nulo é (1) hellip a (p) 0. Se o segundo argumento for um inteiro escalar, k. Realize o mesmo teste em um modelo de auto-agressão linear de ordem k. Isto é, com a t-ésima fileira de x. Sob o nulo, LM tem aproximadamente uma distribuição chisquare com p graus de liberdade e pval é o valor p (1 menos o CDF desta distribuição no LM) do teste. Se nenhum argumento de saída for dado, o valor de p será exibido. Retornar uma simulação do modelo ARMA. O modelo ARMA é definido pelo qual k é o comprimento do vetor a. L é o comprimento do vetor b e e é o ruído branco gaussiano com variância v. A função retorna um vetor de comprimento t. O parâmetro opcional n dá o número de simbólico x (i) usado para inicialização, ou seja, uma sequência de comprimento t n é gerada e x (n 1: t n) é retornado. Se n for omitido, n 100 é usado. Dado uma série de tempo (vetor) y. Devolver uma matriz com os na primeira coluna e os primeiros k valores de lagged de y nas outras colunas. Em outras palavras, para t gt k. 1, y (t -1), hellip, y (t-k) é a t-ésima linha do resultado. A matriz resultante pode ser usada como uma matriz de regressores em autoregressões. Retornar os coeficientes de filtro de uma janela Bartlett (triangular) de comprimento m. Para uma definição da janela Bartlett, veja, por exemplo, A. V. Oppenheim amp. R. W. Schafer, processamento de sinal de tempo discreto. Retornar os coeficientes de filtro de uma janela Blackman de comprimento m. Se o argumento opcional quotperiodicquot for dado, a forma periódica da janela será retornada. Isto é equivalente à janela do comprimento m 1 com o último coeficiente removido. O argumento opcional quotsymmetricquot é equivalente a não especificar um segundo argumento. Para uma definição da janela Blackman, veja, p. Ex. A. V. Oppenheim amp. R. W. Schafer, processamento de sinal de tempo discreto. Se x é um vetor, detrend (x. P) remove o melhor ajuste de um polinômio de ordem p a partir do dado x. Se x é uma matriz, detrend (x. P) faz o mesmo para cada coluna em x. O segundo argumento p é opcional. Se não for especificado, um valor de 1 é assumido. Isso corresponde à remoção de uma tendência linear. A ordem do polinômio também pode ser dada como uma string, caso em que p deve ser quotconstantquot (corresponde a p 0) ou quotlinearquot (corresponde a p 1). Retornar o estimador d para o parâmetro de diferenciação de uma série de tempo integrada. As frequências de 2piat, 2pibT são utilizadas para a estimativa. Se b for omitido, o intervalo 2piT, 2piaT é usado. Se ambos b e a são omitidos, usa-se 0,5 sqrt (T) e b 1,5 sqrt (T), onde T é o tamanho da amostra. Se x for uma matriz, estima-se o parâmetro de diferenciação de cada coluna. Os estimadores para todas as freqüências nos intervalos descritos acima são retornados em dd. O valor de d é simplesmente a média de dd. Referência: P. J. Brockwell ampère R. A. Davis. Série temporal: teoria e métodos. Springer 1987. Execute um passo do algoritmo Durbin-Levinson. O vetor c especifica as autocovariâncias gamma0, hellip, gammat do intervalo 0 a t. Oldphi especifica os coeficientes com base em c (t -1) e oldv especifica o erro correspondente. Se oldphi e oldv forem omitidos, todos os passos de 1 a t do algoritmo são executados. Execute uma mudança do vetor x. Para uso com as funções fft e ifft, para mover a freqüência 0 para o centro do vetor ou matriz. Se x é um vetor de N elementos correspondentes a N amostras de tempo espaçadas por dt. Então fftshift (fft (x)) corresponde a freqüências Se x é uma matriz, o mesmo é válido para linhas e colunas. Se x for uma matriz, então o mesmo se mantém ao longo de cada dimensão. O argumento dim opcional pode ser usado para limitar a dimensão ao longo da qual a permutação ocorre. Desfazer a ação da função fftshift. Para comprimento uniforme x. Fftshift é o seu próprio inverso, mas os comprimentos estranhos diferem ligeiramente. Calcule as diferenças fracionárias (1-L) d x onde L denota o lag-operator e d é maior que -1. Retornar os coeficientes de filtro de uma janela Hamming de comprimento m. Se o argumento opcional quotperiodicquot for dado, a forma periódica da janela será retornada. Isto é equivalente à janela do comprimento m 1 com o último coeficiente removido. O argumento opcional quotsymmetricquot é equivalente a não especificar um segundo argumento. Para uma definição da janela de Hamming, veja, por exemplo, A. V. Oppenheim amp. R. W. Schafer, processamento de sinal de tempo discreto. Retornar os coeficientes de filtro de uma janela Hanning de comprimento m. Se o argumento opcional quotperiodicquot for dado, a forma periódica da janela será retornada. Isto é equivalente à janela do comprimento m 1 com o último coeficiente removido. O argumento opcional quotsymmetricquot é equivalente a não especificar um segundo argumento. Para uma definição da janela Hanning, veja, por exemplo, A. V. Oppenheim amp. R. W. Schafer, processamento de sinal de tempo discreto. Estimar o parâmetro Hurst da amostra x através da estatística do intervalo redimensionado. Se x for uma matriz, o parâmetro é estimado para cada coluna. Retorne o polinômio de interpolação de Hermite Cubic Piecewise (pchip) dos pontos x e y. Se for chamado com dois argumentos, retorne o polinômio por partes pp que pode ser usado com ppval para avaliar o polinômio em pontos específicos. Quando chamado com um terceiro argumento de entrada, pchip avalia o polinômio pchip nos pontos xi. O terceiro formulário de chamada é equivalente a ppval (pchip (x. Y), xi). A variável x deve ser um vetor estritamente monotônico (aumentando ou diminuindo) do comprimento n. Y pode ser um vetor ou matriz. Se y for um vetor, ele deve ter o mesmo comprimento n como x. Se y for uma matriz, então o tamanho de y deve ter o formulário s1. S2. Hellip, sk. N A matriz é remodelada internamente para uma matriz onde a dimensão principal é dada por s1 s2 hellip sk e cada linha dessa matriz é tratada separadamente. Observe que isso é exatamente o oposto do interp1, mas é feito para compatibilidade MATLAB. Retornar o periodograma (Power Spectral Density) de x. As possíveis entradas são: vetor de dados. Se x for real-valorado, um espectro unilateral é estimado. Se x for de valor complexo, ou o quot quot especifica quottwosidedquot. O espectro completo é estimado. Dados de peso da janela. Se a janela estiver vazia ou não especificada, uma janela retangular padrão é usada. Caso contrário, a janela é aplicada ao sinal (x. Win) antes de calcular o periodograma. Os dados da janela devem ser um vetor do mesmo comprimento que x. Número de caixas de freqüência. O padrão é 256 ou a potência superior seguinte de 2 maiores que o comprimento de x (max (256, 2.nextpow2 (comprimento (x)))). Se nfft for maior do que o comprimento da entrada, o x será preenchido com zero ao comprimento de nfft. taxa de amostragem. O padrão é 1. faixa de espectro. Quotonesidedquot calcula o espectro de 0..nfft21. Quottwidesidedquot calcula o espectro de 0..nfft-1. A segunda saída opcional w são as frequências angulares normalizadas. Para um cálculo unilateral w está no intervalo 0, pi se nfft é igual e 0, pi) se nfft for estranho. Da mesma forma, para um cálculo de dois lados w está na faixa de 0, 2pi ou 0, 2pi) dependendo de nfft. Se for especificada uma frequência de amostragem, Fs. Então as freqüências de saída f estarão na faixa 0, Fs 2 ou 0, Fs 2) para cálculos unilaterais. Para cálculos de dois lados, o intervalo será 0, Fs). Quando chamado sem saídas, o periodograma é imediatamente plotado na janela da figura atual. Retornar um sinetone de frequência freq com um comprimento de seg segundos a taxa de amostragem e com amplitude ampl. Os argumentos freq e ampl podem ser vetores de tamanho comum. Os padrões são a taxa 8000, segundo 1 e o ampl 64. Retorna um vetor de elemento m com i-ésimo elemento dado pelo pecado (2 pi (i d -1) n). O valor padrão para d é 0 e o valor padrão para n é m. Retornar o estimador de densidade espectral dado um vetor de autocovariâncias c. Nome da janela ganha. E largura de banda, b. O nome da janela, p. ex. Quottrianglequot ou quotrectanglequot é usado para procurar uma função chamada win lw. Se a vitória for omitida, a janela do triângulo é usada. Se b for omitido, use 1 sqrt (comprimento (x)). Retornar o estimador de densidade espectral dado um vetor de dados x. Nome da janela ganha. E largura de banda, b. O nome da janela, p. ex. Quottrianglequot ou quotrectanglequot é usado para procurar uma função chamada win sw. Se a vitória for omitida, a janela do triângulo é usada. Se b for omitido, use 1 sqrt (comprimento (x)). Retorna a média móvel de 15 pontos de Spencerrsquos de cada coluna de x. Calcule a transformada de Fourier de curto prazo do vetor x com coeficientes numcoef aplicando uma janela de pontos de dados de winsize e um incremento de pontos inc. Antes de computar a transformada de Fourier, uma das seguintes janelas é aplicada: os nomes das janelas podem ser passados ​​como strings ou pelo número wintype. Os padrões a seguir são usados ​​para argumentos não especificados: winsize 80, inc 24, numcoef 64 e wintype 1. y stroft (x. Hellip) retorna os valores absolutos dos coeficientes de Fourier de acordo com as freqüências positivas numcoef. Y. C stft (x. Hellip) retorna toda a matriz STFT y e um vetor de 3 elementos c que contém o tamanho da janela, o incremento e o tipo de janela, que é necessário para a função de síntese. Calcule um sinal de sua transformada de Fourier de curta duração y e um vetor de 3 elementos c especificando tamanho da janela, incremento e tipo de janela. Os valores y e c podem ser derivados por Ajustar um modelo de AR (p) com estimativas de Yule-Walker, dado um vetor c de autocovariâncias gamma0, hellip, gammap. Retorna os coeficientes AR, a. E a variância do ruído branco, v. Weighted Moving Average Em resposta a esta publicação por Luca Delucchi Um método on-line para fazer wma (ou seja, onde você pode adicionar um único valor ao mesmo tempo) com uma janela de expinência (que pesa eventos mais recentes mais fortes do que Menos recentes) funciona assim: xmean (t1) exp (-1tau) (xmean (t) x (t) tau) isso lhe dá uma wma de x (t). Tau é algo como o comprimento da memória (eventos mais atrasados ​​do que o tau não será muito lembrado demais). Am 10.08.2007 um 09:19 schrieb Luca Delucchi: gt Oi, eu posso fazer uma função na Média de Movimento Ponderada onde o valor é gt tomar no modo automático esta minha idéia gt gt yy1, y2, y3, y4, y5 gt função wma (Y) gt (y12y2y3) 4 gt (y22y3y4) 4 gt etc gt etc gt final função gt gt Não consegui repetir a fórmula (y12y2y3) 4 (porque se o vetor longo do gt for diferente eu preciso mudar a função) mas tenho Apenas uma fórmula gt que usa a fórmula para todos os valores de vetor gt gt Espero ter dado uma explicação clara gt gt Luca gt gt Help-octave mailing list gt email escondido gt cae. wisc. edumailmanlistinfohelp-octave Não é uma coisa de oitava, mas Um processo de processamento de sinal. Um filtro FIR (filtro de entrada finito) é determinado pelo vetor dos coeficientes, portanto, se o filtro tiver o comprimento 4, a saída seria algo como: y (t) b (1) x (t) b (2) x ( T-1) b (3) x (t-2) b (4) x (t-3) Então, quando b (1,4) 4, é apenas a média dos últimos quatro elementos. Na oitava, você pode usar a função quotfilterquot para fazer exatamente isso, se x for seu sinal, você pode simplesmente fazer o filtro (b, 1, x) P. S. Isso é quase o mesmo que a sugestão de Sren para usar a convolução (usando a função conv). A única diferença (eu acredito) é que o filtro dará a mesma saída que conv, mas truncado ao comprimento de x. Em 81007, Luca Delucchi enviou email gt escreveu: 2007810, Schirmacher, Rolf ltidden email gt: gt Filtrar com um filtro FIR Os coeficientes seriam gt gt b 1 1 1 1. 4 gt O que é isso Desculpe, mas I39m um novato de octave gt gt ----- Mensagem Original ----- gt gt De: Luca Delucchi mailto: email escondido gt gt Enviado: sexta-feira, 10 de agosto de 2007 9:20 Gt gt Para: octave gt gt Assunto: Média móvel ponderada gt gt gt gt gt gt Oi, eu posso fazer uma função na média de movimento ponderada onde o valor é gt gt tomar no modo automático esta minha idéia gt gt gt gt yy1, y2, Y3, y4, y5 gt gt função wma (y) gt gt (y12y2y3) 4 gt gt (y22y3y4) 4 gt gt etc gt gt etc gt gt função final gt gt gt gt Não consegui repetir a fórmula (y12y2y3) 4 (porque Se o longo vetor gt gt for diferente eu preciso mudar a função), mas tenho apenas uma fórmula gt gt que usa a fórmula para todos os valores do vetor gt gt gt gt Espero ter dado uma explicação clara gt gt gt gt Luca gt gt gt Gt Help-octave mailing list gt gt email escondido gt gt cae. wisc. edumailmanlistinfohelp-octave gt gt gt Em resposta a esta publicação por Luca Delucchi Oi, eu posso fazer uma função na Média de Movimento Ponderada onde o valor é assumido no automat Ic mode this my idea yy1, y2, y3, y4, y5 função wma (y) (y12y2y3) 4 (y22y3y4) 4 etc etc final função Eu não poderia repetir a fórmula (y12y2y3) 4 (porque se o longo do vetor é diferente Eu devo mudar a função), mas tenho apenas uma fórmula que usa a fórmula para todos os valores do vetor. Você deve pensar sobre seus dados de uma maneira diferente se você quiser usar o Matlaboctave de forma eficiente. Os dados são representados como vetores ou matrizes, e você deve fazer todas as operações em todos os dados --- não pense nos elementos y1, y2, etc., mas sim lidar com todo o vetor y. Você terá que misturar elementos vetoriais de diferentes posições, então você precisa construir versões deslocadas do vetor. Por exemplo, y (2: end) é o vetor cujo primeiro elemento é o segundo elemento de y. Quando você faz isso dessa forma, isso obriga você a reconhecer vários problemas que são varridos sob o tapete de outra forma, por exemplo, qual o significado da sua média ponderada para y1, que não tem um ponto de dados anterior. Uma abordagem pode ser duplicar a Primeiro e último ponto: temp y (1) yy (final) média (temp (1: final-2) 2temp (2: final-1) temp (3: fim)) 4 ou desista e admita que você só pode calcular A média em um subconjunto de pontos: média (y (1: fim-2) 2y (2: final-1) y (3: fim)) 4 Existe uma função de oitava chamada filtro () que pode aplicar um filtro linear arbitrário. Bastante complicado porque permite um feedback linear que você não está interessado em que você use uma forma específica de um vetor de feedback b1 0 0 0 0 0. filtro médio (1 2 14,1, y)) Finalmente, a Octave possui alguns filtros criados no meu Favorito é o filtro de conservação Savitzky-Golay de 2 momentos