Ik zoek een 2d-fft routine in basic

[Wimapon]

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

[Kaspace]

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)

[Wimapon]

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.

[Wimapon]

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

[317070]

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?

[Wimapon]

hoi 317070

De beide bestanden bestaan uit 256 getallen.

Ik zou ze hier onder elkaar kunnen zetten.... werkt dat?

Groetjes

Wim

[317070]

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.

[Wimapon]

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

[317070]

Mijn resultaat:



Die FFT2 werkt dus perfect

Maar ik heb nog altijd ernstige twijfels dat dat de hemel voorstelt.

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...

[Wimapon]

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!!!

[Wimapon]

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.....

[317070]

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?

[Wimapon]

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....

[317070]

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.

[Wimapon]

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