MATLAB: Writing this Mathcad equation in Matlab

bessel functionpsf

hi, Im trying to convert an equation from mathcad into matlab but am struggling. the equation is:

Best Answer

Try this:

psi = 0.142;
Io = 1;
f_vir2 = @(r) 4*Io * abs(integral(@(rho)besselj(0, pi*rho*r*psi).*rho, 0, 1)).^2;
r = pi;
test_call = f_vir2(r)

Related Solutions

MATLAB: Double integral of a dot product error

integral2() assumes that the function handle is able to accept matrix inputs. However, the fun() function in your code does not work correctly with matrix inputs. The following codes shows a quick fix using arrayfun

light_dir = [-0.6626;0;-0.7490];
camera_dir = [-1;0;0];
R = sqrt(2/(4*pi));
x_c =@(theta,phi) R.*sin(theta).*cos(phi);
y_c =@(theta,phi) R.*sin(theta).*sin(phi);
z_c =@(theta) R.*cos(theta);
normal =@(theta,phi) [x_c(theta,phi); y_c(theta,phi); z_c(theta)];
fun = @(theta,phi) (dot(light_dir,normal(theta,phi))).*(dot(camera_dir,normal(theta,phi))<0).*R^2.*sin(theta);
intensity = integral2(@(THETA, PHI) arrayfun(@(theta,phi) fun(theta, phi), THETA, PHI) ,0,pi,0,2*pi);

MATLAB: Error while executing integral function

Change the line

fun = @(r) F(r)*b*r;

fun = @(r) F(r)*b(r)*r;

and replace

A(ii) = a * integral(@(r) fun,0,1,"ArrayValued",true);

with

A(ii) = a * integral(fun,0,1,"ArrayValued",true);

Best Answer

Related Solutions

MATLAB: Double integral of a dot product error

MATLAB: Error while executing integral function

Related Question