1.int

适用于能找到被积函数的原函数,定积分和不定积分皆可。

1.1 不定积分
syms x  
f=x^2;  
int(f,x) 

ans =
x^3/3

1.2 定积分
syms x  
f=x^2;  
int(x^2,x,0,1)

ans =
1/3

2.quad

不能找到被积函数的原函数使用int无法得到精确解,如:

syms x 
int(exp(cos(x)),x,0,1)

ans =
int(exp(cos(x)), x, 0, 1)

这时候可以使用积分的数值求解,使用quad命令数值求解。

f=@(x)exp(cos(x)) 
quad(f,0,1) 

ans =
2.3416

那么,如果是带参积分诸如怎么求呢?诸如:

y=∫0xesin(s)sdsy=\int_0^xe^{sin(s)}sdsy=0xesin(s)sds

想要绘制x在[0,1]范围时y-x曲线,这时候使用for循环转换为多次定积分来计算咯~当然上面的函数我们还是用数值求解quad命令:

% 被积函数定义
g=@(s)exp(sin(s)).*s; 
% 预定义x和y
x=linspace(0,1,21); 
y=zeros(1,21); 
% for循环对每个x的值进行定积分
for ii=1:length(x) 
	y(ii)=quad(g,0,x(ii)); 
end 
% 绘制结果
plot(x,y) 
title('y=\int_0^xe^{sin(s)}sds')

在这里插入图片描述

Logo

魔乐社区(Modelers.cn) 是一个中立、公益的人工智能社区,提供人工智能工具、模型、数据的托管、展示与应用协同服务,为人工智能开发及爱好者搭建开放的学习交流平台。社区通过理事会方式运作,由全产业链共同建设、共同运营、共同享有,推动国产AI生态繁荣发展。

更多推荐