Tuesday, 17 October 2017

Movendo média vs baixa passagem


Estou codificando algo no momento em que eu estou tomando um monte de valores ao longo do tempo de uma bússola de hardware. Esta bússola é muito precisa e atualiza com muita freqüência, com o resultado de que se ele jiggles ligeiramente, eu acabar com o valor estranho que é descontroladamente inconsistente com seus vizinhos. Quero suavizar esses valores. Tendo feito alguma leitura ao redor, parece que o que eu quero é um filtro passa-alta, um filtro passa-baixa ou uma média móvel. Movendo a média que eu posso começar com, basta manter um histórico dos últimos 5 valores ou o que quer que, e usar a média desses valores a jusante no meu código onde eu estava uma vez apenas usando o valor mais recente. Isso deve, eu acho, suavizar esses jiggles muito bem, mas parece-me que o seu provavelmente bastante ineficiente, e este é provavelmente um daqueles Problemas Conhecidos aos Programadores Adequados para que theres uma solução de matemática inteligente realmente limpo. Eu sou, entretanto, um daqueles programadores self-taught terríveis sem um shred da instrução formal em qualquer coisa vagamente relacionado a CompSci ou à matemática. Ler em torno de um pouco sugere que este pode ser um filtro de alta ou baixa passagem, mas eu não consigo encontrar nada que explique em termos compreensíveis para um hack como eu o que o efeito desses algoritmos seria sobre uma matriz de valores, muito menos como a matemática trabalho. A resposta dada aqui. Por exemplo, tecnicamente responde à minha pergunta, mas apenas em termos compreensíveis para aqueles que provavelmente já sabem como resolver o problema. Seria uma pessoa muito bonita e inteligente quem poderia explicar o tipo de problema que isso é, e como as soluções funcionam, em termos compreensíveis para um graduado em Artes. Se a sua média móvel tem que ser longa, a fim de alcançar o alisamento necessário, e você realmente não precisa de qualquer forma especial de kernel, então você está melhor se você usar uma média móvel exponencial decadência: onde você Escolha tiny para ser uma constante apropriada (por exemplo, se você escolher minúsculo 1 - 1N, ele terá a mesma quantidade de média como uma janela de tamanho N, mas distribuídos de forma diferente sobre pontos mais antigos). Enfim, uma vez que o próximo valor da média móvel depende apenas do anterior e seus dados, você não tem que manter uma fila ou qualquer coisa. E você pode pensar nisso como fazendo algo como, Bem, eu tenho um novo ponto, mas eu realmente não confio, então eu vou manter 80 da minha antiga estimativa da medição, e só confiar neste novo ponto de dados 20. Isso é Praticamente o mesmo que dizer, Bem, eu só confiar neste novo ponto 20, e eu uso 4 outros pontos que eu confio na mesma quantidade, exceto que em vez de tomar explicitamente os 4 outros pontos, você está supondo que a média que você fez na última vez Foi sensato para que você possa usar seu trabalho anterior. Respondeu Sep 21 10 at 14:27 Ei, eu sei que isso é 5 anos de atraso, mas obrigado por uma resposta incrível. I39m trabalhando em um jogo onde o som muda com base na sua velocidade, mas devido à execução do jogo em um computador lento ass, a velocidade flutuaria descontroladamente, o que era bom para a direção, mas super irritante em termos de som. Esta foi uma solução realmente simples e barata para algo que eu pensei que seria um problema muito complexo. Ndash Adam Mar 16 15 at 20:20 Se você estiver tentando remover o valor ímpar ocasional, um filtro passa-baixa é a melhor das três opções que você identificou. Os filtros de passagem baixa permitem mudanças de baixa velocidade, como as causadas pela rotação de uma bússola à mão, ao mesmo tempo em que rejeitam mudanças de alta velocidade, como as causadas por solavancos na estrada, por exemplo. Uma média móvel provavelmente não será suficiente, uma vez que os efeitos de um único blip em seus dados afetarão vários valores subseqüentes, dependendo do tamanho de sua janela de média móvel. Se os valores ímpares são facilmente detectados, você pode até ser melhor com um algoritmo de remoção de falhas que ignora completamente eles: Aqui está um gráfico guick para ilustrar: O primeiro gráfico é o sinal de entrada, com uma falha desagradável. O segundo gráfico mostra o efeito de uma média móvel de 10 amostras. O gráfico final é uma combinação da média de 10 amostras e do algoritmo de detecção de falhas simples mostrado acima. Quando a falha é detectada, a média de 10 amostras é usada em vez do valor real. Respondeu Sep 21 10 at 13:38 Bem explicado, e pontos de bônus para o gráfico) ndash Henry Cooke Sep 22 10 em 0:50 Uau. Seldomly viu uma resposta tão agradável ndash Muis Jun 4 13 at 9:14 A média móvel é um filtro passa-baixa. Ndash nomen Oct 21 13 at 19:36 Experimente uma medalha runningstreaming vez. Ndash kert Apr 25 14 at 22:09 Movendo média, eu posso descer com. Mas parece-me que o seu provavelmente bastante ineficaz. Não há realmente nenhuma razão uma média móvel deve ser ineficiente. Você mantém o número de pontos de dados desejados em algum buffer (como uma fila circular). Em cada novo ponto de dados, você pop o valor mais antigo e subtraí-lo de uma soma, e empurrar o mais novo e adicioná-lo à soma. Assim, cada novo ponto de dados realmente só implica um poppush, uma adição e uma subtração. Sua média móvel é sempre esta soma deslocada dividida pelo número de valores em seu buffer. Ele fica um pouco mais complicado se você está recebendo dados simultaneamente a partir de vários segmentos, mas desde que seus dados vem de um dispositivo de hardware que parece altamente duvidoso para mim. Oh e também: programadores auto-didáticos horríveis se unem) A média móvel parecia ineficiente para mim, porque você tem que armazenar um buffer de valores - melhor apenas fazer algumas Matemática Inteligente com seu valor de entrada e valor de trabalho atual Eu acho que é como média móvel exponencial trabalho. Uma otimização que eu tenho visto para esse tipo de média móvel envolve o uso de um amplificador de fila de comprimento fixo, um ponteiro para onde você está nessa fila, e apenas envolver o ponteiro ao redor (com ou um if). Voila Nenhum pushpop caro. Poder para os amadores, irmão ndash Henry: Para uma média móvel em linha reta você precisa do buffer simplesmente para que você saiba o valor obtém estourou quando o próximo valor ser empurrado. Dito isto, o amplificador de fila de comprimento fixo, um ponteiro que você está descrevendo, é exatamente o que eu quis dizer com fila quotcircular. Por isso eu estava dizendo que ele não é eficiente. O que você acha que eu quis dizer E se sua resposta é quotan array que muda seus valores de volta em cada remoção indexada (como std :: vector em C). Bem, eu não sei sobre AS3, mas um programador Java tem coleções como CircularQueue em sua disposição (I39m não é um Desenvolvedor Java para I39m certeza há melhores exemplos lá fora that39s apenas o que eu encontrei a partir de uma pesquisa rápida do Google), que implementa precisamente a funcionalidade we39re falando. I39m bastante confiante a maioria das linguagens de médio e baixo nível com bibliotecas padrão têm algo semelhante (por exemplo, no QueueltTgt lá). Enfim, eu era filosofia, então. tudo é perdoado. Ndash Dan Tao Sep 22 10 at 12:44 Uma média móvel decrescente exponencialmente pode ser calculada manualmente com apenas a tendência se você usar os valores adequados. Veja fourmilab. chhackdiete4 para uma idéia sobre como fazer isso rapidamente com uma caneta e papel, se você está procurando exponencialmente suavizada média móvel com 10 suavização. Mas desde que você tem um computador, você provavelmente quer fazer deslocamento binário ao contrário de deslocamento decimal) Desta forma, tudo que você precisa é uma variável para o seu valor atual e um para a média. A média seguinte pode então ser calculada a partir daí. Há uma técnica chamada uma porta de intervalo que funciona bem com amostras espúrias de baixa ocorrência. Assumindo o uso de uma das técnicas de filtro mencionadas acima (média móvel, exponencial), uma vez que você tenha histórico suficiente (uma Constante de Tempo), você pode testar a nova amostra de dados para a razoabilidade, antes de ser adicionada à computação. É necessário algum conhecimento da taxa máxima de mudança razoável do sinal. A amostra bruta é comparada com o valor mais recente suavizado e se o valor absoluto dessa diferença for maior do que o intervalo permitido, essa amostra é descartada (ou substituída por alguma heurística, por exemplo, uma previsão baseada no diferencial de inclinação ou na tendência Valor de predição do duplo suavização exponencial) respondido Abr. 30 16 a 6: 56 Um filtro de média móvel mede um número de amostras de entrada e produz uma única amostra de saída. Esta acção de média remove os componentes de alta frequência presentes no sinal. Filtros de média móvel são normalmente utilizados como filtros de passa-baixa. No algoritmo de filtragem recursiva, amostras de saída anteriores também são tomadas para a média. Um filtro de média móvel mede um número de amostras de entrada e produz uma única amostra de saída. Esta acção de média remove os componentes de alta frequência presentes no sinal. Filtros de média móvel são normalmente utilizados como filtros de passa-baixa. No algoritmo de filtragem recursiva, amostras de saída anteriores também são tomadas para a média. Essa é a razão pela qual sua resposta ao impulso se estende até o infinito. Como usar o programa de exemplo O arquivo. zip contém tanto o código-fonte como o executável. Para compilar e executar o código-fonte que você precisa ter o Visual Basic 6.0 instalado no seu computador. Para executar o executável, você deve baixar e instalar arquivos de tempo de execução do Visual Basic 6.0. Execute movavgfilt. exe e você verá a janela principal. Na janela principal. A parte mais superior é o gerador de função. Que produz diferentes formas de onda para testar o filtro. Podemos interativamente alterar a amplitude, freqüência e forma do sinal gerado. Para testar o programa, primeiro devemos gerar uma forma de onda apropriada. Aqui nós geraremos uma forma de onda complexa que consiste em duas freqüências diferentes. Deixe tudo nas configurações padrão e clique no botão quotgeneratequot. Agora você pode ver um sinal de 10 Hz no gráfico ao lado do gerador de sinal. A figura abaixo mostra a forma de onda. Agora mude a Freqüência para 100 Hz e clique no botão quotgeneratequot novamente. A forma de onda recém-gerada é adicionada à forma de onda existente e a forma de onda resultante parece uma onda sin 10Hz com ruído de 100 Hz. Veja a forma de onda abaixo. Esta forma de onda é mais adequada para testar o filtro, uma vez que contém duas freqüências diferentes. Você pode executar o filtro clicando no botão quotFilterquot. Das opções disponíveis à esquerda para o botão quotFilterquot. Você pode escolher Filtração recursiva, não recursiva ou sem filtragem. A figura abaixo mostra a saída do filtro. Download Moving Average Filtro de código de origemMoving Average Filter kate escreveu: gt Oi, gt gt Estou procurando algum código para um filtro passa-baixa que eu posso aplicar a gt um sinal antes de realizar a análise espectral. Gt gt Eu apoligise para minha ignorância, mas esta é maneira fora de meu campo assim Im gt que não faz realmente nenhum sentido dele. No domínio analógico, as pessoas usam a filtragem de baixa passagem por pelo menos um par de razões que vêm à mente (i) fazer o sinal parecer melhor (pt) Ii) evitar o aliasing durante a conversão Analog-to-Digital, o que resulta em sinais de ruído de alta freqüência sendo alias para baixas freqüências, o que pode corromper os sinais de menor frequência de interesse e aumentar o piso de ruído. Não parece que nenhuma destas considerações se aplique à sua situação (i) você não está olhando o sinal diretamente (você está indo fazer análise espectral) (ii) o seu sinal já está digitalizado. Especificamente, quando você faz análise espectral, o material de alta freqüência vai aparecer no final de alta freqüência e você pode optar por ignorá-lo. Para qualquer técnica linear (isto inclui FFT e a função Matlab filter ()), o conteúdo de alta freqüência não interferirá com a análise espectral do conteúdo de baixa freqüência. A menos que você deseje dizimar seus dados antes de filtrar. Existe uma razão específica que você quer se livrar do conteúdo de alta freqüência antes de análise espectral kate escreveu: gt Oi, gt gt Estou procurando algum código para um filtro passa-baixo que posso aplicar a gt um sinal antes de transportar Análise espectral. Gt gt Eu apoligise para minha ignorância, mas esta é maneira fora de meu campo assim Im gt que não faz realmente nenhum sentido dele. No domínio analógico, as pessoas usam a filtragem de baixa passagem por pelo menos um par de razões que vêm à mente (i) fazer o sinal parecer melhor (pt) Ii) evitar o aliasing durante a conversão Analog-to-Digital, o que resulta em sinais de ruído de alta freqüência sendo alias para baixas freqüências, o que pode corromper os sinais de menor frequência de interesse e aumentar o piso de ruído. Não parece que nenhuma destas considerações se aplique à sua situação (i) você não está olhando o sinal diretamente (você está indo fazer análise espectral) (ii) o seu sinal já está digitalizado. Especificamente, quando você faz análise espectral, o material de alta freqüência vai aparecer no final de alta freqüência e você pode optar por ignorá-lo. Para qualquer técnica linear (isto inclui FFT e a função Matlab filter ()), o conteúdo de alta freqüência não interferirá com a análise espectral do conteúdo de baixa freqüência. A menos que você deseje dizimar seus dados antes de filtrar. Existe uma razão especial que você quer se livrar do conteúdo de alta freqüência antes da análise espectral Para ser honesto, eu não sei por que estou tentando se livrar das altas freqüências. Im basicamente seguindo as instruções em um ISO. Como você pode ter adivinhado, programação de computador e processamento de sinal não é realmente a minha área para que a linguagem utilizada é estranho para mim O que estou fazendo é a seguinte - Im um engenheiro civil e estou tentando analisar um perfil de superfície da estrada. O perfil é basicamente o equivilent de um sinal que varia com a distância (mas desde que a velocidade é constante, isto é o mesmo que variando com o tempo). A formulação exata do ISO é filtros de pré-processamento deve ser usado, por exemplo, butterworth. No entanto, eu pensei que a média móvel pode ser um lugar mais fácil para começar, presumo que a razão Im tentando erradicar as altas freqüências é porque eles seriam insignificantes em termos de danos na superfície da estrada. Eu aprecio muito o seu tempo, Katherine Rajeev escreveu: gt gt gt kate escreveu: gtgt Oi, gtgt gtgt Estou procurando algum código para um filtro passa-baixa que eu posso gt aplicar gtgt um sinal antes de realizar a análise espectral. Eu apoligise para minha ignorância, mas esta é maneira fora de meu campo assim gt Im gtgt que não faz realmente nenhum sentido dele. Gtgt gt gt gtgt gtgt gtgt No domínio analógico, as pessoas usam filtragem de baixa passagem para pelo menos um gt duas razões que vêm à mente (i) fazer o sinal (Ii) evitar alias durante a conversão Analog-to-Digital, que gt resulta em sinais de ruído de alta freqüência sendo alias para baixas freqüências gt, que pode corromper os sinais de freqüência mais baixa de gt interesse gt e aumentar o piso de ruído. Gt gt Não parece que qualquer destas considerações se aplicam a gt sua situação gt (i) você não está olhando para o sinal diretamente (youre gt vai gt fazer análise espectral) (ii) o seu sinal já está digitalizado. Gt gt Especificamente, quando você faz análise espectral, a alta freqüência gt stuff gt aparecerá no final de alta freqüência e você pode optar por ignorar gt-lo. Gt Para qualquer técnica linear (isto inclui FFT e a função gt do filtro Matlab), o conteúdo de alta freqüência não interferirá com a análise espectral gt do conteúdo de baixa freqüência. A menos que você deseja gt decimate seus dados antes de filtrar. Gt gt Existe uma razão específica que você quer se livrar do gt alta freqüência gt conteúdo antes de análise espectral gt gt HTH gt - rajeev - gt gt Katherine escreveu: gt Para ser honesto, eu não sei por que estou tentando se livrar do Altas freqüências gt. Im basicamente seguindo as instruções em um ISO. Gt Como você pode ter adivinhado, programação de computadores e gt de processamento de sinal não é realmente a minha área para que a linguagem utilizada é estranho para mim gt gt O que estou fazendo é a seguinte - Im um engenheiro civil e estou tentando gt analisar um perfil de superfície da estrada. O perfil é basicamente o gt equivilent de um sinal que varia com a distância (mas desde que a velocidade gt é constante, isto é o mesmo que variando com o tempo). O exato gt redacção do ISO é pré-processamento filtros devem ser utilizados para Algumas perguntas vêm à mente. uma. O que o ISO lhe pede que faça após os filtros de pré-processamento b. Como é realizada a análise espectral c. O ISO especifica a freqüência de corte para o filtro. Ie livrar-se de freqüências acima de X gt example butterworth. No entanto, eu pensei que a média móvel gt poderia ser um lugar mais fácil para começar eu tendem a concordar, a média móvel seria mais fácil. Ele também tem uma propriedade que todos os componentes de freqüência são atrasados ​​exatamente pela mesma quantidade, o que significa que a forma de onda é preservada através do filtro (claro que alguns compnents de freqüência serão atenuados, mas eles não serão deslocados por, digamos, 90 graus , Em relação a outras frequências). O filtro de Butterworth (e em vários graus todos os filtros analógicos) não tem essa propriedade, que é conhecida como linear-fase ou fase-linear. Butterworth refere-se a uma classe de filtros analógicos com uma fase particular e resposta de freqüência, que passa a ser fácil de implementar com componentes eletrônicos como resistores, capacitores e indutores. (Minha suposição razoável é que) as pessoas desenvolveram equivalentes digitais para estes e outros filtros analógicos porque estavam familiarizados com suas propriedades. No entanto, um monte de gente hoje perguntar, se você está indo para operar em um sinal digitalizado, por que se preocupar com um análogo-look-alike filtro. Gt Presumo que a razão pela qual estou tentando erradicar as altas freqüências é gt porque seria insignificante em termos de danos na superfície da estrada. Gt gt Eu aprecio muito o seu tempo, gt Katherine Mais uma vez, estou muito grato a você por tomar o tempo que eu tentei responder a sua qs abaixo: gt Algumas perguntas vêm à mente. Gt gt a. O que o ISO lhe pede que faça após os filtros de pré-processamento Após os filtros de pré-processamento pede que eu realize uma FFT que eu acho que é também uma resposta para a sua próxima pergunta. O grande problema de compreensão que estou tendo é que eu gerei o perfil da estrada, especificando que eu queria que as freqüências fossem um mínimo de 0,01cyclesmeter e um máximo de 4cyclesmeter. Por que então eu preciso para filtrar as altas freqüências gt gt b. Como é realizada a análise espectral gt gt c. O ISO especifica a freqüência de corte para o filtro. Ou seja, gt get gt livre de freqüências acima de X Não especifica qualquer freqüência de corte. Gtgt exemplo butterworth. No entanto, eu pensei que a média móvel gtgt poderia ser um lugar mais fácil para começar gt gt eu tendem a concordar, a média móvel seria mais fácil. Ele também tem uma propriedade gt de que todos os componentes de freqüência são retardados exatamente pela mesma quantidade de gt, gt o que significa que a forma de onda é preservada através do filtro gt gt (claro que alguns compnents de freqüência serão atenuados, mas eles gt wt gt Ser deslocada por, digamos, 90 graus, em relação a outras freqüências). Gt O filtro gt Butterworth (e em vários graus todos os filtros analógicos) gt não gt tem essa propriedade, que é conhecida como linear-fase ou fase-linear. Gt gt Butterworth refere-se a uma classe de filtros analógicos com uma particular gt fase gt e resposta de freqüência, que passa a ser fácil de implementar com gt eletrônicos gt componentes como resistências, capacitores e indutores. (Meu gt razoável gt acho que gt é que) as pessoas desenvolveram equivalentes digitais para estes e outros gt gt filtros analógicos porque estavam familiarizados com suas propriedades. No entanto gt um lote gt gt de pessoas hoje perguntaria, se você está indo para operar em um sinal gt digitalizado, gt por que se preocupar com um análogo-look-alike filtro. Eu presumo que a razão que estou tentando erradicar as altas freqüências é gtgt porque seriam insignificantes em termos de danos na superfície da estrada. Gtgt gtgt Eu aprecio muito o seu tempo, gtgt Katherine gt gt lt. Gt gt gt HTH gt - rajeev - Obrigado. Katherine Parece que você pode estar filtrando os dados já da maneira como você está especificando a faixa de freqüência. Qual é sua taxa de amostragem É espacial ou temporal Se você está especificando 4 ciclosmetro para o sistema é muito improvável que seria apenas a amostragem para obter essa taxa (Fs18 metros) sem algum tipo de filtro de média móvel construído dentro O que é o ISO Requisito (padrão ISO, a partir de onde) Um efeito da filtragem é deslocar a energia para as freqüências mais baixas em vez de apenas cortá-lo como você faria no domínio da freqüência. Se o objetivo final é calcular um IRI ou algum tipo de outra métrica de rugosidade do caminho que isso pode ser crítico. Gt gt Após o pré-processamento de filtros que pede que eu realizar uma FFT que gt eu acho que também é uma resposta para a sua próxima pergunta. O grande problema de compreensão do gt que estou tendo é que eu gerei o perfil da estrada gt, especificando que eu queria que as freqüências fossem um mínimo de 0,01cyclesmeter e um máximo de 4cyclesmeter. Por que então gt deve eu preciso para filtrar as altas freqüências gt Charlie, eu sou muito ignorante sobre a terminologia correta neste material e Im não sei o que você quer dizer por taxa de amostragem. Vou te dizer o que estou fazendo. Primeiro eu estou gerando um perfil de estrada aleatória que tem freqüências espaciais variando de 0,01 - 4 ciclosm. A ISO 8608: 1995 tem classificações de estrada e, dependendo disso, dá um valor PSD para cada uma das freqüências entre 0,01 e 4 thats você deseja. Estes valores são então colocados numa equação para a geração de estradas que cria uma estrada com qualquer número de pontos (no meu caso, 8000, ou 400 metros, isto é, cada 0,05 metros). Eu então grafo todos os valores ISO para o PSD contra as freqüências espaciais que eu tinha acima. Eu estou tentando então trabalhar para trás para ver se eu posso gerar esse mesmo gráfico usando o mesmo perfil da estrada, e encontrando o FFT dele e então o PSD. Eu não sei o que você quer dizer com frequência de amostragem Tem medo, talvez esteja lá em cima no que eu descrevi Muito obrigado pelo seu tempo, Estou completamente como um peixe fora da água sobre este Charlie escreveu: gt gt gt Katherine, Gt gt Parece que você pode estar filtrando os dados já a maneira que você está gt especificando gt a faixa de freqüência. Qual é sua taxa de amostragem É espacial ou gt temporal gt Se você está especificando 4 ciclosmetro para o sistema é muito improvável gt que gt seria apenas a amostragem para obter essa taxa (Fs18 metros), sem alguns gt tipo de gt média móvel filtro construído Gt gt O que é o requisito ISO (norma ISO, de onde) gt gt Um efeito da filtragem é deslocar a energia para as freqüências gt inferior, em vez de apenas cortá-lo como você faria em gt o domínio da freqüência gt. Se o objetivo final é calcular um IRI ou algum tipo gt de outra métrica de rugosidade de estrada que isso pode ser crítico. Gt gt gt gtgt gtgt Após o pré-processamento de filtros que pede que eu realizar uma FFT gt que gtgt acho que também é uma resposta para a sua próxima pergunta. O grande problema de compreensão que eu tenho é que eu gerei o perfil gtgt, especificando que eu queria que as freqüências fossem um mínimo gtgt de 0,01cyclesmeter e um máximo de 4cyclesmeter. Por que então gtgt devo precisar filtrar as altas freqüências gtgt gt gt gt Obrigado pela informação sobre ISO 8608: 1995 parece ser uma boa referência para alguns dos meus trabalhos sobre o processamento de perfil de estrada. De volta ao seu projeto. Como eu entendo você está fazendo: 1. Criar perfil de estrada em domínio de freqüência espacial com conteúdo em 0.01-4 ciclos 2. Gerar perfil espacial de 1 usando algumas equações (400 metros de comprimento, dx0,05 m, Freqüência de amostragem espacial1dx20 ciclosm) 3 Gráfico seu PSD de estrada de 1 contra os valores de ISO de ISO 8608 4. Calcule o fft e o PSD de 2 e compare-o a 3 para ver se você pode re-produzi-lo. Se isso é correto e eu entendo o padrão ISO. Eu não acredito que você precisa fazer qualquer filtragem em tudo. Seu perfil de 2 deve ser capaz de gerar dados de freqüência de 0,0025-10 ciclos, mas você não deve ver qualquer conteúdo acima de 4 ciclos. Espero que isso ajude em vez de confunde. Você pode querer olhar para o pequeno livro de perfis em umtri. umich. eduerdroughnessindex. html para mais informações. Katherine ltkatherine. cashellucd. iegt escreveu na mensagem news: ef02d7a.7webx. raydaftYaTP. Gt Charlie, gt Eu sou muito ignorante sobre a terminologia correta neste material e Im gt não sei o que você quer dizer com a taxa de amostragem. Vou te dizer o que estou fazendo. Gt gt gt Primeiro eu estou gerando um perfil de estrada aleatória que tem freqüências espaciais gt variando de 0,01-4 ciclosm. A ISO 8608: 1995 tem gt classificações de estrada e, dependendo disso, dá um valor PSD gt para cada uma das freqüências entre 0,01 e 4 thats você deseja. Estes valores gt são então colocados numa equação para a geração de estrada que gt cria uma estrada com qualquer número de pontos (no meu caso 8000 ou gt 400meters, isto é, cada 0,05 metro). Gt I então grafo todos os valores ISO para o PSD contra as freqüências espaciais gt que eu tinha acima. Gt Eu estou tentando então trabalhar para trás para ver se eu posso gerar esse mesmo gráfico do gt usando o mesmo perfil da estrada, e encontrando o FFT dele e gt então o PSD. Gt eu não sei o que você quer dizer com freqüência de amostragem Im medo, talvez ele é lá em cima no que eu tenho descrito gt gt Muito obrigado pelo seu tempo, estou completamente como um peixe de água gt sobre este gt gt Katherine Gt Obrigado por isso - realmente é útil apenas para ver a terminologia correta sendo usado para as figuras Charlie escreveu: gt gt gt Katherine, gt gt Obrigado pela informação sobre ISO 8608: 1995 parece boa referência gt para alguns gt do meu trabalho No processamento do perfil da estrada. De volta ao seu projeto. Como gt gt entendê-lo você está fazendo: gt gt 1. Criar perfil de estrada em domínio de freqüência espacial com conteúdo em gt 0.01-4 gt cyclesm gt 2. Gerar perfil espacial de 1 usando algumas equações (400 gt metros de comprimento, gt dx0. 05 m, Frequência de amostragem espacial1dx20 cyclesm) gt 3. Grafique o PSD de estrada de 1 contra os valores ISO de ISO gt 8608 gt 4. Calcule o fft e o PSD de 2 e compare-os com 3 para gt veja se gt você é capaz de Re-produzi-lo. Gt gt Se isso é correto e eu entendo o padrão ISO. Eu não gt acredito que você gt necessidade de fazer qualquer filtragem em tudo. Seu perfil de 2 deve ser capaz de gt gerar dados de freqüência de 0,0025-10 ciclos, mas você não deve ver nenhum conteúdo gt acima de 4 ciclos. Espero que isso ajude em vez de confunde. Você pode querer olhar para o gt Little gt livro de perfis em ltumtri. umich. eduerdroughnessindex. html gt gt gt ou mais informações. Gt gt Charlie gt gt Katherine ltkatherine. cashellucd. iegt escreveu na mensagem gt news: ef02d7a.7webx. raydaftYaTP. Gtgt Charlie, gtgt Eu sou muito ignorante sobre a terminologia correta neste material e gt gtg gtgt não tenho certeza o que você quer dizer com a taxa de amostragem. Vou te dizer o que estou fazendo. Gtgt gtgt gtgt Primeiro estou gerando um perfil de estrada aleatória que tem freqüências espaciais gtgt variando de 0,01-4 ciclosm. A ISO 8608: 1995 tem gtgt classificações de estrada e, dependendo disso, dá um PSD gt gtgt valor para cada uma das freqüências entre 0,01 e 4 thats você deseja. Gt Estes valores gtgt são então colocados numa equação para a geração de estrada que gtgt cria uma estrada com qualquer número de pontos (no meu caso 8000 ou gtgt 400meters, isto é, cada 0,05 metro). Gtgt Eu então grafo todos os valores ISO para o PSD contra o gt gtgt espacial freqüências que eu tinha acima. Gtgt Eu estou tentando então trabalhar para trás para ver se eu posso gerar esse mesmo gtgt gt usando o mesmo perfil de estrada, e encontrando o FFT dele gt e gtgt então o PSD. Gtgt eu não sei o que você quer dizer com freqüência de amostragem Im medo, talvez gt que gtgt está lá em cima no que eu tenho descrito gtgt gtgt Muito obrigado pelo seu tempo, estou completamente como um peixe gt de gtgt água sobre este gtgt Gtgt Katherine gtgt gt gt gt O que é uma lista de observação Você pode pensar em sua lista de observação como segmentos que você marcou. Você pode adicionar tags, autores, threads e até mesmo resultados de pesquisa à sua lista de observação. Desta forma, você pode facilmente acompanhar os tópicos que você está interessado polegadas Para ver a sua lista de observação, clique no link quotMas newsreaderquot. Para adicionar itens à sua lista de observação, clique no link quotadd para assistir listquot na parte inferior de qualquer página. Como adicionar um item à minha lista de observação Para adicionar critérios de pesquisa à sua lista de observação, pesquise o termo desejado na caixa de pesquisa. Clique no botão quotAdicionar esta pesquisa ao meu link de listagem de visualizações na página de resultados de pesquisa. Você também pode adicionar uma tag à sua lista de observação procurando a tag com a diretiva quottag: tagnamequot onde tagname é o nome da tag que você gostaria de assistir. Para adicionar um autor à sua lista de observação, acesse a página de perfil dos autores e clique no botão quotAdicionar este autor ao meu link de lista de atalhos na parte superior da página. Você também pode adicionar um autor à sua lista de observação, indo a um tópico que o autor postou e clicando no quotAdicionar este autor ao meu link listquot do relógio. Você será notificado sempre que o autor fizer um post. Para adicionar um tópico à sua lista de observação, vá para a página do tópico e clique no link quotAdicionar este tópico ao meu link de lista de observações na parte superior da página. Sobre Newsgroups, Newsreaders e MATLAB Central O que são newsgroups Os newsgroups são um fórum mundial aberto a todos. Os newsgroups são usados ​​para discutir uma enorme variedade de tópicos, fazer anúncios e trocar arquivos. As discussões são encadeadas ou agrupadas de forma a permitir que você leia uma mensagem postada e todas as suas respostas em ordem cronológica. Isto torna mais fácil seguir o fio da conversa e ver whatrsquos já foi dito antes de postar sua própria resposta ou fazer uma nova postagem. O conteúdo do grupo de notícias é distribuído por servidores hospedados por várias organizações na Internet. As mensagens são trocadas e gerenciadas usando protocolos de padrão aberto. Nenhuma entidade única ldquoownsrdquo os newsgroups. Existem milhares de newsgroups, cada um abordando um único tópico ou área de interesse. O MATLAB Central Newsreader publica e exibe mensagens no newsgroup comp. soft-sys. matlab. Como posso ler ou publicar nos newsgroups Você pode usar o newsreader integrado no site da MATLAB Central para ler e publicar mensagens neste newsgroup. MATLAB Central é hospedado por MathWorks. As mensagens enviadas através do Central Newsreader MATLAB são vistas por todos os grupos de notícias, independentemente de como eles acessam os grupos de notícias. Há várias vantagens em usar o MATLAB Central. Uma conta A sua conta MATLAB Central está ligada à sua conta MathWorks para fácil acesso. Use o endereço de e-mail da sua escolha O MATLAB Central Newsreader permite que você defina um endereço de e-mail alternativo como seu endereço de postagem, evitando a confusão em sua caixa de correio principal e reduzindo o spam. Controle de Spam A maioria do spam do newsgroup é filtrada para fora pelo newsreader central de MATLAB. Marcação As mensagens podem ser marcadas com um rótulo relevante por qualquer usuário conectado. As tags podem ser usadas como palavras-chave para encontrar determinados arquivos de interesse ou como uma maneira de categorizar suas postagens marcadas. Você pode optar por permitir que outras pessoas visualizem suas tags e você pode exibir ou pesquisar outras tags, assim como as da comunidade em geral. Tagging fornece uma maneira de ver tanto as grandes tendências e as menores, mais obscuros idéias e aplicações. Listas de vigilância A configuração de listas de observação permite que você seja notificado das atualizações feitas nos lançamentos selecionados por autor, segmento ou qualquer variável de pesquisa. As notificações da sua lista de observação podem ser enviadas por email (resumo diário ou imediato), exibidas em Meu leitor de notícias ou enviadas via feed RSS. Outras formas de aceder aos newsgroups Utilize um leitor de notícias através da sua escola, entidade patronal ou fornecedor de serviços Internet Pagar o acesso a grupos de notícias de um fornecedor comercial Utilizar Grupos do Google Mathforum. org fornece um leitor de notícias com acesso ao grupo de notícias comp. soft sys. matlab Execute o seu próprio servidor. Para obter instruções típicas, consulte: slyckng. phppage2 Selecione Seu PaísEsses termos de alta, baixa e de banda referem-se a freqüências. No high-pass, você tenta remover freqüências baixas. No low-pass, você tenta remover alta. Na passagem de banda, você só permite que uma faixa de freqüência contínua permaneça. Escolher a frequência de corte depende da sua aplicação. A codificação desses filtros pode ser feita simulando circuitos RC ou jogando com transformadas de Fourier de seus dados baseados no tempo. Veja os artigos da Wikipédia para exemplos de código. Aqui está como você implementar um filtro passa-baixa usando convolução: Observe que o exemplo é extremamente simplificado. Ele não faz verificações de intervalo e não manipula as bordas corretamente. O filtro usado (caixa-carro) é um filtro lowpass particularmente ruim, porque causará muitos artefatos (toque). Leia acima no projeto do filtro. Você também pode implementar os filtros no domínio da freqüência. Aqui está como você implementar um filtro passa-alta usando FFT: Novamente, isso é simplificado, mas você começa a idéia. O código não parece tão complicado quanto a matemática. Respondeu Sep 17 08 at 12:06 Muito legal ter amostras de código. Por que convolução em um caso e FFT no outro ndash dfrankow Mar 13 09 em 19:03 dfrankow Nenhuma razão específica. Só para mostrar como ele se parece nos diferentes domínios. Atualizado o texto para refletir isso. Obrigado. Você tem certeza que a primeira parte de sua resposta está correta, onde você aplica convolução no domínio do tempo usando uma função de retângulo Eu pensei que um filtro passa-baixa no domínio do tempo exigia a convolução de um sinc A filtragem descreve o ato de processar dados de uma forma que aplica diferentes níveis de atenuação para diferentes freqüências dentro dos dados. Um filtro passa-alta aplicará uma atenuação mínima (isto é, níveis de permissão inalterados) para altas frequências, mas aplica atenuação máxima a baixas frequências. Um filtro passa-baixa é o inverso - ele não aplicará atenuação para baixas freqüências por aplica atenuação para altas freqüências. Há uma série de algoritmos de filtragem diferentes que são usados. Os dois mais simples são, provavelmente, o filtro de Resposta de Impulso Finito (também conhecido como filtro FIR) e o filtro de Resposta de Impulso Infinito (também conhecido como filtro IIR). O filtro FIR funciona mantendo uma série de amostras e multiplicando cada uma dessas amostras por um coeficiente fixo (que é baseado na posição na série). Os resultados de cada uma dessas multiplicações são acumulados e é a saída para essa amostra. Isto é referido como Multiply-Accumulate - e no hardware DSP dedicado há uma instrução MAC específica para fazer exatamente isso. Quando a amostra seguinte é tirada, é adicionada ao início da série e a amostra mais antiga da série é removida eo processo é repetido. O comportamento do filtro é fixado pela seleção dos coeficientes do filtro. Um dos filtros mais simples que muitas vezes é fornecido pelo software de processamento de imagem é o filtro de média. Isso pode ser implementado por um filtro FIR, definindo todos os coeficientes de filtro para o mesmo valor. Respondeu Oct 5 08 em 2:29

No comments:

Post a Comment