void LatexInsert3DPlot(klMatrix<double>& mat, ofstream &_tex, string dir,string filename,string title,klHoldOnStatus holdon,const char* markerType)
{
	klMatlabEngineThreadMap klmtm;
	Engine* matlabEngine=klmtm.find(klThread<klMutex>::getCurrentThreadId() );
	char* arg = new char[2048];
	char* evalString = new char[2048];
	sprintf(arg,"%s//%s.eps",dir.c_str(),filename.c_str());

	const char* xAxis=NULL;
	const char* yAxis=NULL;
	const char* zAxis=NULL;
	bool useExtents=true;
	
	klScatterPlot3D(mat,filename.c_str(),title.c_str(), xAxis, yAxis,zAxis, useExtents, holdon, markerType);

	if(holdon==klHoldOnStatus::LastPlot  || holdon==klHoldOnStatus::NoHold )
	{
		sprintf(evalString,"print -r1200 -depsc %s;",arg);
		engEvalString(matlabEngine, evalString);

		engEvalString(matlabEngine, "hold off;close(gcf);");

		sprintf(evalString,"epstopdf   %s",arg);

		system(evalString);
		sprintf(arg,"%s.pdf",filename.c_str());
		_tex<<"\\includegraphics[width=10.0cm,height=10.0cm]{"<<arg<<"}"<<endl<<endl;
	}
	delete arg;
	delete evalString;
}
예제 #2
0
파일: test.c 프로젝트: kristinbranson/JAABA
int main(int argc, const char *argv[])

