Question

如何用matlab求该函数的积分

a=1,b=2,c=3; y=a/x^5/(exp(b/(c*x))-1) x[0,1]

另外，我用以下程序为什么不行
a=1,b=2,c=3;
f=sym('y=a/x^5/(exp(b/(c*x))-1)')
f1=subs(f1);
result=int(f1,'x',0,1)
问题应如何修改

Answer 1

1,
（1）先在命令窗口给s,a赋值
（2）>>
y=@(x)0.65*(s+x)*x/(a+x)
（3）[x,fval]=
fminbnd(y,0,inf)%用于求函数y在x在[0,inf]内的最小值
则在取负即为求最大值。不过本题我试验了，本题好像有问题。手算都没有最大值的，因为对y求导后y'在x>0内恒大于零，故y为增函数，所以应为在x取无穷大时，y最大，为无穷大。
2，同理
>>w=@(x)s+x-0.65*(s+x)*x/(a+x)
>>
[x,fval]=
fminbnd(w,0,inf)
再取负
好像也无解
你再看看题吧！是否有错！希望对你有所帮助！

Answer 2

>> a=1,b=2,c=3;
>> syms x
>> y=a/x.^5/(exp(b/(c.*x))-1)

y =

1/(x^5*(exp(2/(3*x)) - 1))

>> result=int(y,0,1)
Warning: Explicit integral could not be found.

result =

int(1/(x^5*(exp(2/(3*x)) - 1)), x = 0..1)此函数有错不能进行符号积分
>> f=@(x)a./x.^5/(exp(b./(c*x))-1)

f =

@(x)a./x.^5/(exp(b./(c*x))-1)

>> result=quadl(f,0,1)

Answer 3

>> a=1,b=2,c=3;
>> syms x
>> y=a/x.^5/(exp(b/(c.*x))-1)

y =

1/(x^5*(exp(2/(3*x)) - 1))

>> result=int(y,0,1)
Warning: Explicit integral could not be found.

result =

int(1/(x^5*(exp(2/(3*x)) - 1)), x = 0..1)此函数有错不能进行符号积分
>> f=@(x)a./x.^5/(exp(b./(c*x))-1)

f =

@(x)a./x.^5/(exp(b./(c*x))-1)

>> result=quadl(f,0,1)