Hjelp meg å løse dette på lang tid hjemsøkt problem ...

P

powersys

Guest
Dette problemet har vært uforglemmelige meg i flere måneder.Kindly råd.

La oss si at brukeren er gitt et datasett av signal som vist i Figur 1.
Anta at han ikke vet hvor signalet er opprettet (det vil si han ikke kjenner kodene i Part-A).
Han utfører FFT på signal (med del-B-koden) og resultatene han oppnår er gitt i figur 2 og tabell 1 (kun harmoniske 2-10 blir vist).
Clearly, finner han ut det er 4 store harmoniske i signalet, dvs. 1., 3., 7., og 9. oktober.
Selv om det er noen små verdier på 4de, 8., og 10., bestemmer han seg for å anta at størrelsen på disse harmoniske er null.
Han lærte at et signal kan være konstruert av en rekke sinusoids.
Siden han kjenner størrelsen av store harmoniske, tror han at han kan rekonstruere signalet bruker koden i Part-C (antar det er ingen Matlab ifft funksjon).
Men han snart finner ut rekonstruerte signalet (som vist i figur 3) er ikke lik som i figur 1..
Hvorfor?

Spørsmål:
[1] Er koden utført i del-C for å rekonstruere signalet korrekt?
[2] Er det informasjon (særlig Y (2:10) og YPHASE (2:10)) gitt i tabell 1 nyttig signal gjenoppbygging?Hvis ja, hvordan kan vi bruke dem riktig?

Mange takkCode:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% DEL-A: For å lage signal

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

x = (0:1:359);

Y1 = 2 * sin (x * pi/180);

y3 = 1 * sin (3 * x * pi/180 30/180 * pi);

Y7 = 0.5 * sin (7 * x * pi/180);

y9 = 0.1 * sin (9 * x * pi/180);

y = y1 y3 Y7 y9;figure (1);

plot (x, y); xlim ([0 360]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% DEL-B: For å beregne FFT av signalet

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

N = lengde (y);

Y = FFT (y) / N;YMAG = 2 * abs (Y);

YPHASE = vinkel (Y);

realY = real (Y);

imagY = forestille (Y);figure (2);

bar (YMAG (2:20));[Y (2:10) 'YMAG (2:10)' YPHASE (2:10) ']

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% DEL-C: Å rekonstruere signalet

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

rcy1 = YMAG (1 1) * sin (1 * x * pi/180);

rcy3 = YMAG (1 3) * sin (3 * x * pi/180);

rcy7 = YMAG (1 7) * sin (7 * x * pi/180);

rcy9 = YMAG (1 9) * sin (9 * x * pi/180);

rcy = rcy1 rcy3 rcy7 rcy9;figure (3);

plot (x, rcy);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
Hei
Det er noen feil i koden.Du tok ikke i kontoen i fase sikt.Følgende kode vil ta seg av det.Du kan utvide dette konseptet for mer harmoniske.
Håper det vil være nyttig:).
Code:%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% DEL-A: For å lage signal

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

x = (0:1:359);

Y1 = 2 * sin (x * pi/180);

y3 = 1 * sin (3 * x * pi/180 30/180 * pi);

Y7 = 0.5 * sin (7 * x * pi/180);

y9 = 0.1 * sin (9 * x * pi/180);

y = y1 y3 Y7 y9;figure (1);

plot (x, y); xlim ([0 360]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% DEL-B: For å beregne FFT av signalet

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

N = lengde (y);

Y = FFT (y) / N;Ycos = (Y (2: N / 2) Y (N: -1: N / 2 2));

Ysin = j * (Y (2: N / 2)-Y (N: -1: N / 2 2));%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% DEL-C: Å rekonstruere signalet

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

rcy1 = Ycos (1) * cos (1 * x * pi/180) Ysin (1) * sin (1 * x * pi/180);

rcy3 = Ycos (3) * cos (3 * x * pi/180) Ysin (3) * sin (3 * x * pi/180);

rcy7 = Ycos (7) * cos (7 * x * pi/180) Ysin (7) * sin (7 * x * pi/180);

rcy9 = Ycos (9) * cos (9 * x * pi/180) Ysin (9) * sin (9 * x * pi/180);

rcy = rcy1 rcy3 rcy7 rcy9;figure (2);

plot (x, rcy);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 
Hei MHanif, takk for svaret ditt.Jeg testet koden din og den arbeider stor!Men jeg forstår ikke noen deler i koden.For eksempler:

[1] I del B, kan det hende jeg vet hva som er hensikten med Ycos og Ysin?

[2] I del C, kan det hende jeg vet hvorfor cosinus-funksjonen brukes til å contruct signalet i stedet for sinusfunksjon?

[3] Er beregning basert på invers DFT formulering?

Mange takk

 
Hi Power sys (Mest sannsynligvis det ikke navn

<img src="http://www.edaboard.com/images/smiles/icon_biggrin.gif" alt="Very Happy" border="0" />

),
Denne formuleringen er basert på forholdet mellom komplekse eksponential-og cosinus og sinus form for Fourierrekker / Transform.Hvis du vet relasjonene som er tilstede i alle standard tekst og mye av det er også tilgjengelig på internett, så det er lett å følge.
Utgangspunktet for Euler Ligningen brukes for å utvide exp (j * w * t) som er cos (w * t) j * sin (w * t).Omorganisere koeffisientene vil gi koeffisientene av cosinus og sinus treff.Håper dere forstår ovennevnte uttalelser

<img src="http://www.edaboard.com/images/smiles/icon_smile.gif" alt="Smil" border="0" />

.
Endelig blir du bedt om sinus og cosinus vilkår, disse kommer alltid i representasjon av noe periodiske signal.Fasen term i sinus sikt introduserer to semestre dvs. Asin (w * t) Bcos (w * t).
Jeg vet ikke hvor mye er dette svaret nyttig for deg

<img src="http://www.edaboard.com/images/smiles/icon_smile.gif" alt="Smil" border="0" />

.

 
Jeg brukte 'powersys "som brukernavn, fordi på den tiden jeg registrert som medlem av Edaboard var jeg på kurs på kraftsystemet.Ha Ha ...

Alle dine svar ikke gir meg masse innganger.Men fremdeles har jeg ikke egentlig forstå følgende uttalelse:

"Den fasen begrep i sinus sikt introduserer to semestre, dvs Asin (w * t) Bcost (w * t)"

Jeg er ikke i stand til å visualisere hvordan fase i sinus sikt knyttet til de to treff.Min unnskyldning for min stakkars bakgrunn DFT.

Takk.

 
Hei
La oss se denne enkle setningen.
Code:

y3 = 1 * sin (3 * x * pi/180 30/180 * pi);
 
You must log in or register to reply here.

Welcome to EDABoard.com

Sponsor

Back
Top