求计算一重积分的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函数
这个可以用matlab的符号积分或者数值积分解决，下面提供4种方法：
1、直角坐标系符号积分
syms x yint_y = int(sin(pi*(x^2+y^2)),y,-sqrt(1-x^2),sqrt(1-x^2)); % 先对y积分I = vpa(int(int_y,x,-1,1)) I = 2.0
2、极坐标系符号积分
syms rI = 2*pi*int(r*sin(pi*r^2),r,0,1)I = 2
3、dblquad()是数值积分（针对matlab早期版本）
f=@(x,y)sin(pi*(x.^2+y.^2)).*(x.^2+y.^2<=1);I = dblquad(f,-1,1,-1,1)I = 2.0000
4、integral2()数值积分（较新版本matlab）
f=@(x,y)sin(pi*(x.^2+y.^2));y_l = @(x) -sqrt(1-x.^2);y_u = @(x) sqrt(1-x.^2);I = integral2(f,-1,1,y_l,y_u)I = 2.0000
syms x
f=2*x+2;%设函数
n=(b-a)/0.01;%间距设为0.001；
lie=linspace(a,b,n);%横坐标数列
s=1;
for i=1:n-1
h=(subs(f,x,lie(i))+subs(f,x,lie(i+1)))/2;
area(s)=0.001*h;
s=s+1;
end
s
sum(area)
这是本人编的积分程序，间距可调控制精度，不过有个缺点，如果区间太大，电脑慢的要死！！
这是用分解梯形在求和的简单思想做的程序，思想肯定没错。
求解物理中的微分方程:mg-kv=m(dv\/dt),其中mg是一个物体的重力
dv\/dt=a mg-kv=ma 这道题应该是物体下落时，受到的空气阻力f与速度成正比，即kv。整个方程的意思是F合力＝mg-f＝mg-kv 而合力又等于ma，所以mg-kv=ma。因为物体下落时速度不断改变，加速度也不断改变，所以用积分式。
求解转动惯量相关问题
先求A的转动惯量（可直接套公式，也可积分求解）。我从不记这些公式，都是积分，好在积分很简单。在A上任取一小段长度为dr的质量元，设该质元距O点距离为r，该质元的质量dm=mAdr\/L，则转动惯量 J=∫r^2 dm =∫r^2 mAdr\/L=mA\/L ∫r^2 dr =1\/3 mA L^2，其中积分上下限分别为L，0...
看了你的回答,那你觉得从狭义相对论到广义就是一个简单的积分过程...
总的说就是一个积分过程，但是对于相对论的发展来说，并不是从简单微分到积分这样探索过来的，经过了很多曲折。最初是洛伦兹是想通过一些分析，找出麦克尔逊莫雷实验测量不到光速变化的原因。他也做了大量的分析和推导。最后得到了一个 t'＝t×√(1＋v²\/c²) 的变换式，但是因为他从来...
高中物理题?
W=是F对X的积分，平均P是F对一段X的积分（功）\/对应的时间段（肯定）如果从P瞬间=F瞬间V瞬间来看，那么只要P能表示成时间t的函数，就可以用P对T的积分\/对应时间段来计算（应该是）求重力做功的平均功率，那么，求出W重力，再求出时间，时间的话，使用动能定理：1\/2mv2=pt-mgH,得出pt的值，...
微积分问题
阻力与速度成正比，即f=kv,其中k为比例系数。设物体的质量为m,当t=0时，物体初速度为v0.用牛顿第二定律列微分方程：mg-kv=m(dv\/dt)解此方程并带入初始条件得：v(t)=mg\/k-(mg\/k-v0)e^(-kt\/m)
牛顿第二定律是怎么得出的?
现在公认的看法是:牛顿首先发明微积分,莱布尼兹后来也发明了微积分,但他早于牛顿向外公布出来。由于他曾看过牛顿的论文手稿,因此,争议的焦点在于,到底按...把万有音引力定律应用于在地球上做自由落体运动中,那么G=Km1(地球质量)\/r2(地球半径的平方)*m2(物体质量)。K,m1和r2都是常数,算出来就是重力加速度...
数学的发展历史
微分学包括求导数的运算，是一套关于变化率的理论。它使得函数、速度、加速度和曲线的斜率等均可用一套通用的符号进行讨论。积分学，包括求积分的运算，为定义和计算面积、体积等提供一套通用的方法。第四时期 现代数学。现代数学时期，大致从19世纪初开始。数学发展的现代阶段的开端，以其所有的基础---...
数学的发展史是什么?
第三时期：变量数学时期（公元十七世纪初—十九世纪末）变量数学产生于17世纪，经历了两个决定性的重大步骤：第一步是解析几何的产生；第二步是微积分的创立。第四时期：现代数学时期（十九世纪末开始），数学发展的现代阶段的开端，以其所有的基础---代数、几何、分析中的深刻变化为特征。数学需要严谨...
如何推导单摆周期计算公式
采用牛顿第二定律推导：如下图，摆长为l,重物受力为：重力mg和绳子的张力T。取如图所示的二维坐标系，张力T可以分解为垂直和水平方向的二个力。L与垂线的夹角为θ。F=ma，可以列出重物在x和y二个方向上的运动方程：这二个微分方程相当难解，所以只能采用一种“小角度近似”的方法进行处理，解的...
整体牛顿第二定律条件
但由于当时没有规定1个单位的力的大小,于是取k=1,就有∑F=ma,这就是今天我们熟悉的牛顿第二定律的表达式。2.公式 F合=ma (单位:N(牛)或者千克米...联立各方程得到物体的运动学方程,然后依据题目要求积分求出位移、速度等。编辑本段牛顿第二定律的应用 牛顿第二定律是经典力学的基础和核心,是分析、研究和...
壤殷复方： 程序: 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
中方县18713262707： 如何用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为二重积分上下限三重积分
中方县18713262707： 怎样用matlab编程求解积分 - ？
壤殷复方： int函数 例如s = sin(x),求积分int('s'),即可
中方县18713262707： 求助:用matlab求积分,求个程序, - ？
壤殷复方： 求助:用matlab求积分,求个程序,ms x>> S=int(1617/(0.4083*sqrt(2*pi))*exp(-(x-0.8965).^2/(2*(0.4087^2))),15458,15708) S =-(868581423309425*6018094969932590^(1/2)*pi^(1/2)*(erf((129663582076928*6018094969932590^(1/2)...
中方县18713262707： 求高手帮忙看一下Matlab求积分的程序,谢谢 - ？
壤殷复方： clear a=0:0.5:2; n=5; for i=1:n; a(i)*quadl(@(x)exp(-x.^2),0,1) end
中方县18713262707： 求助:用matlab求积分,求个程序, - ？
壤殷复方： >> S=int(1578.394./(0.33407138.^x),0.6493,0.6743) S = (6941850224847487*4096^(6493/10000)*4398046511104^(6493/10000)*6018094969932589^(3257/10000)...
中方县18713262707： 急: matlab中求积分 - ？
壤殷复方： 因为你这里面有多个符号变量,而且积分变量是符号,那只能符号积分了.看样子式子中的符号变量应该有L Ts u t.程序如下:syms L Ts u t %定义符号变量 g=(1/2*L*Ts)*(1-u*(cos(2*pi*t))/(L*Ts)); %你的式子是这样吧?int(g,t,0,t) %符号积分 结果为:ans =-1/4*(-2*t*L*Ts*pi+u*sin(2*pi*t))/pi
中方县18713262707： 用matlab求积分 - ？
壤殷复方： 用quad函数,例如求f(x)=1./(x.^3-2*x-5)从0到2的积分: F = @(x)1./(x.^3-2*x-5); %生成函数句柄 Q = quad(F,0,2);
中方县18713262707： 怎么用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)%三重(闭型)数值积分指令 .
中方县18713262707： MATLAB中积分的计算指令 - ？
壤殷复方： syms x; int(x^3*log10(x))答案:4 4x log(x) x1/4 --------- - 1/16 -------log(10) log(10)如果对您有帮助,请记得采纳为满意答案,谢谢!祝您生活愉快! vaela
