Numerical integration for a vector.
3 views (last 30 days)
Show older comments
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
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.
Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!