请问这个积分在matlab里用数值积分的方法应该如何做？

matlab中怎样对二元函数中的一个变量做数值积分？~

用int函数不行吗？
我给你举一个可以通用的例子：
设z=axy+by^2 ，其中a、b为常量，x、y为变量
下面求a=1,b=2,x=3;时，对y在[0,1]上的积分，步骤如下：
首先，执行如下命令：yy=int('a*x*y+b*y^2','y',0,1)
会显示以下结果：
yy =

1/2*a*x+1/3*b
然后输入以下命令：
a=1;b=2;x=3;
yyyyy=subs(yy)
如下是结果：
yyyyy =

2.1667
不知道，说明白了吗！！！！！！！

1、在matlab中，积分运算有多种方式，为了便于查看不同方式处理异同，以下面这个积分为例：

2、梯形积分法
第一种，采用最简单的方式，以函数trapz为例，z = trapz(x,y) 其中x表示积分区间的离散化向量，y是与x同维数的向量，表示被积函数，z是返回的积分近似值。
clc;clear;
% 梯形积分法
x = -1:0.001:1;
y = exp(-x.^2);
s = trapz(x,y)
% 计算结果： s = 1.4936
3、高精度数值积分（1）
为了克服梯形积分法精度低的问题，可以采用高精度积分方式，第一种可以采用 z = quad(Fun，a,b) 该方式是自适应步长Simpson计分法求得函数Fun在区间[a,b]上定积分，如下：
clc;clear;
% 梯形积分法
s = quad(inline('exp(-x.^2)'),-1,1)
% 计算结果： s = 1.4936
4、高精度数据积分（2）
采用高精度Lobatto积分法，格式： z = quadl(Fun,a,b)
clc;clear;
% 梯形积分法
s = quadl(inline('exp(-x.^2)'),-1,1)
% 计算结果： s = 1.4936
% 注：在编写完代码后，要按如下图红色箭头所指处运行程序才会有输出！

T_0=linspace(300,1500,101);

%T_0=linspace(300,1500,1201);

Eg=1.519-5.405*(10^(-4))*((T_0.^2)./(T_0+204));

%T_0=linspace(300,1500,1201);

Tc=repmat(T_0,100,1);

%Tc=repmat(T_0,1e4,1);

E0=repmat(Eg,100,1);

%E0=repmat(Eg,1e4,1);

k=1.38*10^(-23);

C0=linspace(1,10,100);

%C0=linspace(1,10,1e4);

C=repmat(C0',1,101);

%C=repmat(C0',1,1201);

hv=E0.*C;

R_matrix=hv.^2./(exp(hv./k./T_0)-1);

%R_matrix=hv.^2./(exp(hv./k./T)-1);

R=sum(R_matrix)*10.0250719315895/100;

R1=sum(sum(R_matrix)*10.0250719315895/100)

plot(T_0,R);

---

望城县17176491856： 如何用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

望城县17176491856： 求教大神用matlab数值解积分方程 - ？
褚泽地巴： 这个方程无解.如果 |x| < 1,积分项现复数,而另一项为实数,无解.而当 |x| >= 1时,积分项的值大于0.5,另一项最大为0.5,仍然无解. 可以很容易画出方程左端的函数图像:12 syms x e real ezplot(int(sqrt(x^2-e^2),e,0,1)-1/(x^2+1),[-22]) 由图可见,函数值始终为正(其中-1~1区间为复数的实部),所以方程无解.

望城县17176491856： 请问如何用matlab解这个积分方程组?? - ？
褚泽地巴： 这个积分比较复杂,得不到解析解的表达式 所以用符号运算给不出解 可以考虑用数值办法求解B=15;A=0.2;cv=5;r=0.07;f=@(t,x) 1/B*log(B./(cv+x*exp(r*t))); f1=@(T,x) quad(@(t) f(t,x),0,T)-A; f2=@(T,x) T-1/r*log((B-cv)/x); fun=@(X) [f1(X(1),X(2));f2(...

望城县17176491856： matlab求带参数的积分,符号积分解不出,怎么用数值积分? - ？
褚泽地巴： t=0:0.01:10;a= ;b= ;%自己对ab赋值f=(exp((-(a-t).^2-b^2)./t)+t)./t;y=sum(f)*0.01;%积分后的近似结果

望城县17176491856： Matlab里面怎么求一个分段函数的积分 - ？
褚泽地巴： 图像可分段画,如下语句:x0=-1:0.01:0;y0=zeros(size(x0));x1=0:0.01:1;y1=x1;x2=1:0.01:2;y2=-ones(size(x2));x3=2:0.01:3;y3=zeros(size(x3));plot(x0,y0,x1,y1,x2,y2,x3,y3) 积分可分段积分,如下 f=0+int(x,0,1)-1*2+0

望城县17176491856： 用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]上的定积分

望城县17176491856： 怎么用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)%三重(闭型)数值积分指令 .

望城县17176491856： matlab中怎样对二元函数中的一个变量做数值积分? - ？
褚泽地巴： clear all clc syms va rho tau f=@(va,rho,tau) (2.5-0.1*(va+rho*sin(tau)))./(exp(2.5-0.1*(va+rho*sin(tau)))-1); n = 5000;%当想要结果更精确时,可以把n设置更大 tau = linspace(0,2*pi,n);%自编的简单方法 Tn =@(va,rho) pi/n*(sum(f(va,rho,tau(1:end-...

望城县17176491856： 怎样用Matlab求符号表达式的数值积分? - ？
褚泽地巴： f=@(x)sin(x)+cos(x); %定义被积函数Intf=quad(f,0,1) %对f进行积分,下限为0,上限为1

望城县17176491856： 高分求助:matlab求积分？
褚泽地巴： 是sin(o1^2),还是sin(o1)*sin(o1)啊 可以用数值积分,也可以用符号积分,但数值积分的应用范围比符号的广. 这里用的是数值积分. 在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%...