bubuluandut
2011-07-31, 16:31
各位前辈们，小弟刚接触matlab，要做符号运算，做一个积分。积分如下：
syms rho g A h x y z x0 y0 z0 k kx ky w t theta
kx=k*cos(theta);
ky=k*sin(theta);
p=rho*g*A*(exp(k*(z+h))+exp(-k*(z+h)))/(exp(k*h)+exp(-k*h))*cos(kx*x+ky*y-w*t);
My=int(int(p*(x-x0)*sqrt(diff(y,x)^2+1+diff(y,z)^2),x,(1-cos(45/180*pi))/sin(45/180*pi)*(z-z0)+x0-0.3,(1-cos(45/180*pi))/sin(45/180*pi)*(z-z0)+x0-0.225),z,z0,z0+0.225*sin(45/180*pi))
要做的积分就是上边那个My，求其他积分还可以，但是这个积分运行后即显示busy，没有停止也没有出错，还没有结果，一直busy，请问是为什么
内容是一个求压力的积分，压强函数是p
有劳各位了
syms rho g A h x y z x0 y0 z0 k kx ky w t theta
kx=k*cos(theta);
ky=k*sin(theta);
p=rho*g*A*(exp(k*(z+h))+exp(-k*(z+h)))/(exp(k*h)+exp(-k*h))*cos(kx*x+ky*y-w*t);
My=int(int(p*(x-x0)*sqrt(diff(y,x)^2+1+diff(y,z)^2),x,(1-cos(45/180*pi))/sin(45/180*pi)*(z-z0)+x0-0.3,(1-cos(45/180*pi))/sin(45/180*pi)*(z-z0)+x0-0.225),z,z0,z0+0.225*sin(45/180*pi))
要做的积分就是上边那个My，求其他积分还可以，但是这个积分运行后即显示busy，没有停止也没有出错，还没有结果，一直busy，请问是为什么
内容是一个求压力的积分，压强函数是p
有劳各位了