| 查看: 1008 | 回复: 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) 请问高手为什么会这样，请高手帮忙指点 |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有138人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】matlab求解一个简单的积分（急！！！）
已经有7人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab 积分的问题
已经有5人回复
- 【求助】用matlab编过有限差分法的请进
已经有9人回复
【答案】应助回帖
★ ★ ★
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
|
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