How to write summation in matlab with two variables?

10 Ansichten (letzte 30 Tage)
Alicia
Alicia am 21 Apr. 2021
Bearbeitet: Alicia am 22 Apr. 2021
I am having trouble transferring these formulas onto matlab. My code compiles, however, my variables are all 1x1 and 0. I'm not sure why it is not calculating these values using my data properly. Here is my code so far:
syms u(i) v(i)
ubar = (1/N)*symsum(u,i,1,N);
uprime2 = (1/N)*symsum(u-ubar,i,1,N);
vbar = (1/N)*symsum(v,i,1,N);
vprime2 = (1/N)*symsum(v-vbar,i,1,N);
uprimevprime = (1/N)*symsum((u-ubar)*(v-vbar),i,1,N);
  2 Kommentare
Paul Hoffrichter
Paul Hoffrichter am 21 Apr. 2021
Bearbeitet: Paul Hoffrichter am 21 Apr. 2021
Could you tell me what the last equation represents? I was a little surprised about the subscript i being used for both u and v.
Alicia
Alicia am 22 Apr. 2021
Yes, u represents the 2D matrix of the u-component of the velocity on the vector grid with likely erroneous vectors replaced by NaN with units of m/s, and v represents the 2D matrix of the v-component of the velocity on the vector grid with likely erroneous vectors replaced by NaN with units of m/s

Melden Sie sich an, um zu kommentieren.

Akzeptierte Antwort

Paul Hoffrichter
Paul Hoffrichter am 21 Apr. 2021
I am new to sym, so I thought I would give this a try. Please check carefully, as I did not have time to do so. I hope you are doing this in Live Editor. Will be back in 16 hours to see if any problems.
N=4;
u = (sym('u%d',[1 N])).';
ubar = mean(u);
u_ubar = (u - ubar);
uprime2 = mean( u_ubar.^2);
v = (sym('v%d',[1 N])).';
vbar = mean(v)
v_vbar = (v - vbar)
vprime2 = mean(v_vbar.^2)
uprimevprime = mean(dot(u_ubar, v_vbar))
  1 Kommentar
Alicia
Alicia am 22 Apr. 2021
That gave me much better results, I am not receiving any 0 values anymore. Thank you for your help

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Kategorien

Mehr zu Energy Production finden Sie in Help Center und File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by