f=(x)sin(x)+cos(x);定义被积函数Intf=quad(f,0,1)%对f进行积分,下限为0,上限为1www.mh456.com防采集。
其实,你仔细观察一下会发现,用第二种方法得到百xk表达式之后再subs代入的结果有这样的特点:
syms s f=1/250*(1-s/188.5+1/(2*pi)*sin(4/377*pi*s));a1=vpa(int(f,s,0,1),10) a2=int(f,s);a3=subs(a2,1)-subs(a2,0) a1和a3是基本一样的
1、大部分项的值都很小,量级在10^(-16),可以度看作0;
3、如果使用符号运算,不同版本下可能遇到不同的现象。在6.5版的结果: syms x z a=int(int(exp(-z)/z,((x-10)^2+10)/4e-6,inf),x,-10,10);Warning:Explicit integral could not be found. In D:\\MATLAB6
2、对应于第一种方法得到的非零项的结果是Inf或NaN。
这么复杂的问题。你分开问或许还有人答。数值方法: 解非线性方程(组)用fsolve 解一阶常微分方程(组)用ode系列函数 符号计算方法: 解非线性方程(组)用solve 解常微分方程(组)用dsolve
符号 syms x;int(exp(2*x),x,0,1) ans= exp(2)/2-1/2 数值 f=(x)exp(2*x);quad(f,0,1) ans= 3.1945 符号积分精确度高但速度慢,有时候有些函数没有解析解,就得用数值积分,并且数值积分速度快,但精确度
这有知两个原因: