解方程组(方程组中含有二重积分,并且积分的上下限中含有变量)

问题来源:http://www.ilovematlab.cn/thread-57006-1-3.html

问题描述:http://hi.csdn.net/attachment/200912/20/37536_1261296414KS8i.png

主要问题就是方程组中含有二重积分,并且积分的上下限中含有变量,所以我现在是无从下手,望高手指点

Forcal代码:

(::p,q,m,C1,C2,C3,C4,k,g,Txing,T2)= { p=0.020,q=0.219,m=10369.6,C1=800,C2=2,C3=6,C4=8,k=3,g=4,Txing=8,T2=12 }; t_T2(u:a:p,q,m)= a=exp[-(p+q)*u],m*p*(p+q)^2*a/(p+q*a)^2; t2_T2(t::h_t_T2,T2,p,q,g)= if[t==T2,return(0)], XSLSF::fpqg[h_t_T2,t,T2,1e-6]/[g*(p+q)]; f(t2,T1,y1,y2:a,b:h_t_T2,h_t2_T2,p,q,m,C1,C2,C3,C4,k,g,Txing,T2)= //函数定义 { a=exp[-(p+q)*Txing], b=m*p*(p+q)^2*a/(p+q*a)^2, y1=-C2*XSLSF::fpqg[h_t2_T2,t2,T2,1e-6]+C4*b*[1-k*(p+q)] +C3*{k*(p+q)*XSLSF::fpqg[h_t_T2,T1,Txing,1e-6]+b*k*(p+q)*t2-k*(p+q)*Txing}, y2=C2*T1*T1/[2*g*(p+q)]-k*(p+q)*C3*(t2-T1)-C4*[1-k*(p+q)] }; main(:x,i,t2,T1:h_t_T2,h_t2_T2)= { h_t_T2=HFor("t_T2"),h_t2_T2=HFor("t2_T2"), x=new[rtoi(real_s),rtoi(4),rtoi(EndType),1.0,1.0], i=XSLSF::netn[HFor("f"),1e-6,0.1,0.1,x,100], x.XSLSF::getra[0,&t2,&T1], printff{"/r/n实际迭代次数={1,i}, t2={2,r}, T1={3,r}/r/n",i,t2,T1}, delete[x] };

结果:


实际迭代次数=4, t2=5.6891908532692748, T1=3.4017560817759569

当然,选择其他的初始参数时,还会有其他解。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值