【求助】一个比较复杂的积分运算

magrog

syms h h2 A A2 alfa fai d r
>> h2=h+A*cos(alfa)*sin(fai);
A2=d^2+(1-sin(fai)^2*cos(alfa)^2)*A^2-2*d*A*cos(alfa)*sqrt(1-sin(fai)^2*cos(alfa)^2);
>> f=simple((h2^2*(h2^2+A2^2+r^2))/(pi*((h2^2+A2^2+r^2+2*A2*r)^3/2)*((h2^2+A2^2+r^2-2*A2*r)^3/2)))


>>a=int(f,-pi/2,pi/2)


??? 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 39

r = reshape(maple('map','int',f(:),x),size(f));

兄弟刚开始接触matlab，请各位先学多多帮助，谢谢！

sigma665

表达式太复杂，用数值积分

magrog

我给其中的参量赋值后进行积分，好像还是报错，正在计算中，谢谢ls

magrog

>> f=sqrt(1-0.25*cos(x)^2);
??? 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.sqrt at 7
Y = maple('map','sqrt',X);

magrog

>> clear all;
>> syms A x;
>> f=(1+A*cos(x))^2*((1+A*cos(x))^2+A^2*sin(x)^2+A*cos(x)*sqrt(1-0.25*cos(x)^2))/(((1+A*cos(x))^2+A^2*sin(x)^2+A*cos(x)*sqrt(1-0.25*cos(x)^2))^2-4*(A^2*sin(x)^2+Acos(x)*sqrt(1-0.25*cos(x)^2)))^3/2;
Warning: Could not find an exact (case-sensitive) match for 'Acos'. C:\Program Files\MATLAB704\toolbox\symbolic\@sym\acos.m is a case-insensitive match and will be used instead.  You can improve the performance of your code by using exact name matches and we therefore recommend that you update your usage accordingly.  Alternatively, you can disable this warning using warning('off','MATLAB:dispatcher:InexactMatch').
>> a=int(f,'x',-pi/2,pi/2)
Warning: Explicit integral could not be found.
> In sym.int at 58

a =

int(1/2*(1+A*cos(x))^2*((1+A*cos(x))^2+A^2*sin(x)^2+1/2*A*cos(x)*(4-cos(x)^2)^(1/2))/(((1+A*cos(x))^2+A^2*sin(x)^2+1/2*A*cos(x)*(4-cos(x)^2)^(1/2))^2-4*A^2*sin(x)^2-2*acos(x)*(4-cos(x)^2)^(1/2))^3,x = -1/2*pi .. 1/2*pi)

sigma665

用数值积分

help quad

magrog

??? 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))
我把参数赋值后积分，还是出错。不知道这个报错究竟是什么意思啊？麻烦了，能否解释一下阿？

messenger

不明白，你连数都没算出来呢，为什么着急reshape呢？

alljoyland

这个是说reshape 不能进行下去,reshape函数是一个变形矩阵的函数,
元素的 size 必须相等,
你 尝试 把 所有的 * -> .*  
/-> ./
^ -> .^
看看 是否会好

magrog

不是我要用reshape ，是运行后报错显示中有这个。我也搞不懂了，这个就是函数复杂一点，怎么就积分不出来呢？
f=(h+0.5*A*cos(x))^2*((h+0.5*A*cos(x))^2+d^2+A^2*(1-0.25*cos(x)^2)+2*d*A*sqrt(1-0.25*cos(x)^2)*cos(x)+r^2)/(((h+0.5*A*cos(x))^2+d^2+A^2*(1-0.25*cos(x)^2)+2*d*A*sqrt(1-0.25*cos(x)^2)*cos(x)+r^2)^2-4*(d^2+A^2*(1-0.25*cos(x)^2)+2*d*A*sqrt(1-0.25*cos(x)^2)*cos(x))*r^2)^3/2，就是这个函数，这里面的h=0.1,A=0.05,d=0.1,r=0.05,这样就只有一个x是变量了，可是还是积分不出来，老师出现上面的报错信息。各位，这到底是怎么回事请啊？

[ 本帖最后由 magrog 于 2008-8-26 19:49 编辑 ]

sigma665

说了n篇了，用数值积分

如果不知道什么是数值积分
请查看帮助里的quadl

messenger

试了一下，语法上没问题，就是一直在运算

你这个式子这么复杂，积不出来也正常

还是按版主说的，用数值积分吧

magrog

谢谢了，数值积分是不是就是用quad命令，不用int命令啊？