Springen naar inhoud

Ik zoek een 2d-fft routine in basic


  • Log in om te kunnen reageren

#16

Wimapon

    Wimapon


  • >100 berichten
  • 114 berichten
  • Ervaren gebruiker

Geplaatst op 11 mei 2011 - 22:35

Hoi 317070

Allereerst toen in deze figuur op het internet zette,had ik de kaart er nog niet in gezien.
Als je het kwadrant rechtsboven rechts/links omklapt en dan in het linkeronder kwadrant zet,
en het kwadrant linksboven horizontaal gezien omklapt en dan naar kwadrant rechtsboven
verplaats, krijg je een hemelkaart.
alleen de hoekpunten moet je dan weglaten, want die bevatten geen zinnige punten aan de hemelbol.
Dit doe ik dus voortaan in alle kaartjes die ik presenteer.

Nu je vragen:
1: de pixels linksboven krijg je altijd als je natuurlijke data door een fft haalt. men noemt dit de
gelijkstroom waarde. als de wisselspanning zijn nullijn precies op 0 volt heeft, krijg je deze waarde niet als
uitschieter.

2:De boel is inderdaad 2zijdig symetrisch. Als ik de boel omklap laat ik het symetrische deel ook weg.
Een ster in het Noordwesten komt op dezelfde plek als een ster in het zuidoosten.

3: de buine vraagteken staat in de zg. locatiekaart. Ik heb op de grond een grid met afstanden van 4 meter in
het gras hier uitgezet. Op elk punt van dit grid moet de tweede antenne een keer komen te staan.
De vaste staat links boven op het eerste punt van het grid.
De plekken waar de tweede antenne heeft gestaan heb ik met een kleurtje aangegeven.
Dit kleurtje is zo gekozen dat het de waarde van de op die plek gemeten cross-correlatie aangeeft.
Dus rood is een hoge waarde en blauw is een lage waarde.
De schaal staat er rechts van.
Theoretisch zou je het gehele grid van 16 bij 16 moeten opvullen.. maar daar kom ik niet aan toe want het
zou dan 16 maal 16 dagen meten zijn....
het blijkt dat je met 4 meetpunten
al een kaartje krijgt.. hoemeer punten hoe scherper het kaartje wordt.
Dit zie je hogerop in de pagina.

4: de resolutie kan niet groter zijn dan het meetgrid.(denk ik) Ik koos als eerste test voor een grid van 16 bij 16 met 4 m
tussen de punten.... dat leek me redelijk.... ik wil later nog andere gridspacings gaan proberen.....
in theorie: hoe groter de afstanden tussen de antennes, hoe groter de resolutie.
strakjes kies ik een grid van 32 bij 32 punten... dan meet ik er precies tussen in....
Of ik start opnieuw met antenne afstanden zo rond de 150 meter. ( nu is dat ongeveer max 50 m)
Ik moet dan de grid spacing ook kleiner maken..tis nog even zoeken en testen...
in elk geval kun je met een antenne afstand van 150 meter de positie van een ster op 1 graad nauwkeurig
bepalen.. ik weet nog niet hoe ik dat met een grid zou moeten aanpakken.. ik kan geen grid van
128 bij 128 punten gaan meten,,,, ik denk hier natuurlik wel aan!... een deel van dit grid meten is een
optie, maar ik moet hier nog eens goed over nadenken.

5: ik zou het heel leuk vinden om je een dataset toe te spelen.
ik moet daar eerst een programma voor maken, want ik haal de data automatisch uit alle meetfiles....en dat
is nog al een gedoe... maar ik ga dat zeker voor je doen.. ik ben heel benieuwd wat jij ermee kunt doen.
je krijgt dan een file die in mijn figuur het locatiekaartje geeft.
Ik zal je er eentje geven die bij mij een mooie kaart geeft.
het wordt een rij van 256 getallen... voor elk punt van het grid 1 getal... dus voornamelijk nullen, want daar is niet gemeten.


Lang verhaal

Groetjes
Wim

Veranderd door Wimapon, 11 mei 2011 - 22:49


Dit forum kan gratis blijven vanwege banners als deze. Door te registeren zal de onderstaande banner overigens verdwijnen.

#17

Kaspace

    Kaspace


  • >100 berichten
  • 202 berichten
  • Ervaren gebruiker

Geplaatst op 11 mei 2011 - 22:42

Zoek eens op "BBC Basic" en "FFT".

Er circuleren goede FFT demos in BBC basic.

De oude BBC Basic is geschikt gemaakt voor de moderne PC. Bij meneer R.T. Russell kan je een trial versie van BBC Basic for Windows downloaden of de full version kopen voor weinig geld.

