MATLAB高斯-勒让德积分求解

2024-04-12 23:59:27

MATLAB高斯-勒让德积分求解

积分函数

function [ y ] = CalcuFunctionValue( x )
y=exp(-x^2);%积分函数
end

求积分算法函数

function [ result ] = gslrdjf( a,b,n )
if n2
GaussP=[-0.5773503 0.5773503];
GaussA=[1 1];
end
if n3
GaussP=[-0.7745966 0 0.7745966];
GaussA=[0.5555556 0.888888889 0.5555556];
end
if n4
GaussP=[-0.8611363 -0.3399810 0.3399810 0.8611363];
GaussA=[0.3478548 0.65214515 0.65214515 0.3478548];
end
if n5
GaussP=[-0.9061798 -0.5384593 0 0.5384593 0.9061798];
GaussA=[0.2369268 0.47862867 0.56888889 0.47862867 0.2369268];
end
A=(b-a)/2;
B=(b+a)/2;
result=0.0;
for i=1:n
result=result+GaussA(1,i)CalcuFunctionValue(AGaussP(1,i)+B);
end
result=A*result;
end

调用语句

clear all
close all
clc

a=0;%函数求积分的区间下限
b=1;%函数求积分的区间上限
h=(b-a)/4;%函数区间分份
A=0;
for i=1:4
A=A+gslrdjf( a+(i-1)h,a+ih ,2 );
end
A
A=0;
for i=1:4
A=A+gslrdjf( a+(i-1)h,a+ih ,3 );
end
A
A=0;
for i=1:4
A=A+gslrdjf( a+(i-1)h,a+ih ,4 );
end
A
A=0;
for i=1:4
A=A+gslrdjf( a+(i-1)h,a+ih ,5 );
end
A

运行截图

MATLAB高斯-勒让德积分求解相关推荐

高斯-勒让德积分——MATLABPython混合实现
高斯-勒让德积分前言一种很NB的数值积分.原理我们就不多说了,这里只介绍怎么直接拿来用. 核心问题是如何获取系数表. 参数介绍对于函数 f(x),积分区间为 a 到 b,求在这一区间上的定积分. ...
高斯勒让德(Gauss-legendre)求解多重积分(python,数值积分)
第四十四篇高斯勒让德求解多重积分多重积分在工程分析中,经常需要在一个面积或体积上对函数进行积分.多重积分的解析方法在有限的情况下是可能的,但在这一篇中使用数值积分去求解.一维的函数积分详见重复牛 ...
GSL科学计算库——计算高斯-勒让德积分
相关文章: Windows系统Qt5配置GSL科学计算库 Windows系统下Cygwin+Dev C ++ 配置GSL科学计算库假设计算下列积分: ∫0πexcos(x)dx\int_0^\pi ...
高斯-勒让德积分中不同阶数下最大高斯节点间距的关系
高斯-勒让德计算方法的高斯节点是勒让德多项式的根值,很多教材都把常用的根值和其对应的加权系数表示在如下的表格(P403 SEC. 7.5 GAUSS-LEGENDRE INTEGRATION (OPT ...
数据积分-牛顿科茨法与高斯勒让德法对比及示例
在结点xi上插值f的次数最多是n次的多项式 p(x)=∑f(xi)li(x) 拉格朗日插值多项式: li(x)= ∏(x-xj)/(xi-xj) 基于插值的数值积分: ∫f(x)dx≈∫p(x)dx= ...
matlab for求积分,MATLAB数值分析之数值积分（一）
[实验原理] 一.复合辛普森方法 1)复合辛普森方法基本思路把积分区间分成若干子区间,再在每个子区间上用低阶辛普森求积公式. 2)复合辛普森算法描述二.龙贝格方法 1)基本思路实质上是梯形公式的 ...
高斯勒让德求积公式matlab通用程序,数值分析——Gauss-Legendre 求积分(Matlab实现)...
2020 4.4 题目:编写Gauss求积法计算积分的程序(Gauss点数取1,2,3,4,5即可)并用于计算积分 ∫10sinxxdx∫01sin⁡xxdx\int_{0}^{1} \frac{\s ...
重复高斯勒让德法则(gauss-legendre)求积分(python,数值积分)
第四十一篇高斯勒让德法则牛顿-科特斯法则很方便使用,因为样本点在积分范围内均匀分布,而且权重系数很容易记住(至少到辛普森法则).详情可见重复牛顿-科特斯法则高斯法则采取样本点在积分范围内的最佳间 ...
高斯-勒让德求积分的Python程序
理论知识: 勒让德多项式及性质第四节高斯(Gauss)求积公式高斯-勒让德求积公式及Matlab实现目标函数求积区间为[-1,1]时代码: #高斯-勒让德求积公式 from sympy im ...

最新文章

热门文章