Integral calculation using matlab

3 次查看（过去 30 天）
MINAS KOUROUMPLAKIS
编辑： Bruno Luong 2020-8-4
Hello i wish to calculate the following integral
My difficulty stems from the absolute value which is "outside" of the integral with respect to θ

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

采纳的回答

Bruno Luong
Bruno Luong 2020-8-4
编辑：Bruno Luong 2020-8-4
Feel free to change these hyper-parameters
n = 3
k = 2*pi*rand()
integral(@(phi) arrayfun(@(phi) abs(integral(@(theta) exp(1i*k*cos(theta-phi)).*cos(n*theta), 0, 2*pi)).^2, phi), 0, 2*pi)
  4 个评论
显示 2更早的评论
Matt J
Matt J 2020-8-4
Looks like a square .^2 is missing from the absolute value.
Bruno Luong
Bruno Luong 2020-8-4
Oh heah, thank you for spotting the error Matt. I'll fix it.

请先登录，再进行评论。

更多回答（1 个）

Matt J
Matt J 2020-8-4
编辑：Matt J 2020-8-4
My difficulty stems from the absolute value which is "outside" of the integral
Why is this a difficulty? Because it prevents you from using integral2()? Just apply integral() twice instead.
  5 个评论
显示 3更早的评论
MINAS KOUROUMPLAKIS
MINAS KOUROUMPLAKIS 2020-8-4
This is what i do and i get an error message
fun=@(x,y)exp(-1i*cos(x-y)).*cos(x)
fun = @(x,y)exp(-1i*cos(x-y)).*cos(x)
>> integral(@(x) f(x,y),0,2*pi)
Error using subsindex
Function 'subsindex' is not defined for values of class 'function_handle'.
Error in sym/subsref (line 771)
R_tilde = builtin('subsref',L_tilde,Idx);
Steven Lord
Steven Lord 2020-8-4
Write a function to compute the inner integral. This function should accept phi. One wrinkle here is that since the outer integral call is going to pass a vector of values into it as phi, your inner integral call will need to say that the function is ArrayValued.
Write a function that takes the absolute value of that inner integral function.
Write a function to compute the outer integral of the absolute value function.
If you try this and it doesn't work, show us the code you wrote and we may be able to offer some guidance.

请先登录，再进行评论。

类别

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

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by