{
  Engine *ep;
  char buff[10240];
  int i;

  /* matlab must be in the PATH! */
  if (!(ep = engOpen("matlab -nodisplay"))) {
    fprintf(stderr, "Can't start MATLAB engine\n");
    return -1;
  }
  engOutputBuffer(ep, buff, 10239);

  /* load the mex file */
  if(argc<2){
    fprintf(stderr, "Error. Give full path to the MEX file as input parameter.\n");
    return -1;
  }
  void *handle = dlopen(argv[1], RTLD_NOW);
  if(!handle){
    fprintf(stderr, "Error loading MEX file: %s\n", strerror(errno));
    return -1;
  }

  /* grab mexFunction handle */
  mexFunction_t mexfunction = (mexFunction_t)dlsym(handle, "mexFunction");
  if(!mexfunction){
    fprintf(stderr, "MEX file does not contain mexFunction\n");
    return -1;
  }

  /* load input data - for convenience do that using MATLAB engine */
  /* NOTE: parameters are MEX-file specific, so one has to modify this*/
  /* to fit particular needs */
  engEvalString(ep, "load input.mat");
  mxArray *arg1 = engGetVariable(ep, "im1");
  mxArray *arg2 = engGetVariable(ep, "im2");
  mxArray *arg3 = engGetVariable(ep, "szx");
  mxArray *arg4 = engGetVariable(ep, "szy");
  mxArray *arg5 = engGetVariable(ep, "ngh");
  mxArray *pargout[1] = {0};
  const mxArray *pargin[5] = {arg1, arg2,arg3, arg4,arg5};

  /* execute the mex function */
  mexfunction(1, pargout, 5, pargin);

  /* print the results using MATLAB engine */
  engPutVariable(ep, "result", pargout[0]);
  engEvalString(ep, "result");
  printf("%s\n", buff);

  /* cleanup */
  mxDestroyArray(pargout[0]);
  engEvalString(ep, "clear all;");
  dlclose(handle);
  engClose(ep);

  return 0;
}
void LatexInsertHistogram(klVector<double>& vec, unsigned int numBins,ofstream &_tex, string dir,string filename,string title)
{
	/*vec.setupRange();
	klVector<double> hist=vec.histogram(numBins,vec.y0,vec.y1);
	hist.setupRange();
	hist.setupDomain(vec.y0,vec.y1);*/

	klMatlabEngineThreadMap klmtm;
	Engine* matlabEngine=klmtm.find(klThread<klMutex>::getCurrentThreadId() );
	char* arg = new char[2048];
	char* evalString = new char[2048];
	sprintf(arg,"%s\\%s.eps",dir.c_str(),filename.c_str());
	
	//klPlot1D<double>(hist,arg,title.c_str());

	klPlotHistogram<double>(vec,arg,title.c_str());
		
	sprintf(evalString,"print -r1200 -depsc %s;",arg);
	engEvalString(matlabEngine, evalString);
	engEvalString(matlabEngine, "hold off;close(gcf);");
	sprintf(evalString,"epstopdf   %s",arg);
	system(evalString);
	sprintf(arg,"%s.pdf",filename.c_str());
	_tex<<"\\includegraphics[width=10.0cm,height=10.0cm]{"<<arg<<"}"<<endl<<endl;
	delete arg;
	delete evalString;
}
예제 #4
0
파일: demoeng.cpp 프로젝트: haishanwu/cnpkg
bool Eval(Engine *ep, char *out, const char *cmd, ...) {

    char buf[BUFSIZE];
    va_list args;
    va_start(args, cmd);
    vsprintf(buf, cmd, args);
    va_end(args);

    engEvalString(ep, "clear engerr; lasterr('');");

    engEvalString(ep, buf);
    printf("%s", out);

    engEvalString(ep, "engerr = lasterr;");

    mxArray *err = engGetVariable(ep, "engerr");
    if (err == NULL) {
        printf("unable to retrieve variable 'engerr'\n\n");
        return false;
    }
    bool ok = (mxGetNumberOfElements(err) == 0);
    mxDestroyArray(err);

    return ok;

}
예제 #5
0
int PASCAL WinMain (HINSTANCE hInstance,
                    HINSTANCE hPrevInstance,
                    LPSTR     lpszCmdLine,
                    int       nCmdShow){
	char buffer[BUFSIZE];
	mxArray *app;
	Engine *ep;

	/* 啟動 MATLAB 引擎 */
	if (!(ep = engOpen(NULL))){		// 產生一個 MATLAB 引擎物件
		MessageBox ((HWND)NULL, (LPSTR)"Can't start MATLAB engine",
			(LPSTR)"plotViaMatlab01.c", MB_OK);
		exit(-1);
	}

	/* 切換目錄並執行 plotSine.m */
	_getcwd(buffer, BUFSIZE);		// 將此程式所在目錄存入字串 buffer
	app = mxCreateString(buffer);		// 產生 MATLAB 內部的字串變數 app
	engPutVariable(ep, "appDir", app);	// 將字串變數 app 至入工作空間的變數 appDir
	engEvalString(ep, "cd(appDir)");	// 將 MATLAB 的工作目錄切換至字串 appDir 所指定的目錄
	engEvalString(ep, "plotSine");		// 執行同目錄下的 plotSine.m
	
	/* 取得 MATLAB 輸出訊息 */
	engOutputBuffer(ep, buffer, BUFSIZE);	// 設定 buffer 可以接收 MATLAB 的輸出訊息
	engEvalString(ep, "whos");		// 在 MATLAB 引擎執行 whos 指令
	MessageBox((HWND)NULL, (LPSTR)buffer, (LPSTR) "MATLAB - whos", MB_OK);		// 顯示 buffer 的內容

	engClose(ep);				// 最後關閉 MATLAB 引擎
	return(0);
}
예제 #6
0
    void MatlabAdapter::initModel()
    {
        mEngine = engOpen(mEngineArgs.size() ?  mEngineArgs.toStdString().c_str(): "\0");
        engOutputBuffer(mEngine, mBuffer, mBufferSize);

        if(mEngine == 0)
        {
            TRACE("warning: connection failed!");
        }

        const QString modelDir = "cd '" +  mWorkingDir + "/Projects/MMM" + "'";
        TRACE("working dir is: %s", modelDir.toStdString().c_str());

        int retVal = 0;
        retVal += engEvalString(mEngine, modelDir.toStdString().c_str());
        retVal += engEvalString(mEngine, "initializeModelAdapter");

        // *************** read out dofs ****************************
        retVal += engEvalString(mEngine, "numDOFs = length(Model.q0);");
        mxArray* mxNumDOFs = 0;

        if((mxNumDOFs =engGetVariable(mEngine, "numDOFs")) == 0)
        {
            TRACE("WARNING: can't evaluate degrees of freedom");
            mxDestroyArray(mxNumDOFs);
            return;
        }

        tree.dofs.resize(static_cast<int>(mxGetPr(mxNumDOFs)[0]));
        TRACE("model has %lu degrees of freedom", tree.dofs.size());

        // ************** read out segment / joint information *************
        retVal += engEvalString(mEngine, "numSegs = size(Model.Body,2);");
        retVal += engEvalString(mEngine, "numJoints = size(Model.Joint,2);");
        mxArray *mxNumSegs = 0, *mxNumJoints = 0;

        if((mxNumSegs =engGetVariable(mEngine, "numSegs")) == 0) {
            TRACE("WARNING: can't evaluate number of segments");
            mxDestroyArray(mxNumSegs);
            return;
        }

        if((mxNumJoints =engGetVariable(mEngine, "numJoints")) == 0) {
            TRACE("WARNING: can't evaluate number of joints");
            mxDestroyArray(mxNumJoints);
            return;
        }
        tree.jointNodes.resize(static_cast<int>(mxGetPr(mxNumJoints)[0]));
        tree.partNodes.resize(static_cast<int>(mxGetPr(mxNumSegs)[0]));
        //model.numJoints = static_cast<int>(mxGetPr(mxNumJoints)[0]);
        //model.numSegments = static_cast<int>(mxGetPr(mxNumSegs)[0]);

        mxDestroyArray(mxNumDOFs);
        mxDestroyArray(mxNumJoints);
        mxDestroyArray(mxNumSegs);

    }
