神马品牌网哈尔滨城市通充值地址:急询：用Matlab对phai^(k-1)*cos(phai)的积分程序？感谢
来源：百度文库 编辑：神马品牌网 时间：2024/04/20 07:52:03
急询：用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))
不过积分下限被零除了。