复化梯形公式求二重积分matlab源码

这段代码具有很好的交互性和通用性，将代码复制到matlab编辑器之后，按照提示操作即可。傻瓜式操作，结果一目了然～

%%%%%%%%%%     2020.6.5        %%%%%%%%%
%%%%%%%%%%复化梯形公式求二重积分%%%%%%%%%
%%%%%%%%%%     Liu Deping      %%%%%%%%%
clear all;
%%%被积函数及积分上下限导入；
s=input('请输入函数表达式：f = ','s');
f=inline(s);
a = input('请输入积分变量x左边界a的值:');
b = input('请输入积分变量x右边界b的值:');
c = input('请输入积分变量y左边界c的值:');
d = input('请输入积分变量y右边界d的值:');
h1 = input('请输入积分变量x步长h1的值:');
h2 = input('请输入积分变量y步长h2的值:');
m=round((b-a)/h1);
n=round((d-c)/h2);
%%%系数矩阵T，t[i,j]为复化梯形公式的系数；
TT=zeros(m+1,n+1);
TT(1,1)=1;
TT(m+1,1)=1;
TT(1,n+1)=1;
TT(m+1,n+1)=1;
for i=2:mTT(i,1)=2;TT(i,n+1)=2;
end
for j=2:nTT(1,j)=2;TT(m+1,j)=2;
end
TT(TT==0)=4   %%系数表%%%计算各二维节点函数值，并存放于矩阵F中；
F=zeros(m+1,n+1);
for i=0:mfor j=0:nF(i+1,j+1)=f(a+i*h1,c+j*h2);    end
end %%%结果输出；
format long;
fprintf('各节点函数值；')
F
fprintf('复化梯形公式计算结果；')
Tnm=(b-a)*(d-c)/(4*m*n)*sum(sum(TT.*F))

之后还有辛普森公式求二重积分，希望喜欢的小伙伴持续关注篱落～～成殇～～

复化梯形公式求二重积分matlab源码相关推荐

微分方程一维抛物热传导方程向前向后欧拉C-N格式二阶BDF格式MATLAB源码显式欧拉，隐式欧拉，梯形公式，改进欧拉
微分方程一维抛物热传导方程向前向后欧拉C-N格式二阶BDF格式MATLAB源码显式欧拉,隐式欧拉,梯形公式,改进欧拉五点差分,九点差分差分格式,紧差分格式直拍,只有pdf版方法说明 word版 ...
matlab的数值计算方法,数值计算方法中的一些常用算法的Matlab源码
数值计算方法中的一些常用算法的Matlab源码,这些程序都是原创,传上来仅供大家参考,不足之处请大家指正,切勿做其它用途-- 说明:这些程序都是脚本函数,不可直接运行,需要创建函数m文件,保存时文件名 ...
【故障检测问题】基于matlab免疫算法求解故障检测问题【含Matlab源码 196期】
一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[故障检测问题]基于matlab免疫算法求解故障检测问题[含Matlab源码 196期] 获取代码方式2: 通过订阅紫极神光博客付费专栏,凭 ...
详解非局部均值滤波原理以及用MATLAB源码实现
详解非局部均值滤波原理以及用MATLAB源码实现序言均值滤波.中值滤波.高斯滤波在滤除噪声的过程中,无可避免的使图像的边缘细节和纹理信息所被滤除.针对此问题,Buades[1]等人提出了非局部均值 ...
fama matlab源码_基于优化算法改造的Fama-French三因子模型
基于光大证券金融工程研报<站在巨人的肩膀上,从牛基组合到牛股发现 --FOF 专题研究系列之十六 >中提及的Carhart四因子Alpha优化模型,本文在Fama-French三因子模型上 ...
Kmeans聚类算法及其matlab源码
本文介绍了K-means聚类算法,并注释了部分matlab实现的源码. K-means算法 K-means算法是一种硬聚类算法,根据数据到聚类中心的某种距离来作为判别该数据所属类别.K-means算法 ...
【APF三维路径规划】基于matlab人工势场算法无人机三维路径规划【含Matlab源码 168期】
一.获取代码方式获取代码方式1: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码. 获取代码方式2: 完整代码已上传我的资源:[三维路径规划]基于matlab人工势场算法无人机三维 ...
【RRT三维路径规划】基于matlab RRT算法无人机三维路径规划【含Matlab源码 155期】
一.获取代码方式获取代码方式1: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码. 获取代码方式2: 完整代码已上传我的资源:[三维路径规划]基于matlab RRT算法无人机三维 ...
【A_star二维路径规划】基于matlab A_star算法无人机二维路径规划（起终点障碍物可设置）【含Matlab源码 1321期】
⛄一.获取代码方式获取代码方式1: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码. 获取代码方式2: 完整代码已上传我的资源:[A_star二维路径规划]基于matlab A_s ...
果蝇优化算法（Fruit Fly Optimization Algorithm，FOA）-Matlab源码
获取更多资讯,赶快关注上面的公众号吧! 文章目录果蝇优化算法(Fruit Fly Optimization Algorithm,FOA) 启发初始化食物搜索计算味道浓度判定值适应度评估确定 ...

复化梯形公式求二重积分matlab源码

复化梯形公式求二重积分matlab源码

复化梯形公式求二重积分matlab源码相关推荐

最新文章

热门文章