Deze Basic is helemaal op maat gemaakt voor de PC met flink wat kritische delen in directe X86-assembler geschreven. Er is ontzettend veel mooie documentatie, je kan goed zien dat de BBC Basic een serieuze programmeertaal was/is. De huidige ARM processors waren bedoeld om met BBC - Acorn - Archimedes vanuit Engeland Intel/Microsoft/Apple naar de kroon te steken.

Run analyser.exe
of
Ga naar de website en zoek op "Audio Spectrum Analyser"

(demo is direct uitvoerbaar als .exe)

Veranderd door Kaspace, 11 mei 2011 - 22:52


#18

Wimapon

    Wimapon


  • >100 berichten
  • 114 berichten
  • Ervaren gebruiker

Geplaatst op 11 mei 2011 - 22:51

Ik ga kijken, want ik programmeerde vroeger op een BBC... in de basic ervan en in de assembler ervan...grin

inmiddels heb ik het snel bekeken... tis een 1d-fft.....maar ik kijk nog verder.

Veranderd door Wimapon, 11 mei 2011 - 22:55


#19

Wimapon

    Wimapon


  • >100 berichten
  • 114 berichten
  • Ervaren gebruiker

Geplaatst op 12 mei 2011 - 12:49

Ik heb de set data klaar.... voor zowel de input van de 2d-FFT als de output van de 2d-FFT.

Hoe kan ik ze aanleveren??

Wim

#20

317070

    317070


  • >5k berichten
  • 5613 berichten
  • Moderator

Geplaatst op 14 mei 2011 - 11:26

Ik heb de set data klaar.... voor zowel de input van de 2d-FFT als de output van de 2d-FFT.

Hoe kan ik ze aanleveren??

Dankjewel.

Euhm, hoe groot zijn die bestand(en)?
En hoe groot zijn ze als je ze in 1 bestandje comprimeert met 7-zip?
What it all comes down to, is that I haven't got it all figured out just yet
And I've got one hand in my pocket and the other one is giving the peace sign
-Alanis Morisette-

#21

Wimapon

    Wimapon


  • >100 berichten
  • 114 berichten
  • Ervaren gebruiker

Geplaatst op 14 mei 2011 - 12:52

hoi 317070
De beide bestanden bestaan uit 256 getallen.
Ik zou ze hier onder elkaar kunnen zetten.... werkt dat?

Groetjes
Wim

#22

317070

    317070


  • >5k berichten
  • 5613 berichten
  • Moderator

Geplaatst op 14 mei 2011 - 12:54

hoi 317070
De beide bestanden bestaan uit 256 getallen.
Ik zou ze hier onder elkaar kunnen zetten.... werkt dat?

Ah, wel, als je een antwoord op mijn bericht typt en een beetje naar beneden scrollt, dan zie je een venstertje bijlagen, waarmee je bestanden kunt uploaden en toevoegen in je bericht.
What it all comes down to, is that I haven't got it all figured out just yet
And I've got one hand in my pocket and the other one is giving the peace sign
-Alanis Morisette-

#23

Wimapon

    Wimapon


  • >100 berichten
  • 114 berichten
  • Ervaren gebruiker

Geplaatst op 14 mei 2011 - 14:15

hoi hier is de inputfile voor de 2d-fft.
het onderstaande plaatje geeft alle benodigde info... zowel de locatie van de getallen en het resultaat van de 2d-fft
in kaartvorm.

Ik geloof dat ik de rij getallen er niet als bijlage bij krijg.
In elk geval... het eerste getal staal in de locatiekaart linksboven.
Het tweede getal staat eronder. het derde getal daar weer onder....
HET laatste getal staat rechtsonder...
let op het ene getal dat zwART ( = 0.0000 ) is op het kaartje.. dit kun je als controle gebruiken.


hier het databestand

Hier de bijlage.....als het lukt...


Hier de bijlage.....als het lukt...

ik krijg het niet voor elkaar


hij doet het wel uploaden, maar ik zie het nergens meer....

pffff


Dan maar met de botte bijl:


0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
31591.7246093750
28479.2167968750
-7721.1274414063
-861.8062133789
10353.8027343750
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
-22403.6777343750
2060.5156250000
-2423.5686035156
7632.9648437500
52650.1210937500
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
29484.5800781250
48344.1523437500
45977.5742187500
0.0000000000
-919.7631835938
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
-9737.3720703125
-13119.2724609375
7155.3227539063
439.6993103027
33177.3007812500
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
20780.8867187500
12920.9296875000
24306.5546875000
21226.0625000000
13275.3154296875
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
23929.1503906250
11130.1845703125
15820.0000000000
-2682.6337890625
9610.9472656250
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
-13840.5761718750
-9342.3330078125
21622.7871093750
10353.8027343750
16472.1660156250
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
27780.8066406250
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000
0.0000000000


dit moet hem zijn...

Wim