WaveData SignalProcessTools::waveRead(string filename)
{
	engEvalString(ep, ("[R,Fs]=audioread('"+filename+"');").c_str());
	engEvalString(ep, "R=R(:,1);");
	WaveData data;
	getVector(data.wave, "R");
	mxArray *Fs = engGetVariable(ep, "Fs");
	data.sampleRate = (int)mxGetPr(Fs)[0];
	return data;
}
void SignalProcessTools::plot(vector dataX, vector dataY, \
							  string title, string xLabel, string yLabel)
{
	putVector(dataX, "X");
	putVector(dataY, "Y");
	engEvalString(ep, "plot(X,Y);");
	if (title != "")  engEvalString(ep, ("title('"+title+"');").c_str());
	if (xLabel != "")  engEvalString(ep, ("xlabel('"+xLabel+"');").c_str());
	if (yLabel != "")  engEvalString(ep, ("ylabel('"+yLabel+"');").c_str());
}
void SignalProcessTools::plotSpec(FreqData data)
{
	int n1 = data.freq.length;
	int n2 = data.time.length;
	put2DData(data.real.data, n1, n2, "BReal");
	put2DData(data.imag.data, n1, n2, "BImag");
	string cmd = "B=BReal+i*BImag;imagesc(time,freq,abs(B));";
	engEvalString(ep, cmd.c_str());
	cmd = "axis('xy');xlabel('ʱ¼ä (Ãë)'); ylabel('ƵÂÊ (Hz)');";
	engEvalString(ep, cmd.c_str());
}
예제 #10
0
PyObject * mlabraw_eval(PyObject *, PyObject *args)
{
  //XXX how large should this be?
  const int  BUFSIZE=10000;
  char* fmt = "try, %s; MLABRAW_ERROR_=0; catch, MLABRAW_ERROR_=1; end;";
  char buffer[BUFSIZE];
  char cmd[BUFSIZE];
  char *lStr;
  char *retStr = buffer;
  PyObject *ret;
  PyObject *lHandle;
  if (! PyArg_ParseTuple(args, "Os:eval", &lHandle, &lStr)) return NULL;
  if (! PyCObject_Check(lHandle)) {
    PyErr_SetString(PyExc_TypeError, "Invalid object passed as mlabraw session handle");
    return NULL;
  }
  sprintf(cmd, fmt, lStr); //FIXME check buffer overflow
  // std::cout << "DEBUG: CMD " << cmd << std::endl << std::flush; 
  engOutputBuffer((Engine *)PyCObject_AsVoidPtr(lHandle), retStr, BUFSIZE-1);
  if (engEvalString((Engine *)PyCObject_AsVoidPtr(lHandle), cmd) != 0) {
    PyErr_SetString(mlabraw_error, 
                   "Unable to evaluate string in MATLAB(TM) workspace");
    return NULL;
  }
  {
    mxArray *lArray = NULL; 
    char buffer2[BUFSIZE];
    char *retStr2 = buffer2;
    bool __mlabraw_error;
    if (NULL == (lArray = _getMatlabVar(lHandle, "MLABRAW_ERROR_")) ) {
      PyErr_SetString(mlabraw_error, 
                      "Something VERY BAD happened whilst trying to evaluate string " 
                      "in MATLAB(TM) workspace.");
      return NULL;
    }
    __mlabraw_error = (bool)*mxGetPr(lArray);
    mxDestroyArray(lArray);
    if (__mlabraw_error) {
      engOutputBuffer((Engine *)PyCObject_AsVoidPtr(lHandle), retStr2, BUFSIZE-1);
      if (engEvalString((Engine *)PyCObject_AsVoidPtr(lHandle), 
                        "disp(subsref(lasterror(),struct('type','.','subs','message')))") != 0) {
        PyErr_SetString(mlabraw_error, "THIS SHOULD NOT HAVE HAPPENED!!!");
        return NULL;
      }
      PyErr_SetString(mlabraw_error, retStr2 + ((strncmp(">> ", retStr2, 3) == 0) ?  3 : 0));
      return NULL;
    }
  }
  if (strncmp(">> ", retStr, 3) == 0) { retStr += 3; } //FIXME
  ret = (PyObject *)PyString_FromString(retStr);
  return ret;
}
예제 #11
0
    void MatlabAdapter::initGenCoords()
    {
        engEvalString(mEngine, "coordNames = {Model.GenCoor.Name}");
        engEvalString(mEngine, "coordRanges = [ Model.lb; Model.ub]");

        mxArray *mxNames = engGetVariable(mEngine, "coordNames");
        mxArray *mxRanges = engGetVariable(mEngine, "coordRanges");

        double* values = mxGetPr(mxRanges);
        double min, max, stepSize;

        unsigned size = tree.dofs.size();
        for(unsigned i=0; i < size; i++)
        {
            mxArray *mxCell = mxGetCell(mxNames, i);
            char* name = mxArrayToString(mxCell);

            min = values[i];
            max = values[i+size];
            stepSize = (max-min) / 100.0;
            if(stepSize > 0.05) stepSize = 0.05;

            DegreeOfFreedom dof(min, max, 0, stepSize);
            dof.name = name;

            if(dof.name.size() > 6)
            {
            	if(dof.name.substr(0,6).compare("root_r") == 0)
            	{
            		TRACE("limited");
            		dof.rangeType = DegreeOfFreedom::limitedRange;
            		dof.min = -5;
            		dof.max = 5;
            	} else
            	{
            		TRACE("circular");
            		dof.rangeType = DegreeOfFreedom::circularRange;
            		dof.min = - M_PI;
            		dof.max = + M_PI;
            	}
            }

            tree.dofs.at(i)=dof;
            TRACE("generalized coordinat %s @ %i : min %f to %f", name, i, min, max);
        }


        mxDestroyArray(mxNames);
        mxDestroyArray(mxRanges);
    }
