Discussion :

[clem3900]

Bonjour,

actuellement j'ai un programme que je lis dans un script comme ceci:

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
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
clear all;
clc;
close all;
 
file1='20120731-full-350nm-slit100ShTi0.1Gain5.asc';
S1=load(file1);
file2='20120731-full-450nm-slit100ShTi0.1Gain5.asc';
S2=load(file2);
file3='20120731-full-550nm-slit100ShTi0.1Gain5.asc';
S3=load(file3);
file4='20120731-full-650nm-slit100ShTi0.1Gain5.asc';
S4=load(file4);
file5='20120731-full-750nm-slit100ShTi0.1Gain5.asc';
S5=load(file5);
file6='20120731-full-850nm-slit100ShTi0.1Gain5.asc';
S6=load(file6);
 
l1=size(S1,1);
l2=size(S2,1);
l3=size(S3,1);
l4=size(S4,1);
l5=size(S5,1);
l6=size(S6,1);
 
c1=size(S1,2);
c2=size(S2,2);
c3=size(S3,2);
c4=size(S4,2);
c5=size(S5,2);
c6=size(S6,2);
 
Q1= S1(129:l1-133,1);
Q2= S2(129:l2-133,1);
Q3= S3(129:l3-133,1);
Q4= S4(129:l4-133,1);
Q5= S5(129:l5-133,1);
Q6= S6(129:l6-133,1);
 
 
R1 = sum(S1(129:l1-133,150:c1-150),2);
R2 = sum(S2(129:l2-133,150:c2-150),2);
R3 = sum(S3(129:l3-133,150:c3-150),2);
R4 = sum(S4(129:l4-133,150:c4-150),2);
R5 = sum(S5(129:l5-133,150:c5-150),2);
R6 = sum(S6(129:l6-133,150:c6-150),2);
 
 
% figure(1);
 
% subplot(2,3,1);plot(Q1,R1);title('350nm');
% subplot(2,3,2);plot(Q2,R2);title('450nm');
% subplot(2,3,3);plot(Q3,R3);title('550nm');
% subplot(2,3,4);plot(Q4,R4);title('650nm');
% subplot(2,3,5);plot(Q5,R5);title('750nm');
% subplot(2,3,6);plot(Q6,R6);title('850nm');
% 
% saveName = ('Lines');  
% saveas(gcf, saveName, 'jpg');      
% saveas(gcf, saveName, 'fig');
% 
% figure(2);
% 
% plot(Q1,R1,'b');
% hold on
% plot(Q2,R2,'r');
% hold on
% plot(Q3,R3,'y');
% hold on
% plot(Q4,R4,'g');
% hold on
% plot(Q5,R5,'c');
% hold on
% plot(Q6,R6,'k');
% title('Full spectra at 100% - at the top (21,5mm from the bottom)');
% xlabel('Wavelength in nm');
% 
% hold off
% hold off
% hold off
% hold off
% hold off
% 
% saveName = ('Full spectra at 100% - at the top (21,5mm from the bottom)');  
% saveas(gcf, saveName, 'jpg');      
% saveas(gcf, saveName, 'fig');

Chacun de ces fichiers étant des tableaux de 1024x1024, j'effectue des opérations dessus et je trace les résultats.
Le problème c'est qu'à chaque fois je dois modifier mon script en copiant/collant les noms des fichiers .asc, parce que j'exécute le même programme avec beaucoup d'autres fichiers.

J'aurais donc voulu savoir s'il y a un moyen de créer une fonction qui ferait les même calculs, mais au lieu de copier/coller les noms des fichiers dans le script, avoir 6 sorties dans ma command window où je rentrerais les noms de mes fichiers et les calculs s’exécuteraient dans la fonction. Ou ce qui serait vraiment parfait c'est que je puisse directement coller les noms des 6 fichiers en un seul bloc dans la commande window et enregistrer chaque fichier dans un tableau.

J'espère avoir été assez clair, et merci pour l'aide

[Invité]

Bonjour,

Il suffit de faire une boucle :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
files = { '20120731-full-350nm-slit100ShTi0.1Gain5.asc'
    '20120731-full-450nm-slit100ShTi0.1Gain5.asc'
    '20120731-full-550nm-slit100ShTi0.1Gain5.asc'
    '20120731-full-650nm-slit100ShTi0.1Gain5.asc'
    '20120731-full-750nm-slit100ShTi0.1Gain5.asc'
    '20120731-full-850nm-slit100ShTi0.1Gain5.asc'};
 
for i = 1:length(files)
    S = load(files{i});
    % ou si besoin de stocker :     S{i} = load(files{i});
 
    % ...
 
end

[clem3900]

Merci beaucoup, tu m'as fait économiser une bonne trentaine de lignes

J'ai une autre question du coup...
Au lieu d'écrire:

Code :

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

files = {%copier/coller tous les noms}

Est-ce que c'est possible de créer un dossier avec uniquement les 6 fichiers à lires, se placer dans ce dossier, et demander à Matlab de mettre lui même les fichiers dans un tableau sans qu'il y ait besoin d'écrire les noms des fichiers soi-même?

Merci

[le fab]

salut

Code :

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

1
2
3
4
rep = 'chemin du repertoire';
files = dir(fullfile(rep,'*.asc'));
files = {files.name};

edit : y avait une erreur

[clem3900]

Nickel, merci beaucoup