MATLAB 一重积分
syms x y;
f_1=0.2*exp(-0.5*((x-16)^2)+((y-10)^2)/0.64-((x-16)*(y-10))/2);
xup=y/4+sqrt(6-1.5*y^2);
xlow=y/4-sqrt(6-1.5*y^2);
fs1=int(f_1,'x',xlow,xup);
fs2=int(fs1,'y',-2,2);
fs2 =
int((exp(-(135*y)/4)*exp(675/4)*exp((27*y^2)/16)*(erf((21*2^(1/2))/2 - (3*2^(1/2)*y)/8 + (12 - 3*y^2)^(1/2)/2) + erf((3*2^(1/2)*y)/8 - (21*2^(1/2))/2 + (12 - 3*y^2)^(1/2)/2))*(2*pi)^(1/2))/10, y, -2, 2)
syms x y;
f1=x*y;
f2=int(f1,'y',1,x);
f3=int(f2,'x',1,2);
扩展资料:
运用的符号积分命令int用法:
int(fun):求函数fun的不定积分;
int(fun,var):求函数fun关于变量var的不定积分;
int(fun, var, a,b,):求函数fun的在[a,b]间的定积分或广义积分;
示例
>> clear;syms x y z;
>> int(sin(x*y+z),z)
ans = -cos(x*y+z)
参考资料来源:百度百科—INT函数
syms x y z
f = y * sin(x) + z * cos(x);
int1 = int(f,z, z_1, z_2); % z_1, z_2 分别为z的上下限
int2 = int(int1, y, y_1, y_2); % y_1, y_2 分别为y的上下限
int_f = int(int2, x, x_1, x_2); % x_1, x_2 分别为x的上下限
syms t;
f = 3696.0*t/(222.0-6*t)*(1-11/3*t/(222-6*t));
g =int(f);
这时得到,
g= -8932/9*t+4638172/9/(-37+t)-455840/9*log(-37+t)
因为matlab的这个积分结果只是当-37+t》0时成立,实际上我们考虑的情况是t从零开始的,所以这时t-37一般是小于零,(Ta是不能达到37的,因为222-6t不能在积分过程中为零)。
所以我们需要手工的将log里面的-37+t 改为37-t,再算。即输入
g= -8932/9*t+4638172/9/(-37+t)-455840/9*log(37-t);
然后当t=0时得,
ans1 = limit(g,0.0);
equ = g-ans1-248;
solve(equ)
既得到答案
28.130706728683984639179887364200
然后可以通过matlab数值积分函数quad验证其正确性:
f = @(t)3696.*t./(222-6.*t).*(1-11./3.*t./(222-6*t));
quad(f,0,28.130706728683984639179887364200)
即可得到248,证明正确
玉仇信必: 首先 syms t; f = 3696.0*t/(222.0-6*t)*(1-11/3*t/(222-6*t)); g =int(f); 这时得到, g= -8932/9*t+4638172/9/(-37+t)-455840/9*log(-37+t) 因为matlab的这个积分结果只是当-37+t》0时成立,实际上我们考虑的情况是t从零开始的,所以这时t-37一般是小...
长春市18366521000: 如何用matlab算一重,二重,三重积分.请高人解答,最好举例子.不甚感激.财富值只剩4分了.. - ?
玉仇信必: syms x fx=sin(x);%fx为被积分函数 fx1=int(fx,'x','0','pi')%fx为被积分函数,x,为自变量0,pi为一重积分上下限 fx2=int(int(int(fx,'x','0','pi')),'x','0','2*pi')%x为自变量0,2*pi为二重积分上下限 fx3=int(int(int(fx,'x','0','pi'),'x','0','2*pi'),'x','0','t')%x为自变量0,t为二重积分上下限三重积分
长春市18366521000: matlab求一重和二重积分的近似值 - ?
玉仇信必: syms x y int(exp(-x^2)*cos(x^2),x,0,pi); vpa(ans,5) ans =.68850-.64202e-5*i int(int(cos(x^2-y^2),y,0,sqrt(pi)),x,0,sqrt(pi)); vpa(ans,5) ans =1.2401
长春市18366521000: 如何利用matlab进行一次线性函数积分 比如对y=3x+1在x等于2到5上积分,求详细代码及讲解 - ?
玉仇信必: fun = @(x) 3*x+1 ; % 定义函数 I = quad(fun, 2, 5) % 2~5积分
长春市18366521000: matlab怎么算多重积分 - ?
玉仇信必: 符号积分 >> syms x y z>> int(int(int(x+y+z,y,0,1-x-z),x,0,1-z),z,0,1) ans = 1/8数值积分 >> triplequad(@(x,y,z)(x+y+z).*(x+y+z<=1),0,1,0,1,0,1)ans = 0.1250
长春市18366521000: matlab重积分函数 x允许向量输入 什么意思 - ?
玉仇信必: 意思就是X可以是向量...X=[a1 a2 a3...]
长春市18366521000: 如何用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
长春市18366521000: 大家帮忙用matlab进行数值积分,并列出命令行 - ?
玉仇信必: 用matlab可以如下数值积分法,来求解定积分、二重积分、三重积分的数值解问题. 1、梯形数值积分计算 trapz() X = 0:pi/100:pi; Y = sin(X); Z = pi/100*trapz(Y) 2、自适应辛普森数值积分计算 quad() F = @(x)1./(x.^3-2*x-5); Q = quad(F,0,2); 3、自...
长春市18366521000: 用matlab编写一个复杂的公式,其中包含2重积分,应该怎么编 - ?
玉仇信必: quadl(@(x) arrayfun(@(xx) quadl(@(y) xx*y,sin(xx),cos(xx)),x),1,2)
长春市18366521000: 用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]上的定积分