Numerical integration for a vector.

3 Ansichten (letzte 30 Tage)
mads skibsted
mads skibsted am 17 Apr. 2022
Kommentiert: Star Strider am 17 Apr. 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 Kommentare
mads skibsted
mads skibsted am 17 Apr. 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 am 17 Apr. 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.

Melden Sie sich an, um zu kommentieren.

Antworten (0)

Community Treasure Hunt

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

Start Hunting!

Translated by