Kolozzeum Forum - Sveriges största träningsforum
Gå tillbaka   Kolozzeum Forum - Sveriges största träningsforum > Träning > Allmänt Forum

Svara
 
Ämnesverktyg Visningsalternativ
Gammal 2008-01-11, 13:32   #1
Palm
Registered User
 
Palms avatar
 
Reg.datum: May 2005
Ort: Strax ovanför LHC
Inlägg: 315
Mattehjälp!

Behöver lite hjälp av er matematiker eller signalbehandlare.

Jag har en samplad signal u[k] som med god approximation kan representeras av en annan (kortare) signal r[k], tidsskiftad och linjärt skalad på formen:

u[k] = sum(i=1:P){lambda_i * r[k-t_i]} + E[k]
där r[k] är en lämplig referensfunktion som innehåller en topp och ringningen efter toppen. (P=16 i det här fallet, E = differens mellan fit och samplad signal)

""

Som synes dränks transienten från första excitationen av nästa osv. Det jag vill få fram är bästa* valet av lambda_i, t_i, och "single excitation response" r[k].
Det jag inte vill göra är att försöka filtrera bort "oscillationerna" mellan topparna. Har ganska många signaler av den här typen och frekvenskarakteristiken på ringningen ser olika ut från gång till gång.

Knåpade ihop en MonteCarlo-algoritm som mer eller mindre slumpar sig fram till lämpliga val av r,lambda_i och t_i vilket funkar bra, men det går ganska långsamt och ger inga garantier att man inte landat i ett lokalt, icke-globalt, optimum (även om man kan minimera den risken genom heroiska konvergensvillkor och dynamiska slumpintervall)

Funderar på att kasta in allt i fourierrummet, men det blir rätt grisigt där också...

Tipsa gärna om du sitter på något lämpligt angreppssätt. Systemet antas alltså vara linjärt och jag vill få fram en hyfsad uppskattning av "impulssvaret" (inkommande signaler är iofs för långa för att kallas impulser, men vi kan kalla dem det iaf) ur signalen ovan som är responsen för flera impulser som endast skiljer sig åt i amplitud och starttid.

*Nån minsta-kvadrat duger, även om det vore önskvärt att skala felstraffning med signalamplitud i någon form (viktigare att ha bra fit kring huvudtopparna).

(Nej, det är inte läxa )
Palm är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 13:42   #2
Muskelbyggaren
Löjligt internskämt
 
Muskelbyggarens avatar
 
Reg.datum: Sep 2006
Inlägg: 206
Tyvärr kan jag för lite om transformer och sånt för att kunna hjälpa dig.

Får man fråga om vad det är i för sammanhang du ska göra detta?
Muskelbyggaren är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 13:45   #3
jakke
Registered User
 
Reg.datum: Aug 2002
Ort: Kungsholmen
Inlägg: 3 009
Pinsammt att jag typ glömt allt från mina signalbehandlingskurser. Men en simpelt winerfilter löser inte biffen då?
jakke är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 13:50   #4
The_RobRoy
Living the good life
 
The_RobRoys avatar
 
Reg.datum: Jul 2004
Ort: Hongkong
Inlägg: 1 788
Vadå är inte spikarna tillräckligt tydliga som dom är, det är väl dom du är intresserad av?
__________________
Iron addict for life
The_RobRoy är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 13:54   #5
Palm
Registered User
 
Palms avatar
 
Reg.datum: May 2005
Ort: Strax ovanför LHC
Inlägg: 315
Citat:
Ursprungligen postat av Muskelbyggaren Visa inlägg
Tyvärr kan jag för lite om transformer och sånt för att kunna hjälpa dig.

Får man fråga om vad det är i för sammanhang du ska göra detta?
Oki. Det borde nog gå att lösa med massiv matrisalgebra. Frågan är om det går att få en entydig lösning... Hmmm...

Signalen i bilden är från en current transformer och motsvarar (idealt) intensiteten i en protonstråle över tid några meter innan den brakar in i en jetstråle av kvicksilver. Av någon anledning är det en massa skräp mellan topparna vilket ställer till det när jag ska normera signalerna från sekundärpartikeldetektorerna mot primär intensitet.
Palm är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 13:59   #6
Palm
Registered User
 
Palms avatar
 
