i want to store the result of 191 iterations in one variable and then i want to find the slope of those values.

A = load ('example.txt') size(A) n = 191; for i = 1:191 A (i,:) A([i] , [2 3 4 5 6 7 8 9])

freq(i) = A (i,1) s11mag(i) = A(i,2) s11ang(i) = A(i,3) s12mag(i) = A(i,4) s12ang(i) = A(i,5) s21mag(i) = A(i,6) s21ang(i) = A(i,7) s22mag(i) = A(i,8) s22ang(i) = A(i,9)

% 1deg = 0.0174532925 radians q = 0.0174532925

s11real(i) = s11mag(i) * cos(s11ang(i) * q) s11img(i) = s11mag(i) * sin (s11ang(i) * q) s11(i) = s11mag(i).* exp( j * (s11ang(i) * q))

s12real(i) = s12mag(i) * cos(s12ang(i) * q) s12imag(i) = s12mag(i) * sin (s12ang(i) * q) s12(i) = s12mag(i).* exp( j * (s12ang(i) * q))

s21real(i) = s21mag(i) * cos(s21ang(i) * q) s21imag(i) = s21mag(i) * sin (s21ang(i) * q) s21(i) = s21mag(i).* exp( j * (s21ang(i) * q))

s22real(i) = s22mag(i) * cos(s22ang(i) * q) s22imag(i) = s22mag(i) * sin (s22ang(i) * q) s22(i) = s22mag(i).* exp( j * (s22ang(i) * q))

s = [s11 s12 ; s21 s22] D = ((1+s11(i))*(1+s22(i))-s12(i)*s21(i)) Z0 = 50

%s- parameter to y-parameter

y11(i) = ((1-s11(i))*(1+s22(i))+s12(i)*s21(i)) / (D * Z0) y12(i) = -2*s12(i) / (D * Z0) y21(i) = -2*s21(i) / (D * Z0) y22(i) = ((1+s11(i))*(1-s22(i))+s12(i)*s21(i)) / (D * Z0)

y = [y11 y12 ; y21 y22]

y11real(i) = real(y11(i)) y11imag(i) = imag (y11(i)) y12real(i) = real(y12(i)) y12imag(i) = imag (y12(i)) y21real(i) = real(y21(i)) y21imag(i) = imag (y21(i)) y22real(i) = real(y22(i)) y22imag(i) = imag (y22(i))

%extrinsic parasitic capacitances

Cgdo(i) = -y12imag(i) / (2 * pi * freq(i)) Cgda(i) = [0.01:0.5] * Cgdo(i) Cgdi(i) = 2 * Cgda(i) Cgd(i)= Cgdo(i) - Cgda(i) - Cgdi(i) Cdso(i) = y22imag(i) /(2 * pi * freq(i)) - Cgdo(i) Cpda(i) = [0.01:0.5]* Cdso(i) Cpdi(i) = 3* Cpda(i) Cds(i) = Cdso(i) - (4 * Cpda(i)) Cgs(i) = Cgd(i) Cpga(i) = [0.01:0.5] * Cdso(i)

Cgso(i) = y11imag(i) /(2 * pi * freq(i)) - Cgdo(i) Cpgi(i) = Cgso(i) - Cgd(i) - Cpga(i)

 %plotting capacitances vs frequency

% subplot(221);plot(freq,Cpga);hold on; % subplot(222);plot(freq,Cpgi);hold on; % subplot(223);plot(freq,Cpda);hold on; % subplot(224);plot(freq,Cpdi);hold on;

% de- embedding capacitances Cpga,Cpda and Cgda from Y-parameters

y11int(i) = y11(i) - (j* (2 * pi * freq(i)) (Cpga(i) + Cgda(i))) y12int(i) = y12(i) + (j (2 * pi * freq(i)) * Cgda(i)) y21int(i) = y21(i) + (j* (2 * pi * freq(i)) * Cgda(i)) y22int(i) = y22(i) - (j* (2 * pi * freq(i)) *(Cgda(i) +Cpda(i)))

yint = [y11int y12int ; y21int y22int]

