Borrar filtros
Borrar filtros

Numerical integration for a vector.

2 visualizaciones (últimos 30 días)
mads skibsted
mads skibsted el 17 de Abr. de 2022
Comentada: Star Strider el 17 de Abr. de 2022
Can anybody help me with this? I do have problems to be integrating, when z is also defined as a vector, so I can't be using syms.
  5 comentarios
mads skibsted
mads skibsted el 17 de Abr. de 2022
And at t=0, then the higest value of u_e is equal u_i, where u_i is a constant.
Star Strider
Star Strider el 17 de Abr. de 2022
Try either of these, after supplying , , and t
z = ...;
d = numel(z)/2;
n = linspace(0, 1E4, 2*d);
ftrpz = (ui,cv,t,z) trapz(ui.*sin(pi*n.*z/(2*d))) .* sin(pi*n.*z/(2*d)) .* exp(-n.^2.*pi^2.*cv.*t/(4*d^2));
fintg = (ui,cv,t,z) integral(@(z)ui.*sin(pi*n.*z/(2*d)), 0, inf, 'ArrayValued',1) .* sin(pi*n.*z/(2*d)) .* exp(-n.^2.*pi^2.*cv.*t/(4*d^2))
See the documentation on trapz and integral for details. Then decide which version of ‘f’ works best in your application. Use the sum function complete the expression.
I cannot test this, so I am posting it as UNTESTED CODE. It will be necessary to experiment with it with your data to produce the desired result. The variable does not appear anywhere, so I have no idea what to do with it.
There is not enough information to proceed further.

Iniciar sesión para comentar.

Respuestas (0)

Categorías

Más información sobre Numerical Integration and Differentiation en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by