예제 #12
0
MATLABPLOTTER::MATLABPLOTTER()
{
	//stuff	
	// Start the MATLAB engine 
	
	if (!(ep = engOpen(NULL))) {
		MessageBox ((HWND)NULL, (LPSTR)"Can't start MATLAB engine", 
			(LPSTR) "Matlab.cpp", MB_OK);
		exit(-1);
	}
	engEvalString(ep, "iter=1;");
	engEvalString(ep, "color_codes=['r' 'g' 'b' 'k' 'y'];");
	

}
예제 #13
0
void MATLABPLOTTER::plotfreqscan(std::vector <double> &FREQ_vector, std::vector <double> &DAQ_vector)
{
	
	mxArray *data_freq = NULL;
	mxArray *data_DAQ = NULL;
		//Scan laser


		double *freq_data_ptr = (double*) calloc(FREQ_vector.size(), sizeof(double));
		if (freq_data_ptr==NULL) exit(1);
		double *DAQ_data_ptr = (double*) calloc(DAQ_vector.size(), sizeof(double));
		if (DAQ_data_ptr==NULL) exit(1);

	for(unsigned int i=0;i < DAQ_vector.size(); i++) {
		double stra = FREQ_vector.at(i);
		
		double strb = DAQ_vector.at(i);
		freq_data_ptr[i] = stra;
		DAQ_data_ptr[i] = strb;
		}


	data_freq = mxCreateDoubleMatrix(1, FREQ_vector.size(), mxREAL);
	memcpy((char*) mxGetPr(data_freq), freq_data_ptr, FREQ_vector.size()*sizeof(double));

	data_DAQ = mxCreateDoubleMatrix(1, FREQ_vector.size(), mxREAL);
	memcpy((char*) mxGetPr(data_DAQ), DAQ_data_ptr, FREQ_vector.size()*sizeof(double));

	/*
	 * Place the variable dataFreq & dataDAQ into the MATLAB workspace
	 */

	engPutVariable(ep, "dataFreq", data_freq);
	engPutVariable(ep, "dataDAQ", data_DAQ);

	/* Plot the result
	 */
	//engEvalString(ep, "hold on;");
	engEvalString(ep, "figure(1);");
	engEvalString(ep, "plot(dataFreq(:)',dataDAQ(:)');");
	engEvalString(ep, "title('Rb Scan');");
	engEvalString(ep, "xlabel('Frequency (GHz)');");
	engEvalString(ep, "ylabel('Absorbtion (V)');");

	free(DAQ_data_ptr);
	free(freq_data_ptr);

}
예제 #14
0
void MATLABPLOTTER::plotlockpoints(std::vector <double> &FITFREQ_vector, std::vector <double> &FITDAQ_vector)
{
	
	mxArray *fit_data_freq = NULL;
	mxArray *fit_data_DAQ = NULL;
		//Scan laser


		double *fit_freq_data_ptr = (double*) calloc(FITFREQ_vector.size(), sizeof(double));
		if (fit_freq_data_ptr==NULL) exit(1);
		double *fit_DAQ_data_ptr = (double*) calloc(FITDAQ_vector.size(), sizeof(double));
		if (fit_DAQ_data_ptr==NULL) exit(1);

	for(unsigned int i=0;i < FITDAQ_vector.size(); i++) {
		double stra = FITFREQ_vector.at(i);
		
		double strb = FITDAQ_vector.at(i);
		fit_freq_data_ptr[i] = stra;
		fit_DAQ_data_ptr[i] = strb;
		}


	fit_data_freq = mxCreateDoubleMatrix(1, FITFREQ_vector.size(), mxREAL);
	memcpy((char*) mxGetPr(fit_data_freq), fit_freq_data_ptr, FITFREQ_vector.size()*sizeof(double));

	fit_data_DAQ = mxCreateDoubleMatrix(1, FITFREQ_vector.size(), mxREAL);
	memcpy((char*) mxGetPr(fit_data_DAQ), fit_DAQ_data_ptr, FITFREQ_vector.size()*sizeof(double));

	/*
	 * Place the variable dataFreq & dataDAQ into the MATLAB workspace
	 */

	engPutVariable(ep, "fitdataFreq", fit_data_freq);
	engPutVariable(ep, "fitdataDAQ", fit_data_DAQ);

	/* Plot the result
	 */
	engEvalString(ep, "i=1");
	for(unsigned int j=1; j <= FITFREQ_vector.size(); j++) {
		engEvalString(ep, "hold on;");
		engEvalString(ep, "plot(fitdataFreq(i)',fitdataDAQ(i)', 's','MarkerSize',10, 'MarkerFaceColor','g','MarkerEdgeColor', 'r');");
		engEvalString(ep, "i=i+1;");
	}

	free(fit_DAQ_data_ptr);
	free(fit_freq_data_ptr);

}
예제 #15
0
파일: uav.cpp 프로젝트: mvsframework/mvs
void Uav::initializeMatlabEngine()
{
  matlabEngine = engOpen(NULL);
    
  if( !matlabEngine )
  { qCritical("%s : Cannot open the MATLAB engine.",
                objectName().toStdString().c_str() ); }; 
                
 
  // prepare an empty matrix to hold the ownship's track
  engEvalString(matlabEngine, "ownshipTrack = [];");
  
  // prepare a plot figure to show that track
  engEvalString(matlabEngine, "figure('Name','UAV Track');");
  engEvalString(matlabEngine, "axis equal;");
}
void MatLabEngine::executeCommand( const std::string &command, const bool unchecked ) {

	if ( _matLabLogFileStreamPtr ) {
		std::cerr << command << std::endl;
		(*_matLabLogFileStreamPtr) << "try\n\t" << command << "\nend" << std::endl;
	}

	if (  engEvalString( _ep, command.c_str() )  !=  0  ) {
		_ep = 0;
		throw MatLabUdm::Exception( "MatLab connection lost." );
	}

	_matLabOutput = getBuffer();
	if (  _matLabOutput.find( "???" ) == 0  ) {
#ifdef _DEBUG
	// KMS: is there a better way to check for an error???
		std::string error = "Executing MatLab command '";
		error += command + "' resulted in error '" + _matLabOutput;
		if (!unchecked) {
#ifdef _WIN32
			OutputDebugString(error.c_str());
			DebugBreak();
#endif
			throw MatLabUdm::Exception(error);
		} else {
#ifdef _WIN32
			OutputDebugString(error.c_str());
#endif
		}
#endif
	}

}
예제 #17
0
int
MatlabEvaluator::runAnalysis(const Vector &x)
{	
	
	// Let's just make a direct call since we have the pointer to OpenSees domain
	// This replaces above call to Tcl command; however, in the reset command
	// revertToStart() is also called on theTransientIntegrator -- MHS needs to check
	if (theOpenSeesDomain->revertToStart() != 0) {
		opserr << "ERROR MatlabEvaluator -- error in resetting Domain" << endln;
		return -1;
	}
	
	// Put random variables into the structural domain according to the RandomVariablePositioners
	int rvIndex;
	RandomVariablePositionerIter rvPosIter = theReliabilityDomain->getRandomVariablePositioners();
	RandomVariablePositioner *theRVPos;
	while ((theRVPos = rvPosIter()) != 0) {
		rvIndex = theRVPos->getRvIndex();
		theRVPos->update(x(rvIndex));
	}
	
	// Start a Matlab engine
	Engine *ep;
	ep = engOpen("\0");
	
	// Execute a Matlab function called 'matlabgfun'
	char theMatlabCommand[50];
	sprintf(theMatlabCommand,"matlabgfun");
	engEvalString(ep, theMatlabCommand);
	
	// Shut down the Matlab engine
	engClose(ep);
	
	return 0;
}
예제 #18
0
파일: cmd.c 프로젝트: shanecelis/mEngine
void engcmd(const char* command)
{
	char buffer[BUFSIZE];	//MatLab output buffer
	bool SUCCESS = true;		//success flag

	if (NULL == Eng)	//if not opened yet
	{
		msg("eng::noMLB");	//message 
		SUCCESS = false;
	}
	else
	{
		engOutputBuffer(Eng, buffer, BUFSIZE);	//for return output
		//issue command
		if(engEvalString(Eng, command))	//if unsucessful
		{
			msg("engCmd::erexe");
			SUCCESS = false;
		}
	}

	if(SUCCESS)
		MLPutByteString(stdlink, buffer, strlen(buffer));
	else
		MLPutSymbol(stdlink, "$Failed");

}
void MatLabEngine::executeInteractiveCommands( void ) {

	std::string command;
	while( true ) {

		char *commandString = readline( ">> " );
		if ( !commandString ) {
			continue;
		}
		if ( *commandString ) {
			add_history( commandString );
		}
		command = commandString;
		free( commandString );

		if (  command.substr( 0, 4 ) == "quit"  ) {
			break;
		}

		zeroBuffer();

		if (  engEvalString( _ep, command.c_str() )  !=  0  ) {
			_ep = 0;
			throw MatLabUdm::Exception( "MatLab connection lost." );
		}
		std::string output( getBuffer() );
		output = output.substr( output.find( '\n' ) + 1 );
		std::cout << output << std::endl;
	}
}
FreqData SignalProcessTools::stft(WaveData data, string windowName /* = "hamming" */, 
					int windowN /* = 256 */, int NOverlap /* = 192 */, int NFFT /* = 256 */)
{
	putWaveData(data);
	string cmd = "h=window('"+windowName+"',"+intToString(windowN)+")";
	engEvalString(ep, cmd.c_str());
	cmd = "[B,freq,time]=spectrogram(W,h,"+intToString(NOverlap)+","+intToString(NFFT)+",Fs);";
	engEvalString(ep, cmd.c_str());
	cmd = "BReal=real(B);BImag=imag(B);";
	engEvalString(ep, cmd.c_str());
	FreqData freq;
	getVector(freq.freq, "freq");
	getVector(freq.time, "time");
	getVector(freq.real, "BReal");
	getVector(freq.imag, "BImag");
	return freq;
}
예제 #21
0
void MATLABPLOTTER::plotData(std::vector <double> &timeVector, std::vector <double> &signalVector, bool new_figure)
{
	
	mxArray *data_time = NULL;
	mxArray *data_Signal = NULL;
		//Scan laser


		double *time_data_ptr = (double*) calloc(timeVector.size(), sizeof(double));
		if (time_data_ptr==NULL) exit(1);
		double *Signal_data_ptr = (double*) calloc(signalVector.size(), sizeof(double));
		if (Signal_data_ptr==NULL) exit(1);

	for(unsigned int i=0;i < signalVector.size(); i++) {
		double stra = timeVector.at(i);
		
		double strb = signalVector.at(i);
		time_data_ptr[i] = stra;
		Signal_data_ptr[i] = strb;
		}


	data_time = mxCreateDoubleMatrix(1, timeVector.size(), mxREAL);
	memcpy((char*) mxGetPr(data_time), time_data_ptr, timeVector.size()*sizeof(double));

	data_Signal = mxCreateDoubleMatrix(1, timeVector.size(), mxREAL);
	memcpy((char*) mxGetPr(data_Signal), Signal_data_ptr, timeVector.size()*sizeof(double));

	
	 // Place the variable dataTime & dataSignal into the MATLAB workspace
	 

	engPutVariable(ep, "dataTime", data_time);
	engPutVariable(ep, "dataSignal", data_Signal);

	// Plot the result
	 
	
	if (new_figure) {
		engEvalString(ep, "close;");
		engEvalString(ep, "figure;");
	}
	else {
	engEvalString(ep, "hold on;");
	}
	
	engEvalString(ep, "plot(dataTime(:)',dataSignal(:)', color_codes(round(rand*4+1)) );");
	engEvalString(ep, "title('Scope Trace');");
	engEvalString(ep, "xlabel('Time (s)');");
	engEvalString(ep, "ylabel('Volts');");
	engEvalString(ep, "iter=iter+1;");

	free(Signal_data_ptr);
	free(time_data_ptr);

}
void LatexInsertLegend(string legend)
{
	klMatlabEngineThreadMap klmtm;
	Engine* matlabEngine=klmtm.find(klThread<klMutex>::getCurrentThreadId() );
	stringstream evalString;
	evalString<<"legend("<<legend<<");";
	engEvalString(matlabEngine, evalString.str().c_str());
	
}
예제 #23
0
파일: Matlab.cpp 프로젝트: jasonhogan/sti
void MATLABPLOTTER::saveVector(std::string filename)
{
		
	std::string path = "\\\\epsrv1\\EP\\data\\ThermocoupleAutoSave\\";
	std::string save_command = "dlmwrite('" + path + filename + "', inputData, '-append');";

	
	engEvalString(ep, save_command.c_str());

}
예제 #24
0
파일: uav.cpp 프로젝트: mvsframework/mvs
void Uav::useMatlabEngine()
{
  Q_ASSERT(matlabEngine);
  
  /* 
   * Create a variable for our data
   */
  mxArray* pos2d = mxCreateDoubleMatrix(1, 2, mxREAL);
  
  
  //NOTE: the UAV's (x,y,z) coordinates are given wrt. to the 
  // North/x-East/y-Down/y frame of reference.
  // Matlab's plot function assume the standart mathematical 
  // X/East-Y/North-Z/Up frame of reference...
  
  memcpy(static_cast<void*>(mxGetPr(pos2d)),
         static_cast<const void*>( &(ownship().position2D().y()) ),
         sizeof(ownship().position2D().x())) ;
         
  memcpy(static_cast<void*>(mxGetPr(pos2d)+1 ),
         static_cast<const void*>( &(ownship().position2D().x()) ),
         sizeof(ownship().position2D().y())) ;         
        
  /*
   * Place the variable "pos2d" into the MATLAB workspace
   */
  engPutVariable(matlabEngine, "pos2d", pos2d);
  
  // append the current position to the track
  engEvalString(matlabEngine, "ownshipTrack = [ownshipTrack; pos2d];");
  
  /*
   * Plot the result
   */
  engEvalString(matlabEngine, "plot(ownshipTrack(:,1),ownshipTrack(:,2),'-b',pos2d(1),pos2d(2),'+r');");
  
  /*
   * Clean up
   */
  mxDestroyArray(pos2d);
}
void LatexInsertHeatMap(klMatrix<double>& mat, ofstream &_tex, string dir,string filename,string title)
{
	klMatlabEngineThreadMap klmtm;

	Engine* matlabEngine=klmtm.find(klThread<klMutex>::getCurrentThreadId() );

	char* arg = new char[2048];
	char* evalString = new char[2048];
	sprintf(arg,"%s//%s.eps",dir.c_str(),filename.c_str());
	//klMatrix<TYPE>  c,const char* filename,  const char* title=NULL,const char* xAxis=NULL,const char* yAxis=NULL,const char* zAxis=NULL,bool useExtents=true,bool holdOn=false,const char* marker=NULL
	klHeatMapPlot<double>(mat,arg,title.c_str());
	sprintf(evalString,"print -r1200 -depsc %s;",arg);
	engEvalString(matlabEngine, evalString);
	engEvalString(matlabEngine, "hold off;close(gcf);");
	sprintf(evalString,"epstopdf   %s",arg);
	system(evalString);
	sprintf(arg,"%s.pdf",filename.c_str());
	_tex<<"\\includegraphics[width=10.0cm,height=10.0cm]{"<<arg<<"}"<<endl<<endl;
	delete arg;
	delete evalString;
}
예제 #26
0
파일: LSI_2.cpp 프로젝트: shiyqw/LSI-System
void Search_Synonym_LSI() //这里主要是把相关度前二十的诗词打包成向量投射到词汇空间

