yhkn.net
当前位置:首页 >> 用mAtlAB计算定积分, 谢谢!急求! >>

用mAtlAB计算定积分, 谢谢!急求!

可以用符号运算来求:>> J=int(exp(x)*(x^2+sqrt(x)),0,1) J = 2*exp(1)+1/2*i*pi^(1/2)*erf(i)-2>> double(J)ans = 1.9739>> vpa(J,30) ans = 1.97391191101090886191652635584得到的结果是解析表达式,要想得到具体的数,可以用double或vpa函数(其中vpa可以得到任意精度的结果,而double得到的是双精度浮点数).当然,也可以用数值积分:>> quadl(inline('exp(x).*(x.^2+sqrt(x))'),0,1)ans = 1.9739

第一步 将function y = myfun(x) y = 1./(x.*(x.^2 - 3*x + 2 ).^(1/3));保存为M文件.(文件→新建→M文件)第二步将warning off allQ = quad(@myfun,2,100000)Q = quad(@myfun,2,10000000)Q = quad(@myfun,2,1000000000)warning on 粘贴至(命令窗口)结果显示,当上限为100000,10000000,1000000000时,定积分的值为x=1.4389,1.4396,1.4396.因此,可以将1.4396作为此定积分的值.

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) % 积分值

你是SYMBOL数学还是其他,如果是数值解的话用ODE命令,常用ODE45 LUNGEN KUTTA 法.

>> syms x f1 f2f1= ((sin(x)).^3 - (sin(x)).^5).^(1/2);f2=int(f1,0,pi);simplify(f2)ans = 0.8000000000

syms x a;b=int(1/(0.000411*(209.4*a^3*x^2+3117856727*a^4*x^3-0.033*x^2-215.4)),x,20,40)基本不管多复杂的式子,MATLAB都能解决,只要你输入正确,解出来只是时间问题,我给你一个解题的参考格式.设f为你的被积函数,只能自己输,别输错;x为积分变量,a,b分别为积分下限和上限.MATLAB中计算函数f的不定积分的格式如下:syms x;int(f,x)计算定积分的格式如下:syms a b x;int(f,x,a,b).例子: syms x a; b=int(a*exp(-a*x),x,0,10)结果b =1 - 1/exp(10*a)

1, 30);y = (k-t)*sin(t);Q = int(int(y, t, 0, k), ksyms k t

syms xs=int(exp(0.35*x)-exp(0.69)*cos(0.35*sqrt(1-x^2)),'x',0.24,0.97);vps(s,5)计算得到,s=-.49861

syms xint(abs(x)*exp(x),x,-1,2)ans =exp(2) - 2/exp(1) + 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]上的定积分

网站首页 | 网站地图
All rights reserved Powered by www.yhkn.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com