Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Storing Iteration Results in Matlab

  1. Sep 2, 2012 #1
    Dear All,


    Please help me on this little problem.

    I have the code below.

    for k=40:10:80
    T(k)=273.15+k;
    z=[0.2 0.2 0.2 0.4];
    W_PR=0.245;
    C=4;
    omega=[0.344 0.467 0.578 0.789];
    Tc=[600 700 500 570];
    Pc=[50 70 58 76];
    for c=1:C

    x_PR(1,c)=z(c)/(1+W_PR*(K_PR(c)-1));
    x_PR(2,c)=K_PR(c)*x_PR(1,c);
    end
    for c=1:C
    kappa_PR=0.37464+1.54226.*omega(c)-0.26992.*omega(c).^2;
    alpha_PR=(1+kappa_PR.*(1-sqrt(T(k)./Tc(c)))).^2;
    a_PR(c,c)=0.45724.*R.^2.*Tc(c).^2./Pc(c).*alpha_PR;
    b_PR(c)=0.07780*R.*Tc(c)./Pc(c);
    end
    for c=2:C
    for n=1:(c-1)
    a_PR(c,n)=sqrt(a_PR(c,c).*a_PR(n,n));
    a_PR(n,c)=a_PR(c,n);
    end
    end
    for c=1:C
    A_PR(c,c)=a_PR(c,c).*P./(R.*T(k)).^2;
    B_PR(c)=b_PR(c).*P./(R.*T(k));
    end
    for c=1:C
    Z(c,c)=A_PR(c,c)./5;
    V(c)=B_PR(c).*6;
    end
    end

    QUESTION: Each time I run the code, I want result for Z and V at each T(k). The code as it is only gives result for one T value though I want it to run the loop and give result for Z and V for all T(k).

    Thanks

    adeeyo
     
  2. jcsd
  3. Sep 2, 2012 #2

    Simon Bridge

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

    [edited with apologies]
    you have managed to make a vector for T, and store the kth value in the kth position in the vector.
    So what is stopping you doing the same for Z and V.

    You'd have to put Z and V inside the k loop - encoding Z(k,c,c) and V(k,c)
     
    Last edited: Sep 2, 2012
  4. Sep 13, 2012 #3
    Hi Simon,
    I put Z and V inside k loop, but is giving me 80 by 4 by 4 matrix for Z instead of 5 by 4 by 4 matrix. Since my k=40:10:80; and T(k)=273.15+k;
    Any advice?

    Thanks
     
  5. Sep 13, 2012 #4

    Simon Bridge

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

    Well yes - you will have to change how these things get calculated too.
    Clearly you have missed a step or terminated the cycle too soon.
     
  6. Sep 13, 2012 #5
    Hi Simon,
    I have checked the code over and over again I could not figure out where I got it wrong.
    Please help me to check and verify the code.

    Thanks
     
  7. Sep 13, 2012 #6

    Simon Bridge

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

    flow-chart it.

    you can hire me to do it...
     
  8. Sep 13, 2012 #7
    Please assist. I am a student.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Storing Iteration Results in Matlab
  1. Iterating on Matlab (Replies: 0)

Loading...