Problèmes multiples de programmation

Version imprimable

05/06/2007, 15h04
VanessaDu67

Problèmes multiples de programmation

Bonjour,
Je bosse sur Matlab.
J'ai un fichier CalculConstante.m qui me calcule des constantes.
Pour plus de clarté, je fais un autre programme qui doit utiliser le calcul des constantes précédemment utilisées.
Comment faire pour appeler le premier fichier CalculConstante.m.................?
Merci d'avance

Configuration: Windows XP
Firefox 2.0.0.4
05/06/2007, 15h12
coolzy

Peux-tu preciser ce qu'il y a dans calculconstante.m? une fonction?
en regle generale tu as une fonction calcule, ecrite dans le fichier calcule.m, et pour l'appeler quelquepart il suffit de citer son nom, calcul, matlab ouvrira automatiquement le fichier du meme nom en '.m' et executra ce qu'il y a dedans.
si tu as une fonction du type par exemple:

Code:

[output 1,output2] = function calcule(arg1,arg2,arg3)

ecrite dans ton fichier calcule.m, alors dans ton autre fichier, pour appeler cette fonction, il te suffit d'ecrire

Code:

[a,b]=calcule(argument1,argument2,argument3)

avec les arguments voulus.
05/06/2007, 15h31
Lili42

