（有偿）有人做过matlab c混编吗，matlab调用c程序，c传出三维矩阵

猪猪虾

下面是C的程序，在matlab里面调用的时候，也不报错，就是直接闪退，我查了网上说可能是数组溢出了，断点调试了，但是也没找到是哪有问题

1. #include "mex.h"
   
2. #include "math.h"
   
3. void Recon(double *dims, double *m, double *Pro, double *P)
   
4. {
   
5.         int i, j, k;
   
6.         int u, v  ;
   
7.         double pu, pv, lamda, x, y, z, dx, dy, dz;
   
8.         double inter_x, inter_y ,sin_w,cos_w ,s;
   
9.         double Nx, Ny, Nz, nu , nv ,light_to_center;
   
10. 
    
11.         Nx = dims[0];
    
12.         Ny = dims[1];
    
13.         Nz = dims[2];
    
14.         dx = dims[3];
    
15.         dy = dims[4];
    
16.         dz = dims[5];
    
17.         nu = dims[6];
    
18.         nv = dims[7];
    
19.     sin_w = dims[8];
    
20.     cos_w = dims[9];
    
21.     light_to_center = dims[10];
    
22.    
    
23.         for (k = 1; k < Nz+1; k++)
    
24.         {
    
25.                 for (i = 1; i < Nx+1; i++)
    
26.                 {                       
    
27.                         for (j = 1; j < Ny+1; j++)
    
28.                         {
    
29.                                 x = (i-0.5)*dx - dx*(Nx / 2 + 0.5);                       
    
30.                                 y = (j-0.5)*dy - dy*(Ny / 2 + 0.5);
    
31.                                 z = (k-0.5)*dz - dz*(Nz / 2 + 0.5);
    
32.                                
    
33.                                 //加权参数
    
34.                                 //s = - x * sin(deg2rad(beta)) + y * cos(deg2rad(beta));
    
35.                 s = -x * sin_w + y * cos_w;
    
36.                                 pu = m[0] * x + m[1] * y + m[2] * z + m[3];
    
37.                                 pv = m[4] * x + m[5] * y + m[6] * z + m[7];
    
38.                                 lamda = m[8] * x + m[9] * y + m[10] * z +1;
    
39.                                
    
40.                                 inter_x = pu /lamda ;
    
41.                                 inter_y = pv /lamda ;
    
42.                                
    
43.                                 u = int(inter_x); // exchange u,v
    
44.                                 v = int(inter_y);
    
45.                                
    
46.                                 if ((inter_x >= 1 && inter_x <= (nu-1 )) && (inter_y >= 1 && inter_y <= (nv-1)))
    
47.                                 {
    
48.                                         //在投影图中插值获得当前点的像素
    
49.                                         P[int(j-1 + Ny*((i-1) + Nx*(k-1)))] = Pro[int( inter_x - 1 + nv* (inter_y - 1))] * (light_to_center * light_to_center)/ ((light_to_center - s )* (light_to_center - s ));                                       
    
50.                                 }
    
51.                                 else
    
52.                                 {
    
53.                                          P[int(j-1 + Ny*((i-1) + Nx*(k-1)))]  = 0;
    
54.                                 }
    
55. 
    
56. //                                if ((u >= 1 && u <= (Ni-1 )) && (v >= 1 && v <= (Nj-1)))
    
57. //                                {
    
58. //                                P[u - 1 + Nj* (v - 1)] = P[u - 1 + Nj* (v - 1)] + invoxel[j-1 + Ny*((i-1) + Nx*(k-1))] * (u - inter_x + 1)*(v - inter_y + 1);
    
59. //                                P[u - 1 + Nj* v]       = P[u - 1 + Nj* v]       + invoxel[j-1 + Ny*((i-1) + Nx*(k-1))] * (u - inter_x + 1)*(inter_y - v);
    
60. //                                P[u + Nj* (v - 1)]     = P[u + Nj* (v - 1)]     + invoxel[j-1 + Ny*((i-1) + Nx*(k-1))] * (inter_x - u)*(v - inter_y + 1);
    
61. //                                P[u + Nj* v]           = P[u + Nj* v]           + invoxel[j-1 + Ny*((i-1) + Nx*(k-1))] * (inter_x - u)*(inter_y - v);
    
62. //                                }
    
63.                         }
    
64.                 }
    
65.         }
    
66. 
    
67. }
    
68. 
    
69. 
    
70. //接口
    
71. //nlhs 为输出参数的数目，
    
72. //plhs是输出参数数组，指针数组
    
73. //nrhs是为输入参数的数目
    
74. //prhs是输入参数数组吗，是一个指针数组
    
75. void mexFunction(int nlhs, mxArray *plhs[],int nrhs, const mxArray *prhs[])
    
76. {
    
77.         double *dims;
    
78.         double *m, *Pro;
    
79.         double *P;
    
80.         if (nrhs != 3)
    
81.                 mexErrMsgTxt("3 inputs required.");
    
82.         if (nlhs != 1)
    
83.                 mexErrMsgTxt("One output required.");
    
84. 
    
85.         dims = mxGetPr(prhs[0]);
    
86.         m = mxGetPr(prhs[1]);
    
87.         Pro = mxGetPr(prhs[2]);
    
88. 
    
89.         /*  set the output pointer to the output matrix */
    
90. //         plhs[0] = mxCreateDoubleMatrix(dims[3],dims[4], mxREAL);
    
91. 
    
92.     const mwSize ds[3] = {*dims,*(dims +1),*(dims +2)};
    
93.     plhs[0] = mxCreateNumericArray(3,ds,mxSINGLE_CLASS,mxREAL);
    
94.         /*  create a C pointer to a copy of the output matrix */
    
95.         P = mxGetPr(plhs[0]);   //获得矩阵第一个元素的指针
    
96.         /*  call the C subroutine */
    
97.         Recon(dims, m, Pro, P);
    
98. }
    
99. 
    
100. //plhs：指向输出参数的指针
     
101. //nrhs：输入参数数目

复制代码