y11intreal(i) = real(y11int(i)) y11intimag(i) = imag (y11int(i)) y12intreal(i) = real(y12int(i)) y12intimag(i) = imag (y12int(i)) y21intreal(i) = real(y21int(i)) y21intimag(i) = imag (y21int(i)) y22intreal(i) = real(y22int(i)) y22intimag(i) = imag (y22int(i))

%converting Yint-para to z-parameters

z11(i) = y22int(i) / (y11int(i) * y22int(i) - y12int(i) * y21int(i)) z12(i) = -y12int(i) / (y11int(i) * y22int(i) - y12int(i) * y21int(i)) z21(i) = -y21int(i) / (y11int(i) * y22int(i) - y12int(i) * y21int(i)) z22(i) = y11int(i) / (y11int(i) * y22int(i) - y12int(i) * y21int(i))

z = [z11 z12 ; z21 z22]

w(i) =(2* pi* freq(i))

z11real(i) = real (z11(i)) z11imag(i) = imag (z11(i)) z12real(i) = real (z12(i)) z12imag(i) = imag (z12(i)) z21real(i) = real (z21(i)) z21imag(i) = imag (z21(i)) z22real(i) = real (z22(i)) z22imag(i) = imag (z22(i))

m(i) = w(i) * (z11imag(i) - z12imag(i)) n(i) = w(i) * (z22imag(i) - z12imag(i)) o(i) = w(i) * z12imag(i)

subplot(221);plot(w(i)^2, m(i)); hold on; subplot(222);plot(w(i)^2, n(i)); hold on; subplot(223);plot(w(i)^2, o(i)); hold on; subplot(224);plot(w(i)^2,m(i),’b’, w(i)^2,n(i),’r’, w(i)^2,o(i),’g’); hold on ;

end

[Matt Fig]

Two things:

1. FORMAT YOUR CODE!
2. Ask a question.

[ria]

        A = load ('example.txt')
        size(A) n = 191; 

        for i = 1:191 
        A (i,:) 
        A([i] , [2 3 4 5 6 7 8 9])

        freq(i) = A (i,1)
        s11mag(i) = A(i,2) 
        s11ang(i) = A(i,3) 
        s12mag(i) =  A(i,4)
        s12ang(i) = A(i,5) 
        s21mag(i) = A(i,6) 
        s21ang(i) = A(i,7) 
        s22mag(i) = A(i,8) 
        s22ang(i) = A(i,9)

        % 1deg = 0.0174532925 radians q = 0.0174532925

        s11real(i) = s11mag(i) * cos(s11ang(i) * q) 
        s11img(i) = s11mag(i) * sin (s11ang(i) * q) 
        s11(i) = s11mag(i).* exp( j * (s11ang(i) * q))

        s12real(i) = s12mag(i) * cos(s12ang(i) * q) 
        s12imag(i) = s12mag(i) * sin (s12ang(i) * q)
        s12(i) = s12mag(i).* exp( j * (s12ang(i) * q))

        s21real(i) = s21mag(i) * cos(s21ang(i) * q) 
      s21imag(i) = s21mag(i) * sin (s21ang(i) * q) 
      s21(i) = s21mag(i).* exp( j * (s21ang(i) * q))

        s22real(i) = s22mag(i) * cos(s22ang(i) * q) 
      s22imag(i) = s22mag(i) * sin (s22ang(i) * q) 
      s22(i) = s22mag(i).* exp( j * (s22ang(i) * q))

        s = [s11 s12 ; s21 s22] 

      D = ((1+s11(i))*(1+s22(i))-s12(i)*s21(i)) 
      Z0 = 50

        %s- parameter to y-parameter

        y11(i) = ((1-s11(i))*(1+s22(i))+s12(i)*s21(i)) / (D * Z0)

       y12(i) = -2*s12(i) / (D * Z0) 

      y21(i) = -2*s21(i) / (D * Z0)

       y22(i) = ((1+s11(i))*(1-s22(i))+s12(i)*s21(i)) / (D * Z0)

        y = [y11 y12 ; y21 y22]

        y11real(i) = real(y11(i)) 
      y11imag(i) = imag (y11(i))
       y12real(i) = real(y12(i))
       y12imag(i) = imag (y12(i)) 
      y21real(i) = real(y21(i)) 

      y21imag(i) = imag (y21(i)) 
      y22real(i) = real(y22(i))
       y22imag(i) = imag (y22(i))

        %extrinsic parasitic capacitances

        Cgdo(i) = -y12imag(i) / (2 * pi * freq(i)) 
      Cgda(i) = [0.01:0.5] * Cgdo(i) Cgdi(i) = 2 * Cgda(i) Cgd(i)= Cgdo(i) - Cgda(i) - Cgdi(i) Cdso(i) = y22imag(i) /(2 * pi * freq(i)) - Cgdo(i) Cpda(i) = [0.01:0.5]* Cdso(i) Cpdi(i) = 3* Cpda(i) Cds(i) = Cdso(i) - (4 * Cpda(i)) Cgs(i) = Cgd(i) Cpga(i) = [0.01:0.5] * Cdso(i)

        Cgso(i) = y11imag(i) /(2 * pi * freq(i)) - Cgdo(i) Cpgi(i) = Cgso(i) - Cgd(i) - Cpga(i)

         %plotting capacitances vs frequency

        % subplot(221);plot(freq,Cpga);hold on; % subplot(222);plot(freq,Cpgi);hold on; % subplot(223);plot(freq,Cpda);hold on; % subplot(224);plot(freq,Cpdi);hold on;

        % de- embedding capacitances Cpga,Cpda and Cgda from Y-parameters

        y11int(i) = y11(i) - (j* (2 * pi * freq(i)) (Cpga(i) + Cgda(i))) y12int(i) = y12(i) + (j (2 * pi * freq(i)) * Cgda(i)) y21int(i) = y21(i) + (j* (2 * pi * freq(i)) * Cgda(i)) y22int(i) = y22(i) - (j* (2 * pi * freq(i)) *(Cgda(i) +Cpda(i)))

        yint = [y11int y12int ; y21int y22int]

        y11intreal(i) = real(y11int(i))
 y11intimag(i) = imag (y11int(i))
 y12intreal(i) = real(y12int(i)) 
