这个式子帮算一算

szl

syms x t;
f=(-2*x^2+1)/(2*x^2-3*x+1)^2;
I=simple(int(f,x,cos(t),exp(-2*t)))

此为《高等应用数学问题的MATLAB求解（第2版）》（薛定宇,陈阳泉）例3-17
结果为“Warning: Explicit integral could not be found. ”
式子哪写错了？

winner245

直接用你的代码求解，很慢，最后得到的是一个piecewise function，因为它要判断上限比下限大。为了避免这个判断，可以先求不定积分，再用 subs：

1. syms x t;
   
2. f=(-2*x^2+1)/(2*x^2-3*x+1)^2;
   
3. y = int(f)
   
4. I = subs(f,exp(-2*t)) - subs(f, cos(t))

复制代码

y =

1/(x - 1) - 1/(2*x - 1)


I =

(2*cos(t)^2 - 1)/(2*cos(t)^2 - 3*cos(t) + 1)^2 - (2*exp(-4*t) - 1)/(2*exp(-4*t) - 3*exp(-2*t) + 1)^2

nwcwww

结果为“Warning: Explicit integral could not be found. ”
式子哪写错了？
-------------------

式子没错，不过作者成书期间MATLAB的符号计算引擎还是MAPLE, 近年来改成MUPAD了。
二者的差异还是不少的，很可能就是这个原因。

szl

谢谢二位高人。二楼的式子我要写在书边。
nwcwww 版主：我用的是2013版，下次装个2008版试试。

szl

2006b版下程序运行正常。儿子才学高数，让他接触MatLab 以期培养兴趣。入门级，书与程序匹配最重要，否则就一砖头拍死了。:lol