MATLAB 定积分上限是变量a,如何求解目标函数关于变量a的最大值?谢谢

作者&投稿:左丘甄 (若有异议请与网页底部的电邮联系)
有关上限是变量的定积分函数及其作图,Matlab~

syms x Ef Et y;
Ei=0.5426;
Ef=0.842;
Et=linspace(1.1,1.4,100);
%%%%好像是9.11e-31
n=((4*pi*(2*1.18*9.11*10^31)^(3/2)/(4.14*10^(-15))).*(sqrt(x-1.1)/(1+2.71.^((x-Ef)./0.026))));
for i=1:length(Et)
i
y(i)=double(int(n,x,1.1,Et(i)));
end
plot(log(y),Et);

按你的方程积分不出来吧
syms x z
int(pi^2/(sqrt(1-x^2)*sqrt(1-(x-z)^2)),x,-1,z)
Warning: Explicit integral could not be found.

无约束极值问题,求原积分的最大值等效为求其相反数的极小值

>> fun=@(x)0.82*0.49*exp(6.72-47.7./x).*exp(-0.05*x); % 被积函数
>> f = @(a) integral(fun,0.01,a); % 定积分
>> obj = @(x) -f(x); % 定积分的相反数作为要求极小值的目标函数
>> [a,val] = fminsearch(obj,1) % 无约束极值

a =

  938.5500


val =

 -7.430055148537375e+02

最后结果表明在a=938.55时,f(a)取最大值743.0055148537375

实际上,根据被积函数的图像可知,这个定积分即下图与坐标轴围成的面积是递增的,只不过当a达到一定值后,变化较小了

例如,当a=938和a=1000时,f(a)的差异微乎其微:

>> f(968)-f(1000)

ans =

     1.298076313105412e-09



兴宁区17090795996: matlab 求定积分,积分的上限是变量 -
安泡小牛: 按你的方程积分不出来吧 syms x z int(pi^2/(sqrt(1-x^2)*sqrt(1-(x-z)^2)),x,-1,z) Warning: Explicit integral could not be found.

兴宁区17090795996: matlab 求定积分,积分的上限是变量,函数带有其他参数积分表达式为:y=1/(t+2(r - sqrt(r^2 - (m/2 - x)^2))),积分区间为(0,m),其中0 -
安泡小牛:[答案] >> syms t r m x >> y=1/(t+2*(r-sqrt(r^2-(m/2-x)^2))); >> f=int(y,0,m); f=simple(f); f 1/2*(-i*t*atanh((4*r^2+(-t*(t+4*r))^(1/2)*m)/((t+2*r)^2)^(1/2)/(-m^2+4*r^2)^(1/2))+2*t*atan(m/(t^2+4*t*r)^(1/2))-i*t*atanh((-4*r^2+(-t*(t+4*r))^(1/2)*m)/((t+2*r)^2)^(1/2)/(-m^2+4*r^2)^...

兴宁区17090795996: 在matlab中,如何积分上下限为变量 -
安泡小牛: >> syms x a r >> int('x^3*cos(a)','x')>> int(int('x^3*cos(a)','x','r-3','r+3'),'a','-0.21/r','0.21/r')ans =12.*r*(r^2+9.)*sin(.21/r)>> %直接使用符号积分,把b用r替换就可以了.

兴宁区17090795996: matlab怎么求参数为a式子的积分,Rr都已经定义 -
安泡小牛: 运行结果显示如下:Warning: Length of lower bounds is > length(x); ignoring extra bounds.> In checkbounds at 27 In optim\private\lsqncommon at 39 In lsqcurvefit at 165 In myfun_nihe_lianxi at 11??? Error using ==> optim\private\lsqncommon ...

兴宁区17090795996: matlab中quad积分,其中积分的上下限可以是变量吗? -
安泡小牛: 此乃正解>> syms x t; int(x*x+2*x,0,t) ans =(t^2*(t + 3))/3>> 此外 quad注释里面明确说明 Numerically evaluate integral, adaptive Simpson quadrature 只用于数值计算,不能用于符号积分 int就可以用符号计算

兴宁区17090795996: Matlab怎么计算定积分 -
安泡小牛: int(f,v,a,b) f是函数,v是变量(如是x可省略),a下限,b是上限.如syms x f=exp(x) int(f,0,1)

兴宁区17090795996: matlab积分问题,上限为变量. -
安泡小牛: k= 0.9092;syms x af=1./sqrt(1-k.^2*(sin(x).^2));int=int(f,x,a,pi/2);b=0:pi/100:pi/2;for n=1:length(b)c(n)=double(subs(int,a,b(n)));endplot(b,c)

兴宁区17090795996: 用matlab求解:已知定积分的值,求其上限 -
安泡小牛: syms x Q solve(int(1/100*(x/100)^4*exp(-x/100)/24,x,0,Q)-4/5,Q) ans =672.09787874865558455510741713836

兴宁区17090795996: matlab中:不定积分怎么写 -
安泡小牛: matlab中的不定积分(定积分)用 int()函数来表示. 使用格式: int(expr,var) %计算不定积分 int(expr,var,a,b) %计算定积分 expr——被积函数 var——积分变量 a——积分下限 b——积分上限 使用实例: 实例1, syms x int(-2*x/(1 + x^2)^2) 实例2, syms x z int(x/(1 + z^2), x) int(x/(1 + z^2), z) 实例3, syms a t int([exp(t), exp(a*t); sin(t), cos(t)])

本站内容来自于网友发表,不代表本站立场,仅表示其个人看法,不对其真实性、正确性、有效性作任何的担保
相关事宜请发邮件给我们
© 星空见康网