Bijgevoegde miniaturen

  • scherm.jpg

Veranderd door Wimapon, 14 mei 2011 - 14:17


#24

317070

    317070


  • >5k berichten
  • 5613 berichten
  • Moderator

Geplaatst op 14 mei 2011 - 20:03

Mijn resultaat:
Geplaatste afbeelding

Die FFT2 werkt dus perfect ;)
Maar ik heb nog altijd ernstige twijfels dat dat de hemel voorstelt. :P
Mag ik vragen welke berekening je doet om die kaart te maken? Of waar je je algoritme vandaan haalt om van je antennedata naar de hemelkaart te gaan?

In ieder geval, mij lijkt dat je met 2 antennes voldoende hebt om een hemelkaart te maken, zonder dat je die moet verplaatsen...
What it all comes down to, is that I haven't got it all figured out just yet
And I've got one hand in my pocket and the other one is giving the peace sign
-Alanis Morisette-

#25

Wimapon

    Wimapon


  • >100 berichten
  • 114 berichten
  • Ervaren gebruiker

Geplaatst op 15 mei 2011 - 13:56

Hoi 317070
Jeetje wat goed dat jij zo maar even dezelfde soort kaartjes maakt als ik na weken programmeren!

Even niet alles tegelijk. Ik zal nu een kaartje maken van deze data zonder het resjuvvelen van de
blokjes tot een hemelkaart... dan kunnen we goed zien of we echt hetzelfde krijgen..

Ik zal je andere vragen daarna beantwoorden. Dat wordt een hele interessante discussie, want ik
ben op zoek naar een methode om met 2 niet verplaatsbare antennes een hemelkaart te maken...
Bij de profs noemen ze dat aperture synthese..ze gebruiken de draaing van de aarde voor de
verplaatsing...


spannend!!!

#26

Wimapon

    Wimapon


  • >100 berichten
  • 114 berichten
  • Ervaren gebruiker

Geplaatst op 15 mei 2011 - 15:41

Er gaat toch iets mis...

ziehier mijn presentatie van de uitkomst van mijn 2dfft
Dit ziet er toch anders uit dan het jouwe....

de grafiek linksboven zou gelijk aan de jouwe moeten zijn.....

Bijgevoegde miniaturen

  • zonderkaartmaaktruuk.jpg

#27

317070

    317070


  • >5k berichten
  • 5613 berichten
  • Moderator

Geplaatst op 15 mei 2011 - 16:11

Hmm, ik heb de nulfrequenties wel in het midden van mijn grafiek geplaatst, zou het kunnen dat jij die nog ergens anders staan hebt? Waar staat de DC-waarde op je linkergrafiek?
What it all comes down to, is that I haven't got it all figured out just yet
And I've got one hand in my pocket and the other one is giving the peace sign
-Alanis Morisette-

#28

Wimapon

    Wimapon


  • >100 berichten
  • 114 berichten
  • Ervaren gebruiker

Geplaatst op 15 mei 2011 - 17:18

De DC-waarde staat linksbovenaan. Dit is toch gelijk aan de 0 frequentie.....?
Maar ik ben mijn houwvast door al die tellertjes wat kwijtgeraakt....
Pas toen ik zag dat de diverse azimuths en elevaties van gesimuleerde bronnen zich begrijpelijk
over het kaartje bewogen, begon ik er in te geloven en besloot ik nergens meer aan te komen....

#29

317070

    317070


  • >5k berichten
  • 5613 berichten
  • Moderator

Geplaatst op 16 mei 2011 - 17:20

Geplaatste afbeelding

Dit is zonder de nulfrequentie in het midden te plaatsen.
Volgens mij heb je dus je 2 assen omgewisseld en zijn er nog een paar dingen omgewisseld. En voor mij gaat dit inderdaad een stuk sneller eenvoudiger, doordat ik gebruik kan maken van een licentie van de universiteit voor matlab, die anders stukken van mensen kost.
What it all comes down to, is that I haven't got it all figured out just yet
And I've got one hand in my pocket and the other one is giving the peace sign
-Alanis Morisette-

#30

Wimapon

    Wimapon


  • >100 berichten
  • 114 berichten
  • Ervaren gebruiker

Geplaatst op 16 mei 2011 - 17:56

Hoi 317070
Als ik mijn plaatje 90 graden linksom draai, krijg ik zo goed als hetzelfde plaatje als jij.

Mag ik dus concluderen dat mijn 2d-fft routine goed is??? ( ook al is hij wat geroteerd ofzoiets , want dat roteren maak
ik verderop in mijn programme weer goed.)

Wim





0 gebruiker(s) lezen dit onderwerp

0 leden, 0 bezoekers, 0 anonieme gebruikers

Ook adverteren op onze website? Lees hier meer!

Gesponsorde vacatures

Vacatures