衡阳邮政编码:急询:用Matlab对phai^(k-1)*cos(phai)的积分程序?感谢!
来源:百度文库 编辑:神马品牌网 时间:2024/08/26 18:41:51
急询:用Matlab对phai^(k-1)*cos(phai)的积分程序?感谢
syms phai k t
y=int(phai^(k-1)*cos(phai),phai,0,t)
结果如下:
Warning: Explicit integral could not be found.
> In D:\EDS\MATLAB6p5\toolbox\symbolic\@sym\int.m at line 58
y =
int(phai^(-1+k)*cos(phai),phai = 0 .. t)
怎么没有具体的运行结果啊?
请问高手,是不是程序有问题,请给出新程序。
谢谢
syms phai k t
y=int(phai^(k-1)*cos(phai),phai,0,t)
结果如下:
Warning: Explicit integral could not be found.
> In D:\EDS\MATLAB6p5\toolbox\symbolic\@sym\int.m at line 58
y =
int(phai^(-1+k)*cos(phai),phai = 0 .. t)
怎么没有具体的运行结果啊?
请问高手,是不是程序有问题,请给出新程序。
谢谢
遇到这种情况,先做不定积分,再将积分上下限代入,如果还不行,就用maple了。
>> y=int(phai^(k-1)*cos(phai),phai)
结果如下:
y =
2^(k-1)*pi^(1/2)*(3/pi^(1/2)/k*phai^(k-1)*2^(-k)*(4+2*k+2*phai^2)/(6+3*k)*sin(phai)+2^(1-k)/pi^(1/2)/k*phai^(k-1)*(cos(phai)*phai-sin(phai))-3/pi^(1/2)/k*phai^(1/2)*2^(1-k)/(6+3*k)*LommelS1(k+3/2,3/2,phai)*sin(phai)-1/pi^(1/2)/k/phai^(1/2)*2^(1-k)*(cos(phai)*phai-sin(phai))*LommelS1(k+1/2,1/2,phai))
不过积分下限被零除了。