hudu11111 发表于 2009-11-4 11:44

关于积分

请教各位
一个有十个变量的复杂表达式,对其中的一个变量x10求积分,最后表达式含有其他九个变量,怎么做,int(f,x10,0,1.5)这样可以吗?会不会求不出,因为我等了很长时间,仍然busy,所以怕是变量太多解不出,各位提提意见
谢谢了

VibrationMaster 发表于 2009-11-4 19:31

不好办。第10个变量和前面9个变量的耦合有没有规律。如果有规律自己先凑出点,然后再积。否则真的计算机要累死

hudu11111 发表于 2009-11-5 10:33

回复 沙发 VibrationMaster 的帖子

为什么这么做有错误?

syms x1 x2 x3 x4 x5 x6 x7 x8 x9 x10;
a=(-x1)+x3*cos((4*pi/3)*x10);
b=x2+x3*sin((4*pi/3)*x10);
e=2*atan((1-((x4^2+a^2+b^2-x5^2)/(2*x4*(a^2+b^2)^(1/2)))*((-a)/((a^2+b^2)^(1/2)))-((1-((x4^2+a^2+b^2-x5^2)/(2*x4*(a^2+b^2)^(1/2)))^2)^(1/2))*(b/(a^2+b^2)^(1/2)))/(((1-((x4^2+a^2+b^2-x5^2)/(2*x4*(a^2+b^2)^(1/2)))^2)^(1/2))*((-a)/((a^2+b^2)^(1/2)))-((x4^2+a^2+b^2-x5^2)/(2*x4*(a^2+b^2)^(1/2)))*(b/(a^2+b^2)^(1/2))));
c=a+x4*cos(e);
d=b+x4*sin(e);
m=a+x7*cos(e+x9);
n=b+x7*sin(e+x9);
q=asin((m+x8)/x6);
p=n+x6*cos(q);
f=abs(diff(p,x10));
>> s=int(f,x10,0,1.5)
??? Error using ==> reshape
To RESHAPE the number of elements must not change.

Error in ==> sym.maple at 94
      result = reshape(result,size(varargin{3}));

Error in ==> sym.int at 51
   r = reshape(maple('map','int',f(:),),size(f));
程序及错误如上

我还发现一个小问题,如果写一个简单的a^2+b^2这样的积分,本来能求出,可是运行完上面的程序,出现错误以后,clc clear 然后再运行这个简单的积分,他也会出现上面的错误,这是为什么,对了,这是用7.1解得,7.8运行一天解不出[
这是对一个表达式求导得到的,可以求导,不可以积分?

beyondhxf 发表于 2009-11-5 11:29

http://forum.vibunion.com/forum/thread-87767-1-1.html

hudu11111 发表于 2009-11-7 14:59

回复 地板 beyondhxf 的帖子

看了看,好像解决不了问题呀?请多指教

hudu11111 发表于 2009-11-13 09:38

自己顶!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
页: [1]
查看完整版本: 关于积分