Homec4science
Diffusion yearII_reports (master) Browse

rSTICAZZI/physnum/rap5/Convergence.m
5fb2fcbcde85master

Last Change
Raw File

Convergence.m
View Options

% Ce script Matlab automatise la production de resultats
% lorsqu'on doit faire une serie de simulations en
% variant un des parametres d'entree.
%
% Il utilise les arguments du programme (voir ConfigFile.h)
% pour remplacer la valeur d'un parametre du fichier d'input
% par la valeur scannee.
%
%% Parametres %%
%%%%%%%%%%%%%%%%
repertoire = './'; % Chemin d'acces au code compile (NB: enlever le ./ sous Windows)
executable = 'Exercice5'; % Nom de l'executable (NB: ajouter .exe sous Windows)
input = 'configuration.in'; % Nom du fichier d'entree de base
nsimul = 150; % Nombre de simulations a faire
dt = logspace(-5,-3, nsimul);
N = 40;
%outfile = sprintf('outs/conv%d.dat', N);
outfile = 'outs/conv.dat';
paramstr = 'dt'; % Nom du parametre a scanner
param = dt; % Valeurs du parametre a scanner
point = [0.625; 0.425];
op = [0; 0]; % output point, ex1
op(1) = point(1) * N;
op(2) = point(2) * N;
outdir = 'm_outs/';
%% Simulations %%
%%%%%%%%%%%%%%%%%
output = cell(1, nsimul); % Tableau de cellules contenant le nom des fichiers de sortie
for i = 1:nsimul
output{i} = [outdir, paramstr, '=', num2str(param(i))];
% Execution du programme en lui envoyant la valeur a scanner en argument
cmd = sprintf('%s%s %s %s=%.15g N=%d \"%s=(%d, %d)\" output=%s', repertoire, executable, input, paramstr, param(i), N, 'outpoint', op(1), op(2), output{i});
disp(cmd)
system(cmd);
end
%% Analyse %%
%%%%%%%%%%%%%
% Parcours des resultats de toutes les simulations
Tp = zeros(1,nsimul);
for i = 1:nsimul % Parcours des resultats de toutes les simulations
data = load([output{i} '_T.out']);
s = size(data);
Tp(i) = data(s(1), 2);
end
% get minimal value
first = Tp(1);
for i = 1:nsimul
Tp(i) = abs(Tp(i) - first);
end
%% Linear regression %%
%%%%%%%%%%%%%%%%%%%%%%%
m = polyfit(dt, Tp, 1);
linreg = sprintf('%.8g * x + %.8g', m(1), m(2));
disp(linreg);
%% File writing %%
%%%%%%%%%%%%%
% print data on file
out = fopen(outfile, 'w');
fprintf(out, 'dt T\n');
for i = 1:nsimul
fprintf(out, '%.15g %.15g\n', dt(i), Tp(i));
end
fclose(out);
%% Figures %%
%%%%%%%%%%%%%
% make figure with data
figure
loglog(dt,Tp,'k+')
xlabel('\Delta t [s]')
ylabel(sprintf('T(%0.2f,%0.2f) [°C]',point(1),point(2)))
grid on
c4science · Help