源程序清单

2) 再选择菜单：Project->Settings 在对话框中选择 Link，在 Object\Library Modules 中添加：libmx.lib,libmat.lib,libeng.lib 三项。此操作每次打开编 译器都要修改一次。
engEvalString(ep,"r=a*t;"); engEvalString(ep,"FA=compose(FA,r);"); engEvalString(ep, " syms t;"); engEvalString(ep,"u=u+diff(FA,t)"); } cout<<"please input du/dt(x,y,z,0)"<<endl; 输 入 du/dt(x,y,z,0) 的 解 析 表 达式，处理过程与处理 u(x,y,z,0)类似。 int legs=SaveList(s); // if(legs){ double sai[40]; b=s; for(i=0;i<legs;i++){ sai[i]=b->data;cout<<i<<"\t"<<sai[i]<<(char)sai[i]<<endl; b=b->next; } Delete(s); SAI=mxCreateDoubleMatrix(1,legs,mxREAL); memcpy((char*)mxGetPr(SAI),(char*)sai,legs*sizeof(double)); engPutVariable(ep,"SAI",SAI); engEvalString(ep,"SA=trilinear(SAI);"); engEvalString(ep,"r=a*t;"); engEvalString(ep,"S=compose(SA,r);"); engEvalString(ep, " u=u+S;"); } cout<<"please input f(x,y,z,t)"<<endl;输入非齐次项 int lege=SaveList(e); / 将非齐次项 f(x,y,z)按照上面类似方式处理存入 MATLAB 计算引擎。 if(lege){ double ef[40]; b=e; for(i=0;i<lege;i++){ ef[i]=b->data; b=b->next; } Delete(e); EF=mxCreateDoubleMatrix(1,lege,mxREAL);
1)
在 开 始 编 译 之 前 ， 要 先 选 择 菜 单 ： Toos->Options 在 对 话 框 中 选 择 Directories,在 include files 项中添加：C:\...\matlab\extern\include(请读 者 自 行 修 改 为 自 己 的 MATLAB 安 装 路 径 ) ， 再 在 library files 中 添 加 C:\...\matlab\extern\lib\win32\microsoft\msvc60(请读者自行修改为自己 的 MATLAB 安装路径),本操作对每个项目修改一次即可。
2、
MATLAB 下的 M 文件部分：
说明：这里由于时间原因，我们没能给出三维非齐次波动方程的求 解 M 文件。
function Q=trilinear(T) ascii_T=double(T); p=char(ascii_T);//将数组 T的 ASCII标准码转换成字符形式的函数解析 表达式。 syms r m n x y z a t X Y Z;//设置符号变量 X=x+r*sin(m)*cos(n);//球坐标代换。 Y=y+#include<iostream> using namespace std; struct List{//创建链表数据结构。 double data; struct List *next; }; typedef List L,*LI; extern void Init(LI &l){//新建链表。 l=NULL; } extern void Save(LI &head,L *p){// 将 输 入 的 字 符 所 转 化 的 ASCII 标 准 代码存入链表。 l=head; if(!head){ head=p; } else{ while(l->next)l=l->next; l->next=p; } } extern void Print(LI head){//输出链表以便测试检查。 L *l;l=head; while(l){ cout<<"l="<<l->data<<"address"<<l<<"\t"<<l->next<<endl; l=l->next; } } extern void Delete(L *head){//销毁链表以保护内存。 L *l; l=head; while(head){ head=head->next; delete l;
memcpy((char*)mxGetPr(EF),(char*)ef,lege*sizeof(double)); engPutVariable(ep,"EF",EF); engEvalString(ep,"EF=triunlinear(EF);"); 与前两个条件处理方式不 同的是这一次调用 triunlinear 的 M 文件求解。但是由于时间愿意，这 个 M 文件还没有完成。 engEvalString(ep, " u=u+EF;"); 将 计 算 所 得 结 果 与 齐 次 方 程 的 解 相加。 } if(legf)mxDestroyArray(FAI);释放所有内存空间。 if(legs)mxDestroyArray(SAI); if(lege)mxDestroyArray(EF); int ch; cout<<"close engine?"<<endl; cin>>ch; if(ch) engClose(ep); //关闭 MATLAB 引擎。 }
int i=0; LI f,s,e,b; Engine *ep; //打开 MATLAB 计算引擎，赋指针为 EP mxArray *FAI=NULL,*SAI=NULL,*EF=NULL; if (!(ep=engOpen(NULL))) //测试引擎是否打开。 { cout <<"Can't start Matlab engine!" <<endl; exit(1); } engEvalString(ep, " syms X Y Z x y z m n a t;"); engEvalString(ep, " u='0';"); cout<<"please input the type of the qutation:"<<endl;// 输 入 方 程 的 初 值 条 件 u(x,y,z,0)，以 ’0’字符结束。 cout<<"please input u(x,y,z,0)"<<endl; int legf=SaveList(f); //生成解析表达式的字符链表并返回字符个数。 if(legf){//将解析表达式的 ASCII 代码存入数组 fai[40] b=f; double fai[40]; for(i=0;i<legf;i++){ fai[i]=b->data; b=b->next; } Delete(f);清除链表。 FAI=mxCreateDoubleMatrix(1,legf,mxREAL); 在 MATLAB 引 擎 里 创建一个一维矩阵，长度为所输入的解析表达式的字符个数 memcpy((char*)mxGetPr(FAI),(char*)fai,legf*sizeof(double)); ASCII 标准代码数组复制进入 MATLAB 引擎。 engPutVariable(ep,"FAI",FAI);为 MATLAB 引擎赋予变量。 engEvalString(ep,"FA=trilinear(FAI);"); 调 用 trilinear 文 件 将 ASCII 标准代码转换为字符的解析函数表达式并积分求解。 将
l=head; } delete head; } extern int SaveList(LI &a)// 输 入 解 析 表 达 式 ， 并 将 其 以 字 符 的 ASCII 标准码形式存入链表。 { Init(a); char le; int length=0; cout<<"please input"<<endl; cin>>le; while(le!='0'){ L *temp=new L; temp->data=(double)le; temp->next=NULL; Save(a,temp); length++; cin>>le; } return length; } 文件 3：三维波动方程求解计算核心模块 #include<iostream> #include"engine.h" #include<string.h> #include<cmath> using namespace std; struct List{ double data; struct List *next; }; typedef List L,*LI; extern int SaveList(LI &a);/对调用的外部函数进行声明。 extern void Delete(L *head); extern void Print(LI head); extern void Save(LI &head,L *p); extern void Print(LI head); extern void TRI(){//三维方程求解计算核心。