clear all close all clc %%%%% %% %Vitesse du vent m/s wind_speed_2 = [ 0.00200000000000000 0.118000000000000 0.632000000000000 0.805000000000000 0.999000000000000 1.19000000000000 1.21900000000000 1.44500000000000 1.59600000000000 1.78400000000000 1.86500000000000 2.07100000000000 2.10700000000000 1.76000000000000 2.51300000000000 2.42900000000000 2.71900000000000 2.86300000000000 3.01300000000000 3.32900000000000 3.69200000000000 3.96100000000000 4.19500000000000 5.10800000000000 5.28000000000000 5.95300000000000 7.11600000000000]; %Vitesse du feuillage (cm/s) Foliage_velocity = [0 0 0 0 1.96508278500587 2.09086791224201 2.10720024656275 2.44840766015495 2.48309431881800 2.68287974524885 2.76434370644185 3.16378587612085 2.94795086646812 3.10859959240870 3.56365355443517 3.48920114147074 3.83884402843078 4.04247846932406 4.42206131685584 4.73589966837583 4.56677359976928 3.51960622653381 4.02499187591986 3.82804823183098 4.65476326886195 4.49574241775635 5.52948465445743]; %Vitesse Large scale Vitesse_l=[0 0 0 0 0 0 0 0 0.180368886434187 0.478981285093914 0.660807899310191 0.710310680086996 0.577569779839189 0.443398019670609 1.03430059546277 2.34842490558222 0.522732325996892 0.933163522342116 1.71855035884597 2.51359703543046 2.14473897578224 1.24246265273439 2.46681252005416 3.13978966639822 3.72028254479681 3.28057796284205 4.56685772224079]; %Vitesse Small scale Vitesse_s= [0 0 0 0 1.96508278500587 2.09086791224201 2.10720024656275 2.44840766015495 2.47653476877525 2.63977662994359 2.68419988964197 3.08301797070600 2.89081778404081 3.07681485014548 3.41025643233229 2.58058618698966 3.80308748124159 3.93329861763872 4.07445836328501 4.01379815292168 4.03181297151917 3.29300995237121 3.18047097017484 2.18990276405215 2.79755587540741 3.07400522389420 3.11753285286456]; %velocity nomega=0; %displacement %nomega=-1; facteur=1.8; %nomega=1; figure('Position',[100 100 400 400], 'PaperPosition', [0 0 3.5 3.5]) loglog(wind_speed_2 ,.01*Vitesse_l.*(2*pi*1)^nomega, 'ok', 'Markersize', 6); hold on; loglog(wind_speed_2 ,facteur*0.002*wind_speed_2.^(11/6), '-k', 'LineWidth',1) loglog([1.2 1.2],facteur*[0.0028 0.008],'-k') loglog([1.2 2.13],facteur*[0.008 0.008],'-k') text(1.5,facteur*0.01, '11') xlim([.8, 10]) ylim([.001,.1]) %xlabel('Wind mean velocity (m/s)') %ylabel('Foliage RMS velocity (m/s)') xlabel('Wind') ylabel('Fol') axis square ax=gca; ax.XTick=[1 10]; ax.XTickLabel={'1', '10'}; ax.YTick=[0.001 0.01 0.1]; ax.YTickLabel={'a', 'b', 'c'}; ax.XLabel.Position=[2.8284 0.0007 -1.0000]; ax.YLabel.Position=[0.5396 0.0100 -1.0000]; %saveas(gcf, '/Users/loictadrist/Documents/Ecole Polytechnique/Article2017/Latex-rsos/Fig/Branch_scaling.eps', 'eps2') %%check %plot(Foliage_velocity, (Vitesse_l.^2+Vitesse_s.^2).^.5) %% repartition %lx= (Vitesse_l./Vitesse_s); %plot(wind_speed_2,lx,'ok') % lx= (Vitesse_s./ Foliage_velocity'); %plot(wind_speed_2,lx.^2,'ok') %hold on ; % lx= (Vitesse_l./ Foliage_velocity'); % plot(wind_speed_2,lx.^2,'ok','Markersize', 6, 'MarkerFaceColor', 'k'); % hold off;