清空筛选条件
清空筛选条件

How to take an integral symbolically and then convert it to type double?

2 次查看（过去 30 天）
L'O.G.
L'O.G. 2023-5-16
评论： Walter Roberson 2023-5-16
How do I take the following integral symbolically? And convert it to double? I have a 1000 x 1 numeric vector of type double for A(t).
The following takes a minute or two to run and I'm not sure how to plot it to check.
syms t omega real
int(A*sin(omega*t),0,2*pi/omega)
  6 个评论
显示 4更早的评论
L'O.G.
L'O.G. 2023-5-16
Thank you. I think the symbolic aspect is confusing me, so I am trying to do it in a more straightforward numerical approach.
Walter Roberson
Walter Roberson 2023-5-16
Ran in:
format long g
syms t omega real
A = rand(5,1)
A = 5×1
0.988717788679954 0.961075790428598 0.548371988302065 0.728110260321502 0.293761135174945
Asym = sym(A);
symbolic_result = int(Asym*sin(omega*t),0,2*pi/omega)
symbolic_result = 
omega = abs(randn())
omega =
1.29628039951883
numeric_result = integral(@(t) A .* sin(omega.*t), 0, 2*pi./omega, 'Arrayvalued', true)
numeric_result = 5×1
1.0e+00 * -2.20309881449055e-16 4.33680868994202e-17 -5.55111512312578e-17 -2.91433543964104e-16 -6.11490025281825e-17
Zero to within round-off error.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

VBBV
VBBV 2023-5-16
编辑：VBBV 2023-5-16
Ran in:
syms t omega real
A = rand(1000,1);
% specify the integration variable t as the argument
I = double(vpa(int(A.*sin(omega*t),t,0,2*pi/omega),2))
I = 1000×1
0 0 0 0 0 0 0 0 0 0

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Assumptions 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by