Question

MATLAB求积分？

这个函数怎么在MATLAB里求不定积分，k^2=0.0521

Answer 1

题主给出的积分问题，应该是一个特殊函数的定积分，而不是不定积分问题。我们先假定φ的变化范围从0到π/2，则可以用quad（）函数求得其数值解，然后用plot（）函数绘制φ—f（φ）关系曲线图。

实现过程：

k2=0.0521;

fun=@(phi)1./sqrt(1+k2*sin(phi).^2);

phi0=0:pi/10:pi/2;

for i=1:length(phi0)

f=quad(fun,0,phi0(i))

f0(i)=f;

end

plot(phi0,f0)

xlabel('φ'),ylabel('f（φ）')

追问

我想求不定积分，将v用关于phi的函数表示出来

)

追答

>> syms a k rb phi
>> v=-int(a/sqrt(1+k^2*sin(phi)^2)-rb,phi)
v =phi*rb - a*ellipticF(phi, -k^2) %ellipticF（） 第一类不完全椭圆积分

Answer 2

这个表达式里面被积函数应该是φ，你应该在Matlab里把它定义成符号变量。然后再使用积分命令进行符号运算。