How to write summation in matlab with two variables?

29 visualizaciones (últimos 30 días)
Alicia
Alicia el 21 de Abr. de 2021
Editada: Alicia el 22 de Abr. de 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 comentarios
Paul Hoffrichter
Paul Hoffrichter el 21 de Abr. de 2021
Editada: Paul Hoffrichter el 21 de Abr. de 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 el 22 de Abr. de 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

Iniciar sesión para comentar.

Respuesta aceptada

Paul Hoffrichter
Paul Hoffrichter el 21 de Abr. de 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 comentario
Alicia
Alicia el 22 de Abr. de 2021
That gave me much better results, I am not receiving any 0 values anymore. Thank you for your help

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Numbers and Precision en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by