|
楼主 |
发表于 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(:),[x.s '=(' a.s ')..(' b.s ')']),size(f));
程序及错误如上
我还发现一个小问题,如果写一个简单的a^2+b^2这样的积分,本来能求出,可是运行完上面的程序,出现错误以后,clc clear 然后再运行这个简单的积分,他也会出现上面的错误,这是为什么,对了,这是用7.1解得,7.8运行一天解不出[
这是对一个表达式求导得到的,可以求导,不可以积分? |
|