- 在线时间
- 2 小时
- 最后登录
- 2016-1-31
- 注册时间
- 2016-1-22
- 听众数
- 12
- 收听数
- 1
- 能力
- 0 分
- 体力
- 17 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 8
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 10
- 主题
- 4
- 精华
- 0
- 分享
- 0
- 好友
- 3
升级   3.16% TA的每日心情 | 奋斗 2016-1-31 16:59 |
---|
签到天数: 2 天 [LV.1]初来乍到
 |
20体力
参数积分问题!求带有几个未知参数函数的积分,因为不能用数值积分的方法,但是int符号积分又给不出显示的解析式,积分结果还是用int、erfi表示的。但是我试过给几个参数先赋值再int积分是能得到积分值的,现在我想得到带有参数(alpha)的积分表达式。求助呀~~~谢谢了!!
以下是我完整的代码,麻烦大家帮我看一下,可以直接复制到matlab运算。式子有些复杂,不知道是不是积分不存在。
syms r phi z y alpha
k=0.012347^2;
t=239;
T0=38+273.5;
T1=70+273.5;
rho=1;
f1=(exp(-(sqrt(r^2+z^2)-y).^2./(4*alpha^2*t))-exp(-(sqrt(r^2+z^2)+y).^2./(4*alpha^2*t)))*(T0-T1);
f2=(exp(-(r*cos(phi)-y).^2./(4*k*t))+exp(-(r*cos(phi)+y).^2./(4*k*t)))*(T0-310.5*y);
intf1=int(f1,y,0,inf);
intf2=int(f2,y,0,inf);
T_1=1./(sqrt(4*pi*t)*alpha).*intf1+T1;
T_2=1./sqrt(4*pi*k*t).*intf2+310.5*r*cos(phi);
f3=(T_1+T_2-T0)*r;
T_3=int(int(int(f3,r,0,1),phi,0,pi/2),z,-0.35,0.35)
mT=rho*T_3
|
zan
|