Discussion :

[khalil.ajmi]

bonjour tout le monde ;
en fait mon probleme est le suivant ,j'ai crée une fonction en c qui se nomme calcul.c et maintenant je n'arrive pas a faire appel de cette fonction depuis un code matlab radio.m ,j'ai essayé de voir des tutorial sur mexfile mais j'ai pa su comment l'utiliser .
est ce que quelqu'un peut m'aider ?
prière de me répondre
merci bien a tous.

[Jerome Briot]

Montre nous le code de cette fonction calcul.c

Voici un exemple très simpliste :

Code C :

Code C :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
#include <stdio.h>
 
double somme(double a, double b) {
    return a+b;
}
 
int main(void) {
    double a, b, c;
 
    a = 2.5;
    b = 4.8;
 
    c = somme(a, b);
 
    printf("\nLa somme de %3.1f et %3.1f vaut %3.1f\n", a, b, c);
 
    return 0;
 
}

Equivalent C-MEX :

Code C :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
#include "mex.h"
 
double somme(double *a, double *b) {
    return a[0]+b[0];
}
 
void mexFunction( int nlhs, mxArray *plhs[],
        int nrhs, const mxArray *prhs[]) {
 
    double *a, *b, c;
 
    a = mxGetPr(prhs[0]);
    b = mxGetPr(prhs[1]);
 
    c = somme(a, b);
 
    mexPrintf("\nLa somme de %3.1f et %3.1f vaut %3.1f\n",a[0], b[0], c);
 
}

[khalil.ajmi]

merci bien pour votre réponse
en fait j'ai essaye d'appliquer ca sur ma fonction calcul.c et voila le code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
#include "mex.h"
#include "matrix.h"
#include "fintrf.h"
void calcul (double *x,  double *z,int n)
{
 int i,j,k;
 double temp;
 double *b;
 b=x;
for (i=1;i++;i<=n){
 
z[i]=0;
for (j=1;j++;j<=n){
 
     z[i]=z[i]+x[j]*b[j];
}
 
temp=b[n];
 
    for (k=1;k++;k<=n-1){
    b[n-k+1]=b[n-k];
    }
 
 b[1]=temp;
 }}
 
/* The gateway function */
void mexFunction( int nlhs, mxArray *plhs[],
                  int nrhs, const mxArray *prhs[])
{
/* variable declarations here */
 
     double *inMatrix;       /* 1xN input matrix */
     int ncols;           /* size of matrix */
     double *outMatrix;      /* output matrix */
 
/* code here */
    if(nrhs!=1) {
    mexErrMsgIdAndTxt("MyToolbox:arrayProduct:nrhs",
                      "one inputs required.");
}
if(nlhs!=1) {
    mexErrMsgIdAndTxt("MyToolbox:arrayProduct:nlhs",
                      "one output required.");
}
 
 
/* create a pointer to the real data in the input matrix  */
inMatrix = mxGetPr(prhs[0]);
 
/* get dimensions of the input matrix */
ncols = mxGetN(prhs[0]);
/* create the output matrix */
plhs[0] = mxCreateDoubleMatrix(1,ncols,mxREAL);
 
/* get a pointer to the real data in the output matrix */
outMatrix = mxGetPr(plhs[0]);
calcul(inMatrix,outMatrix,ncols);
 
}

je complique sur matlab avec la commande "mex calcul.c" et puis j 'essaie de le tester sur un vecteur de taille 3 mais ca me donne toujours une erreur .je me demande est ce que j'ai une erreur dans mon code.
merci bien.
prière de me rependre

[Jean Dumoncel]

Bonjour,

quel est le message d'erreur?

[khalil.ajmi]

bonjour
lorsque je met par exemple un vecteur

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

a=[1 2 3]

et je lance sur matlab ma fonction "calcul(a)" ca me donne cette erreur

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
>> calcul(a)
??? one output required.

en fait cette fonction doit retourner un vecteur de même taille que a mais la ca marche pas
merci pour votre réponse

[khalil.ajmi]

j'ai passé des heure a tester mais toujours le même problème.
voila ma fonction calcul.c

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
void calcul (double *x,  double *z,int n)
{
 int i,j,k;
 double temp;
 double *b;
 b=x;
for (i=1;i++;i<=n){
 
z[i]=0;
for (j=1;j++;j<=n){
 
     z[i]=z[i]+x[j]*b[j];
}
 
temp=b[n];
 
    for (k=1;k++;k<=n-1){
    b[n-k+1]=b[n-k];
    }
 
 b[1]=temp;
  }}

est ce que quelqu'un pourrait m'envoyer le code C-mex de cette fonction ?
j'en ai vraiment besoin ,merci d'avance

[Jerome Briot]

As-tu au moins vérifié que ton code C fonctionne ?

Code C :

Sélectionner tout - Visualiser dans une fenêtre à part

b=x;

Je pense que tu souhaites copier le contenu du tableau x dans le tableau b... si c'est le cas, la ligne que tu as codée ne copie pas le contenu mais l'adresse => http://c.developpez.com/faq/?page=ta..._tableau_copie

Code C :

Sélectionner tout - Visualiser dans une fenêtre à part

for (i=1;i++;i<=n){

Tu devrais revoir la syntaxe de la boucle FOR en C => http://melem.developpez.com/tutoriel...sions#LIII-F-6

[khalil.ajmi]

bonjour;
merci bien pour votre réponse ca la fonction marche bien avec des valeurs réelles mais si jamais j'introduis des variable complexe j'obtiens des beugs sur matlab
je crois que l'erreur est dans la déclaration de mes variable ,
comment déclarer une variable complexe sur c?
merci bien .

[Jerome Briot]

Pourquoi utiliser un fichier C-MEX dans ce cas ?
Pourquoi ne pas directement coder ton algorithme sous MATLAB ?

[khalil.ajmi]

bonjour;
en fait j'ai implémenter ca sur matlab mais ca prend énormément de temps pour la simulation c'est pour ca j'ai voulu implémenter en c parce que je vais gagner énormément en temps de simulation mais j'ai pas su comment declarer des variables complexe en c ?