matlab黎曼和定义求解定积分的代码

七十二五

已于 2024-08-29 08:58:59 修改

阅读量140

点赞数 2

分类专栏： matlab代码文章标签： matlab 经验分享算法青少年编程矩阵

于 2024-08-29 08:45:00 首次发布

本文链接：https://blog.csdn.net/byllalala/article/details/141644467

版权

matlab代码专栏收录该内容

21 篇文章 0 订阅

订阅专栏

定积分值的定义及代码：

$\int_{a}^{b}f(x)dx=\frac{1}{n}\sum_{i=1}^{n}f\left ( i\frac{{b-a}}{n}\right )=F$

clear;clc;
n=1000;      %n为整数   可以修改
a=1;   %定积分下限为1
b=2;   %定积分上限为2    可以修改
hx=(b-a)/(n)
x=a:hx:b;

f=@(X)(X*X);      %以x*x为例   求x*x的定积分结果    

F=0;              %叠加定积分的初值为0

for i=1:n
F=F+f(x(i))/n;
end

%F即为x从a到b上函数f的定积分结果，提高n的值可以提高精度

变上限的积分的定义及代码：

$\int_{a}^{T}f(x)dx=\frac{1}{n}\sum_{i=1}^{n}f\left ( i\frac{T-a}{n} \right )=F(T)$

clear;clc;
n=1000;      %n为整数   可以修改
a=1;   %定积分下限为1
b=2;   %定积分上限为2    可以修改
hx=(b-a)/(n)
x=a:hx:b;
T=x;

f=@(X)(X*X);      %以x*x为例   求x*x的定积分结果    

F1=0;              %叠加定积分的初值为0
F=zeros(n+1,1)
for i=1:n
F1=F1+f(x(i))/n;
F(i+1,1)=F1;
end

%F即为x从a到T上函数f的定积分结果F(T)，提高n的值可以提高精度