作业帮 > 数学 > 作业

matlab 一表达式求积分出来含有EllipticE(1/2*(-1/z^2)^(1/2))和csgn(z),

来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/09/14 22:39:05
matlab 一表达式求积分出来含有EllipticE(1/2*(-1/z^2)^(1/2))和csgn(z),
求r/(r^2+z^2)^(5/2)的积分过程如下:
>> syms sita r z d
>> S=r/(r^2+z^2)^(5/2);
>> s1=int(S,r,0,d*sin(sita))
s1 =
1/3*(-(z^2)^(3/2)+(d^2*sin(sita)^2+z^2)^(3/2))/(d^2*sin(sita)^2+z^2)^(3/2)/(z^2)^(3/2)
>> int(s1,sita,0,pi)
ans =
-1/3*(-pi*(d^2+z^2)^(3/2)+2*csgn(z)*z^3*((d^2+z^2)/z^2)^(1/2)*EllipticE((-d^2/z^2)^(1/2)))/(z^2)^(3/2)/(d^2+z^2)^(3/2)
那位好心人帮我看看其中的csgn(z)和EllipticE((-d^2/z^2)^(1/2)))这两个函数什么意识?有谁能指点我下,对工程进行简化!
其中的Z值是大于零的数!
对于"第二类椭圆积分 EllipticE()"的问题能麻烦您详细的说下吗!
同时其中Z肯定是大于零的 我如何知道EllipticE()具体的数值,
csgn( ) 称为"复符号函数",在实数中也有一个"符号函数" sgn( )
sgn( ) 的( )里是一个实数
csgn( ) 的( )里是一个复数
但sgn( )和csgn( ) 的返回值只有两个-1和1
这就是它们为什么都叫"符号"函数,也就是说,如果一个数你不知道它们的符号是正还是负,但是这个数的正负可以从sgn( )和csgn( ) 的返回值是+1还是-1就可以确定这个数的正负性.
sgn(x)的定义为:
当x>=0,sgn(x) = 1 ;
当x 0 或 ( Re(x) = 0 且 Im(x) > 0 )时,csgn (x)= 1
当Re(x) < 0 或 ( Re(x) = 0 且 Im(x)