Numerical integration for a vector.

2 ビュー (過去 30 日間)
mads skibsted
mads skibsted 2022 年 4 月 17 日
コメント済み: Star Strider 2022 年 4 月 17 日
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 件のコメント
mads skibsted
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
Star Strider 2022 年 4 月 17 日
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.

サインインしてコメントする。

回答 (0 件)

カテゴリ

Help Center および File ExchangeNumerical Integration and Differentiation についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by