Double Integral using Integral2 Error

1 次查看(过去 30 天)
I am trying to write code for the attached equation in matlab but I am getting errors.
My code is as under;
phi=0.0:0.01:90;
P=sqrt((pi^2*Ef*df^3*rounot*S*(1+n))/(4)+(pi^2*Ef*df^3*Gd)/(2));
p_phi=sin(phi);
p_z=2/Lf;
m1=P*exp(f*phi); %multiplier,inside the integral part of eq#8.
m2=(4*Vf)/(pi*df^2); %multiplier, outside the integral part of eq#8.
myfun_eq8=@(z,phi) m1.*p_phi.*p_z;
I= integral2(myfun_eq8,0,1,0,pi/2);
sigma_b=m2*I;
I am getting the following error, can someone help me correct it.
Error using integral2Calc>integral2t/tensor (line 241) Integrand output size does not match the input size.
Error in integral2Calc>integral2t (line 55) [Qsub,esub] = tensor(thetaL,thetaR,phiB,phiT);
Error in integral2Calc (line 9) [q,errbnd] = integral2t(fun,xmin,xmax,ymin,ymax,optionstruct);
Error in integral2 (line 106) Q = integral2Calc(fun,xmin,xmax,yminfun,ymaxfun,opstruct);
Error in myeqsetlin (line 160) I= integral2(myfun_eq8,0,1,0,pi/2);
  6 个评论
Andrew Newell
Andrew Newell 2017-4-18
I notice that p(z) is a constant, so integrating p(z)dz gives simply cos(phi). Thus, you really have a one-dimensional integral to solve numerically.
jack carter
jack carter 2017-4-18
Thanks Andrew for your reply. I was writing a series of equations in a function file to run a simulation after completion. That is why there seem to had variation in yours and mine results. I wrote delta as S in my script.
I re-ran the above mentioned script separately and I have received the same error as you mentioned in your first reply for m1. If I define as a range of numeric values (S=0:0.01:30), the function still does not work. But if I select S=1; then I get the same set of integral errors.
Your simplification of the equation i.e. p(z) as a constant is helpful

请先登录,再进行评论。

采纳的回答

Andrew Newell
Andrew Newell 2017-4-18
Based on the above discussion, the critical part of code you need is as follows:
myfun_eq8=@(phi) exp(f*phi).*cos(phi).*sin(phi);
I= integral(myfun_eq8,0,pi/2);
All the other terms can stay outside of the integral.

更多回答(2 个)

David Goodmanson
David Goodmanson 2017-4-17
编辑:David Goodmanson 2017-4-18
Hi jack, I believe this is occurring because your definition of myfun_eqn8 uses the two variables z,phi, but the integrand m1.*p_phi.*p_z is not explicitly a function of either of those variables. Matlab does not know how to create the integrand for arbitrary z,phi. You need to define m1, p_phi,p_z explicitly in myfun_eqn8.
Once this gets working the first and third lines of code will probably become superfluous anyway, but right now you are taking the sine of angles in degrees. You can convert to radians or use the sind function there.
  5 个评论
David Goodmanson
David Goodmanson 2017-4-18
编辑:David Goodmanson 2017-4-18
ok then it does appear to reduce to a one-dimensional integral in the variable phi, and the answer for that only depends on the constant f.

请先登录,再进行评论。


Shahid Hasnain
Shahid Hasnain 2018-11-15
I have following expression with x=[1 2 3 4 5], y=[1 2 3 4 5],
L = 100; % Length of the box
p=1;
q=1;
Ui = @(x,y) sin(pi*x).*sin(pi*y); % An initial distribution for u, i.e. the u(x,y,0)
Vi = @(x,y) sin(pi*x).*sin(pi*y); % An initial distribution for v, i.e. the v(x,y,0)
Wi = Ui(x,y) + Vi(x,y); % Initial condition for w, i.e. w(x,y,0)
% Computation of the b-coefficients
if p < Ntrunc
if q < Ntrunc
my_integrand = @(x,y) 4*Wi.*sin(p.*x).*cos(q.*y)/L^2;
bd(p,q) = integral2(my_integrand(x,y),0,L,0,L);
q = q+1;
end
p = p+1;
end
I got following error,
Error using integral2Calc>integral2t/tensor (line 231)
Input function must return 'double' or 'single' values. Found 'function_handle'.
Error in integral2Calc>integral2t (line 55)
[Qsub,esub] = tensor(thetaL,thetaR,phiB,phiT);
Error in integral2Calc (line 9)
[q,errbnd] = integral2t(fun,xmin,xmax,ymin,ymax,optionstruct);
Error in integral2 (line 106)
Q = integral2Calc(fun,xmin,xmax,yminfun,ymaxfun,opstruct);
Correction will be highly appreciated.

类别

Help CenterFile Exchange 中查找有关 Matrix Computations 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by