如何用matlab求积分？

~ 于无法求得exp(x^2)的原函数，我们只能用数值算法来求解，可以用复化梯形公式、Romberg公式、Gauss公式等，有好多种。我用Matlab编了一个用Gauss公式求解积分的函数。

function S=GaussIntegrate()
%运用Gauss求积公式计算数值积分

%f为被积函数，Rho为权函数，二者均为符号函数
x=sym('x');
f=exp(x^2);
Rho=1;

%a,b分别为求积区间的左界和右界
a=1;
b=2;

%n表示求积结点的个数，是一正整数
n=8;

%本程序利用线性变换将区间[a,b]变换到[-1.1]，
%同时令g=f*Rho为被积函数，然后利用
%古典的Gauss求积公式进行计算，此时直交多项式即为Legendre多项式

if n=0||n~=floor(n)
error('错误，n必须是一个非负整数！');
end;

if a>b
error('错误，区间的左界a一定不大于右界b！');
end;

%计算n次Legendre多项式
syms x;
P=1/(2^n*factorial(n))*diff((x^2-1)^n,n);
w=roots(sym2poly(P));

%计算数值积分
A=zeros(1,n);
S=0;
for k=1:n
A(k)=2/((1-w(k)^2)*(subs(diff(P),w(k))^2));

t=a+(b-a)/2*(w(k)+1);
g=(b-a)/2*subs(f*Rho,t);
S=S+A(k)*g;
end;

--------------------------------
我取了8个结点，计算精度就已经达到了小数点后8位，效率还是很高的。
注意：由于Matlab调用Maple的符号计算工具箱，第一次运行时会加载一小会，耐心等待。
以后再运行速度就很快了。

一、符号积分
符号积分由函数int来实现。该函数的一般调用格式为：
int(s)：没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分；
int(s,v)：以v为自变量，对被积函数或符号表达式s求不定积分；
int(s,v,a,b)：求定积分运算。a,b分别表示定积分的下限和上限。该函数求被积函数在区间[a,b]上的定积分。a和b可以是两个具体的数，也可以是一个符号表达式，还可以是无穷(inf)。当函数f关于变量x在闭区间[a,b]上可积时，函数返回一个定积分结果。当a,b中有一个是inf时，函数返回一个广义积分。当a,b中有一个符号表达式时，函数返回一个符号函数。

2.数值积分的实现方法
基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数的调用格式为：
[I,n]=quad('fname',a,b,tol,trace)
基于变步长、牛顿－柯特斯(Newton-Cotes)法，MATLAB给出了quadl函数来求定积分。该函数的调用格式为：
[I,n]=quadl('fname',a,b,tol,trace)
其中fname是被积函数名。a和b分别是定积分的下限和上限。tol用来控制积分精度，缺省时取tol=0.001。trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。返回参数I即定积分值，n为被积函数的调用次数。

---

兴庆区19783037645： 用MATLAB算积分 - ？
南费十味： 使用MATLAB软件,可以用int()计算不定积分或定积分.计算方法如下:syms x int((log(x))^2/x) %这里 ln(x) 用 log(x) 来表示 计算结果 log(x)^3/3 %ln³(x)/3

兴庆区19783037645： 如何用matlab求积分的程序 - ？
南费十味： 程序: fun=sin(0.5*pi*x./y);%% a=int(int(fun,y,sqrt(x),x),x,1,2); b=simple(a) %化简 I=vpa(b,4) %得到4位近似解,也可以任意N位解 结果: I = 0.2719

兴庆区19783037645： 在matlab中怎样求带有参数的定积分 - ？
南费十味：[答案] 设被积函数为f=ax+b,积分区间为(c,d) MATLAB代码如下: syms a b c d x; f=a*x+b; answer=int(f,x,c,d).

兴庆区19783037645： 用matlab定积分怎么求 - ？
南费十味： int(s,v,a,b):以v为自变量,对被积函数s在区间[a,b]上的定积分.a和b可以是两个具体的数,也可以是一个符号表达式,还可以是无穷(inf).例:syms x %定义符号变量 f=int(x^3,x,1,2) % 求x^3在区间[1,2]上的定积分

兴庆区19783037645： 怎么用MATLAB求积分 - ？
南费十味： 基本不管多复杂的式子,MATLAB都能解决,只要你输入正确,解出来只是时间问题,我给你一个解题的参考格式.设f为你的被积函数,只能自己输,别输错;x为积分变量,a,b分别为积分下限和上限.MATLAB中计算函数f的不定积分的格式如下: syms x;int(f,x) 计算定积分的格式如下: syms a b x;int(f,x,a,b).

兴庆区19783037645： 怎样用matlab编程求解积分 - ？
南费十味： int函数 例如s = sin(x),求积分int('s'),即可

兴庆区19783037645： matlab中怎么求函数的积分 - ？
南费十味： 用 quad 命令: 具体用法如下:1 2y = @(x)(sin(x)); quad(y,0,pi)

兴庆区19783037645： matlab 如何求积分？
南费十味： quad('函数名',a,b)我这里的a和b指的是积分下限和上限

兴庆区19783037645： 怎么用matlab算广义积分? - ？
南费十味： 用matlab算广义积分方法:%符号积分 int(f,v) int(f,v,a,b)%数值积分 trapz(x,y)%梯形法沿列方向求函数Y关于自变量X的积分 cumtrapz(x,y)%梯形法沿列方向求函数Y关于自变量X的累计积分 quad(fun,a,b,tol)%采用递推自适应Simpson法计算积分 quad1(fun,a,b,tol)%采用递推自适应Lobatto法求数值积分 dbquad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)%二重(闭型)数值积分指令 triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)%三重(闭型)数值积分指令 .