convolutie in T domein en vermenigvuldiging in F-domein

J

J_expoler2

Guest
hoi
i studeren in DSP maar ik heb het probleem?
- Convolutie in de tijd domein en vermenigvuldiging in frequentiedomeinen
Ex 1
a = [1 4 3 2];
b = [0,5 2 1 3];
ans = conv (a, b);
ans = 2,0000 8,5000 7,5000 20,0000 10,0000 11,0000 6,0000en Ex2
a = [1 4 3 2];
b = [0,5 2 1 3];
A = fft (a);
B = fft (b);
ans = ifft (A * B);
ans = 12,0000 19,5000 13,5000 20,0000

ex1 waarom niet gelijk Ex2?
Met vriendelijke groet.

 
hai,

post toegevoegd ...door mezelf ...
FFT voor convolutie vereist nul opvulling zoals gezegd door de onder 2 leden ...i mistook ur questiion ...en beantwoordt niet snel een antwoord ...te laat om te beantwoorden ..maakt niet uit ..i couldnot verwijder dit bericht tooo, dus ik heb gevoegd ...goede uitleg hieronder te volgen ...

groeten,
Arun.
Laatst gewijzigd door cedance op 24 februari 2004 19:54; bewerkt 1 keer in totaal

 
Wanneer convoluting a en b de lengte van de gecompliceerde signaal is:
lengte (a) lengte (b) - 1
Het originele signaal en het filter een impuls respons b niet periodieke signalen.

Bij de omvorming van A en B bij frequentiedomeinen je moeten weten wat je doet.De Fourier coëfficiënten vertegenwoordigen perodic tijd signalen.Bij de omvorming van een = [1 4 3 2] om Freqeuncy domein wordt automatisch betekent dat je knippen een periode van een moment als dit signaal: een '= [1 4 3 2 4 1 3 2 4 1 2 3 ..... ..].Als u wilt dat de juiste weergave van het signaal dat u de tijd hebt tot nul toevoegen.Zoals: een''= [1 4 3 2 0 0 0 0 0 0 ....] aan een aantal van 2 ^ n waarden.Door dit u bovendien een toename van de nauwkeurigheid van de frequentie resolutie.

Bij de omvorming van dit een''te frequentiedomeinen de amplitude zal veranderen.Dat
is de reden waarom u te delen door het aantal monsters omgezet.

Probeer het opnieuw met nul voegen.

 
Aoxomox is goed, u moet toevoegen nullen te krijgen waar de berekening,
als een minimale lengte van a en b, moet gelijk zijn aan de lengte van de convolutie.

lengte van convolutie = lengte (a) lengte (b) -1

Code:

a = [1 4 3 2 0 0 0];

b = [0,5 2 1 3 0 0 0];

A = fft (a);

B = fft (b);

ans = ifft (A * B)
 

Welcome to EDABoard.com

Sponsor

Back
Top