Hello everyone, I have a problem and need its solution. I have applied loop (say chi angle iteration) on input data that was loaded as .dat file. Now I want to see results (graphically and numerically) against each rotation angle, how can I find

Hello, This code does not have any trouble. Problem is to move next step. I will try to clear my question. See, I have three input columns of P wave velocity (Vp), S wave velocity (Vs) and density (Rho). Each column (Vp, Vs, Rho) has 202 elements (numerical values) and by using these elements I have computed EEI (a fourth column). But EEI is also dependent on the angle (Chi angle or rotation angle). So I used iteration/loop and try to vary chi angles between -90 to +90. In other words, by using three columns of Vp, Vs, Rho I need to compute 181 EEI columns (between -90 to +90). But I am not clear how to plot and show these 181 graphs of EEI further

%% The Input data comprizes on seismic velocities and density...

load VelRhob.dat -ASCII % load input file in ascii format

size(VelRhob)

     Depth = VelRhob(:,1);     % Depth in meters
     Vp = VelRhob(:,2);        % P wave velocity in km/s
     Vs = VelRhob(:,3);        % S wave velocity in km/s
     Rho = VelRhob(:,4);       % Density in gm/cm3

     Vp2 = power(Vp,2);        % Element-wise power/sq. of input P velocity
     Vs2 = power(Vs,2);        % Element-wise power/sq. of input S velocity

     Vp0 = mean(Vp);           % Mean of P velocity
     Vs0 = mean(Vs);           % Mean of S velocity
     Rho0 = mean(Rho);         % Mean of density
     k = (Vs0/Vp0)^2;          % Constant for EEI-- Extended elastic imp.
     a = Vp0*Rho0;
     min_angle = -90;
     max_angle = 90;
     int_angle = 5;

     angle = min_angle:int_angle:max_angle;

for ii=1:length(angle)

      Chi(ii) = angle(ii)*pi/180;

      p = cos(Chi(ii))+sin(Chi(ii));
      q = -8*k*sin(Chi(ii));
      r = cos(Chi(ii))-4*k*sin(Chi(ii));
      eei1 = rdivide(Vp,Vp0);
      eei2 = rdivide(Vs,Vs0);
      eei3 = rdivide(Vs,Vs0);
      eei11(ii) = power(eei1(ii),p);
      eei12(ii) = power(eei2(ii),q);
      eei13(ii) = power(eei3(ii),r);
      eei(ii) = times(eei11(ii),eei12(ii));
      eeii(ii) = times(eei(ii),eei13(ii));
      EEI(ii) = times(a,eeii(ii));
   I need to find out EEI as function of Chi angle on coulmns data

[Jan]

The question is not clear. Does the shown code cause any troubles? If not, how does it help to understand the problem? Please explain "see results (graphically and numerically)" with any details.

[Image Analyst]

"I want to see results (graphically..."

Are you just looking for someone to tell you that there is a "plot()" function? Or area() or bar() or other functions?

[Ahmed]

Dear Simon,

This code does not have any trouble. Problem is to move next step. I will try to clear my question.

See, I have three input columns of P wave velocity (Vp), S wave velocity (Vs) and density (Rho). Each column (Vp, Vs, Rho) has 202 elements (numerical values) and by using these elements I have computed EEI (a fourth column). But EEI is also dependent on the angle (Chi angle or rotation angle). So I used iteration/loop and try to vary chi angles between -90 to +90. In other words, by using three columns of Vp, Vs, Rho I need to compute 181 EEI columns (between -90 to +90). But I am not clear how to plot and show these 181 graphs of EEI further