From c0f1862d5172c67d0f4c4735c97b82591366b5f1 Mon Sep 17 00:00:00 2001 From: Jean-Denis Vauguet Date: Mon, 1 Jun 2009 01:33:12 +0200 Subject: [PATCH] plop --- .../MATLAB/vectorial/data_processing.m | 11 ++++-- .../IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m | 12 ++++++ .../IPGP Saint Maur/MATLAB/vectorial/wind_effect.m | 46 ++++++++++++---------- 3 files changed, 44 insertions(+), 25 deletions(-) diff --git a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/data_processing.m b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/data_processing.m index 059da1b..dfcc379 100644 --- a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/data_processing.m +++ b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/data_processing.m @@ -14,7 +14,7 @@ data = load(thefile, '-ascii'); binf = 1; res = 1; -bsup = 1; +bsup = size(data,1); if optargin > 0 % real part, anyone? @@ -36,15 +36,15 @@ if optargin > 0 end if ~isempty(varargin{2}) - binf = varargin{2}; + binf = varargin{2} end if ~isempty(varargin{3}) - res = varargin{3}; + res = varargin{3} end if ~isempty(varargin{4}) - bsup = varargin{4}; + bsup = varargin{4} end else % on demand… @@ -54,6 +54,9 @@ else end end +disp '---' +size(w) +disp '---' w = w(binf:res:bsup); u = u(binf:res:bsup); p = p(binf:res:bsup); diff --git a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m index f88089e..8b4e499 100644 --- a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m +++ b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/modesT_vect.m @@ -159,6 +159,18 @@ nt = max(size(wn)); load BestView2 +% petits contrôles… + +pause +disp '+++' +binf +res +bsup +T +L +disp '+++' +pause + %---------------------------------- % boucle sur les nombres d'onde (sur les longeurs d'onde), ici trois. diff --git a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/wind_effect.m b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/wind_effect.m index 3029b2c..fccce9a 100644 --- a/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/wind_effect.m +++ b/docs/fac/stage/IPGP Saint Maur/MATLAB/vectorial/wind_effect.m @@ -31,7 +31,7 @@ end % these are obtained by running modesT_vect setting winds to 0 m/s prefix = 'data/wind/'; -res = 100; % FIXME faire gaffe ! +res = 1; % attention ! ne pas faire de tests à résolution inf. à 10, le code devient nécessairement instable binf = 1; bsup = 9999; @@ -47,7 +47,7 @@ for j = 1:length(L) name = [prefix 'data_wind0to0_L' num2str(L(j)) '_period' num2str(T(k)) '.txt']; % FIXME j'ai ajouté L, mais ça donne une dimension de plus dans les tableaux, à faire ! wind_ref(:,:,j,k) = data_processing(name, 'real', binf, res, bsup); - omega = 2.*pi./T; + %omega = 2.*pi./T; end end @@ -55,21 +55,21 @@ size(wind_ref) %pause % check -for j = 1:length(L) - for k = 1:length(T) - figure(k); - set(k, 'position', [10 1 750 600]); - set(gcf, 'Color', 'w'); - - plot(wind_ref(:,1,j,k), Z, 'r', 'linewidth', 2); - v = axis; - axis([v(1) v(2) Zmin Zmax]); - xlabel('Vertical V (m/s)', 'fontsize', 12); - ylabel('Altitude (km)', 'fontsize', 16); - title(['Period = ', num2str(T), ' mn' ' (omega = ', num2str(omega), ' rad/s)'], 'fontsize', 12); - pause - end -end +%for j = 1:length(L) + %for k = 1:length(T) + %figure(k); + %set(k, 'position', [10 1 750 600]); + %set(gcf, 'Color', 'w'); + + %plot(wind_ref(:,1,j,k), Z, 'r', 'linewidth', 2); + %v = axis; + %axis([v(1) v(2) Zmin Zmax]); + %xlabel('Vertical V (m/s)', 'fontsize', 12); + %ylabel('Altitude (km)', 'fontsize', 16); + %title(['Period = ', num2str(T), ' mn' ' (omega = ', num2str(omega), ' rad/s)'], 'fontsize', 12); + %pause + %end +%end % }}} @@ -80,8 +80,10 @@ modesT_vect(u0b, u0t, T, L, res); % outputs in data_wind[u0b]to[u0t]_period[T]. data = []; for j = 1:length(L) for k = 1:length(T) - name = [prefix 'data_wind' num2str(u0b) 'to' num2str(u0t) '_period' num2str(T(k)) '.txt']; - data(:,:,j,k) = data_processing(name, 'real', binf, res, bsup) + name = [prefix 'data_wind' num2str(u0b) 'to' num2str(u0t) '_L' num2str(L(j)) '_period' num2str(T(k)) '.txt']; + data(:,:,j,k) = data_processing(name, 'real', binf, 1, bsup) + % TODO FIXME là le modesT_vect a fait nawak, du 10e80 sur w, doit y avoir une connerie avec les arguments ! + % donc, les faire s'afficher, avec des pauses, ainsi que K etc. pause end end @@ -95,7 +97,6 @@ end % idées : % faire les deux, ie. pour chaque L, le max en fonction de T, puis pour chaque T, le max en fonction de L, % et essayer un pcolor pour avoir le max pour tous les couples L/T (mais ça doit plus ou moins revenir au même que le pcolor du script principal : ses pcolor plottent l'amplification relative par rapport au sol, ici ce serait, intégré, par rapport à un cas de référence à vent nul, donc c'est pas exactement la même chose quand même) -% FIXME ne pas faire de abs() en fait, sinon on ne voit pas si ça amplifie ou atténue ! % ne pas oublier de plotter le cas de référence à vent nul, d'ailleurs % un truc qui serait chanmé : un pcolor en couche, avec une croix indiquant la position du cas de référence, une autre indiquant la position du cas étudié là, le vecteur vent moyen, en gros. À réfléchir… @@ -134,11 +135,14 @@ amp_max = amp_max' % L in columns (one max per line for a given T), T in rows ( % ça permet de plotter direct versus L et T, pour un T ou un L donné respectivement ! % TODO plot plot ! +% FIXME pas vraiment sûr que ça fonctionne, les valeurs de amp_max sont cheloues +% tout vérifier pas à pas, ie. quelles valeurs ça manipule par rapport aux fichiers .txt +% (faire des load manuels et vérifier les amplitudes max etc.) hold on; figure(100); disp 'plot pour T = ' T(1) -plot(amp_max(:,1), L); +plot(L, amp_max(:,1)); %plot(T, amp_ref, '--b'); %title('maximum w amplitude for each period'); hold off; -- 2.11.4.GIT