Matlab 积分区域为任意四边形
来源:学生作业帮 编辑:搜狗做题网作业帮 分类:数学作业 时间:2024/04/27 11:13:51
Matlab 积分区域为任意四边形
我求一个函数的二重积分,积分函数很简单.但是积分区域为平面上给定的四个节点.例如1、(x1,y1).4、(x4,y4),感觉处理上很麻烦,
积分函数为:z=2(x-y^2)^3
我想求出这个积分的解,积分区域为任意的四边形
我求一个函数的二重积分,积分函数很简单.但是积分区域为平面上给定的四个节点.例如1、(x1,y1).4、(x4,y4),感觉处理上很麻烦,
积分函数为:z=2(x-y^2)^3
我想求出这个积分的解,积分区域为任意的四边形
计算速度较慢,请耐心等待
clear;clc
pos=[0 -2;2 0;0 2;-2 0];%凸多边形顶点,自己改
f=@(x,y)2*(x-y.^2).^3;%积分函数
X=pos(:,1);Y=pos(:,2);
k=convhull(X,Y);
X=X(k);Y=Y(k);
fill(X,Y,'g');title('积分区域');
xmin=min(X);xmax=max(X);
ymin=min(Y);ymax=max(Y);
g=@(x,y)f(x,y).*inpolygon(x,repmat(y,size(x)),X,Y);
val=quadl(@(y)arrayfun(@(y)quadl(@(x)g(x,y),xmin,xmax),y),ymin,ymax)%这就是你要的结果
再问: 朋友,你好! 谢谢你的回答,可不可以问下程序设计的基本思路?后面几行不是很明白! 谢谢啦!
再答: 基本思路很简单,quadl命令求解一维数值积分,于是多次利用quadl可以求解矩形区域的多维数值积分,那么怎么求解非矩形区域的积分呢?只要让指定区域外的函数值等于0就可以了,倒数第二行就是干这个事情,最后一行就是利用两次quald命令求解二维积分。
clear;clc
pos=[0 -2;2 0;0 2;-2 0];%凸多边形顶点,自己改
f=@(x,y)2*(x-y.^2).^3;%积分函数
X=pos(:,1);Y=pos(:,2);
k=convhull(X,Y);
X=X(k);Y=Y(k);
fill(X,Y,'g');title('积分区域');
xmin=min(X);xmax=max(X);
ymin=min(Y);ymax=max(Y);
g=@(x,y)f(x,y).*inpolygon(x,repmat(y,size(x)),X,Y);
val=quadl(@(y)arrayfun(@(y)quadl(@(x)g(x,y),xmin,xmax),y),ymin,ymax)%这就是你要的结果
再问: 朋友,你好! 谢谢你的回答,可不可以问下程序设计的基本思路?后面几行不是很明白! 谢谢啦!
再答: 基本思路很简单,quadl命令求解一维数值积分,于是多次利用quadl可以求解矩形区域的多维数值积分,那么怎么求解非矩形区域的积分呢?只要让指定区域外的函数值等于0就可以了,倒数第二行就是干这个事情,最后一行就是利用两次quald命令求解二维积分。
关于积分区域Ω为椭球的三重积分
matlab 解系数为矩阵的积分
请问matlab中计算数值积分(三角形区域)有哪几种函数?
考研数学题:积分区域为球体的三重积分.利用极坐标系.
广义积分 ∫ln(1-x^2)dx收敛于________(积分区域为0-1)
matlab 区域膨胀
matlab 空白区域搜索
把f(x,y) 形成的二次积分化为极坐标形式的二次积分,其中积分区域D为
如何用matlab求任意一函数的定积分,相关的程序有吗?
对定积分∫(x-1)∧2dx(积分区域为0到2)用换元法t=(x-1)∧2时,t的积分区域
Matlab的Monte carlo法求定积分 区域(-π/2,π/2),y=cosx所围面积
求问 如何用matlab 对图像中的一块四边形区域进行处理,比如进行傅里叶变化?