Fehlersuche mex-File

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Achi

Forum-Guru


	Beiträge: 250

	Anmeldedatum: 14.04.10

	Wohnort: ---

	Version: R2007a

Verfasst am: 02.12.2010, 17:00 Titel: Fehlersuche mex-File

Ich finde schon einen ganzen Tag den Fehler in meinem mex-File nicht

Es werden Strahlen berechnet und wenn diese nicht auf den Empfänger treffen sollen diese gar nicht aufgenommen werden.
1. Sie gehen nicht ganz auf den Empfänger
2. Es kommen auch Straheln heraus, die nicht auf den Empfänger treffen.
So hier der interessante Ausschnitt:

Code:

/* Reflexion 1.Ordnung */

pointerwinkel=&Winkelmatrix[0];
if (n>0 && winkel[n]>0 && winkel[n]!=winkel[n-1]){ // Bedignung nur hinreichend wenn Kanten nicht direkt || aufeinander folgen
//-------------------------------------------------------------------------
/* Koordinaten berechnen des reflektierten Strahls */

for(m=0;;m++){
//-------------------------------------------------------------------------

/* Berechnung der Größe des Ausgabearrays*/
werte=0;
for (s=0;s<=counter;s++)
{werte=werte+mArray[s];}

mexPrintf("werte= %i \n",werte);

/* Koordinatenberechnung der Reflexionen */

*(xArrayReflex+werte+m)=*(xArray+n)+m*Step*-cosd(reflexionswinkel);
*(yArrayReflex+werte+m)=*(yArray+n)+m*Step* sind(reflexionswinkel);
/*Abbruchbedingungen der Berechnungen*/

if (*(xArrayReflex+werte+m)<=0.5 ||
*(yArrayReflex+werte+m)<=0.5 ||
*(xArrayReflex+werte+m)>=maxsensorradius*2 ||
*(yArrayReflex+werte+m)>=maxsensorradius*2 ||
isnan(*(xArrayReflex+werte+m)) ||
isnan(*(yArrayReflex+werte+m)) )
{break;}
//-------------------------------------------------------------------------
/* Dynamische Speicher Allokierung */

if (m==speicherreflex) {

speicherreflex*=2; // Doppelte Länge

xArrayReflex=mxRealloc(xArrayReflex,speicherreflex*sizeof(double)); // mxRealloc: Berechnete Werte werden umkopiert
yArrayReflex=mxRealloc(yArrayReflex,speicherreflex*sizeof(double)); // (Speicher wurde vorher mit mxMalloc erzeugt)
linearerind=mxRealloc(linearerind,speicherreflex*sizeof(int));
mArray=mxRealloc(mArray,speicherreflex*sizeof(int));
}
}
/* Überprüfung ob ein reflektierter Strahl auf den Empfänger getroffen ist */

schalterreflex=0;
speicher=floor(10*(matrixgroesse/step)*strahlen+0.5);
abstandsarrayreflex=mxMalloc(speicher*sizeof(double));

z=0;
for(r=0;r<=m;r++) {
for(u=0;u<=strahlen;u++){
//-------------------------------------------------------------------------
/* Abstandsberechnung reflektierter Strahl */

abstandsarrayreflex[z]=(sqrt(pow(xArrayReflex[werte+r]-empfreflexx[u],2)+(pow(yArrayReflex[werte+r]-empfreflexy[u],2))));
if (abstandsarrayreflex[z]<=empfaengerbreitereflexion)
{
schalterreflex=1;
m=r;
break;
}

}
//-------------------------------------------------------------------------
/* Dynamische Speichererweiterung */
if (z==speicher) {

speicher*=2;
abstandsarrayreflex=mxRealloc(abstandsarrayreflex,speicher*sizeof(double));
}
//-------------------------------------------------------------------------
z++;
}
//-------------------------------------------------------------------------

if (schalterreflex==1) {

mArray[counter]=m; // Externer Schleifendzähler um Anzahl der Schleifen die beendet wurden zu zählen
counter+=1;
}
}

Funktion ohne Link?

Ausgabe erfolgt dann so:

Code:

/*Ausgabe der Koordinaten reflektierte Strahlen und Intensität*/

/* Berechnete Koordinaten der Reflexionen müssen aufsummiertwerden,
* um Speicherbedarf zu ermitteln */