y12intimag(i) = imag (y12int(i)) 
y21intreal(i) = real(y21int(i)) 
y21intimag(i) = imag (y21int(i))
 y22intreal(i) = real(y22int(i))
 y22intimag(i) = imag (y22int(i))

        %converting Yint-para to z-parameters

        z11(i) = y22int(i) / (y11int(i) * y22int(i) - y12int(i) * y21int(i)) z12(i) = -y12int(i) / (y11int(i) * y22int(i) - y12int(i) * y21int(i)) z21(i) = -y21int(i) / (y11int(i) * y22int(i) - y12int(i) * y21int(i)) z22(i) = y11int(i) / (y11int(i) * y22int(i) - y12int(i) * y21int(i))

        z = [z11 z12 ; z21 z22]

        w(i) =(2* pi* freq(i))

        z11real(i) = real (z11(i)) 
      z11imag(i) = imag (z11(i)) 
      z12real(i) = real (z12(i))
       z12imag(i) = imag (z12(i))
       z21real(i) = real (z21(i)) 
      z21imag(i) = imag (z21(i)) 
      z22real(i) = real (z22(i))
       z22imag(i) = imag (z22(i))

        m(i) = w(i) * (z11imag(i) - z12imag(i))
       n(i) = w(i) * (z22imag(i) - z12imag(i))
       o(i) = w(i) * z12imag(i)

        subplot(221);plot(w(i)^2, m(i)); hold on;
       subplot(222);plot(w(i)^2, n(i)); hold on; 
      subplot(223);plot(w(i)^2, o(i)); hold on; subplot(224);plot(w(i)^2,m(i),’b’, w(i)^2,n(i),’r’, w(i)^2,o(i),’g’); hold on ;

        end

[Image Analyst]

A tip: in the editor, type control-A and then control-I and it will make your code look nice - fix all the indenting.

Answer 1

[John Petersen]

gradient(y)

will give you the slope at every point of the variable y.

[ria]

slope is between two variables X and Y .then i can i give the command as gradient(y)..please help..

[John Petersen]

The default spacing used by gradient() is sequential points. The actual spacing of those points is defined by the user. You CAN specify a different spacing if you want. Here's an example:

t = 0:.01:10;
ct = cos(t);
slope_of_ct = gradient(ct);
figure;plot(t,ct,t,slope_of_ct); grid on;
legend('signal','slope of signal');