19.01.2013 Views

VBN - Aalborg Universitet

VBN - Aalborg Universitet

VBN - Aalborg Universitet

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Appendiks A - Transformation og anblok<br />

Indholdet i residualvektoren er residualet mellem x- og y-koordinaterne for punkterne i Model A og<br />

Model B efter transformationen:<br />

Side | 20<br />

[ ]<br />

r = rX rY rX rY rX rY rX rY<br />

1 1 2 2 3 3 4 4<br />

[ -0,22 -0,02 0,14 -0,04 0,24 0,01 -0,17 0,05] T<br />

r =<br />

Ovennævnte residualer er identiske med residualerne fra den lineære metode.<br />

Til 2D transformation med den ulineære metode har projektgruppen ligeledes udarbejdet et script i<br />

MATLAB. Dette script hedder D2_trans_sincos.m. Denne fil er på Bilags-CD’en i mappen Appendiks<br />

B. I scriptet er ovennævnte procedure foretaget, dog er der nogle forhold der her skal gøres opmærksom<br />

på i forbindelse med gennemløb af scriptet. Som med 2D transformation med den lineære<br />

metode hentes de to modeller ind i scriptet fra txt-filer ved navn modelA.txt og modelB.txt. Disse<br />

filer indeholder matricer med tre søjler, hvor første søjle er punktnummer, mens de øvrige to er<br />

henholdsvis x- og y-koordinater. Søjlerne i txt-filerne er adskilt af mellemrum. Disse txt-filer skal<br />

som tidligere beskrevet ligge i samme mappe som scriptet køres fra. For at genskabe eksemplet<br />

som er gennemgået i dette afsnit skal filerne modelA.txt og modelB.txt hentes fra mappen 2D transformation<br />

under mappen 2D koordinatfiler og placeres direkte under mappen Appendiks B på Bilags-CD’en<br />

inden scriptet gennemløbes.<br />

Rækkefølgen af punkterne i de to filer, modelA.txt og modelB.txt, skal, som med den forrige metode,<br />

være den samme.<br />

Inden scriptet til 2D transformation med den ulineære metode, D2_trans_sincos.m, gennemløbes<br />

skal scriptet med 2D transformation med den lineære metode, D2_trans_ab.m, gennemløbes. Dette<br />

skyldes at den foreløbige værdi for drejningen hentes fra x_for_p.txt, som genereres ved gennemløb<br />

af den lineære metode, D2_trans_ab.m.<br />

Rækkefølgen af rækkerne i A-matricen, b-vektoren og r-vektoren er lidt anderledes i scriptet end i<br />

den ovenfor gennemgåede teori. Den ændrede rækkefølge skyldes, at det programmeringsmæssigt<br />

er lettere. Rækkefølgen af de ændrede matricer/vektorer er identiske med rækkefølgen som i scriptet<br />

til 2D transformation med den lineære metode, som er beskrevet i afsnit 4.1 Lineær metode.<br />

En anden måde at beregne de partielt afledede værdier til A-matricen er at anvende Peter Cederholms<br />

script numafl.m. Scriptet foretager en numerisk approksimation af de partielt afledede af en<br />

funktion. Inputtet til scriptet er en funktion, hvortil de partielt afledede ønskes, samt foreløbige<br />

værdier til de variable i funktionen. Outputtet fra numafl.m er de afledede udtryk af transformationsparametrene<br />

og 0’te ordens leddet.<br />

Scriptet approksimerer værdierne af 0. og 1. ordens afledede. De 0. ordens afledede anvendes i bvektoren<br />

ligesom i det foregående eksempel. Mens de 1. ordens afledede anvendes i A-matricen,<br />

ligeledes som i det foregående eksempel.<br />

T

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!