{
	double rvpoems[POEMSNUM];
	memset(rvpoems , 0 , sizeof(rvpoems));
	for (int i = 0 ; i <= 19 ; ++i)
	{
		rvpoems[ vresult[i] -1] = 1;
	}
	mxArray * vpoems_m = mxCreateDoubleMatrix(1,POEMSNUM, mxREAL);  
	memcpy(mxGetPr(vpoems_m), rvpoems, POEMSNUM * sizeof(double));  
	engPutVariable(ep, "R",vpoems_m);
	engEvalString(ep, "BS=zeros(7474,1);");
	engEvalString(ep, "AS = R*(V)*(inv(U));");
	engEvalString(ep, "k = 1;");
	for (int i=1 ;i<= WORDNUM ; ++i)
	{
		engEvalString(ep, "BS(k,1) = dot(AS,S(k,:)) / (norm(AS)*norm(S(k,:)) );");
		engEvalString(ep, "k = k + 1;");
	}
	engEvalString(ep, "[LS,INDS]=sort(BS,'descend');" );	
	mxArray * vsyn_m = engGetVariable(ep, "INDS" );
	mxArray * srelation_m = engGetVariable(ep , "LS");
	memcpy(rvsyn, mxGetPr(vsyn_m) , WORDNUM * sizeof(double));
	memcpy(srelation , mxGetPr(srelation_m) , WORDNUM * sizeof(double));
	mxDestroyArray(vpoems_m);
	mxDestroyArray(vsyn_m);
}
예제 #27
0
int object_to_string(void *addr, char *value, int size)
{
	mxArray **my = (mxArray**)addr;
	if (*my!=NULL)
	{
		engPutVariable(engine,"object",*my);
		engEvalString(engine,"disp(char(object));");
		output[strlen(output)-1]='\0'; // trim CR off end
		return _snprintf(value,size-1,"[[%s]]",output);
	}
	else
		return 0;
}
예제 #28
0
static mxArray* matlab_exec(MATLABLINK *matlab, char *format, ...)
{
	char cmd[4096];
	va_list ptr;
	va_start(ptr,format);
	vsprintf(cmd,format,ptr);
	va_end(ptr);
	engEvalString(matlab->engine,cmd);
	if ( matlab->output_buffer && strcmp(matlab->output_buffer,"")!=0 )
		gl_verbose( "%s", matlab->output_buffer );
	mxArray *ans = engGetVariable(matlab->engine,"ans");
	return ans; // failed
}
/* Record data for MATLAB processing. */
void recordData() {
	if (dataCounter<numberOfSamples) {
		dataSetVisual[dataCounter] = shift+((sizeOrig*size)/2);
		print_time();
	} else if (dataCounter == numberOfSamples) {
		engPutVariable(ep, "cppVisualSamples", cppVisualSamples); 
		engPutVariable(ep, "cppTimingSamples", cppTimingSamples);
		engEvalString(ep, "save visualdata.mat");
		fprintf(stderr, "Recording of %i samples is complete.  Data saved as \"visualdata.mat\".\n", numberOfSamples);
		recordMode = 0;
		dataCounter++;
	}
}
예제 #30
0
int main(){
	Engine *ep;
	mxArray *mX = NULL, *mY = NULL;
	const size_t N = 1024;
	
	if (!(ep = engOpen(NULL))) {
		fprintf(stderr, "\nCan't start MATLAB engine\n");
		return EXIT_FAILURE;
	}

	mX = mxCreateDoubleMatrix(1, N, mxREAL);
	mY = mxCreateDoubleMatrix(1, N, mxREAL);

	double *x = (double *)mxGetPr(mX);
	double *y = (double *)mxGetPr(mY);
	for (size_t i=0; i<N; ++i){
		x[i] = 2.0*M_PI*i/N;
		y[i] = sin(x[i]);
	}
	
	engPutVariable(ep, "x", mX);
	engPutVariable(ep, "y", mY);

	engEvalString(ep, "h = figure;");
	engEvalString(ep, "plot(x, y); hold on;");
	engEvalString(ep, "plot([0 0],ylim,'k'); plot(xlim,[0 0],'k');");
	engEvalString(ep, "title('y = sin(x)');");
	engEvalString(ep, "xlabel('x');");
	engEvalString(ep, "ylabel('y');");
	engEvalString(ep, "grid on; grid minor; box off;");
	engEvalString(ep, "print(h,'-dpng','sin.png');");
	engEvalString(ep, "close all;");
	
	mxDestroyArray(mX);
	mxDestroyArray(mY);

	printf("Ok. Press ENTER to exit\n");
	getchar();

	engClose(ep);
	
	return EXIT_SUCCESS;
}