VBN - Aalborg Universitet
VBN - Aalborg Universitet
VBN - Aalborg Universitet
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