matlab求解定积分 定积分在matlab中怎么表示

matlab说他也不会做.如果你只是想求解这个积分的值,可按如下方法做1. 首先看看被积函数是什么形态,有没有奇点 ezplot('(x^2*cos(x)^2 + 1)^(1/2)'); 你这个函数是连续的,所以可积,然后2. 用简单的数值积分 step = 1e-2; x = 0:step:1; y = (x.^2.*cos(x).^2 + 1).^(1/2); step*sum(y) % 积分值

分成两段积 ∫[0,1]cos(wt)t dt+∫[1,2]cos(wt)(2-t) dt syms x w int(cos(w*x)*x,x,0,1)+int(cos(w*x)*(2-x),x,1,2) ans =(cos(w)+w*sin(w)-1)/w^2-(cos(2*w)+w*sin(w)-cos(w))/w^2

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]上的定积分

int(函数f(x),a,b) 计算函数f(x)在区间[a,b]的定积分 quadl:该函数使用了称为洛巴托求积(Lobatto Quadrature)的算法,对于高精度和光滑函数效率更高 I=quadI(func,a,b,tol); func是被积函数,a,b是积分限,tot是期望的绝对误差(如果不提供,默认为1e-6) eps是一个浮点误差值,如果不带参数,它的值为1到下一个能表示的比1大的浮点数之间差距的,它用来表示浮点相对误差精度 理解为机器最小数即可 I=quadl(@(t)(t-3*t.^2+2*t.^3).^(-1/3),eps,1/2)表示积分限在[0,0.5]

一个语句搞定(符号运算):>> int('x/sqrt(16+x^2)',0,3) ans = 1 或者用数值方法计算:>> quad(inline('x./sqrt(16+x.^2)'),0,3) ans = 0.999999994947221

一、符号积分 符号积分由函数int来实现.该函数的一般调用格式为:int(s):没有指定积分变量和积分阶数时,系统按findsym函数指示的默认变量对被积函数或符号表达式.

用matlab的int()函数可以快速算相对简单的定积分.使用格式:int(S,v,a,b) S——被积函数 v——积分变量 a,b——积分区间 syms x F = int(cos(x)/sqrt(1 + x^2), x, 0, 10); vpa(F, 5) ans =0.37571

syms t s k=exp(-t/21)/(21*(1-exp(365/21)))*(1.1*365^(5/2)*(365-t)^(-3/2)-1.1*365^(5/2)*(365-t+s)^(-3/2)-1.65*s); J=int(k,t,0,365)

clearsyms t2y2=-3.109*10^20* t2.^9+7.242*10^18* t2.^8-5.806*10^16* t2.^7+1.097*10^14* t2.^6+1.01*10^12* t2.^5-6.515*10^9 *t2.^4+1.447*10^7* t2.^3-13370 *t2.^2+8.83*t2+0.00003053;int(y2,t2,0,0.006126)

for N=10:10:100 x=0:2*pi/(N-1):2*pi; y=sin(x); z(N/10)=0.5*sum((y(1:end-1)+y(2:end)).*2*pi/(N-1)); end N=10:10:100; plot(N,z,'-o',N,N*0,'*')