matlab积分

liuxb

clc;
clear all;
syms w z r a  x k
a=17.5;
r=17.5;
%b=sqrt(r^2-(r-a)^2);
x=66;
k=0.886;
h=2*sqrt(r^2-w^2);
f2=sqrt(2*h/(3.14*z)*tan(3.14*z/(2*h)))*(0.923+0.199*(1-sin(3.14*z/(2*h)))^4)/cos(3.14*z/(2*h));
f11=(1.122-0.561*(z/h)+0.085*((z/h)^2)+0.018*((z/h)^3))/sqrt(1-(z/h));
%被积函数
y=32*x^2*(r^2-w^2)*z*f2^2+2*r^4*k^4*z*f11^2;
%积分上限
fh=sqrt(r^2-w^2)-(r-a);
%积分结果
F1=int(y,z,0,fh)
F2=eval(int(F1,w,-r,r))
我要用matalb将二重积分换为这样的一重积分计算，程序如上面，但是我的二重积分能够计算，这个一重积分却算不了，提示Error using mupadmex
Error in MuPAD command: Out of memory.
Error in sym/int (line 124)
   rSym = mupadmex('symobj::intdef',f.s,x.s,a.s,b.s,options);
Error in Untitled3 (line 17)
F1=int(y,z,0,fh)

求各位高手指导，这样的二重积分是不是不能改为一重积分？

winner245

int 求解的是符号积分，也就是闭式解（closed-form，或解析解）。通常用 int 函数求不出解有以下几种情形：

1. 积分有闭式解，但 mupad 的符号计算并不完善而求不出解。此时，可以考虑使用更加强大的符号计算软件，如 maple 或 mathematica
2. 积分本身无闭式解，那么，无论用什么科学计算软件都难以求出。此时，可以考虑数值积分求解，利用 quad2d 或 integral2 函数，一般都能得到数值解

shanshui281

新年快乐！！