/* * MATLAB Compiler: 3.0 * Date: Wed Jul 06 02:27:15 2005 * Arguments: "-B" "macro_default" "-O" "all" "-O" "fold_scalar_mxarrays:on" * "-O" "fold_non_scalar_mxarrays:on" "-O" "optimize_integer_for_loops:on" "-O" * "array_indexing:on" "-O" "optimize_conditionals:on" "-m" "-W" "main" "-L" * "C" "-t" "-T" "link:exe" "-h" "libmmfile.mlib" "-B" "sgl" "-m" "-W" "main" * "-L" "C" "-t" "-T" "link:exe" "-h" "libmmfile.mlib" "-W" "mainhg" * "libmwsglm.mlib" "micampo.m" "matrices1lazo.m" "leeActivo.m" "leePasivo.m" * "matActivo.m" */ #include "micampo.h" #include "leeactivo.h" #include "leepasivo.h" #include "libmatlbm.h" #include "libmmfile.h" #include "libmwsglm.h" #include "matactivo.h" #include "matrices1lazo.h" #include "title.h" #include "xlabel.h" #include "ylabel.h" static mxChar _array1_[3] = { 'a', 'l', 'l' }; static mxArray * _mxarray0_; static mxChar _array3_[8] = { 'p', 'a', 't', 'h', '.', 't', 'x', 't' }; static mxArray * _mxarray2_; static mxChar _array5_[10] = { 0x005c, 'f', 'a', 's', 'e', 's', '.', 't', 'x', 't' }; static mxArray * _mxarray4_; static mxArray * _mxarray6_; static mxChar _array8_[1] = { ';' }; static mxArray * _mxarray7_; static mxArray * _mxarray9_; static mxChar _array11_[1] = { ',' }; static mxArray * _mxarray10_; static mxChar _array13_[1] = { '.' }; static mxArray * _mxarray12_; static mxArray * _mxarray14_; static mxArray * _mxarray15_; static mxArray * _mxarray16_; static mxChar _array18_[10] = { 0x005c, 'd', 'a', 't', 'o', 's', '.', 't', 'x', 't' }; static mxArray * _mxarray17_; static mxChar _array20_[11] = { 0x005c, 'p', 'a', 's', 'i', 'v', 'o', '.', 't', 'x', 't' }; static mxArray * _mxarray19_; static mxArray * _mxarray21_; static mxChar _array23_[11] = { 0x005c, 'a', 'c', 't', 'i', 'v', 'o', '.', 't', 'x', 't' }; static mxArray * _mxarray22_; static mxArray * _mxarray24_; static mxArray * _mxarray25_; static mxArray * _mxarray26_; static mxArray * _mxarray27_; static mxArray * _mxarray28_; static mxArray * _mxarray29_; static mxChar _array31_[2] = { 'o', 'n' }; static mxArray * _mxarray30_; static mxChar _array33_[1] = { 'r' }; static mxArray * _mxarray32_; static mxChar _array35_[1] = { 'g' }; static mxArray * _mxarray34_; static mxChar _array37_[19] = { 'C', 'a', 'm', 'p', 'o', ' ', 'r', 'e', 's', 'u', 'l', 't', 'a', 'n', 't', 'e', ' ', 'B', 'r' }; static mxArray * _mxarray36_; static mxChar _array39_[5] = { 'x', ' ', '(', 'm', ')' }; static mxArray * _mxarray38_; static mxChar _array41_[7] = { 'B', 'r', ' ', '(', 'u', 'T', ')' }; static mxArray * _mxarray40_; static mxChar _array43_[8] = { 'B', 'r', ' ', 'l', 'i', 'n', 'e', 'a' }; static mxArray * _mxarray42_; static mxChar _array45_[7] = { 'B', 'r', ' ', 'l', 'a', 'z', 'o' }; static mxArray * _mxarray44_; static mxChar _array47_[7] = { 'B', 'r', ' ', 'n', 'e', 't', 'o' }; static mxArray * _mxarray46_; static mxChar _array49_[3] = { 'o', 'f', 'f' }; static mxArray * _mxarray48_; static mxChar _array51_[21] = { 'F', 'a', 'c', 't', 'o', 'r', ' ', 'd', 'e', ' ', 'A', 'p', 'a', 'n', 't', 'a', 'l', 'l', 'a', 'd', 'o' }; static mxArray * _mxarray50_; static mxChar _array53_[7] = { 'B', 'r', 'l', 'i', 'n', 'e', 'a' }; static mxArray * _mxarray52_; static mxChar _array55_[10] = { 0x005c, 'c', 'a', 'm', 'p', 'o', '.', 't', 'x', 't' }; static mxArray * _mxarray54_; static mxChar _array57_[1] = { 'w' }; static mxArray * _mxarray56_; static mxChar _array59_[29] = { '%', '6', '.', '2', 'f', ',', ' ', '%', '6', '.', '2', 'f', ',', ' ', '%', '6', '.', '2', 'f', ',', ' ', '%', '6', '.', '2', 'f', ' ', 0x005c, 'r' }; static mxArray * _mxarray58_; static mxChar _array61_[7] = { 0x005c, 'f', 'a', '.', 't', 'x', 't' }; static mxArray * _mxarray60_; static mxChar _array63_[15] = { '%', '6', '.', '2', 'f', ',', ' ', '%', '6', '.', '4', 'f', ' ', 0x005c, 'r' }; static mxArray * _mxarray62_; void InitializeModule_micampo(void) { _mxarray0_ = mclInitializeString(3, _array1_); _mxarray2_ = mclInitializeString(8, _array3_); _mxarray4_ = mclInitializeString(10, _array5_); _mxarray6_ = mclInitializeDouble(1.0); _mxarray7_ = mclInitializeString(1, _array8_); _mxarray9_ = mclInitializeDoubleVector(0, 0, (double *)NULL); _mxarray10_ = mclInitializeString(1, _array11_); _mxarray12_ = mclInitializeString(1, _array13_); _mxarray14_ = mclInitializeDouble(2.0); _mxarray15_ = mclInitializeDouble(3.0); _mxarray16_ = mclInitializeDouble(4.0); _mxarray17_ = mclInitializeString(10, _array18_); _mxarray19_ = mclInitializeString(11, _array20_); _mxarray21_ = mclInitializeDouble(-1.0); _mxarray22_ = mclInitializeString(11, _array23_); _mxarray24_ = mclInitializeDouble(0.0); _mxarray25_ = mclInitializeDouble(3.141592653589793); _mxarray26_ = mclInitializeDouble(180.0); _mxarray27_ = mclInitializeDouble(10.0); _mxarray28_ = mclInitializeDouble(5.0); _mxarray29_ = mclInitializeDouble(6.0); _mxarray30_ = mclInitializeString(2, _array31_); _mxarray32_ = mclInitializeString(1, _array33_); _mxarray34_ = mclInitializeString(1, _array35_); _mxarray36_ = mclInitializeString(19, _array37_); _mxarray38_ = mclInitializeString(5, _array39_); _mxarray40_ = mclInitializeString(7, _array41_); _mxarray42_ = mclInitializeString(8, _array43_); _mxarray44_ = mclInitializeString(7, _array45_); _mxarray46_ = mclInitializeString(7, _array47_); _mxarray48_ = mclInitializeString(3, _array49_); _mxarray50_ = mclInitializeString(21, _array51_); _mxarray52_ = mclInitializeString(7, _array53_); _mxarray54_ = mclInitializeString(10, _array55_); _mxarray56_ = mclInitializeString(1, _array57_); _mxarray58_ = mclInitializeString(29, _array59_); _mxarray60_ = mclInitializeString(7, _array61_); _mxarray62_ = mclInitializeString(15, _array63_); } void TerminateModule_micampo(void) { mxDestroyArray(_mxarray62_); mxDestroyArray(_mxarray60_); mxDestroyArray(_mxarray58_); mxDestroyArray(_mxarray56_); mxDestroyArray(_mxarray54_); mxDestroyArray(_mxarray52_); mxDestroyArray(_mxarray50_); mxDestroyArray(_mxarray48_); mxDestroyArray(_mxarray46_); mxDestroyArray(_mxarray44_); mxDestroyArray(_mxarray42_); mxDestroyArray(_mxarray40_); mxDestroyArray(_mxarray38_); mxDestroyArray(_mxarray36_); mxDestroyArray(_mxarray34_); mxDestroyArray(_mxarray32_); mxDestroyArray(_mxarray30_); mxDestroyArray(_mxarray29_); mxDestroyArray(_mxarray28_); mxDestroyArray(_mxarray27_); mxDestroyArray(_mxarray26_); mxDestroyArray(_mxarray25_); mxDestroyArray(_mxarray24_); mxDestroyArray(_mxarray22_); mxDestroyArray(_mxarray21_); mxDestroyArray(_mxarray19_); mxDestroyArray(_mxarray17_); mxDestroyArray(_mxarray16_); mxDestroyArray(_mxarray15_); mxDestroyArray(_mxarray14_); mxDestroyArray(_mxarray12_); mxDestroyArray(_mxarray10_); mxDestroyArray(_mxarray9_); mxDestroyArray(_mxarray7_); mxDestroyArray(_mxarray6_); mxDestroyArray(_mxarray4_); mxDestroyArray(_mxarray2_); mxDestroyArray(_mxarray0_); } static void Mmicampo(void); _mexLocalFunctionTable _local_function_table_micampo = { 0, (mexFunctionTableEntry *)NULL }; /* * The function "mlfMicampo" contains the normal interface for the "micampo" * M-function from file "c:\matlab6p5\work\micampo.m" (lines 1-290). This * function processes any input arguments and passes them to the implementation * version of the function, appearing above. */ void mlfMicampo(void) { mlfEnterNewContext(0, 0); Mmicampo(); mlfRestorePreviousContext(0, 0); } /* * The function "mlxMicampo" contains the feval interface for the "micampo" * M-function from file "c:\matlab6p5\work\micampo.m" (lines 1-290). The feval * function calls the implementation version of micampo through this function. * This function processes any input arguments and passes them to the * implementation version of the function, appearing above. */ void mlxMicampo(int nlhs, mxArray * plhs[], int nrhs, mxArray * prhs[]) { if (nlhs > 0) { mlfError( mxCreateString( "Run-time Error: File: micampo Line: 4 Column: " "1 The function \"micampo\" was called with mor" "e than the declared number of outputs (0)."), NULL); } if (nrhs > 0) { mlfError( mxCreateString( "Run-time Error: File: micampo Line: 4 Column:" " 1 The function \"micampo\" was called with m" "ore than the declared number of inputs (0)."), NULL); } mlfEnterNewContext(0, 0); Mmicampo(); mlfRestorePreviousContext(0, 0); } /* * The function "Mmicampo" is the implementation version of the "micampo" * M-function from file "c:\matlab6p5\work\micampo.m" (lines 1-290). It * contains the actual compiled code for that M-function. It is a static * function and must only be called from one of the interface functions, * appearing below. */ /* * % Modulo de calculo de campo * % * * function micampo */ static void Mmicampo(void) { mexLocalFunctionTable save_local_function_table_ = mclSetCurrentLocalFunctionTable(&_local_function_table_micampo); mxArray * u = NULL; mxArray * n = NULL; mxArray * FA = NULL; mxArray * ju = NULL; mxArray * b = NULL; mxArray * a = NULL; mxArray * Brlazo = NULL; mxArray * Bylazo = NULL; mxArray * Bxlazo = NULL; mxArray * campototal = NULL; mxArray * Btotal = NULL; mxArray * Bytotal = NULL; mxArray * Bxtotal = NULL; mxArray * Byitotal = NULL; mxArray * Byrtotal = NULL; mxArray * Bxitotal = NULL; mxArray * Bxrtotal = NULL; mxArray * campolazo = NULL; mxArray * Byilazo = NULL; mxArray * Byrlazo = NULL; mxArray * Bxilazo = NULL; mxArray * Bxrlazo = NULL; mxArray * campo = NULL; mxArray * Br = NULL; mxArray * By = NULL; mxArray * Bx = NULL; mxArray * m = NULL; mxArray * Byi = NULL; mxArray * Byr = NULL; mxArray * Bxi = NULL; mxArray * Bxr = NULL; mxArray * lazosAct = NULL; mxArray * fidAct = NULL; mxArray * anguloLazo = NULL; mxArray * Ilazo = NULL; mxArray * mislazos = NULL; mxArray * I = NULL; mxArray * lazosPas = NULL; mxArray * hayLazosActivos = NULL; mxArray * hayLazosPasivos = NULL; mxArray * fidPas = NULL; mxArray * incremento = NULL; mxArray * xmax = NULL; mxArray * xmin = NULL; mxArray * altura = NULL; mxArray * f = NULL; mxArray * fases = NULL; mxArray * desf = NULL; mxArray * intens = NULL; mxArray * y = NULL; mxArray * x = NULL; mxArray * ii = NULL; mxArray * indices = NULL; mxArray * cont = NULL; mxArray * fdatos = NULL; mxArray * cadena = NULL; mxArray * path = NULL; mxArray * fid = NULL; mxArray * ans = NULL; /* * * % Limpiamos la memoria. * clear */ mlfClear( &Br, &Brlazo, &Btotal, &Bx, &Bxi, &Bxilazo, &Bxitotal, &Bxlazo, &Bxr, &Bxrlazo, &Bxrtotal, &Bxtotal, &By, &Byi, &Byilazo, &Byitotal, &Bylazo, &Byr, &Byrlazo, &Byrtotal, &Bytotal, &FA, &I, &Ilazo, &a, &altura, &anguloLazo, &ans, &b, &cadena, &campo, &campolazo, &campototal, &cont, &desf, &f, &fases, &fdatos, &fid, &fidAct, &fidPas, &hayLazosActivos, &hayLazosPasivos, &ii, &incremento, &indices, &intens, &ju, &lazosAct, &lazosPas, &m, &mislazos, &n, &path, &u, &x, &xmax, &xmin, &y, NULL); /* * close all */ mclPrintAns(&ans, mlfNClose(0, _mxarray0_, NULL)); /* * * % Leemos datos de los ficheros. * * % Fichero path * fid=fopen('path.txt'); */ mlfAssign(&fid, mlfFopen(NULL, NULL, _mxarray2_, NULL, NULL)); /* * path=fgets(fid); */ mlfAssign(&path, mlfFgets(NULL, mclVv(fid, "fid"), NULL)); /* * fclose(fid); */ mclAssignAns(&ans, mlfFclose(mclVv(fid, "fid"))); /* * * % Fichero datos * cadena=strcat(path,'\fases.txt'); */ mlfAssign(&cadena, mlfStrcat(mclVv(path, "path"), _mxarray4_, NULL)); /* * fdatos=fopen(cadena); */ mlfAssign( &fdatos, mlfFopen(NULL, NULL, mclVv(cadena, "cadena"), NULL, NULL)); /* * * cont=1; */ mlfAssign(&cont, _mxarray6_); /* * while(~feof(fdatos)) */ while (mclNotBool(mlfFeof(mclVv(fdatos, "fdatos")))) { /* * * cadena=fgets(fdatos); % Leemos la fase actual */ mlfAssign(&cadena, mlfFgets(NULL, mclVv(fdatos, "fdatos"), NULL)); /* * indices=strfind(cadena,';'); % Aqui tenemos todas las ocurrencias de ";" */ mlfAssign(&indices, mlfStrfind(mclVv(cadena, "cadena"), _mxarray7_)); /* * * for ii=1:length(cadena) % Si el simbolo decimal es ',' lo cambiamos por '.' */ { int v_ = mclForIntStart(1); int e_ = mclLengthInt(mclVv(cadena, "cadena")); if (v_ > e_) { mlfAssign(&ii, _mxarray9_); } else { /* * if cadena(ii)==',' * cadena(ii)='.'; * end * end */ for (; ; ) { if (mclEqBool( mclIntArrayRef1(mclVv(cadena, "cadena"), v_), _mxarray10_)) { mclIntArrayAssign1(&cadena, _mxarray12_, v_); } if (v_ == e_) { break; } ++v_; } mlfAssign(&ii, mlfScalar(v_)); } } /* * * x=cadena(1:indices(1)-1); */ mlfAssign( &x, mclArrayRef1( mclVv(cadena, "cadena"), mlfColon( _mxarray6_, mclMinus( mclIntArrayRef1(mclVv(indices, "indices"), 1), _mxarray6_), NULL))); /* * x=str2double(x); */ mlfAssign(&x, mlfStr2double(mclVv(x, "x"))); /* * * y=cadena(indices(1)+1:indices(2)-1); */ mlfAssign( &y, mclArrayRef1( mclVv(cadena, "cadena"), mlfColon( mclPlus( mclIntArrayRef1(mclVv(indices, "indices"), 1), _mxarray6_), mclMinus( mclIntArrayRef1(mclVv(indices, "indices"), 2), _mxarray6_), NULL))); /* * y=str2double(y); */ mlfAssign(&y, mlfStr2double(mclVv(y, "y"))); /* * * intens=cadena(indices(2)+1:indices(3)-1); */ mlfAssign( &intens, mclArrayRef1( mclVv(cadena, "cadena"), mlfColon( mclPlus( mclIntArrayRef1(mclVv(indices, "indices"), 2), _mxarray6_), mclMinus( mclIntArrayRef1(mclVv(indices, "indices"), 3), _mxarray6_), NULL))); /* * intens=str2double(intens); */ mlfAssign(&intens, mlfStr2double(mclVv(intens, "intens"))); /* * * desf=cadena(indices(3)+1:indices(4)-1); */ mlfAssign( &desf, mclArrayRef1( mclVv(cadena, "cadena"), mlfColon( mclPlus( mclIntArrayRef1(mclVv(indices, "indices"), 3), _mxarray6_), mclMinus( mclIntArrayRef1(mclVv(indices, "indices"), 4), _mxarray6_), NULL))); /* * desf=str2double(desf); */ mlfAssign(&desf, mlfStr2double(mclVv(desf, "desf"))); /* * * fases(cont,1)=x; */ mclArrayAssign2(&fases, mclVv(x, "x"), mclVv(cont, "cont"), _mxarray6_); /* * fases(cont,2)=y; */ mclArrayAssign2( &fases, mclVv(y, "y"), mclVv(cont, "cont"), _mxarray14_); /* * fases(cont,3)=intens; */ mclArrayAssign2( &fases, mclVv(intens, "intens"), mclVv(cont, "cont"), _mxarray15_); /* * fases(cont,4)=desf; */ mclArrayAssign2( &fases, mclVv(desf, "desf"), mclVv(cont, "cont"), _mxarray16_); /* * cont=cont+1; */ mlfAssign(&cont, mclPlus(mclVv(cont, "cont"), _mxarray6_)); /* * * end */ } /* * * fclose(fdatos); */ mclAssignAns(&ans, mlfFclose(mclVv(fdatos, "fdatos"))); /* * * % Leemos datos generales * cadena=strcat(path,'\datos.txt'); */ mlfAssign(&cadena, mlfStrcat(mclVv(path, "path"), _mxarray17_, NULL)); /* * fid=fopen(cadena); */ mlfAssign(&fid, mlfFopen(NULL, NULL, mclVv(cadena, "cadena"), NULL, NULL)); /* * * cadena=fgets(fid); */ mlfAssign(&cadena, mlfFgets(NULL, mclVv(fid, "fid"), NULL)); /* * indices=strfind(cadena,';'); */ mlfAssign(&indices, mlfStrfind(mclVv(cadena, "cadena"), _mxarray7_)); /* * * f=str2double(cadena(1:indices(1)-1)); */ mlfAssign( &f, mlfStr2double( mclArrayRef1( mclVv(cadena, "cadena"), mlfColon( _mxarray6_, mclMinus(mclIntArrayRef1(mclVv(indices, "indices"), 1), _mxarray6_), NULL)))); /* * altura=str2double(cadena(indices(1)+1:indices(2)-1)); */ mlfAssign( &altura, mlfStr2double( mclArrayRef1( mclVv(cadena, "cadena"), mlfColon( mclPlus(mclIntArrayRef1(mclVv(indices, "indices"), 1), _mxarray6_), mclMinus(mclIntArrayRef1(mclVv(indices, "indices"), 2), _mxarray6_), NULL)))); /* * xmin=str2double(cadena(indices(2)+1:indices(3)-1)); */ mlfAssign( &xmin, mlfStr2double( mclArrayRef1( mclVv(cadena, "cadena"), mlfColon( mclPlus(mclIntArrayRef1(mclVv(indices, "indices"), 2), _mxarray6_), mclMinus(mclIntArrayRef1(mclVv(indices, "indices"), 3), _mxarray6_), NULL)))); /* * xmax=str2double(cadena(indices(3)+1:indices(4)-1)); */ mlfAssign( &xmax, mlfStr2double( mclArrayRef1( mclVv(cadena, "cadena"), mlfColon( mclPlus(mclIntArrayRef1(mclVv(indices, "indices"), 3), _mxarray6_), mclMinus(mclIntArrayRef1(mclVv(indices, "indices"), 4), _mxarray6_), NULL)))); /* * incremento = str2double(cadena(indices(4)+1:indices(5)-1)); */ mlfAssign( &incremento, mlfStr2double( mclArrayRef1( mclVv(cadena, "cadena"), mlfColon( mclPlus(mclIntArrayRef1(mclVv(indices, "indices"), 4), _mxarray6_), mclMinus(mclIntArrayRef1(mclVv(indices, "indices"), 5), _mxarray6_), NULL)))); /* * * fclose(fid); */ mclAssignAns(&ans, mlfFclose(mclVv(fid, "fid"))); /* * * % Solo habra un tipo de lazos, asi que si no podemos abrir 'pasivo.txt' es * % porque nos han mandado un lazo activo o porque no hay lazo. * cadena=strcat(path,'\pasivo.txt'); */ mlfAssign(&cadena, mlfStrcat(mclVv(path, "path"), _mxarray19_, NULL)); /* * fidPas=fopen(cadena); % Abrimos el fichero que contiene el/los lazos. */ mlfAssign( &fidPas, mlfFopen(NULL, NULL, mclVv(cadena, "cadena"), NULL, NULL)); /* * * if fidPas~=-1 */ if (mclNeBool(mclVv(fidPas, "fidPas"), _mxarray21_)) { /* * * hayLazosPasivos=true; */ mlfAssign(&hayLazosPasivos, mlfTrue(NULL)); /* * hayLazosActivos=false; */ mlfAssign(&hayLazosActivos, mlfFalse(NULL)); /* * lazosPas=leePasivo(fidPas); % Obtenemos los lazos del fichero. */ mlfAssign(&lazosPas, mlfLeepasivo(mclVv(fidPas, "fidPas"), NULL)); /* * fclose(fidPas); */ mclAssignAns(&ans, mlfFclose(mclVv(fidPas, "fidPas"))); /* * [I,mislazos]=matrices1lazo(f,lazosPas,fases); % Obtenemos los lazos simples equivalentes. */ mlfAssign( &I, mlfMatrices1lazo( &mislazos, mclVv(f, "f"), mclVv(lazosPas, "lazosPas"), mclVv(fases, "fases"))); /* * Ilazo=abs(I); % Modulo de las intensidades de lazos simples. */ mlfAssign(&Ilazo, mlfAbs(mclVv(I, "I"))); /* * anguloLazo=angle(I); % Desfases de las intensidades de lazos simples. */ mlfAssign(&anguloLazo, mlfAngle(mclVv(I, "I"))); /* * * else */ } else { /* * * hayLazosPasivos=false; */ mlfAssign(&hayLazosPasivos, mlfFalse(NULL)); /* * hayLazosActivos=false; */ mlfAssign(&hayLazosActivos, mlfFalse(NULL)); /* * cadena=strcat(path,'\activo.txt'); */ mlfAssign(&cadena, mlfStrcat(mclVv(path, "path"), _mxarray22_, NULL)); /* * fidAct=fopen(cadena); % Intentamos abrir el activo. */ mlfAssign( &fidAct, mlfFopen(NULL, NULL, mclVv(cadena, "cadena"), NULL, NULL)); /* * if fidAct~=-1 */ if (mclNeBool(mclVv(fidAct, "fidAct"), _mxarray21_)) { /* * hayLazosActivos=true; */ mlfAssign(&hayLazosActivos, mlfTrue(NULL)); /* * lazosAct=leeActivo(fidAct); */ mlfAssign(&lazosAct, mlfLeeactivo(mclVv(fidAct, "fidAct"))); /* * fclose(fidAct); */ mclAssignAns(&ans, mlfFclose(mclVv(fidAct, "fidAct"))); /* * end */ } /* * * end */ } /* * * if hayLazosPasivos==true */ if (mclEqBool(mclVv(hayLazosPasivos, "hayLazosPasivos"), mlfTrue(NULL))) { /* * * Ilazo=abs(I); */ mlfAssign(&Ilazo, mlfAbs(mclVv(I, "I"))); /* * anguloLazo=angle(I); */ mlfAssign(&anguloLazo, mlfAngle(mclVv(I, "I"))); /* * * elseif hayLazosActivos==true */ } else if (mclEqBool( mclVv(hayLazosActivos, "hayLazosActivos"), mlfTrue(NULL))) { /* * * [I,mislazos]=matActivo(lazosAct,0); */ mlfAssign( &I, mlfMatactivo(&mislazos, mclVv(lazosAct, "lazosAct"), _mxarray24_)); /* * Ilazo=abs(I); */ mlfAssign(&Ilazo, mlfAbs(mclVv(I, "I"))); /* * anguloLazo=angle(I); */ mlfAssign(&anguloLazo, mlfAngle(mclVv(I, "I"))); /* * * end */ } /* * * % Contador para llevar el numero de puntos * cont=1; */ mlfAssign(&cont, _mxarray6_); /* * * % Inicializacion x * x=xmin; */ mlfAssign(&x, mclVv(xmin, "xmin")); /* * * while(x<=xmax) */ while (mclLeBool(mclVv(x, "x"), mclVv(xmax, "xmax"))) { /* * * % Ponemos las variables a cero para evitar interferencias de valores de * % una iteracion anterior. * Bxr=0; */ mlfAssign(&Bxr, _mxarray24_); /* * Bxi=0; */ mlfAssign(&Bxi, _mxarray24_); /* * Byr=0; */ mlfAssign(&Byr, _mxarray24_); /* * Byi=0; */ mlfAssign(&Byi, _mxarray24_); /* * * % Campo creado por la linea * * for m=1:size(fases,1) */ { int v_ = mclForIntStart(1); int e_ = mclForIntEnd( mlfSize( mclValueVarargout(), mclVv(fases, "fases"), _mxarray6_)); if (v_ > e_) { mlfAssign(&m, _mxarray9_); } else { /* * * % Contribucion fase i * Bxr=Bxr+2*(altura-fases(m,2))*fases(m,3)*cos(fases(m,4)*pi/180)/((altura-fases(m,2))^2+(x-fases(m,1))^2); * Bxi=Bxi+2*(altura-fases(m,2))*fases(m,3)*sin(fases(m,4)*pi/180)/((altura-fases(m,2))^2+(x-fases(m,1))^2); * Byr=Byr+2*(x-fases(m,1))*fases(m,3)*cos(fases(m,4)*pi/180)/((altura-fases(m,2))^2+(x-fases(m,1))^2); * Byi=Byi+2*(x-fases(m,1))*fases(m,3)*sin(fases(m,4)*pi/180)/((altura-fases(m,2))^2+(x-fases(m,1))^2); * * end */ for (; ; ) { mlfAssign( &Bxr, mclPlus( mclVv(Bxr, "Bxr"), mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(fases, "fases"), v_, 2))), mclIntArrayRef2(mclVv(fases, "fases"), v_, 3)), mlfCos( mclMrdivide( mclMtimes( mclIntArrayRef2(mclVv(fases, "fases"), v_, 4), _mxarray25_), _mxarray26_))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2(mclVv(fases, "fases"), v_, 2)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2(mclVv(fases, "fases"), v_, 1)), _mxarray14_))))); mlfAssign( &Bxi, mclPlus( mclVv(Bxi, "Bxi"), mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(fases, "fases"), v_, 2))), mclIntArrayRef2(mclVv(fases, "fases"), v_, 3)), mlfSin( mclMrdivide( mclMtimes( mclIntArrayRef2(mclVv(fases, "fases"), v_, 4), _mxarray25_), _mxarray26_))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2(mclVv(fases, "fases"), v_, 2)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2(mclVv(fases, "fases"), v_, 1)), _mxarray14_))))); mlfAssign( &Byr, mclPlus( mclVv(Byr, "Byr"), mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(fases, "fases"), v_, 1))), mclIntArrayRef2(mclVv(fases, "fases"), v_, 3)), mlfCos( mclMrdivide( mclMtimes( mclIntArrayRef2(mclVv(fases, "fases"), v_, 4), _mxarray25_), _mxarray26_))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2(mclVv(fases, "fases"), v_, 2)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2(mclVv(fases, "fases"), v_, 1)), _mxarray14_))))); mlfAssign( &Byi, mclPlus( mclVv(Byi, "Byi"), mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(fases, "fases"), v_, 1))), mclIntArrayRef2(mclVv(fases, "fases"), v_, 3)), mlfSin( mclMrdivide( mclMtimes( mclIntArrayRef2(mclVv(fases, "fases"), v_, 4), _mxarray25_), _mxarray26_))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2(mclVv(fases, "fases"), v_, 2)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2(mclVv(fases, "fases"), v_, 1)), _mxarray14_))))); if (v_ == e_) { break; } ++v_; } mlfAssign(&m, mlfScalar(v_)); } } /* * * % Pasamos a microteslas (dividimos entre 10) * Bxr=Bxr/10; */ mlfAssign(&Bxr, mclMrdivide(mclVv(Bxr, "Bxr"), _mxarray27_)); /* * Bxi=Bxi/10; */ mlfAssign(&Bxi, mclMrdivide(mclVv(Bxi, "Bxi"), _mxarray27_)); /* * Byr=Byr/10; */ mlfAssign(&Byr, mclMrdivide(mclVv(Byr, "Byr"), _mxarray27_)); /* * Byi=Byi/10; */ mlfAssign(&Byi, mclMrdivide(mclVv(Byi, "Byi"), _mxarray27_)); /* * * % Campo equivalente * Bx=sqrt(Bxr^2+Bxi^2); */ mlfAssign( &Bx, mlfSqrt( mclPlus( mclMpower(mclVv(Bxr, "Bxr"), _mxarray14_), mclMpower(mclVv(Bxi, "Bxi"), _mxarray14_)))); /* * By=sqrt(Byr^2+Byi^2); */ mlfAssign( &By, mlfSqrt( mclPlus( mclMpower(mclVv(Byr, "Byr"), _mxarray14_), mclMpower(mclVv(Byi, "Byi"), _mxarray14_)))); /* * Br=sqrt(Bx^2+By^2); */ mlfAssign( &Br, mlfSqrt( mclPlus( mclMpower(mclVv(Bx, "Bx"), _mxarray14_), mclMpower(mclVv(By, "By"), _mxarray14_)))); /* * * % Guardamos valores en el vector * campo(cont,1)=x; */ mclArrayAssign2(&campo, mclVv(x, "x"), mclVv(cont, "cont"), _mxarray6_); /* * campo(cont,2)=Bxr; */ mclArrayAssign2( &campo, mclVv(Bxr, "Bxr"), mclVv(cont, "cont"), _mxarray14_); /* * campo(cont,3)=Bxi; */ mclArrayAssign2( &campo, mclVv(Bxi, "Bxi"), mclVv(cont, "cont"), _mxarray15_); /* * campo(cont,4)=Byr; */ mclArrayAssign2( &campo, mclVv(Byr, "Byr"), mclVv(cont, "cont"), _mxarray16_); /* * campo(cont,5)=Byi; */ mclArrayAssign2( &campo, mclVv(Byi, "Byi"), mclVv(cont, "cont"), _mxarray28_); /* * campo(cont,6)=Br; */ mclArrayAssign2( &campo, mclVv(Br, "Br"), mclVv(cont, "cont"), _mxarray29_); /* * * % Ahora ańadimos el valor del campo generado por los lazos. * * if (hayLazosPasivos==true | hayLazosActivos==true) */ { mxArray * a_ = mclInitialize( mclEq( mclVv(hayLazosPasivos, "hayLazosPasivos"), mlfTrue(NULL))); if (mlfTobool(a_) || mlfTobool( mclOr( a_, mclEq( mclVv(hayLazosActivos, "hayLazosActivos"), mlfTrue(NULL))))) { mxDestroyArray(a_); /* * * for m=1:size(mislazos,1) */ { int v_ = mclForIntStart(1); int e_ = mclForIntEnd( mlfSize( mclValueVarargout(), mclVv(mislazos, "mislazos"), _mxarray6_)); if (v_ > e_) { mlfAssign(&m, _mxarray9_); } else { /* * * % Conductor 1 del lazo i. * Bxrlazo=2*(altura-mislazos(m,3))*Ilazo(m)*cos(anguloLazo(m))/((altura-mislazos(m,3))^2+(x-mislazos(m,2))^2); * Bxilazo=2*(altura-mislazos(m,3))*Ilazo(m)*sin(anguloLazo(m))/((altura-mislazos(m,3))^2+(x-mislazos(m,2))^2); * Byrlazo=2*(x-mislazos(m,2))*Ilazo(m)*cos(anguloLazo(m))/((altura-mislazos(m,3))^2+(x-mislazos(m,2))^2); * Byilazo=2*(x-mislazos(m,2))*Ilazo(m)*sin(anguloLazo(m))/((altura-mislazos(m,3))^2+(x-mislazos(m,2))^2); * * % Conductor 2 del lazo i. * Bxrlazo=Bxrlazo+2*(altura-mislazos(m,5))*Ilazo(m)*cos(pi+anguloLazo(m))/((altura-mislazos(m,5))^2+(x-mislazos(m,4))^2); * Bxilazo=Bxilazo+2*(altura-mislazos(m,5))*Ilazo(m)*sin(pi+anguloLazo(m))/((altura-mislazos(m,5))^2+(x-mislazos(m,4))^2); * Byrlazo=Byrlazo+2*(x-mislazos(m,4))*Ilazo(m)*cos(pi+anguloLazo(m))/((altura-mislazos(m,5))^2+(x-mislazos(m,4))^2); * Byilazo=Byilazo+2*(x-mislazos(m,4))*Ilazo(m)*sin(pi+anguloLazo(m))/((altura-mislazos(m,5))^2+(x-mislazos(m,4))^2); * * end */ for (; ; ) { mlfAssign( &Bxrlazo, mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 3))), mclIntArrayRef1(mclVv(Ilazo, "Ilazo"), v_)), mlfCos( mclIntArrayRef1( mclVv(anguloLazo, "anguloLazo"), v_))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 3)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 2)), _mxarray14_)))); mlfAssign( &Bxilazo, mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 3))), mclIntArrayRef1(mclVv(Ilazo, "Ilazo"), v_)), mlfSin( mclIntArrayRef1( mclVv(anguloLazo, "anguloLazo"), v_))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 3)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 2)), _mxarray14_)))); mlfAssign( &Byrlazo, mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 2))), mclIntArrayRef1(mclVv(Ilazo, "Ilazo"), v_)), mlfCos( mclIntArrayRef1( mclVv(anguloLazo, "anguloLazo"), v_))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 3)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 2)), _mxarray14_)))); mlfAssign( &Byilazo, mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 2))), mclIntArrayRef1(mclVv(Ilazo, "Ilazo"), v_)), mlfSin( mclIntArrayRef1( mclVv(anguloLazo, "anguloLazo"), v_))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 3)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 2)), _mxarray14_)))); mlfAssign( &Bxrlazo, mclPlus( mclVv(Bxrlazo, "Bxrlazo"), mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 5))), mclIntArrayRef1( mclVv(Ilazo, "Ilazo"), v_)), mlfCos( mclPlus( _mxarray25_, mclIntArrayRef1( mclVv(anguloLazo, "anguloLazo"), v_)))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 5)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 4)), _mxarray14_))))); mlfAssign( &Bxilazo, mclPlus( mclVv(Bxilazo, "Bxilazo"), mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 5))), mclIntArrayRef1( mclVv(Ilazo, "Ilazo"), v_)), mlfSin( mclPlus( _mxarray25_, mclIntArrayRef1( mclVv(anguloLazo, "anguloLazo"), v_)))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 5)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 4)), _mxarray14_))))); mlfAssign( &Byrlazo, mclPlus( mclVv(Byrlazo, "Byrlazo"), mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 4))), mclIntArrayRef1( mclVv(Ilazo, "Ilazo"), v_)), mlfCos( mclPlus( _mxarray25_, mclIntArrayRef1( mclVv(anguloLazo, "anguloLazo"), v_)))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 5)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 4)), _mxarray14_))))); mlfAssign( &Byilazo, mclPlus( mclVv(Byilazo, "Byilazo"), mclMrdivide( mclMtimes( mclMtimes( mclMtimes( _mxarray14_, mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 4))), mclIntArrayRef1( mclVv(Ilazo, "Ilazo"), v_)), mlfSin( mclPlus( _mxarray25_, mclIntArrayRef1( mclVv(anguloLazo, "anguloLazo"), v_)))), mclPlus( mclMpower( mclMinus( mclVv(altura, "altura"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 5)), _mxarray14_), mclMpower( mclMinus( mclVv(x, "x"), mclIntArrayRef2( mclVv(mislazos, "mislazos"), v_, 4)), _mxarray14_))))); if (v_ == e_) { break; } ++v_; } mlfAssign(&m, mlfScalar(v_)); } } /* * * else */ } else { mxDestroyArray(a_); /* * * Bxrlazo=0; */ mlfAssign(&Bxrlazo, _mxarray24_); /* * Bxilazo=0; */ mlfAssign(&Bxilazo, _mxarray24_); /* * Byrlazo=0; */ mlfAssign(&Byrlazo, _mxarray24_); /* * Byilazo=0; */ mlfAssign(&Byilazo, _mxarray24_); } /* * * end */ } /* * * % Expresamos en microTeslas el campo * Bxrlazo=Bxrlazo/10; */ mlfAssign( &Bxrlazo, mclMrdivide(mclVv(Bxrlazo, "Bxrlazo"), _mxarray27_)); /* * Bxilazo=Bxilazo/10; */ mlfAssign( &Bxilazo, mclMrdivide(mclVv(Bxilazo, "Bxilazo"), _mxarray27_)); /* * Byrlazo=Byrlazo/10; */ mlfAssign( &Byrlazo, mclMrdivide(mclVv(Byrlazo, "Byrlazo"), _mxarray27_)); /* * Byilazo=Byilazo/10; */ mlfAssign( &Byilazo, mclMrdivide(mclVv(Byilazo, "Byilazo"), _mxarray27_)); /* * * % Guardamos valores en el vector * campolazo(cont,1)=x; */ mclArrayAssign2( &campolazo, mclVv(x, "x"), mclVv(cont, "cont"), _mxarray6_); /* * campolazo(cont,2)=Bxrlazo; */ mclArrayAssign2( &campolazo, mclVv(Bxrlazo, "Bxrlazo"), mclVv(cont, "cont"), _mxarray14_); /* * campolazo(cont,3)=Bxilazo; */ mclArrayAssign2( &campolazo, mclVv(Bxilazo, "Bxilazo"), mclVv(cont, "cont"), _mxarray15_); /* * campolazo(cont,4)=Byrlazo; */ mclArrayAssign2( &campolazo, mclVv(Byrlazo, "Byrlazo"), mclVv(cont, "cont"), _mxarray16_); /* * campolazo(cont,5)=Byilazo; */ mclArrayAssign2( &campolazo, mclVv(Byilazo, "Byilazo"), mclVv(cont, "cont"), _mxarray28_); /* * campolazo(cont,6)=sqrt(Bxrlazo^2+Bxilazo^2+Byrlazo^2+Byilazo^2); */ mclArrayAssign2( &campolazo, mlfSqrt( mclPlus( mclPlus( mclPlus( mclMpower(mclVv(Bxrlazo, "Bxrlazo"), _mxarray14_), mclMpower(mclVv(Bxilazo, "Bxilazo"), _mxarray14_)), mclMpower(mclVv(Byrlazo, "Byrlazo"), _mxarray14_)), mclMpower(mclVv(Byilazo, "Byilazo"), _mxarray14_))), mclVv(cont, "cont"), _mxarray29_); /* * * % Calculamos el campo resultante de la suma del creado por los * % conductores de fase y el lazo, ya vienen en microTeslas. * Bxrtotal=Bxr+Bxrlazo; */ mlfAssign( &Bxrtotal, mclPlus(mclVv(Bxr, "Bxr"), mclVv(Bxrlazo, "Bxrlazo"))); /* * Bxitotal=Bxi+Bxilazo; */ mlfAssign( &Bxitotal, mclPlus(mclVv(Bxi, "Bxi"), mclVv(Bxilazo, "Bxilazo"))); /* * Byrtotal=Byr+Byrlazo; */ mlfAssign( &Byrtotal, mclPlus(mclVv(Byr, "Byr"), mclVv(Byrlazo, "Byrlazo"))); /* * Byitotal=Byi+Byilazo; */ mlfAssign( &Byitotal, mclPlus(mclVv(Byi, "Byi"), mclVv(Byilazo, "Byilazo"))); /* * * Bxtotal=sqrt(Bxrtotal^2 + Bxitotal^2); */ mlfAssign( &Bxtotal, mlfSqrt( mclPlus( mclMpower(mclVv(Bxrtotal, "Bxrtotal"), _mxarray14_), mclMpower(mclVv(Bxitotal, "Bxitotal"), _mxarray14_)))); /* * Bytotal=sqrt(Byrtotal^2 + Byitotal^2); */ mlfAssign( &Bytotal, mlfSqrt( mclPlus( mclMpower(mclVv(Byrtotal, "Byrtotal"), _mxarray14_), mclMpower(mclVv(Byitotal, "Byitotal"), _mxarray14_)))); /* * Btotal=sqrt(Bxtotal^2+Bytotal^2); */ mlfAssign( &Btotal, mlfSqrt( mclPlus( mclMpower(mclVv(Bxtotal, "Bxtotal"), _mxarray14_), mclMpower(mclVv(Bytotal, "Bytotal"), _mxarray14_)))); /* * * campototal(cont,1)=x; */ mclArrayAssign2( &campototal, mclVv(x, "x"), mclVv(cont, "cont"), _mxarray6_); /* * campototal(cont,2)=Bxrtotal; */ mclArrayAssign2( &campototal, mclVv(Bxrtotal, "Bxrtotal"), mclVv(cont, "cont"), _mxarray14_); /* * campototal(cont,3)=Bxitotal; */ mclArrayAssign2( &campototal, mclVv(Bxitotal, "Bxitotal"), mclVv(cont, "cont"), _mxarray15_); /* * campototal(cont,4)=Byrtotal; */ mclArrayAssign2( &campototal, mclVv(Byrtotal, "Byrtotal"), mclVv(cont, "cont"), _mxarray16_); /* * campototal(cont,5)=Byitotal; */ mclArrayAssign2( &campototal, mclVv(Byitotal, "Byitotal"), mclVv(cont, "cont"), _mxarray28_); /* * campototal(cont,6)=Btotal; */ mclArrayAssign2( &campototal, mclVv(Btotal, "Btotal"), mclVv(cont, "cont"), _mxarray29_); /* * * % Calculamos campo equivalente * Bxlazo=sqrt(Bxrlazo^2+Bxilazo^2); */ mlfAssign( &Bxlazo, mlfSqrt( mclPlus( mclMpower(mclVv(Bxrlazo, "Bxrlazo"), _mxarray14_), mclMpower(mclVv(Bxilazo, "Bxilazo"), _mxarray14_)))); /* * Bylazo=sqrt(Byrlazo^2+Byilazo^2); */ mlfAssign( &Bylazo, mlfSqrt( mclPlus( mclMpower(mclVv(Byrlazo, "Byrlazo"), _mxarray14_), mclMpower(mclVv(Byilazo, "Byilazo"), _mxarray14_)))); /* * Brlazo=sqrt(Bxlazo^2+Bylazo^2); */ mlfAssign( &Brlazo, mlfSqrt( mclPlus( mclMpower(mclVv(Bxlazo, "Bxlazo"), _mxarray14_), mclMpower(mclVv(Bylazo, "Bylazo"), _mxarray14_)))); /* * * cont=cont+1; */ mlfAssign(&cont, mclPlus(mclVv(cont, "cont"), _mxarray6_)); /* * x=x+incremento; */ mlfAssign(&x, mclPlus(mclVv(x, "x"), mclVv(incremento, "incremento"))); /* * * end */ } /* * * figure */ mclPrintAns(&ans, mlfNFigure(0, NULL)); /* * plot(campo(:,1),campo(:,6)) */ mclPrintAns( &ans, mlfNPlot( 0, mclArrayRef2(mclVv(campo, "campo"), mlfCreateColonIndex(), _mxarray6_), mclArrayRef2(mclVv(campo, "campo"), mlfCreateColonIndex(), _mxarray29_), NULL)); /* * hold on */ mlfHold(_mxarray30_); /* * if(hayLazosPasivos==true | hayLazosActivos==true) */ { mxArray * a_ = mclInitialize( mclEq(mclVv(hayLazosPasivos, "hayLazosPasivos"), mlfTrue(NULL))); if (mlfTobool(a_) || mlfTobool( mclOr( a_, mclEq( mclVv(hayLazosActivos, "hayLazosActivos"), mlfTrue(NULL))))) { mxDestroyArray(a_); /* * * plot(campolazo(:,1),campolazo(:,6),'r') */ mclPrintAns( &ans, mlfNPlot( 0, mclArrayRef2( mclVv(campolazo, "campolazo"), mlfCreateColonIndex(), _mxarray6_), mclArrayRef2( mclVv(campolazo, "campolazo"), mlfCreateColonIndex(), _mxarray29_), _mxarray32_, NULL)); /* * plot(campototal(:,1),campototal(:,6),'g') */ mclPrintAns( &ans, mlfNPlot( 0, mclArrayRef2( mclVv(campototal, "campototal"), mlfCreateColonIndex(), _mxarray6_), mclArrayRef2( mclVv(campototal, "campototal"), mlfCreateColonIndex(), _mxarray29_), _mxarray34_, NULL)); } else { mxDestroyArray(a_); } /* * * end */ } /* * * title('Campo resultante Br') */ mclPrintAns(&ans, mlfNTitle(0, _mxarray36_, NULL)); /* * xlabel('x (m)') */ mclPrintAns(&ans, mlfNXlabel(0, _mxarray38_, NULL)); /* * ylabel('Br (uT)') */ mclPrintAns(&ans, mlfNYlabel(0, _mxarray40_, NULL)); /* * * if(hayLazosPasivos==true | hayLazosActivos==true) */ { mxArray * a_ = mclInitialize( mclEq(mclVv(hayLazosPasivos, "hayLazosPasivos"), mlfTrue(NULL))); if (mlfTobool(a_) || mlfTobool( mclOr( a_, mclEq( mclVv(hayLazosActivos, "hayLazosActivos"), mlfTrue(NULL))))) { mxDestroyArray(a_); /* * * legend('Br linea','Br lazo','Br neto',0) */ mclPrintAns( &ans, mlfNLegend( 0, NULL, _mxarray42_, _mxarray44_, _mxarray46_, _mxarray24_, NULL)); /* * [a,b]=size(campo); */ mlfSize(mlfVarargout(&a, &b, NULL), mclVv(campo, "campo"), NULL); /* * for ju=1:a */ { int v_ = mclForIntStart(1); int e_ = mclForIntEnd(mclVv(a, "a")); if (v_ > e_) { mlfAssign(&ju, _mxarray9_); } else { /* * * FA(ju,1)=campo(ju,6)/campototal(ju,6); * * end */ for (; ; ) { mclIntArrayAssign2( &FA, mclMrdivide( mclIntArrayRef2(mclVv(campo, "campo"), v_, 6), mclIntArrayRef2( mclVv(campototal, "campototal"), v_, 6)), v_, 1); if (v_ == e_) { break; } ++v_; } mlfAssign(&ju, mlfScalar(v_)); } } /* * * hold off */ mlfHold(_mxarray48_); /* * figure */ mclPrintAns(&ans, mlfNFigure(0, NULL)); /* * plot(campo(:,1),FA(:,1)) */ mclPrintAns( &ans, mlfNPlot( 0, mclArrayRef2( mclVv(campo, "campo"), mlfCreateColonIndex(), _mxarray6_), mclArrayRef2( mclVv(FA, "FA"), mlfCreateColonIndex(), _mxarray6_), NULL)); /* * xlabel('x (m)') */ mclPrintAns(&ans, mlfNXlabel(0, _mxarray38_, NULL)); /* * ylabel('Factor de Apantallado') */ mclPrintAns(&ans, mlfNYlabel(0, _mxarray50_, NULL)); /* * * else */ } else { mxDestroyArray(a_); /* * * legend('Brlinea'); */ mclAssignAns(&ans, mlfNLegend(0, NULL, _mxarray52_, NULL)); } /* * * end */ } /* * * % Por si se necesitan luego para Excel, se graban las matrices de campo y FA a * % fichero * * fid=fopen(strcat(path,'\campo.txt'),'w'); */ mlfAssign( &fid, mlfFopen( NULL, NULL, mlfStrcat(mclVv(path, "path"), _mxarray54_, NULL), _mxarray56_, NULL)); /* * [m,n]=size(campo); */ mlfSize(mlfVarargout(&m, &n, NULL), mclVv(campo, "campo"), NULL); /* * for u=1:m */ { int v_ = mclForIntStart(1); int e_ = mclForIntEnd(mclVv(m, "m")); if (v_ > e_) { mlfAssign(&u, _mxarray9_); } else { /* * * fprintf(fid,'%6.2f, %6.2f, %6.2f, %6.2f \r',campo(u,1),campo(u,6),campolazo(u,6),campototal(u,6)); * * end */ for (; ; ) { mclAssignAns( &ans, mlfNFprintf( 0, mclVv(fid, "fid"), _mxarray58_, mclIntArrayRef2(mclVv(campo, "campo"), v_, 1), mclIntArrayRef2(mclVv(campo, "campo"), v_, 6), mclIntArrayRef2(mclVv(campolazo, "campolazo"), v_, 6), mclIntArrayRef2(mclVv(campototal, "campototal"), v_, 6), NULL)); if (v_ == e_) { break; } ++v_; } mlfAssign(&u, mlfScalar(v_)); } } /* * fclose(fid); */ mclAssignAns(&ans, mlfFclose(mclVv(fid, "fid"))); /* * * if(hayLazosPasivos==true | hayLazosActivos==true) */ { mxArray * a_ = mclInitialize( mclEq(mclVv(hayLazosPasivos, "hayLazosPasivos"), mlfTrue(NULL))); if (mlfTobool(a_) || mlfTobool( mclOr( a_, mclEq( mclVv(hayLazosActivos, "hayLazosActivos"), mlfTrue(NULL))))) { mxDestroyArray(a_); /* * * fid=fopen(strcat(path,'\fa.txt'),'w'); */ mlfAssign( &fid, mlfFopen( NULL, NULL, mlfStrcat(mclVv(path, "path"), _mxarray60_, NULL), _mxarray56_, NULL)); /* * [m,n]=size(campo); */ mlfSize(mlfVarargout(&m, &n, NULL), mclVv(campo, "campo"), NULL); /* * for u=1:m */ { int v_ = mclForIntStart(1); int e_ = mclForIntEnd(mclVv(m, "m")); if (v_ > e_) { mlfAssign(&u, _mxarray9_); } else { /* * * fprintf(fid,'%6.2f, %6.4f \r',campo(u,1),FA(u,1)); * * end */ for (; ; ) { mclAssignAns( &ans, mlfNFprintf( 0, mclVv(fid, "fid"), _mxarray62_, mclIntArrayRef2(mclVv(campo, "campo"), v_, 1), mclIntArrayRef2(mclVv(FA, "FA"), v_, 1), NULL)); if (v_ == e_) { break; } ++v_; } mlfAssign(&u, mlfScalar(v_)); } } /* * fclose(fid); */ mclAssignAns(&ans, mlfFclose(mclVv(fid, "fid"))); } else { mxDestroyArray(a_); } /* * end */ } mxDestroyArray(ans); mxDestroyArray(fid); mxDestroyArray(path); mxDestroyArray(cadena); mxDestroyArray(fdatos); mxDestroyArray(cont); mxDestroyArray(indices); mxDestroyArray(ii); mxDestroyArray(x); mxDestroyArray(y); mxDestroyArray(intens); mxDestroyArray(desf); mxDestroyArray(fases); mxDestroyArray(f); mxDestroyArray(altura); mxDestroyArray(xmin); mxDestroyArray(xmax); mxDestroyArray(incremento); mxDestroyArray(fidPas); mxDestroyArray(hayLazosPasivos); mxDestroyArray(hayLazosActivos); mxDestroyArray(lazosPas); mxDestroyArray(I); mxDestroyArray(mislazos); mxDestroyArray(Ilazo); mxDestroyArray(anguloLazo); mxDestroyArray(fidAct); mxDestroyArray(lazosAct); mxDestroyArray(Bxr); mxDestroyArray(Bxi); mxDestroyArray(Byr); mxDestroyArray(Byi); mxDestroyArray(m); mxDestroyArray(Bx); mxDestroyArray(By); mxDestroyArray(Br); mxDestroyArray(campo); mxDestroyArray(Bxrlazo); mxDestroyArray(Bxilazo); mxDestroyArray(Byrlazo); mxDestroyArray(Byilazo); mxDestroyArray(campolazo); mxDestroyArray(Bxrtotal); mxDestroyArray(Bxitotal); mxDestroyArray(Byrtotal); mxDestroyArray(Byitotal); mxDestroyArray(Bxtotal); mxDestroyArray(Bytotal); mxDestroyArray(Btotal); mxDestroyArray(campototal); mxDestroyArray(Bxlazo); mxDestroyArray(Bylazo); mxDestroyArray(Brlazo); mxDestroyArray(a); mxDestroyArray(b); mxDestroyArray(ju); mxDestroyArray(FA); mxDestroyArray(n); mxDestroyArray(u); mclSetCurrentLocalFunctionTable(save_local_function_table_); /* * * */ } e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.


PROGRAMA DE DISEÑO DE LAZOS ACTIVOS/PASIVOS PARA AMORTIGUAMIENTO DEL CAMPO MAGNÉTICO EN LÍNEAS ELÉCTRICAS AÉREAS DE ALTA TENSIÓN

: CASAS PÉREZ, JOSÉ EDUARDO DE LAS
: Ingeniería Industrial
Contenido del proyecto:
Directorio raíz  >  Programa  >  Codigo  >  calculos  >  micampo.c