求计算一重积分的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