解决三重积分的功能,上传它的价值每一步

1视图(30天)

显示旧的评论

路易吉Stragapede 2020年7月14日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/564833-solve-triple-integral-of-a-function-that-upload-its-value-every-step

评论道: 沃尔特·罗伯森 2020年7月17日

我必须解决这个三重积分:

P_4 = integral3 (fun_r R1、R2 0,τ,z3、z4);

为了解决这个积分,它是必要的使用一个函数处理。

所以:

fun_r = @ (R,θ,z) abs (Jr4) ^ 2 * R;

问题是,fun_r取决于Jr4但是Jr4上传它的价值在每个迭代中,总结所有先前的值(它是一个总和)。

如果我考虑Jr4(仅为当前迭代计算),我把它立即fun_r括号内,它的工作原理。

但是我需要求和,为了做到这一点,我需要定义外Jr4函数处理。但Jr4开掘R,θ和z。所以我想这样做:

                         信谊zθR
                        
                         Xr4 = A4d。* exp (gamma4。* z) + B4d。* exp (-gamma4。* z);
                        
                         Jr4 =实际(Jr4 +(真实((Xr4。* (1. / R)。* besselj (np, alfak。* R)。* exp(1我。* np。* (theta-OMEGA。* t))))));
                        
                         fun_r = @ (R,θ,z) abs (Jr4) ^ 2 * R;
                        
                         P_4 = integral3 (fun_r R1、R2 0,τ,z3、z4);

错误的是:

                         错误使用integral2Calc > integral2t /张量(第231行)
                        
                         输入函数必须返回“双重”或“单一”的价值观。发现“符号”。
                        
                         错误在integral2Calc > integral2t(55行)
                        
                         [Qsub, esub] =张量(thetaL、thetaR phiB, phiT);
                        
                         错误在integral2Calc(第9行)
                        
                         [q, errbnd] = integral2t(乐趣、xmin xmax、ymin ymax, optionstruct);
                        
                         错误在integral3 / innerintegral(第137行)
                        
                         Q1 = integral2Calc (…
                        
                         错误在integralCalc / iterateScalarValued(第314行)
                        
                         fx =乐趣(t);
                        
                         错误在integralCalc / vadapt(第132行)
                        
                         [q, errbnd] = iterateScalarValued (u tinterval pathlen);
                        
                         错误在integralCalc(第75行)
                        
                         [q, errbnd] = vadapt (@AtoBInvTransform间隔);
                        
                         错误在integral3(第121行)
                        
                         Q = integralCalc (@innerintegral xmin、xmax integralOptions);
                        
                         错误在Untitled5(第91行)
                        
                         P_4 = integral3 (fun_r R1、R2 0,τ,z3、z4);