Reg.datum: May 2005
Ort: Strax ovanför LHC
Inlägg: 315
Citat:
Ursprungligen postat av The_RobRoy Visa inlägg
Vadå är inte spikarna tillräckligt tydliga som dom är, det är väl dom du är intresserad av?
Njae... Jag är rätt kräsen
I det här fallet t.ex. så hamnar andra toppen direkt över ett positivt rippel från den första => Överskattning av lambda_2. Vill ha en precision på sisådär 2 procent. Ibland är rippelamplituden så hög som 1/3 av huvudtopparna. Dessutom tenderar hela signalen att successivt bli mer och mer negativ pga transienterna vilket omöjliggör en enkel intervallintegrering.
Palm är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 14:00   #7
Palm
Registered User
 
Palms avatar
 
Reg.datum: May 2005
Ort: Strax ovanför LHC
Inlägg: 315
Citat:
Ursprungligen postat av jakke Visa inlägg
Pinsammt att jag typ glömt allt från mina signalbehandlingskurser. Men en simpelt winerfilter löser inte biffen då?
Tack, ska kika på det. Har aldrig använd ett sånt
Palm är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 14:02   #8
Snafu
Registered User
 
Snafus avatar
 
Reg.datum: Mar 2003
Ort: Newcastle upon Tyne, UK
Inlägg: 453
Goddamn hade jag varit brud och tänt på geek-talk hade jag nog självexploderat efter den här tråden.
__________________
Sy sy sy sy sy sy, tom rulle!
Snafu är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 14:06   #9
jakke
Registered User
 
Reg.datum: Aug 2002
Ort: Kungsholmen
Inlägg: 3 009
Citat:
Ursprungligen postat av Snafu Visa inlägg
Goddamn hade jag varit brud och tänt på geek-talk hade jag nog självexploderat efter den här tråden.
Finns det sådanna brudar i (R)england?
jakke är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 14:10   #10
Snafu
Registered User
 
Snafus avatar
 
Reg.datum: Mar 2003
Ort: Newcastle upon Tyne, UK
Inlägg: 453
Citat:
Ursprungligen postat av jakke Visa inlägg
Finns det sådanna brudar i (R)england?
Nej, dåligt med det tror jag...här ska man ha humor över allt annat. Gärna en del pengar också. Tur att jag fallerar på båda punkterna.
__________________
Sy sy sy sy sy sy, tom rulle!
Snafu är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 14:14   #11
Palm
Registered User
 
Palms avatar
 
Reg.datum: May 2005
Ort: Strax ovanför LHC
Inlägg: 315
Citat:
Ursprungligen postat av Snafu Visa inlägg
Goddamn hade jag varit brud och tänt på geek-talk hade jag nog självexploderat efter den här tråden.
"Talk geek to me!"
"Who's my little singular value decomposer?"
"More"
"I'm gonna convolute you... hard... all the way through Hilbert space"
Palm är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 14:38   #12
Scratch89
Klädsamt knubbig
 
Scratch89s avatar
 
Reg.datum: Dec 2004
Ort: Strömmensberg, Göteborg
Inlägg: 19 901
Citat:
Ursprungligen postat av Palm Visa inlägg
"Talk geek to me!"
"Who's my little singular value decomposer?"
"More"
"I'm gonna convolute you... hard... all the way through Hilbert space"
__________________
Leg. fysioterapeut.
Skribent på Styrkelabbet
Scratch89 är inte uppkopplad   Svara med citat
Gammal 2008-01-11, 16:12   #13
stridis
stridis- avatartextmästare
 
stridiss avatar
 
Reg.datum: Jan 2003
Ort: lindome
Inlägg: 6 559
*host* om det är någon tjej här som gillar nördsnack, då helst mattetermer alltifrån mångfalder till liealgebror, så pm:a mig så kan vi nog bestämma en träff.
__________________
Stridis invictus.

I won't just gaze into the abyss, I'll maintain uncomfortable eye contact with it.
stridis är inte uppkopplad   Svara med citat


Svara


Regler för att posta
Du får inte posta nya ämnen
Du får inte posta svar
Du får inte posta bifogade filer
Du får inte redigera dina inlägg

BB-kod är
Smilies är
[IMG]-kod är
HTML-kod är av

Forumhopp


Alla tider är GMT +1. Klockan är nu 23:23.


Powered by vBulletin® Version 3.8.8
Copyright ©2000 - 2025, Jelsoft Enterprises Ltd.
Svensk översättning av: Anders Pettersson
Copyright 2015, Kolozzeum (TM)