matrixlib.cpp

MATLAB混合编程想把matlab强大的功能用在自己的工程实践中吗这是本不错的混合编程的书

// MatrixLIB.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#include <matlib.h>

int main(int argc, char* argv[])
{
	initM(MATCOM_VERSION);     //初始化类库
	
	Mm a;  //初始化矩阵变量
	Mm b;
	Mm x;
	a=zeros(3,3);    //给矩阵变量赋初值　　　　
	a.r(1,1)=3;  a.r(1,2)=1;  a.r(1,3)=7;
	a.r(2,1)=11; a.r(2,2)=4;  a.r(2,3)=6;
	a.r(3,1)=9;  a.r(3,2)=2;  a.r(3,3)=0;
	b=(BR(7),semi, 3, semi, -20);
	
	display(a); 
	printf("\n");
	display(b);	
	
	x=mldivide(a,b);   //求解线性系统方程
	//	display(x);
	//	printf("\n");
	
	for (int i=1;i<=x.rows();i++)    //显示解向量
	{
		for (int j=1;j<=x.cols();j++)
		{
			printf("x(%d,%d) = %14.6g",i,j,x.r(i,j));
		}
		printf("\n");
	}
	
	double *a_v=(double* ) malloc(5*sizeof(double));
	Mm a2=(BR(1),1, 3, 5.4, 0, -1);   //假设a是一个512×1的矩阵
	memcpy(a_v, a2.addr(), 512*sizeof(double)); 
	
	//以下是注意dMm和Mm定义变量的细微差别
	dMm (y); 
	y=1.5;
	display(y);
	
	Mm (y2);
	y2=1.5;
	display(y2);
	
	//将字符串转换为字符串矩阵
	Mm x2;  
	x2=TM("This sentence is false");
	x2.setstr(1);
	disp(x2);
	
	//库常量用法
	Mm a3=(BR(1), 0, -1, semi, i, pi, 5, semi, 0.3, -4, 0);
	display(a3);
	int *w=a3.getdims();
    disp(*w);
	Mm b2=a3(c_p,2);;
	display(b2);
	
	//绘图
	Mm t2; t2=linspace(0, 2*pi);
	plot((CL(t2), sin(t2), TM("*")));
	
	exitM(); //退出Matrix<LIB>类库
	
//	return 0;
	
	return 0;
}