Moving Average Filter Kate skrev: gt Hej, gt gt Jag letar efter någon kod för ett lågpassfilter som jag kan ansöka om att gt en signal innan spektralanalys utförs. gt Jag avskyr mig för min okunnighet, men det här är långt utanför mitt fält, så jag känner inte riktigt av det. Vad är de ingångar som behövs för andra än själva signalen gt gt Tack, gt Kate I den analoga domänen använder folk lågpassfiltrering av åtminstone några orsaker som kommer att tänka på (i) få signalen att se bättre ut ( ii) Undvik aliasing under analog-till-digital konvertering, vilket resulterar i att högfrekventa ljudsignaler aliaseras till låga frekvenser, vilket kan korrumpera de lägre frekvenssignalerna av intresse och öka ljudgolvet. Det verkar inte som om någon av dessa överväganden gäller din situation (i) du tittar inte direkt på signalen (du ska göra spektralanalys) (ii) din signal är redan digitaliserad. Specifikt, när du gör spektralanalys kommer högfrekvensen att dyka upp i högfrekvensen och du kan välja att ignorera den. För alla linjära tekniker (detta inkluderar FFT och Matlab filter () - funktionen) kommer högfrekvensinnehållet inte att störa spektralanalysen av lågfrekvensinnehållet. Om du inte vill avkoda dina uppgifter innan du filtrerar. Finns det en särskild anledning att du vill bli av med högfrekvensinnehållet före spektralanalysen skrev kate: gt Hej, gt gt Jag letar efter någon kod för ett lågpassfilter som jag kan ansöka om att gt en signal före bärande ut spektralanalys. gt Jag avskyr mig för min okunnighet, men det här är långt utanför mitt fält, så jag känner inte riktigt av det. Vad är de ingångar som behövs för andra än själva signalen gt gt Tack, gt Kate I den analoga domänen använder folk lågpassfiltrering av åtminstone några orsaker som kommer att tänka på (i) få signalen att se bättre ut ( ii) Undvik aliasing under analog-till-digital konvertering, vilket resulterar i att högfrekventa ljudsignaler aliaseras till låga frekvenser, vilket kan korrumpera de lägre frekvenssignalerna av intresse och öka ljudgolvet. Det verkar inte som om någon av dessa överväganden gäller din situation (i) du tittar inte direkt på signalen (du ska göra spektralanalys) (ii) din signal är redan digitaliserad. Specifikt, när du gör spektralanalys kommer högfrekvensen att dyka upp i högfrekvensen och du kan välja att ignorera den. För alla linjära tekniker (detta inkluderar FFT och Matlab filter () - funktionen) kommer högfrekvensinnehållet inte att störa spektralanalysen av lågfrekvensinnehållet. Om du inte vill avkoda dina uppgifter innan du filtrerar. Finns det en särskild anledning att du vill bli av med högfrekvensinnehållet före spektralanalys För att vara ärlig vet jag inte varför jag försöker bli av med högfrekvenserna. Jag följer i princip instruktionerna i en ISO. Som du kanske har gissat är datorprogrammering och signalbehandling verkligen inte mitt område så språket som används är främmande för mig. Vad jag gör är som följer - Jag är civilingenjör och jag försöker analysera en vägprofil. Profilen är i grunden likvärdig med en signal som varierar med avstånd (men eftersom hastigheten är konstant, är detta samma som varierande med tiden). Den exakta formuleringen av ISO är förbehandlingsfilter bör användas till exempel butterworth. Men jag trodde att det rörliga genomsnittet skulle kunna vara en lättare plats att börja jag antar anledningen till att jag försöker utrota högfrekvenser är att de skulle vara försumbar när det gäller vägarbetsskador. Jag uppskattar mycket din tid, skrev Katherine Rajeev: gt gt gt kate skrev: gtgt Hej, gtgt gtgt Jag letar efter någon kod för ett lågpassfilter som jag kan applicera på gtgt en signal innan spektralanalys utförs. gtgt gtgt Jag ber om ursäkt för min okunnighet, men det här är långt utanför mitt fält, så det gör jag inte riktigt. Vad är de ingångar som behövs andra än själva signalen gtgt Tackgt, gtgt Kate gt gt I den analoga domänen använder folk lågpassfiltrering åtminstone ett antal orsaker som kommer att tänka på (i) göra signalen ser bättre ut gt (ii) undviker aliasing under analog-till-digital konvertering, vilket resulterar i att högfrekventa ljudsignaler aliaseras till låg-gt-frekvenser, vilket kan korrumpera de lägre frekvenssignalerna med gt interest gt och öka ljudgolvet. gt) Det ser inte ut som att någon av dessa överväganden gäller för din gt-situation (i) du tittar inte direkt på signalen (du kommer att göra spektralanalys) (ii) din signal är redan digitaliserad. Gt gt När du gör spektralanalys kommer de högfrekventa gt-grecken gt att dyka upp i högfrekvensänden och du kan välja att ignorera det. gt För någon linjär teknik (detta inkluderar FFT och Matlab filter () gt-funktionen) kommer högfrekvensinnehållet inte att störa GT spektralanalysen av lågfrekvensinnehållet. Om du inte vill att gt decimerar dina data innan du filtrerar. Gt gt Finns det en särskild anledning att du vill bli av med gt högfrekvent gt-innehåll innan spektralanalys gt gt HTH gt - rajeev - gt gt Katherine skrev: gt För att vara ärlig vet jag inte varför jag försöker bli av med hög gt frekvenser. Jag följer i princip instruktionerna i en ISO. gt Som du kanske har gissat är dataprogrammering och signalbehandling gt egentligen inte mitt område så det språk som används är främmande för mig. Det jag gör är som följer - Jag är civilingenjör och jag försöker att analysera en vägprofil. Profilen är i princip gt equivilent av en signal som varierar med avstånd (men eftersom hastigheten gt är konstant, är det samma som varierande med tiden). Den exakta GT formuleringen av ISO är förbehandlingsfilter bör användas för Vissa frågor kommer i åtanke. en. Vad frågar ISO dig efter förbehandlingsfiltret b. Hur implementeras spektralanalysen c. Anger ISO huruvida cutoff frekvensen för filtret. dvs bli av med frekvenser över X gt exempel butterworth. Men jag trodde att det rörliga genomsnittet gt kan vara en lättare plats att börja jag tenderar att hålla med, glidande medelvärde skulle vara enklare. Det har också en egenskap att alla frekvenskomponenter är fördröjda med exakt samma mängd, vilket innebär att vågformsformen bevaras genom filtret (givetvis kommer vissa frekvenskomponenter att dämpas men de kommer inte att förskjutas med exempelvis 90 grader , i förhållande till andra frekvenser). Butterworth-filtret (och i varierande grad alla analoga filter) har inte denna egenskap, som är känd som linjärfas eller faslinjär. Butterworth hänvisar till en klass av analoga filter med ett visst fas - och frekvenssvar, som är lätt att implementera med elektroniska komponenter som motstånd, kondensatorer och induktorer. (Mitt rimliga gissning är det) folk utvecklade digitala ekvivalenter till dessa och andra analoga filter eftersom de var bekanta med deras egenskaper. Men många människor idag frågar, om du ska arbeta på en digitaliserad signal, varför stör med ett analogt utseende-liknande filter. gt Jag antar anledningen till att jag försöker utrota högfrekvenser är gt eftersom de skulle vara försumbara när det gäller skador på vägarna. gt Jag uppskattar mycket din tid, gt Katherine Igen, jag är mycket skuldsatt för dig för att ta den tid jag har försökt att svara på dina qs nedan: gt Några frågor kommer i åtanke. gt gt a. Vad frågar ISO dig efter förbehandlingsfiltret Efter förbehandlingsfiltret frågar jag att jag utför en FFT som jag antar är också ett svar på din nästa fråga. Det stora förståelseproblemet som jag har är att jag genererade vägprofilen själv och specificerade att jag ville att frekvenserna skulle vara minst 0,01cykelsmeter och en max 4cyclesmeter. Varför då borde jag behöva filtrera ut högfrekvenser gt gt b. Hur implementeras spektralanalysen gt gt c. Specificerar ISO ISO-gränsen för filtret. dvs gt får bort frekvensen över X Det anger inte någon cutoff-frekvens. gtgt exempel butterworth. Men jag trodde att det rörliga genomsnittet gtgt kan vara en lättare plats att starta gt gt. Jag tenderar att hålla med, glidande medelvärde skulle vara enklare. Det har också en gt property gt att alla frekvenskomponenter är fördröjda med exakt samma gtmängd, gt vilket betyder att vågformsformen bevaras genom gt filter gt (naturligtvis kommer vissa frekvenskomponenter att dämpas men de kommer inte att vara gt skiftas med 90 grader i förhållande till andra frekvenser). Gt The Butterworth filter (och i varierande grad alla analoga filter) har inte denna egenskap, som är känd som linjärfas eller faslinjär. gt gt Butterworth hänvisar till en klass av analoga filter med ett specifikt gtfas gt - och frekvenssvar, som råkar vara lätt att implementera med gt electronic gt-komponenter som motstånd, kondensatorer och induktorer. (Jag är rimlig), det var folk som utvecklade digitala ekvivalenter till dessa och andra gt analoga gtfilter eftersom de var bekanta med deras egenskaper. Men om en gt lot gt av folk idag skulle fråga, om du kommer att fungera på en digitaliserad gt-signal, varför bry sig med ett analogt utseende-filter. gt gtgt Jag antar anledningen till att jag försöker utrota högfrekvenser är gtgt eftersom de skulle vara försumbar när det gäller skadan på vägarna. gtgt gtgt Jag uppskattar din tid, gtgt Katherine gt gt lt. gt gt gt HTH gt - rajeev - Tack. Katherine Låter som att du kanske filtrerar data redan som du anger frekvensområdet. Vad är samplingsfrekvensen Är det rumsligt eller tidsmässigt Om du specificerar 4 cyklermätare till systemet är det mycket osannolikt att det bara skulle vara provtagning för att få den frekvensen (Fs18 meter) utan någon form av glidande medelfilter som byggdes in. Vad är ISO Krav (ISO-standard, varifrån) En effekt av filtreringen är att flytta energin till lägre frekvenser istället för att bara haka av den som du skulle göra i frekvensdomänen. Om slutmålet är att beräkna en IRI eller någon annan väggruvhet metrisk än detta kan vara kritisk. gt gt Efter förbehandlingsfiltret frågar jag att jag utför en FFT som jag antar också är ett svar på din nästa fråga. Det stora GT-begripsproblemet som jag har är att jag genererade GT gt-profilen själv och specificerade att jag ville att frekvenserna skulle vara minst 0,01cykelsmeter och en max 4cyclesmeter. Varför då gt borde jag behöva filtrera ut högfrekvenser gt Charlie, jag är väldigt okunnig om rätt terminologi i dessa saker och jag är inte säker på vad du menar med samplingsfrekvens. Jag säger bara vad jag gör. Först genererar jag en slumpmässig vägprofil som har rumsfrekvenser varierande från 0,01 till 4 cyclesm. ISO 8608: 1995 har klassificering av väg och beroende på detta, ger det ett PSD-värde för varje frekvens mellan 0,01 och 4 som du vill ha. Dessa värden läggs sedan in i en ekvation för väggenerering som skapar en väg med några punkter (i mitt fall 8000 eller 400 meter, dvs varje 0,05 meter). Jag graver sedan alla ISO-värden för PSD mot de rumsliga frekvenserna som jag hade ovan. Jag försöker sedan arbeta bakåt för att se om jag kan generera samma graf genom att använda samma vägprofil och hitta FFT av den och sedan PSD. Jag vet inte vad du menar med samplingsfrekvens Jag är rädd, kanske är den där uppe i det jag har beskrivit Tack så mycket för din tid, jag är helt som en fisk ur vatten på den här Charlie skrev: gt gt gt Katherine, gt gt Det låter som om du kan filtrera data redan när du anger gt frekvensintervallet. Vad är samplingsfrekvensen Är det rumsligt eller gt temporal gt Om du specificerar 4 cyklarsmätare till systemet är det mycket osannolikt att det bara skulle vara provtagning för att få den frekvensen (Fs18 meter) utan att någon gt av gt moving average filter built I. Gt gt Vad är ISO-kravet (ISO-standard, varifrån) gt gt En effekt av filtreringen är att flytta energin till lägre gt frekvenser istället för att bara haka av den som du skulle göra i gt gt-domänen. Om slutmålet är att beräkna en IRI eller någon gt slags gt annan väggrävhet metrisk än detta kan vara kritisk. Gt gt Charlgt gt gtgt gtgt Efter förbehandlingsfiltret frågar det att jag utför en FFT gt som gtgt jag antar är också ett svar på din nästa fråga. Det stora problemet med problematiken som jag har är att jag genererade GTgt-profilen själv och angav att jag ville att frekvenserna skulle vara minst 0,01cykelsmeter och en max 4cyclesmeter. Varför då gtgt ska jag behöva filtrera ut högfrekvenser gtgt gt gt Tack för informationen om ISO 8608: 1995 ser det ut som en bra referens för några av mitt arbete med bearbetning av vägprofil. Tillbaka till ditt projekt. Som jag förstår det gör du: 1. Skapa vägprofil i rymdfrekvensdomänen med innehåll i 0,01-4 cyclesm 2. Generera rumslig profil från 1 med hjälp av vissa ekvationer (400 meter lång, dx0.05 m, rumslig samplingsfrekvens1dx20 cyclesm) 3 Gradera din PSD-väg från 1 mot ISO-värdena från ISO 8608 4. Beräkna fft och PSD från 2 och jämföra den med 3 för att se om du kan återproducera den. Om detta är korrekt och jag förstår ISO-standarden. Jag tror inte att du behöver göra någon filtrering alls. Din profil från 2 ska kunna generera frekvensdata från 0.0025-10 cyclesm, men du får inte se något innehåll över 4 cyclesm. Hoppas det hjälper till snarare än att förvirra. Du kanske vill titta på The Little Book of Profiling på umtri. umich. eduerdroughnessindex. html för mer info. Katherine ltkatherine. cashellucd. iegt skrev i meddelande nyheter: ef02d7a.7webx. raydaftYaTP. gt Charlie, gt Jag är mycket okunnig på rätt terminologi i dessa saker och jag är inte säker på vad du menar med samplingsfrekvens. Jag säger bara vad jag gör. gt gt gt Först genererar jag en slumpmässig vägprofil som har rumsliga gt frekvenser varierande från 0,01 - 4 cyclesm. ISO 8608: 1995 har gt-klassificeringar av vägen och beroende på detta ger det ett PSD-värde gt för var och en av frekvenserna mellan 0,01 och 4 som du vill ha. Dessa gt-värden läggs sedan in i en ekvation för väggenerering som gt skapar en väg med några punkter (i mitt fall 8000 eller 400meters, dvs varje 0,05 meter). gt Jag graver sedan alla ISO-värden för PSD mot de rumsliga gt frekvenserna som jag hade ovan. gt Jag försöker sedan arbeta bakåt för att se om jag kan generera samma gt graf genom att använda samma vägprofil, och hitta FFT av den och gt sedan PSD. gt jag vet inte vad du menar med samplingsfrekvens Jag är rädd, kanske är det där uppe i det jag har beskrivit gt gt Tack så mycket för din tid, jag är helt som en fisk ute av detta vatten gt Katherine gt Tack för det - är verkligen användbart bara för att se rätt terminologi som används för siffrorna Charlie skrev: gt gt gt Katherine, gt gt Tack för informationen om ISO 8608: 1995 ser det ut som en bra referens för några av mitt arbete på vägprofil bearbetning. Tillbaka till ditt projekt. När jag förstår det gör du: gt gt 1. Skapa vägprofil i rumsfrekvensdomänen med innehåll i gt 0.01-4 gt cyclesm gt 2. Generera rumslig profil från 1 med hjälp av vissa ekvationer (400 gt meter lång, gt dx0. 05 m, Spatial samplingsfrekvens1dx20 cyclesm) gt 3. Grafik din väg PSD från 1 mot ISO-värdena från ISO gt 8608 gt 4. Beräkna fft och PSD från 2 och jämföra den med 3 till gt se om du kan Återproducera det. gt gt Om detta är korrekt och jag förstår ISO-standarden. Jag tror inte att du måste göra någon filtrering alls. Din profil från 2 bör gt kunna generera frekvensdata från 0.0025-10 cyclesm, men du får inte se något gt-innehåll över 4 cyclesm. Gt Hoppas detta hjälper snarare än att förvirra. Du kanske vill titta på The gt Little gt bok av profiling på ltumtri. umich. eduerdroughnessindex. html gt gt gt eller mer info. Gt gt Charlie gt gt Katherine ltkatherine. cashellucd. iegt skrev i meddelande gt news: ef02d7a.7webx. raydaftYaTP. Gtgt Charlie, gtgt Jag är mycket okunnig om rätt terminologi i dessa saker och jag är inte säker på vad du menar med samplingsfrekvens. Jag säger bara vad jag gör. Gtgt gtgt gtgt Först genererar jag en slumpmässig vägprofil som har rumsliga gtgtfrekvenser varierande från 0,01 - 4 cyclesm. ISO 8608: 1995 gt har gtgt klassificeringar av vägen och beroende på detta ger det ett PSD gt värde gtgt för var och en av frekvenserna mellan 0,01 och 4 som du vill ha. gt Dessa gtgt-värden läggs sedan in i en ekvation för väggenerering som gtgt skapar en väg med några punkter (i mitt fall 8000, eller gtgt 400meters, dvs varje 0,05 meter). Gtgt Jag grafar sedan alla ISO-värden för PSD mot gt-rumsliga gtgt-frekvenserna som jag hade ovan. Gtgt Jag försöker sedan arbeta bakåt för att se huruvida jag kan generera den gt samma gtgtgrafen genom att använda samma vägprofil och hitta FFT av det gt och gtgt sedan PSD. gtgt jag vet inte vad du menar med samplingsfrekvens Jag är rädd, kanske är det gtgt där uppe i det jag har beskrivit gtgt gtgt Tack så mycket för din tid, jag är helt som en fisk ute av gtgt vatten på den här gtgt Gtgt Katherine gtgt gt gt gt Vad är en klocklista Du kan tänka på din klocklista som trådar som du har bokmärkt. Du kan lägga till taggar, författare, trådar och till och med sökresultat till din bevakningslista. På så sätt kan du enkelt hålla reda på ämnen som du är intresserad av. För att se din tittellista, klicka på quotMy Newsreaderquot-länken. Om du vill lägga till objekt i din bevakningslista klickar du på citadeln för att titta på listotiklänk längst ner på en sida. Hur lägger jag till ett objekt i min vaktlista För att lägga till sökkriterier i din vaktlista, sök efter önskad term i sökrutan. Klicka på citatetLägg till den här sökningen i min klocklistor på länken på sökresultatsidan. Du kan också lägga till en tagg i din bevakningslista genom att leta efter taggen med direktivet quottag: tagnamequot där tagname är namnet på den tagg du vill titta på. För att lägga till en författare till din bevakningslista, gå till författarens profilsida och klicka på quotAddjär den här författaren till min watch listquot-länk högst upp på sidan. Du kan också lägga till en författare till din bevakningslista genom att gå till en tråd som författaren har publicerat och klicka på quotAdd denna författare till min watch listquot link. Du kommer att få besked när författaren gör ett inlägg. Om du vill lägga till en tråd i din bevakningslista går du till trådsidan och klickar på citatetLägg till den här tråden i min larmlista-länk högst upp på sidan. Om nyhetsgrupper, nyhetsläsare och MATLAB Central Vad är nyhetsgrupper Nyhetsgrupperna är ett globalt forum som är öppet för alla. Nyhetsgrupper används för att diskutera ett stort antal ämnen, göra meddelanden och handelsfiler. Diskussionerna är gängade eller grupperade på ett sätt som låter dig läsa ett upplagd meddelande och alla dess svar i kronologisk ordning. Detta gör det enkelt att följa tråden i samtalet och för att se vad du redan har sagt innan du skickar ditt eget svar eller gör ett nytt inlägg. Nyhetsgruppens innehåll distribueras av servrar som är värd av olika organisationer på Internet. Meddelanden utbyts och hanteras med hjälp av öppna standardprotokoll. Ingen enskild enhet ldquoownsrdquo newsgroups. Det finns tusentals nyhetsgrupper som varje adresserar ett enda ämne eller intresseområde. MATLAB Central Newsreader postar och visar meddelanden i comp. soft-sys. matlab-nyhetsgruppen. Hur läser jag eller postar till nyhetsgrupperna Du kan använda den integrerade nyhetsläsaren på MATLAB Central webbplats för att läsa och skicka meddelanden i den här nyhetsgruppen. MATLAB Central är värd MathWorks. Meddelanden som skickas via MATLAB Central Newsreader ses av alla som använder nyhetsgrupper, oavsett hur de kommer åt nyhetsgrupperna. Det finns flera fördelar med att använda MATLAB Central. Ett konto Ditt MATLAB Central-konto är knutet till ditt MathWorks-konto för enkel åtkomst. Använd E-postadressen till ditt val MATLAB Central Newsreader gör att du kan definiera en alternativ e-postadress som din postadress, för att undvika röran i din primära brevlåda och minska spam. Spamkontroll De flesta nyhetsgruppspamfiler filtreras ut av MATLAB Central Newsreader. Taggningsmeddelanden kan märkas med en relevant etikett av någon inloggad användare. Taggar kan användas som nyckelord för att hitta specifika filer av intresse eller som ett sätt att kategorisera dina bokmärkta inlägg. Du kan välja att låta andra se dina taggar, och du kan visa eller söka otherrsquo-taggar såväl som de i samhället som helhet. Tagging ger ett sätt att se både de stora trenderna och de mindre, mer dunkla idéerna och applikationerna. Vaktlistor Med inställning av vallistor kan du få meddelande om uppdateringar gjorda till inlägg som valts av författare, tråd eller någon sökvariabel. Meddelandena om bevakningslistan kan skickas via e-post (dagligen digest eller omedelbar), visas i My Newsreader, eller skickas via RSS-flöde. Andra sätt att komma åt nyhetsgrupperna Använd en nyhetsläsare via din skola, arbetsgivare eller internetleverantör Betala för nyhetsgruppen tillgång från en kommersiell leverantör Använd Google Groups Mathforum. org ger en nyhetsläsare med tillgång till comp. soft sys. matlab-nyhetsgruppen Kör din egen server. För typiska instruktioner, se: slyckng. phppage2 Välj ditt CountryMoving Average Filter (MA filter) Laddar. Det rörliga genomsnittliga filtret är ett enkelt filter med lågt pass FIR (Finite Impulse Response) som vanligtvis används för att utjämna en rad samplade datasignaler. Det tar M prover av ingång i taget och tar medeltalet av de M-proverna och producerar en enda utgångspunkt. Det är en väldigt enkel LPF-struktur (Low Pass Filter) som är användbar för forskare och ingenjörer för att filtrera oönskade bullriga komponenter från de avsedda data. När filterlängden ökar (parametern M) ökar utjämnets jämnhet, medan de skarpa övergångarna i data görs alltmer stumma. Detta innebär att detta filter har utmärkt tidsdomänsvar men ett dåligt frekvenssvar. MA-filtret utför tre viktiga funktioner: 1) Det tar M-ingångspunkter, beräknar medelvärdet av de M-punkterna och producerar en enda utgångspunkt 2) På grund av beräknade beräkningskalkyler. Filtret introducerar en bestämd mängd fördröjning 3) Filtret fungerar som ett lågpassfilter (med dåligt frekvensdomänsvar och ett bra domänsvar). Matlab-kod: Efter matlab-kod simuleras tidsdomänsvaret för ett M-punkts rörande medelfilter och avbildar även frekvensresponsen för olika filterlängder. Tid Domain Response: På den första tomten har vi inmatningen som går in i det glidande medelfiltret. Inmatningen är bullriga och vårt mål är att minska bruset. Nästa figur är utgångsvaret för ett 3-punkts rörande medelfilter. Det kan härledas från figuren att 3-punkts rörande medelfilter inte har gjort mycket för att filtrera ut bruset. Vi ökar filterkranarna till 51 punkter och vi kan se att bruset i utmatningen har minskat mycket, vilket avbildas i nästa bild. Vi ökar kranarna vidare till 101 och 501 och vi kan observera att även om bullret är nästan noll övergår övergångarna drastiskt (observera lutningen på vardera sidan av signalen och jämföra dem med den ideala tegelväggsövergången i vår ingång). Frekvensrespons: Från frekvenssvaret kan man hävda att avrullningen är väldigt långsam och stoppbandets dämpning inte är bra. Med tanke på detta stoppbanddämpning kan det glidande medelfiltret helt klart inte separera ett frekvensband från en annan. Som vi vet att en bra prestanda i tidsdomänen leder till dålig prestanda i frekvensdomänen och vice versa. Kortfattat är det glidande medelvärdet ett exceptionellt bra utjämningsfilter (åtgärden i tidsdomänen), men ett exceptionellt dåligt lågpassfilter (åtgärden i frekvensdomänen) Externa länkar: Rekommenderade böcker: Primär sidofältMovande medelfilter (MA filter ) Läser in. Det rörliga genomsnittliga filtret är ett enkelt filter med lågt pass FIR (Finite Impulse Response) som vanligtvis används för att utjämna en rad samplade datasignaler. Det tar M prover av ingång i taget och tar medeltalet av de M-proverna och producerar en enda utgångspunkt. Det är en väldigt enkel LPF-struktur (Low Pass Filter) som är användbar för forskare och ingenjörer för att filtrera oönskade bullriga komponenter från de avsedda data. När filterlängden ökar (parametern M) ökar utjämnets jämnhet, medan de skarpa övergångarna i data görs alltmer stumma. Detta innebär att detta filter har utmärkt tidsdomänsvar men ett dåligt frekvenssvar. MA-filtret utför tre viktiga funktioner: 1) Det tar M-ingångspunkter, beräknar medelvärdet av de M-punkterna och producerar en enda utgångspunkt 2) På grund av beräknade beräkningskalkyler. Filtret introducerar en bestämd mängd fördröjning 3) Filtret fungerar som ett lågpassfilter (med dåligt frekvensdomänsvar och ett bra domänsvar). Matlab-kod: Efter matlab-kod simuleras tidsdomänsvaret för ett M-punkts rörande medelfilter och avbildar även frekvensresponsen för olika filterlängder. Tid Domain Response: På den första tomten har vi inmatningen som går in i det glidande medelfiltret. Inmatningen är bullriga och vårt mål är att minska bruset. Nästa figur är utgångsvaret för ett 3-punkts rörande medelfilter. Det kan härledas från figuren att 3-punkts rörande medelfilter inte har gjort mycket för att filtrera ut bruset. Vi ökar filterkranarna till 51 punkter och vi kan se att bruset i utmatningen har minskat mycket, vilket avbildas i nästa bild. Vi ökar kranarna vidare till 101 och 501 och vi kan observera att även om bullret är nästan noll övergår övergångarna drastiskt (observera lutningen på vardera sidan av signalen och jämföra dem med den ideala tegelväggsövergången i vår ingång). Frekvensrespons: Från frekvenssvaret kan man hävda att avrullningen är väldigt långsam och stoppbandets dämpning inte är bra. Med tanke på detta stoppbanddämpning kan det glidande medelfiltret helt klart inte separera ett frekvensband från en annan. Som vi vet att en bra prestanda i tidsdomänen leder till dålig prestanda i frekvensdomänen och vice versa. Kortfattat är det rörliga genomsnittet ett exceptionellt bra utjämningsfilter (åtgärden i tidsdomänen), men ett exceptionellt dåligt lågpassfilter (åtgärden i frekvensdomänen) Externa länkar: Rekommenderade böcker: Primär sidofältBetagen onsdagen den 08 oktober 2008 20:04 Senast uppdaterad den 14 mars 2013 01:29 Skriven av Batuhan Osmanoglu Hits: 41505 Flyttande medelvärde I Matlab Ofta befinner jag mig själv att behöva medelvärda data som jag måste minska bullret lite. Jag skrev några funktioner för att göra exakt vad jag vill, men matlabs inbyggda filterfunktion fungerar ganska bra också. Här skriver jag om 1D och 2D-medelvärde för data. 1D-filter kan realiseras med hjälp av filterfunktionen. Filterfunktionen kräver minst tre ingångsparametrar: täljarkoefficienten för filtret (b), nämnarkoefficienten för filtret (a) och data (X) förstås. Ett löpande medelfilter kan enkelt definieras av: För 2D-data kan vi använda Matlabs filter2-funktionen. För mer information om hur filtret fungerar kan du skriva: Här är en snabb och smutsig implementering av ett 16 med 16 glidande medelfilter. Först måste vi definiera filtret. Eftersom allt vi vill ha är lika stort bidrag från alla grannar kan vi bara använda de funktionerna. Vi dela allt med 256 (1616) eftersom vi inte vill ändra signalens generella nivå (amplitud). För att applicera filtret kan vi helt enkelt säga följande Nedan visas resultaten för fas av ett SAR-interferogram. I detta fall är området i Y-axeln och Azimuth är mappad på X-axeln. Filtret var 4 pixlar brett i räckvidd och 16 pixlar vid Azimuth.
No comments:
Post a Comment