我需要解决以下三重积分

1视图(30天)

显示旧的评论

路易吉Stragapede 2020年7月13日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/564224-i-need-to-solve-the-following-triple-integral

关闭:MATLAB回答机器人 2021年8月20日

我需要解决以下三重积分:

输入数据:

                         np = 10;
                        
                         sigmac = 16 e6;
                        
                         τ= 0.9;
                        
                         = 10 e - 3;
                        
                         b = 10 e - 3;
                        
                         c = 5 e - 3;
                        
                         d = 5 e - 3;
                        
                         e = 8 e - 3;
                        
                         咱=一个;
                        
                         zb =咱+ b;
                        
                         佐= zb + c;
                        
                         zd =佐+ d;
                        
                         泽= zd + e;
                        
                         R1 = 25 e - 3;
                        
                         R2 = 65 e - 3;
                        
                         R3 = 90 e - 3;
                        
                         A4d = 1.62 e - 04 - 9.53 - e - 04 -我;
                        
                         B4d = 1.03 e + 12 - 2.13 - e + 12我;
                        
                         gamma4 e = 5.90 + 2 + 2.14 e + 02年我;
                        
                         alfak = 550.0686;
                        
                         ω= 1;
                        
                         t = 10;

和:

我试图使用这个代码,但它给了我一些错误:

                         Xr4 = A4d * exp (gamma4 * z) + B4d * exp (-gamma4 * z);
                        
                         Jr4 =实际((Xr4 * (1 / R) * besselj (np, alfak * R) * exp(我* np * (theta-OMEGA * t))));
                        
                         Jtheta4 =实际((1 * (1 / np) * Xr4 * (alfak * besselj (np-1 alfak * R) - (np / R) * besselj (np, alfak * R)) * exp(我* np * (theta-OMEGAt))));
                        
                         fun_r = @ (R,θ,z) abs (Jr4) ^ 2 * R;
                        
                         fun_theta = @ (R,θ,z) abs (Jtheta4) ^ 2 * R;
                        
                         P_4 = (1 / sigmac) * integral3 (fun_r、R1、R2 0τ,z3、z4) + integral3 (fun_theta、R1、R2 0τ,z3、z4);
                        
                         P_5 = (1 / sigmac) * integral3 (fun_r 0 R1, 0,τ,z3、z4) + integral3 (fun_theta 0 R1, 0,τ,z3、z4);
                        
                         P_6 = (1 / sigmac) * integral3 (fun_r, R2, R3, 0,τ,z3、z4) + integral3 (fun_theta, R2, R3, 0,τ,z3、z4);

5个评论
显示4年长的评论隐藏4年长的评论

路易吉Stragapede 2020年7月13日

编辑:路易吉Stragapede 2020年7月13日

说实话,如果我考虑前面的代码我有这个错误:

                               未定义的函数或变量“z”。
                              
                               错误在Untitled5(37)行
                              
                               Xr4 = A4d * exp (gamma4 * z) + B4d * exp (-gamma4 * z);

因为这个原因我已经把Jr4 Jtheta4直接在函数处理fun_r和fun_theta。所以新代码是:

                               fun_r = @ (R,θ,z) abs(真实的(((A4d * exp (gamma4 * z) + B4d * exp (-gamma4 * z)) * (1 / R) * besselj (np, alfak * R) * exp(我* np * (theta-OMEGA * t))))) ^ 2 * R;
                              
                               fun_theta = @ (R,θ,z) abs(真实((我* (1 / np) * (A4d * exp (gamma4 * z) + B4d * exp (-gamma4 * z)) * (alfak * besselj (np-1 alfak * R) - (np / R) * besselj (np, alfak * R)) * exp(我* np * (theta-OMEGA * t))))) ^ 2 * R;
                              
                               P_4 = (1 / sigmac) * integral3 (fun_r、R1、R2 0τ,z3、z4) + integral3 (fun_theta、R1、R2 0τ,z3、z4);
                              
                               P_5 = (1 / sigmac) * integral3 (fun_r 0 R1, 0,τ,z3、z4) + integral3 (fun_theta 0 R1, 0,τ,z3、z4);
                              
                               P_6 = (1 / sigmac) * integral3 (fun_r, R2, R3, 0,τ,z3、z4) + integral3 (fun_theta, R2, R3, 0,τ,z3、z4);

现在的错误是:

                               错误使用/
                              
                               矩阵尺寸必须同意。
                              
                               错误在Untitled5 > @ (R,θ,z) abs(真实的(((A4d * exp (gamma4 * z) + B4d * exp (-gamma4 * z)) * (1 / R) * besselj (np, alfak * R) * exp(我* np * (theta-OMEGA * t))))) ^ 2 * R
                              
                               错误在integral3 > @ (y, z)有趣(x(1) *(大小(z)), y, z)(第138行)
                              
                               @ (y, z)有趣(x(1) *(大小(z)), y, z),…
                              
                               错误在integral2Calc > integral2t /张量(第228行)
                              
                               Z =乐趣(X, Y);NFE = NFE + 1;
                              
                               错误在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(48)行
                              
                               P_4 = (1 / sigmac) * integral3 (fun_r、R1、R2 0τ,z3、z4) + integral3 (fun_theta、R1、R2 0τ,z3、z4);