1 function [s v d]=mysvd(a)
2     [s v d]=svd(a);
3 end

 1 #include <iostream>
 2 #include "mclmcr.h"
 3 #include "matrix.h"
 4 #include "mclcppclass.h"
 5 #include "libmysvd.h"
 6 using namespace std;
 7 int main()
 8 {
 9     if(!libmysvdInitialize())    //初始化，在libmysvd.h中
10     {
11         cout << "Could not initialize libmysvd!" <<endl;
12         exit(0);
13     }
14     double data[9]={1,2,3,4,5,6,7,8,9};
15     mwArray a(3,3,mxDOUBLE_CLASS,mxREAL);
16 
17     a.SetData(data,9);
18     cout<<a<<endl<<endl;    //显然，mwArry类对<<是有运算符重载的
19 
20     mwArray s(3,3,mxDOUBLE_CLASS,mxREAL);
21     mwArray v(3,3,mxDOUBLE_CLASS,mxREAL);
22     mwArray d(3,3,mxDOUBLE_CLASS,mxREAL);
23 
24     mysvd(3,s,v,d,a);        //调用自己用matlab写的函数，对a做奇异值分解
25 
26     cout<<s<<endl<<endl;    
27     cout<<v<<endl<<endl;
28     
29     for (int i=1;i<4;i++)    //逐个显示元素
30     {
31         for (int j=1;j<4;j++)
32         {
33             cout<<d.Get(i,j)<<"  ";
34         }
35         cout<<endl;
36     }
37 
38     cin.get();
39     return 0;
40 }