if (counter>0) {
plhs[4]=mxCreateDoubleMatrix(1,(werte+m),mxREAL);
plhs[5]=mxCreateDoubleMatrix(1,(werte+m),mxREAL);
// plhs[6]=mxCreateDoubleMatrix(1,arraylaenge,mxREAL);

memcpy(mxGetPr(plhs[4]),xArrayReflex,(werte+m)*sizeof(double));
memcpy(mxGetPr(plhs[5]),yArrayReflex,(werte+m)*sizeof(double));
// memcpy(mxGetPr(plhs[6]),intensitreflex,arraylaenge*sizeof(double));
}

else {
plhs[4]=mxCreateDoubleMatrix(0,0,mxREAL);
plhs[5]=mxCreateDoubleMatrix(0,0,mxREAL);
// plhs[6]=mxCreateDoubleMatrix(0,0,mxREAL);
}

Funktion ohne Link?

Ich weiss es ist viel, würde mich aber über jede Hilfe freuen!!! Rolling Eyes

Achi

Themenstarter

Forum-Guru


	Beiträge: 250

	Anmeldedatum: 14.04.10

	Wohnort: ---

	Version: R2007a

Verfasst am: 02.12.2010, 17:13 Titel:

An der Stelle:

Code:

z=0;
for(r=0;r<=m;r++) {
for(u=0;u<=strahlen;u++){
//-------------------------------------------------------------------------
/* Abstandsberechnung reflektierter Strahl */

abstandsarrayreflex[z]=(sqrt(pow(xArrayReflex[werte+r]-empfreflexx[u],2)+(pow(yArrayReflex[werte+r]-empfreflexy[u],2))));
if (abstandsarrayreflex[z]<=empfaengerbreitereflexion)
{
schalterreflex=1;
m=r;
break;
}

}

Funktion ohne Link?

wird jeder Punkt der Reflexion (xArrayReflex) mit den Empfängerpositionen verglichen (empfreflexx) ist einer dieser kleine genug, hat er den Empfänger getroffen und soll im Array bleiben. Ansonsten soll er überschrieben werden beim nächsten Durchlauf.
Falls es einen Treffer gab:

Code:

if (schalterreflex==1) {

mArray[counter]=m; // Externer Schleifendzähler um Anzahl der Schleifen die beendet wurden zu zählen
counter+=1;
}
}
wird die Anzahl in mArray[counter] geschrieben.
So sit die Idee, produziert aber noch Fehler...

Funktion ohne Link?

Achi

Themenstarter

Forum-Guru


	Beiträge: 250

	Anmeldedatum: 14.04.10

	Wohnort: ---

	Version: R2007a

Verfasst am: 02.12.2010, 17:49 Titel:

Ok ein Fehler wurde gefunden:
ich musste falls es keinen Treffe gab m auf 0 setzen, weil in der Ausgabe ja:

Code:

if (counter>0) {
plhs[4]=mxCreateDoubleMatrix(1,(werte+m),mxREAL);
plhs[5]=mxCreateDoubleMatrix(1,(werte+m),mxREAL);
// plhs[6]=mxCreateDoubleMatrix(1,arraylaenge,mxREAL);

memcpy(mxGetPr(plhs[4]),xArrayReflex,(werte+m)*sizeof(double));
memcpy(mxGetPr(plhs[5]),yArrayReflex,(werte+m)*sizeof(double));
// memcpy(mxGetPr(plhs[6]),intensitreflex,arraylaenge*sizeof(double))

Funktion ohne Link?

das m steht.
Also:

Code:

if (schalterreflex==0){
mArray[counter]=0;
counter=counter;
m=0;

Funktion ohne Link?

Es gibt aber immer noch Strahlen, die nicht auf den Empfänger treffen und dennoch angezeigt werden...???

Achi

Themenstarter

Forum-Guru


	Beiträge: 250

	Anmeldedatum: 14.04.10

	Wohnort: ---

	Version: R2007a

Verfasst am: 02.12.2010, 18:24 Titel:

Hier ein Beispiel wie es aussieht.
Der Würfel ist das Objekt, der Strahlen an Grenzflächen zurückwirft.
Rot Sender und Empfänger. Blau die Reflexionen ohne weitere Brechung (bis jetzt). Strahlen die der Sender aussendet werden nicht angezeigt.
Es sollen nur Strahlen angezeigt werden, die den Empfänger berühren.

Beispiel.jpg

Beschreibung:

Download

Dateiname:

Beispiel.jpg

Dateigröße:

19.54 KB

Heruntergeladen:

448 mal


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.