mdsk.net
当前位置:首页 >> mAtlAB求二元定积分 >>

mAtlAB求二元定积分

可以用复合的形式,比如说:syms a x y; f=x^2+y^2; I=int(int(f,x,(y-a),y),y,a,3*a) %先对x积分,积分限为y-a到y,再对y积分,积分限为a到3a.执行结果:I =14*a^4

x=0.6时,y=68.581359956511627574834151414947 close allclear,clcx=0:0.05:1;y=zeros(1,length(x));syms y yex1=int(int('exp(-(t-50)^2/200)',-inf,y),0,y);ex2=int(int('exp(-(t-50)^2/200)',-inf,0.6*y),0,0.6*y);ex1=ex1/(75*sqrt(2)*pi^(1/2));ex2=ex2/(75*sqrt

方法一:elta_x=100;delta_y=100;V=f(x,y)*delta_x*delta_y;sum(V(:))方法二:xint=trapz(x,f);yint=trapz(y,xint);

必须给出函数f(r+x,c+y)的显式表达式,才能对它进行数值或者符号积分 比如:f=(r+x)^2+(c+y)^2;>> syms w L x y;>> I=int(int('(r+x)^2+(c+y)^2',x,-w/2,w/2),y,-L/2,L/2) I =1/12*w*L^3+1/12*w^3*L+r^2*w*L+c^2*w*L 另外符号积分用的是INT函数,积分上下限为数值才用QUAD

R=int(s,v) %对符号表达式s中指定的符号变量v计算不定积分.表达式R只是表达式函数s的一个原函数,后面没有带任意常数C.R=int(s) %对符号表达式s中确定的符号变量计算计算不定积分.R=int(s,a,b) %符号表达式s的定积分,a,b分别为积分的上、下限R=int(s,x,a,b) %符号表达式s关于变量x的定积分,a,b分别为积分的上、下限trapz(x,y) 梯形积分法,x时表示积分区间的离散化向量,y是与x同维数的向量,表示被积函数,z返回积分值.fblquad('fun',a,b,c,d) 矩形区域二重数值积分,fun表示被积函数的M函数名,a,b分别为x的上、下限,c,d分别为y的上、下限.

用符号积分求解 clc; syms x y; %定义两个符号变量 a=int(int(x*y,x,y*y,y+2),y,-1,2); %积分限x:y*y,y+2 ,y:,-1,2 b=simple(a); %化简 c=vpa(b,4) %得到4位近似解,也可以任意n位解

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]上的定积分

符号计算 syms x y; %定义两个符号变量 a=int(int(x^y,x,0,1),y,1,2) %积分x,0,1 ,y,1,2 b=simple(a) %化简 c=vpa(b,4) %得到4位近似解,也可以任意N位解 数值计算%%二重积分f= @(x,y)exp(sin(x))*ln(y),y从5*x积分到x^2,x从10积分到201 (7.X后版本

这种问题最好给出具体的函数,而不是泛泛地问,因为根据函数的复杂程度,可能用解析法,也可能用数值法.解析法适用于函数积分能够求出解析表达式的情况,例如:1234 syms x y f=x^2*y^3; J=int(f,x,-5,5); ezplot(J) 如果上述使用int求积分得不出解析表达式,一般就应该考虑使用数值方法求解了.如果需要,还是具体结合你的函数来说吧.

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

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