matlab求积分方程
来源:学生作业帮 编辑:神马作文网作业帮 分类:综合作业 时间:2024/11/10 14:19:23
matlab求积分方程
已知一个积分方程,积分下限为a,上限为3,积分号内部的函数为:(e的a次方减a再减去1)的负二分之一次方,右边为一个常数,假设它为1
那么,这个积分方程的球阀是不是应该是
>> f=@(a)(2.72.^a-a-1).^-0.5;%被积函数
>> b=3;%上限
>> val=1;%定积分值
>> g=@(a)quadl(f,a,b)-val;
>> a=fzero(g,b)
可是这么算的话,会有
Warning:Minimum step size reached; singularity possible.
> In quadl at 95
In @(a)quadl(f,a,b)-val
In fzero at 282
的警告
并且只能把右边的常数最大设为2.7009,2.71就会没有答案
但是如果用定积分反推回去的话,到2.9999都是有答案的,这个是怎么回事啊?
1楼的答案,我想问的是怎么样能正确的求出a的数值啊,就是那个积分方程的a要得到一个具体的数的啊?
已知一个积分方程,积分下限为a,上限为3,积分号内部的函数为:(e的a次方减a再减去1)的负二分之一次方,右边为一个常数,假设它为1
那么,这个积分方程的球阀是不是应该是
>> f=@(a)(2.72.^a-a-1).^-0.5;%被积函数
>> b=3;%上限
>> val=1;%定积分值
>> g=@(a)quadl(f,a,b)-val;
>> a=fzero(g,b)
可是这么算的话,会有
Warning:Minimum step size reached; singularity possible.
> In quadl at 95
In @(a)quadl(f,a,b)-val
In fzero at 282
的警告
并且只能把右边的常数最大设为2.7009,2.71就会没有答案
但是如果用定积分反推回去的话,到2.9999都是有答案的,这个是怎么回事啊?
1楼的答案,我想问的是怎么样能正确的求出a的数值啊,就是那个积分方程的a要得到一个具体的数的啊?
syms a
f=((2.72.^a-a-1).^(-1/2)-1);
g=int(f,a,a,3)
结果:
g =
int(1/((68/25)^a-a-1)^(1/2)-1,a = a ..3)
f=((2.72.^a-a-1).^(-1/2)-1);
g=int(f,a,a,3)
结果:
g =
int(1/((68/25)^a-a-1)^(1/2)-1,a = a ..3)