Numerical integration for a vector.

mads skibsted

2022 4 17

0 个回答

4 次查看（30 天）

0 个投票

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 个评论
显示 3更早的评论隐藏 3更早的评论

mads skibsted 2022-4-17

And at t=0, then the higest value of u_e is equal u_i, where u_i is a constant.

Star Strider 2022-4-17

在 MATLAB Online 中打开

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