求大神指教如何用matlab求下面的积分。公式如图片所示,rt, rh为积分上下限。r为积分变量,rm,c2am,k 等都是常数。 QQ图片20170425112126.png 返回小木虫查看更多
n omega 是否是常数,rt、rh是符号还是数值。
使用MATLAB的定积分就可以了,给你假设了一定参数作了试算尝试了一下: clear all rt=5,rh=19,omg=80,k=50,rm=3,c2am=100,n=1; syms r f2=int(sqrt(c2am^2+2*omg*k*(r^(1-n)-rm^(1-n))+(1-n)*k^2/n*(1/r^(2*n)-1/rm^(2*n)))*2*pi*r,r,rt,rh) 结果为33600*pi。 你要计算其他参数时更换相应参数值即可。
直接用Mathematica吧,把公式敲进去就出来结果了
n omega 是否是常数,rt、rh是符号还是数值。
n omega, rt, rh 都是常数
,
rt,rh到最后都是确定的数值
使用MATLAB的定积分就可以了,给你假设了一定参数作了试算尝试了一下:
clear all
rt=5,rh=19,omg=80,k=50,rm=3,c2am=100,n=1;
syms r
f2=int(sqrt(c2am^2+2*omg*k*(r^(1-n)-rm^(1-n))+(1-n)*k^2/n*(1/r^(2*n)-1/rm^(2*n)))*2*pi*r,r,rt,rh)
结果为33600*pi。
你要计算其他参数时更换相应参数值即可。
直接用Mathematica吧,把公式敲进去就出来结果了