55 lines
945 B
Markdown
55 lines
945 B
Markdown
|
---
|
|||
|
|
title: Matlab使用Simpson辛普森法计算数值积分
|
|||
|
|
date: 2021-11-12 00:00:00
|
|||
|
|
tags:
|
|||
|
|
- Matlab
|
|||
|
|
- 编程
|
|||
|
|
category: 技术
|
|||
|
|
---
|
|||
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|
<!-- more -->
|
|||
|
|
|
|||
|
|
直接附代码
|
|||
|
|
|
|||
|
|
```matlab
|
|||
|
|
simpson.m
|
|||
|
|
function [out] = simpson(f,n,a,b)
|
|||
|
|
%SIMPSON 此处显示有关此函数的摘要
|
|||
|
|
% 此处显示详细说明
|
|||
|
|
|
|||
|
|
h=(b-a)/n;
|
|||
|
|
|
|||
|
|
s1=fi(f,1,h,a)+fi(f,n+1,h,a);
|
|||
|
|
|
|||
|
|
s2=0;
|
|||
|
|
for index=3:2:n-1
|
|||
|
|
s2=s2+2*fi(f,index,h,a);
|
|||
|
|
end
|
|||
|
|
|
|||
|
|
s4=0;
|
|||
|
|
for index=2:2:n
|
|||
|
|
s4=s4+4*fi(f,index,h,a);
|
|||
|
|
end
|
|||
|
|
|
|||
|
|
out=h/3*(s1+s2+s4);
|
|||
|
|
end
|
|||
|
|
|
|||
|
|
|
|||
|
|
function [out] = fi(f,i,h,a)
|
|||
|
|
|
|||
|
|
out=f(a+(i-1)*h);
|
|||
|
|
end
|
|||
|
|
```
|
|||
|
|
|
|||
|
|
调用：
|
|||
|
|
|
|||
|
|
```matlab
|
|||
|
|
simpson(<函数句柄>, <n>, <下限>, <上限>)
|
|||
|
|
```
|
|||
|
|
|
|||
|
|
> 参考文章
|
|||
|
|
>
|
|||
|
|
> [matlab用辛普森公式求积分_积分近似计算之辛普森公式_weixin_39726131的博客-CSDN博客](https://blog.csdn.net/weixin_39726131/article/details/110226598)
|