复化辛普森公式求二重积分matlab源码及例题

复化辛普森公式求二重积分matlab源码

直接拷贝到matlab编辑器，傻瓜式操作。具体算法自行探究，网上都有，小编只提供代码。用的好的请加个关注，篱落～～成殇～～再次先行谢过。

%%%%%%%%%%     2020.6.5        %%%%%%%%%
%%%%%%%%%%复化Simpson公式求二重积分%%%%%%%%%
%%%%%%%%%%     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]为复化梯形公式的系数；
T=zeros(m+1,n+1);
R=[1,4,1;4,16,4;1,4,1];
for p=0:m/2-1for q=0:n/2-1for i=1:3for j=1:3T(i+2*p,j+2*q)=T(i+2*p,j+2*q)+R(i,j);endendend
end
fprintf('系数表：');
T
%%%计算各二维节点函数值，并存放于矩阵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('复化Simpson公式计算结果；')
Snm=(b-a)*(d-c)/(9*m*n)*sum(sum(T.*F))

下面附上一个例题吧，感受一下代码的强大。

只看第(2)问，利用辛普森公式计算积分。(关于第一问的代码可以去关注上一篇文章〈复化梯形公式公式求二重积分matlab源码〉)
步骤:
1.复制代码到matlab编辑器，点击运行;
2.按照提示输入:

3.得到结果(包括系数、节点函数值和积分结果):

如果代码对你有帮助，请加一下关注，如果有问题可以私聊！最后，再次感谢喜欢我的小伙伴。

复化辛普森公式求二重积分matlab源码及例题相关推荐

复化辛普森公式的误差matlab,求两个matlab的程序，分别利用复化梯形公式和辛普森公式求解误差函数erf(x)。...
main.m g=@(t)exp(-t.^2); x=linspace(0,5,100); y1=zeros(1,100); y2=zeros(1,100); for i = 1:100 y1(i)= ...
【数值分析实验MATLAB】数值积分：梯形公式、辛普森公式、复化梯形公式、复化辛普森公式、龙贝格算法、自适应求积方法
(1)Matlab代码: f1: if x==0y=1; else%y=exp(-x^2);y=(sin(x))/x; end 梯形公式: a=0; b=1; tx=(b-a)/2*(f1(a)+f1 ...
matlab的数值计算方法,数值计算方法中的一些常用算法的Matlab源码
数值计算方法中的一些常用算法的Matlab源码,这些程序都是原创,传上来仅供大家参考,不足之处请大家指正,切勿做其它用途-- 说明:这些程序都是脚本函数,不可直接运行,需要创建函数m文件,保存时文件名 ...
详解非局部均值滤波原理以及用MATLAB源码实现
详解非局部均值滤波原理以及用MATLAB源码实现序言均值滤波.中值滤波.高斯滤波在滤除噪声的过程中,无可避免的使图像的边缘细节和纹理信息所被滤除.针对此问题,Buades[1]等人提出了非局部均值 ...
果蝇优化算法（Fruit Fly Optimization Algorithm，FOA）-Matlab源码
获取更多资讯,赶快关注上面的公众号吧! 文章目录果蝇优化算法(Fruit Fly Optimization Algorithm,FOA) 启发初始化食物搜索计算味道浓度判定值适应度评估确定 ...
【故障检测问题】基于matlab免疫算法求解故障检测问题【含Matlab源码 196期】
一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[故障检测问题]基于matlab免疫算法求解故障检测问题[含Matlab源码 196期] 获取代码方式2: 通过订阅紫极神光博客付费专栏,凭 ...
fama matlab源码_基于优化算法改造的Fama-French三因子模型
基于光大证券金融工程研报<站在巨人的肩膀上,从牛基组合到牛股发现 --FOF 专题研究系列之十六 >中提及的Carhart四因子Alpha优化模型,本文在Fama-French三因子模型上 ...
蚁狮优化算法( Ant Lion Optimizer,ALO)-Matlab源码
获取更多资讯,赶快关注上面的公众号吧! 文章目录蚁狮优化器Ant Lion Optimizer (ALO) 启发 ALO算子随机游走困于蚁狮坑构造陷阱蚂蚁滑向蚁狮捕捉猎物和重建坑精英 A ...
Kmeans聚类算法及其matlab源码
本文介绍了K-means聚类算法,并注释了部分matlab实现的源码. K-means算法 K-means算法是一种硬聚类算法,根据数据到聚类中心的某种距离来作为判别该数据所属类别.K-means算法 ...
【APF三维路径规划】基于matlab人工势场算法无人机三维路径规划【含Matlab源码 168期】
一.获取代码方式获取代码方式1: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码. 获取代码方式2: 完整代码已上传我的资源:[三维路径规划]基于matlab人工势场算法无人机三维 ...

复化辛普森公式求二重积分matlab源码及例题

复化辛普森公式求二重积分matlab源码

复化辛普森公式求二重积分matlab源码及例题相关推荐

最新文章

热门文章