Ou simplement (si c'est un script sans fonction):
eval('CalculConstante.m');
qui permet d'exécuter ton script dans un autre...
05/06/2007, 15h35
Jerome Briot

Citation:

Envoyé par Lili42

eval('CalculConstante.m');

NON ! Ne pas utiliser EVAL dans ce cas. Utilisez juste le nom du script sans l'extension.
05/06/2007, 15h36
Lili42

j'utilise tous mes codes avec eval ????!
ça semble marcher sans problème...
05/06/2007, 15h41
Jerome Briot

Un conseil si tu veux continuer à coder proprement avec Matlab, il va falloir penser à se passer rapidement d'EVAL.

A lire : http://www.developpez.net/forums/sho...d.php?t=264608 (voir tous les liens mentionés)
05/06/2007, 15h50
VanessaDu67

Ok merci...
Toujours pour le meme programme me permettant de calculer des constantes de propagation pour une fibre optique. Ce calcul est très long, il peut prendre entre 2 et 3 jours, me "ralentissant" ainsi mon PC lorsque je l'utilise la journée........
C'est pour cela que j'aimerai faire tourner ce programme uniquement la nuit quand je pars du boulot, et des le matin quand j'arrive, j'aimerais l'interrompre pour pouvoir le faire reprendre le soir quand je m'en vais.....
Avez vous une idée de la fonction Matlab à utiliser pour vouvoir faire ce genre de chose...?
05/06/2007, 15h53
Jerome Briot

La taille de tes données justifie-t-elle la durée du calcul ?

Peut-on avoir un aperçu de ton programme ?
05/06/2007, 16h04
VanessaDu67

En fait ce qui est long, c'est de rechercher les zéros des constantes de propagation.
En fait j'ai deux programmes.
Un premier programme qui me fait un calcul approximatif de ces constantes. La limite de validité de ce programme est de 80%. Les 20% restant j'utilise l'autre programme qui me calcule exactement mes denrière constantes de propagation.
Donc par exemple, si ma fibre contient 1000 modes, je vais lancer le premier programme jsuque les 800 premiers modes, et les 200 modes qui resteront seront calculés par le deuxième programme,et là ca prend minimum 6 jours....
05/06/2007, 16h48
sango85

Quand Dutmatlab te demande :

Citation:

Envoyé par Dutmatlab

Peut-on avoir un aperçu de ton programme ?

Il a en tête non pas une remise en question de l'algorithme, mais une optimisation de code...
Juste pour info, Je devais relancer une routine Matlab sur 8700 fichiers... J'étais heureux lorsque j'ai réussi à réduire le temp de calcul sur chaque fichier de 40 à 2 secondes... Juste en gérant différemment mes variables... :king:

Tu es sûre de l'optimitation de ton code ? :?

C mon patorn qui a fait le code mais celui ci est le suivant :

Code:

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
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
%------------------ Constants and initialisation ------------------%
 
warning off MATLAB:divideByZero;
% warning off MATLAB:fzero:UndeterminedSyntax;
options = optimset('TolFun',1e-12);
Lamb_mean = 630e-9;
 
%************ Calcul de la frequence, du nombre de mode et du groupe de mode******************%
 
V = (2 * pi / Lamb_mean) * (d_co / 2) * sqrt(n_co^2 - n_cl^2)
N_modes = 0.5 * V^2
M_modes = sqrt(N_modes)
 
k0 = k0_mean;
 
%******************************************************************%
 
 
%************ Sauvegarde du fichier FibreData******************%
repertoire = strcat ('..\Results\',date,'\FibreData_', num2str(V),'.mat');
repertoire2 = strcat ('..\Results\',date,'\FibreData_', num2str(V),'.txt');
save  (repertoire, 'd_co', 'd_cl', 'd_ct', 'ON', 'n_co', 'n_cl', 'Lamb_mean', 'k0_mean');
save (repertoire2, 'd_co', 'd_cl', 'd_ct', 'ON', 'n_co', 'n_cl', 'Lamb_mean', 'k0_mean', '-ASCII', '-DOUBLE','-TABS'); 
%******************************************************************%
 
%------------------ Search for cut off normalised frequencies -----%
 
dBeta = (n_co - n_cl) * k0 / (20 * floor(N_modes * 1.2));
Interv = [];
Beta_array = [];
save(['Intervals_',num2str(V),'.mat'],'Interv'); 
 
BetaMin = k0 * n_cl
BetaMax = k0 * n_co
BetaLength = BetaMax - BetaMin;
nu = 0;
ZeroFound = 1;
Beta = BetaMin;
h0 = waitbar(0,['% of orders processed']);
BarPosit = get(h0,'Position');
while ZeroFound == 1
    i = 1;
    ZeroFound = 0;
    waitbar(nu/M_modes,h0);
    h1 = waitbar(0,['Order ',num2str(nu),' : search of root segments in progress'],'Position',(BarPosit + [100 -100 0 0]));
    while Beta <= BetaMax
        waitbar((Beta-BetaMin)/BetaLength,h1);
        Disp1 = gDispersion(Beta,nu,k0,n_co,n_cl,d_co,J_bess,K_bess);
        Disp2 = gDispersion(Beta + dBeta,nu,k0,n_co,n_cl,d_co,J_bess,K_bess);
        Disp3 = gDispersion((Beta + 2 * dBeta),nu,k0,n_co,n_cl,d_co,J_bess,K_bess);
        if (isfinite(Disp1) & isfinite(Disp2) & isfinite(Disp3)) & (sign(Disp2) ~= sign(Disp1)) & (sign(Disp3 - Disp2) == sign(Disp2 - Disp1))
            Interv((nu+1),i,1) = Beta;
            Interv((nu+1),i,2) = Beta + dBeta;
            i = i + 1;
            ZeroFound = 1;
        end
        Beta = Beta + dBeta;
    end
    close(h1);
    Beta = BetaMin;
    if ZeroFound == 1
        save(['Intervals_',num2str(V),'.mat'],'Interv','-APPEND');
        nu = nu + 1;
    end
end
close(h0);
 
nu
 
h = waitbar(0,'Roots search in progress');
Betas = [];
Size2 = length(Interv(1,:,1))
for i = 1:nu
    waitbar(i/(nu + 1));
    for j = 1:nnz(Interv(i,:,1))
        Betas(i,(Size2 - j + 1)) = fzero(@gDispersion,[Interv(i,j,1) Interv(i,j,2)],options,(i-1),k0,n_co,n_cl,d_co,J_bess,K_bess);
    end
end
close(h);
 
Betas

[Matlab] programme de mise en attente

Toujours pour le meme programme me permettant de calculer des constantes de propagation pour une fibre optique. Ce calcul est très long, il peut prendre entre 2 et 3 jours, me "ralentissant" ainsi mon PC lorsque je l'utilise la journée........
C'est pour cela que j'aimerai faire tourner ce programme uniquement la nuit quand je pars du boulot, et des le matin quand j'arrive, j'aimerais l'interrompre pour pouvoir le faire reprendre le soir quand je m'en vais.....Ca serait comme une mise en attente..
Avez vous une idée de la fonction Matlab à utiliser pour vouvoir faire ce genre de chose...?
Hier on m'a demandé de donner le code de mon programme pour voir si une optimisation pouvait se faire....mais aucune réponse ne m'a été donné et là cela devient très urgent que j'ai cette réponse car il faut que je génère mes tables de matrices des constantes assez rapidement...

Le code est le suivant :

Code:

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
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
%------------------ Constants and initialisation ------------------%
 
warning off MATLAB:divideByZero;
% warning off MATLAB:fzero:UndeterminedSyntax;
options = optimset('TolFun',1e-12);
Lamb_mean = 630e-9;
 
%************ Calcul de la frequence, du nombre de mode et du groupe de mode******************%
 
V = (2 * pi / Lamb_mean) * (d_co / 2) * sqrt(n_co^2 - n_cl^2)
N_modes = 0.5 * V^2
M_modes = sqrt(N_modes)
 
k0 = k0_mean;
 
%******************************************************************%
 
 
%************ Sauvegarde du fichier FibreData******************%
repertoire = strcat ('..\Results\',date,'\FibreData_', num2str(V),'.mat');
repertoire2 = strcat ('..\Results\',date,'\FibreData_', num2str(V),'.txt');
save  (repertoire, 'd_co', 'd_cl', 'd_ct', 'ON', 'n_co', 'n_cl', 'Lamb_mean', 'k0_mean');
save (repertoire2, 'd_co', 'd_cl', 'd_ct', 'ON', 'n_co', 'n_cl', 'Lamb_mean', 'k0_mean', '-ASCII', '-DOUBLE','-TABS'); 
%******************************************************************%
 
%------------------ Search for cut off normalised frequencies -----%
 
dBeta = (n_co - n_cl) * k0 / (20 * floor(N_modes * 1.2));
Interv = [];
Beta_array = [];
save(['Intervals_',num2str(V),'.mat'],'Interv'); 
 
BetaMin = k0 * n_cl
BetaMax = k0 * n_co
BetaLength = BetaMax - BetaMin;
nu = 0;
ZeroFound = 1;
Beta = BetaMin;
h0 = waitbar(0,['% of orders processed']);
BarPosit = get(h0,'Position');
while ZeroFound == 1
    i = 1;
    ZeroFound = 0;
    waitbar(nu/M_modes,h0);
    h1 = waitbar(0,['Order ',num2str(nu),' : search of root segments in progress'],'Position',(BarPosit + [100 -100 0 0]));
    while Beta <= BetaMax
        waitbar((Beta-BetaMin)/BetaLength,h1);
        Disp1 = gDispersion(Beta,nu,k0,n_co,n_cl,d_co,J_bess,K_bess);
        Disp2 = gDispersion(Beta + dBeta,nu,k0,n_co,n_cl,d_co,J_bess,K_bess);
        Disp3 = gDispersion((Beta + 2 * dBeta),nu,k0,n_co,n_cl,d_co,J_bess,K_bess);
        if (isfinite(Disp1) & isfinite(Disp2) & isfinite(Disp3)) & (sign(Disp2) ~= sign(Disp1)) & (sign(Disp3 - Disp2) == sign(Disp2 - Disp1))
            Interv((nu+1),i,1) = Beta;
            Interv((nu+1),i,2) = Beta + dBeta;
            i = i + 1;
            ZeroFound = 1;
        end
        Beta = Beta + dBeta;
    end
    close(h1);
    Beta = BetaMin;
    if ZeroFound == 1
        save(['Intervals_',num2str(V),'.mat'],'Interv','-APPEND');
        nu = nu + 1;
    end
end
close(h0);
 
nu
 
h = waitbar(0,'Roots search in progress');
Betas = [];
Size2 = length(Interv(1,:,1))
for i = 1:nu
    waitbar(i/(nu + 1));
    for j = 1:nnz(Interv(i,:,1))
        Betas(i,(Size2 - j + 1)) = fzero(@gDispersion,[Interv(i,j,1) Interv(i,j,2)],options,(i-1),k0,n_co,n_cl,d_co,J_bess,K_bess);
    end
end
close(h);
 
Betas
 
%************ Sauvegarde des valeurs de Beta***********************%
beta = strcat ('..\Results\',date,'\Betas_', num2str(V),'.mat');
beta2 = strcat ('..\Results\',date,'\Betas_', num2str(V),'.txt');
save  (beta, 'Betas');
save (beta2, 'Betas','-ASCII', '-DOUBLE','-TABS'); 
%******************************************************************%

En vous disant un grand merci..................

Vaness

06/06/2007, 12h18
Jerome Briot
Ce que tu cherches à faire est quasiment impossible.

Pour optimiser ton code :
- utilise la préallocation de mémoire
- évite d'utiliser la fonction WAITBAR
Peux-tu nous montrer la fonction gDispersion (avec l'autorisation de ton responsable bien sur) ?

La fonction gDispersion est la suivante :
Code:

1 2 3 4 5 6 7 8 function dispersion = gDispersion(beta_nu,nu,k0,n_co,n_cl,d_co,J_bess,K_bess) kappa_nu = sqrt((n_co * k0).^2 - beta_nu.^2); gamma_nu = sqrt(beta_nu.^2 - (n_cl * k0).^2); dispersion = kappa_nu * (d_co / 2) .* (feval(J_bess,(kappa_nu * (d_co / 2)),(nu+1)) ./ feval(J_bess,(kappa_nu * (d_co / 2)),nu))... - gamma_nu * (d_co / 2) .* (feval(K_bess,(gamma_nu * (d_co / 2)),(nu+1)) ./ feval(K_bess,(gamma_nu * (d_co / 2)),nu));
Je ne connaissais pas la fonction WaitBar avant d'utiliser ce programme. Cette question va vous paraitre bete mais pourquoi est-il deconsseillé d'utiliser cette fonction?

06/06/2007, 13h14
Jerome Briot

Je pense que c'est la fonction FEVAL qui ralentit le code. Il faudrait essayer de s'en passer... mais ça peut être assez difficile.

La fonction waitbar sert juste à afficher une barre d'avancement lors du calcul et elle n'est pas très efficace.