| 查看: 1459 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助matlab积分的问题
|
|||
|
我想对函数(f*v + (G*a*b)/(b*x + 1) + (M*c*k*m0*x)/((k*x*(c - 1) + 1)*(x - 1)^2))/(r0 - x)积分，于是编了如下的m文件： clear all; clc; syms a b G M m0 c k r0 f v x; I=(G*a*b/(b*x+1)+f*v+M*m0*c*k*x/(((x-1)^2)*(k*x*(c-1)+1)))/(r0-x) J=int(I,x) 但是结果出来和要积分的公式是一致的。结果是 I = (f*v + (G*a*b)/(b*x + 1) + (M*c*k*m0*x)/((k*x*(c - 1) + 1)*(x - 1)^2))/(r0 - x) 请问高手为什么会这样，请高手帮忙指点 |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有267人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】matlab求解一个简单的积分（急！！！）
已经有7人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab 积分的问题
已经有5人回复
- 【求助】用matlab编过有限差分法的请进
已经有9人回复
|
J = G*a*b/(b*r0+1)*log(b*x+1)-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v-M*m0*c*k/(1+k*c-k)/(-1+r0)/(x-1)-2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^2*k-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*r0*f*v*b+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*r0*f*v*k-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*r0*f*v*k*c+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*r0*f*v-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*M*m0*c*k*r0-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^2+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*k*c*r0+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*k*r0+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^2*k*c-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^3*k*c+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^3*k-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*M*m0*c*k*r0^2*b-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^2-2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^2*k+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*b*r0^2+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*k-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*b+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*b*r0^2*k-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*b*r0^2*k*c+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^2*k*c-2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*k*b-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*k*c+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*k*c*b-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^4*k*c*b+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^4*k*b-M*m0*c*k^3/(2*k^2*r0-k*r0-4*k^2*c*r0+1+k^2-2*k+2*k^2*c^2*r0-k^3*r0+2*k*c+k^2*c^2+3*k^3*c*r0-3*k^3*c^2*r0-2*k^2*c+k*c*r0+k^3*c^3*r0)/(k*c-k)*log(k*x*c-k*x+1)-M*m0*c^3*k^3/(2*k^2*r0-k*r0-4*k^2*c*r0+1+k^2-2*k+2*k^2*c^2*r0-k^3*r0+2*k*c+k^2*c^2+3*k^3*c*r0-3*k^3*c^2*r0-2*k^2*c+k*c*r0+k^3*c^3*r0)/(k*c-k)*log(k*x*c-k*x+1)+2*M*m0*c^2*k^3/(2*k^2*r0-k*r0-4*k^2*c*r0+1+k^2-2*k+2*k^2*c^2*r0-k^3*r0+2*k*c+k^2*c^2+3*k^3*c*r0-3*k^3*c^2*r0-2*k^2*c+k*c*r0+k^3*c^3*r0)/(k*c-k)*log(k*x*c-k*x+1)+M*m0*c^2*k^2/(1+k^2-2*k+2*k*c+k^2*c^2-2*k^2*c)/(-1+r0)^2*log(x-1)-M*m0*c*k^2/(1+k^2-2*k+2*k*c+k^2*c^2-2*k^2*c)/(-1+r0)^2*log(x-1)+M*m0*c*k/(1+k^2-2*k+2*k*c+k^2*c^2-2*k^2*c)/(-1+r0)^2*log(x-1)*r0  |
5楼2012-11-08 16:16:38
【答案】应助回帖
★ ★ ★
likuihei(金币+5): 2011-06-17 12:48:59
sunyang1988(金币+3): 谢谢交流 2011-06-17 18:27:31
likuihei(金币+5): 2011-06-17 12:48:59
sunyang1988(金币+3): 谢谢交流 2011-06-17 18:27:31
|
I=(G*a*b/(b*x+1)+f*v+M*m0*c*k*x/(((x-1)^2)*(k*x*(c-1)+1)))/(r0-x) 但是积分完成后是“J=***”， J =G*a*b/(b*r0+1)*log(b*x+1)-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^2*k*c-M*m0*c*k/(1+k*c-k)/(-1+r0)/(x-1)-M*m0*c*k^3/(2*k^2*r0-k*r0+1-2*k+2*k*c+2*k^2*c^2*r0+k^2-k^3*r0-4*k^2*c*r0-2*k^2*c+k^2*c^2+3*k^3*c*r0+k*c*r0+k^3*c^3*r0-3*k^3*c^2*r0)/(k*c-k)*log(k*x*c-k*x+1)+2*M*m0*c^2*k^3/(2*k^2*r0-k*r0+1-2*k+2*k*c+2*k^2*c^2*r0+k^2-k^3*r0-4*k^2*c*r0-2*k^2*c+k^2*c^2+3*k^3*c*r0+k*c*r0+k^3*c^3*r0-3*k^3*c^2*r0)/(k*c-k)*log(k*x*c-k*x+1)-M*m0*c^3*k^3/(2*k^2*r0-k*r0+1-2*k+2*k*c+2*k^2*c^2*r0+k^2-k^3*r0-4*k^2*c*r0-2*k^2*c+k^2*c^2+3*k^3*c*r0+k*c*r0+k^3*c^3*r0-3*k^3*c^2*r0)/(k*c-k)*log(k*x*c-k*x+1)+M*m0*c^2*k^2/(1-2*k+2*k*c+k^2-2*k^2*c+k^2*c^2)/(-1+r0)^2*log(x-1)-M*m0*c*k^2/(1-2*k+2*k*c+k^2-2*k^2*c+k^2*c^2)/(-1+r0)^2*log(x-1)+M*m0*c*k/(1-2*k+2*k*c+k^2-2*k^2*c+k^2*c^2)/(-1+r0)^2*log(x-1)*r0-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*r0*f*v*b+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*r0*f*v*k-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*r0*f*v*k*c+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*r0*f*v-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*M*m0*c*k*r0-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^2-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*M*m0*c*k*r0^2*b-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^3*k*c-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*k*c*r0+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*k*r0-2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^2*k+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^2*k*c+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0^3*k+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*G*a*b*r0+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*b*r0^2-2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^2*k-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*k*c+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^4*k*b-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^2-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*b*r0^2*k*c+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*b*r0^2*k+2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*k*c*b+1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*k-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*b-2/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^3*k*b-1/(b*r0+1)/(-k*r0+1+k*c*r0)/(-1+r0)^2*log(-r0+x)*f*v*r0^4*k*c*b |
2楼2011-06-16 10:00:03
3楼2011-06-17 12:49:39
★ ★ ★
sunyang1988(金币+3): 谢谢交流 2011-07-21 00:14:34
sunyang1988(金币+3): 谢谢交流 2011-07-21 00:14:34
|
J = -1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v+G*a*b/(b*r0+1)*log(b*x+1)-M*m0*c*k/(1-k+k*c)/(r0-1)/(x-1)+2*M*m0*c^2*k^3/(2*k^2*r0+k^2*c^2-k^3*r0-2*k-k*r0+2*k^2*c^2*r0+2*k*c-3*k^3*c^2*r0-4*k^2*c*r0-2*k^2*c+k^2+k^3*c^3*r0+k*c*r0+1+3*k^3*c*r0)/(k*c-k)*log(k*x*c-k*x+1)-M*m0*c^3*k^3/(2*k^2*r0+k^2*c^2-k^3*r0-2*k-k*r0+2*k^2*c^2*r0+2*k*c-3*k^3*c^2*r0-4*k^2*c*r0-2*k^2*c+k^2+k^3*c^3*r0+k*c*r0+1+3*k^3*c*r0)/(k*c-k)*log(k*x*c-k*x+1)-M*m0*c*k^3/(2*k^2*r0+k^2*c^2-k^3*r0-2*k-k*r0+2*k^2*c^2*r0+2*k*c-3*k^3*c^2*r0-4*k^2*c*r0-2*k^2*c+k^2+k^3*c^3*r0+k*c*r0+1+3*k^3*c*r0)/(k*c-k)*log(k*x*c-k*x+1)+M*m0*c^2*k^2/(k^2*c^2-2*k+2*k*c-2*k^2*c+k^2+1)/(r0-1)^2*log(x-1)-M*m0*c*k^2/(k^2*c^2-2*k+2*k*c-2*k^2*c+k^2+1)/(r0-1)^2*log(x-1)+M*m0*c*k/(k^2*c^2-2*k+2*k*c-2*k^2*c+k^2+1)/(r0-1)^2*log(x-1)*r0+2/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*b*r0^2-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*G*a*b-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*r0*f*v*k*c-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*r0*f*v*b+1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*r0*f*v*k+2/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*r0*f*v-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*M*m0*c*k*r0-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*G*a*b*r0^2+2/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*G*a*b*r0-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*G*a*b*k*c*r0+1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*G*a*b*k*r0-2/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*G*a*b*r0^2*k+2/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*G*a*b*r0^2*k*c-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*G*a*b*r0^3*k*c+1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*G*a*b*r0^3*k-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*M*m0*c*k*r0^2*b-2/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^2*k-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^3*b+1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^3*k-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^2-2/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^3*k*b+1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^4*b*k+1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*b*r0^2*k-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*b*r0^2*k*c+2/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^2*k*c+2/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^3*k*c*b-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^3*k*c-1/(-k*r0+1+k*c*r0)/(b*r0+1)/(r0-1)^2*log(-r0+x)*f*v*r0^4*b*k*c |
4楼2011-07-16 20:38:11
