Clear Filters
Clear Filters

Numerical integration for a vector.

3 views (last 30 days)
mads skibsted
mads skibsted on 17 Apr 2022
Commented: Star Strider on 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 Comments
mads skibsted
mads skibsted on 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 on 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.

Sign in to comment.

Answers (0)

Community Treasure Hunt

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

Start Hunting!