当前位置：文档之家 › matlab代码--数值积分

matlab代码--数值积分

  • 格式：pdf
  • 大小：159.53 KB
  • 文档页数：12

数值积分

1.CombineTraprl复合梯形公式求积分 function [I,step] = CombineTraprl(f,a,b,eps)

%f 被积函数

%a,b 积分上下限

%eps 精度

%I 积分结果 %step 积分的子区间数

if(nargin ==3)

eps=1.0e-4;

end

n=1;

h=(b-a)/2;

I1=0;

I2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h;

while abs(I2-I1)>eps n=n+1;

h=(b-a)/n;

I1=I2;

I2=0;

for i=0:n-1 x=a+h*i;

x1=x+h;

I2=I2+(h/2)*(subs(sym(f),findsym(sym(f)),x)+subs(sym(f),findsym(sym(f)),x1));

end end

I=I2;

step=n;

2.IntSimpson用辛普森系列公式求积分

function [I,step] = IntSimpson(f,a,b,type,eps)

%type = 1 辛普森公式

%type = 2 辛普森3/8公式 %type = 3 复合辛普森公式

if(type==3 && nargin==4)

eps=1.0e-4; %缺省精度为0.0001

end

I=0;

switch type

case 1, I=((b-a)/6)*(subs(sym(f),findsym(sym(f)),a)+...

4*subs(sym(f),findsym(sym(f)),(a+b)/2)+...

subs(sym(f),findsym(sym(f)),b));

step=1;

case 2,

I=((b-a)/8)*(subs(sym(f),findsym(sym(f)),a)+...

3*subs(sym(f),findsym(sym(f)),(2*a+b)/3)+ ...

3*subs(sym(f),findsym(sym(f)),(a+2*b)/3)+subs(sym(f),findsym(sym(f)),b)); step=1;

case 3,

n=2;

h=(b-a)/2; I1=0;

I2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h;

while abs(I2-I1)>eps

n=n+1;

h=(b-a)/n; I1=I2;

I2=0;

for i=0:n-1

x=a+h*i;

x1=x+h; I2=I2+(h/6)*(subs(sym(f),findsym(sym(f)),x)+...

4*subs(sym(f),findsym(sym(f)),(x+x1)/2)+...

subs(sym(f),findsym(sym(f)),x1));

end

end I=I2;

step=n;

end

3.NewtonCotes用牛顿－科茨系列公式求积分

function I = NewtonCotes(f,a,b,type)

%type = 1 科茨公式

%type = 2 牛顿－科茨六点公式 %type = 3 牛顿－科茨七点公式

I=0;

switch type case 1,

I=((b-a)/90)*(7*subs(sym(f),findsym(sym(f)),a)+...

32*subs(sym(f),findsym(sym(f)),(3*a+b)/4)+...

12*subs(sym(f),findsym(sym(f)),(a+b)/2)+...

32*subs(sym(f),findsym(sym(f)),(a+3*b)/4)+7*subs(sym(f),findsym(sym(f)),b));

case 2,

I=((b-a)/288)*(19*subs(sym(f),findsym(sym(f)),a)+...

75*subs(sym(f),findsym(sym(f)),(4*a+b)/5)+... 50*subs(sym(f),findsym(sym(f)),(3*a+2*b)/5)+...

50*subs(sym(f),findsym(sym(f)),(2*a+3*b)/5)+...

75*subs(sym(f),findsym(sym(f)),(a+4*b)/5)+19*subs(sym(f),findsym(sym(f)),b));

case 3,

I=((b-a)/840)*(41*subs(sym(f),findsym(sym(f)),a)+...

216*subs(sym(f),findsym(sym(f)),(5*a+b)/6)+...

27*subs(sym(f),findsym(sym(f)),(2*a+b)/3)+...

272*subs(sym(f),findsym(sym(f)),(a+b)/2)+... 27*subs(sym(f),findsym(sym(f)),(a+2*b)/3)+...

216*subs(sym(f),findsym(sym(f)),(a+5*b)/6)+41*subs(sym(f),findsym(sym(f)),b));

end

4.IntGauss用高斯公式求积分 function I = IntGauss(f,a,b,n,AK,XK)

if(n<5 && nargin == 4)

AK = 0;

XK = 0;

else XK1=((b-a)/2)*XK+((a+b)/2);

I=((b-a)/2)*sum(AK.*subs(sym(f),findsym(f),XK1));

end

ta = (b-a)/2; tb = (a+b)/2;

switch n

case 0,

I=2*ta*subs(sym(f),findsym(sym(f)),tb);

case 1,

I=ta*(subs(sym(f),findsym(sym(f)),ta*0.5773503+tb)+...

subs(sym(f),findsym(sym(f)),-ta*0.5773503+tb));

case 2,

I=ta*(0.55555556*subs(sym(f),findsym(sym(f)),ta*0.7745967+tb)+...

0.55555556*subs(sym(f),findsym(sym(f)),-ta*0.7745967+tb)+...

0.88888889*subs(sym(f),findsym(sym(f)),tb));

case 3,

I=ta*(0.3478548*subs(sym(f),findsym(sym(f)),ta*0.8611363+tb)+...

0.3478548*subs(sym(f),findsym(sym(f)),-ta*0.8611363+tb)+...

0.6521452*subs(sym(f),findsym(sym(f)),ta*0.3398810+tb)...

+0.6521452*subs(sym(f),findsym(sym(f)),-ta*0.3398810+tb));

case 4,

I=ta*(0.2369269*subs(sym(f),findsym(sym(f)),ta*0.9061793+tb)+...

0.2369269*subs(sym(f),findsym(sym(f)),-ta*0.9061793+tb)+... 0.4786287*subs(sym(f),findsym(sym(f)),ta*0.5384693+tb)...

+0.4786287*subs(sym(f),findsym(sym(f)),-ta*0.5384693+tb)+...

0.5688889*subs(sym(f),findsym(sym(f)),tb));

end

5.IntGaussLada用高斯拉道公式求积分

function I = IntGaussLada(f,a,b,n,AK,XK)

if(n<6 && nargin == 4) AK = 0;

XK = 0;

else

XK1=((b-a)/2)*XK+((a+b)/2);

I=((b-a)/2)*((2/n/n)*subs(sym(f),findsym(sym(f)),a)+... sum(AK.*subs(sym(f),findsym(sym(f)),XK1)));

end

ta = (b-a)/2;

tb = (a+b)/2; switch n

case 2,

I=ta*0.5*subs(sym(f),findsym(sym(f)),ta*(-1)+tb)+...

1.5*ta*subs(sym(f),findsym(sym(f)),ta*(1/3)+tb);

case 3,

合集下载

计算方法-数值积分实验

计算方法-数值积分实验

实验二 数值积分实验

一. 实验目的

（1）熟悉数值积分与数值微分方法的基本思想，加深对数值积分与数值微分方法的理解。

（2）熟悉Matlab编程环境，利用Matlab实现具体的数值积分与数值微分方法。

二. 实验要求

用Matlab软件实现复化梯形方法、复化辛甫生方法、龙贝格方法和高斯公式的相应算法，并用实例在计算机上计算。

三. 实验内容

1. 实验题目

已知xexfx4sin1)(的数据表

分别编写用复化梯形法、复化辛甫生公法、龙贝格法、三点高斯法求积分10)(dxxfI近似值的计算机程序。

A. 复化梯形法:

a． 编写文件Trapezoid.m,代码如下所示：

x 0.00 0.25 0.50 0.75 1.00

f(x) 1.00000 1.65534 1.55152 1.06666 0.72159 b． 编写文件f2.m:

c.运行：

B．复化辛甫生公法

a． 编写文件FSimpson.m,代码如下所示：

b．编写文件f2.m:

function f=f2(x)

f=1+exp(-x).*sin(4*x);

c.运行：

C．龙贝格法

a． 编写文件Romberg.m,代码如下所示：

b.运行：

D.三点高斯法

a.编写文件TGauss.m文件，如下所示： b． 运行：

2. 设计思想

要求针对上述题目，详细分析每种算法的设计思想。

总体的思想是化复杂为简单的重复

A． 复化梯形法使用直接法，通过递归，缩减规模；

B． 复化辛甫生也是使用直接法，根据公式直接进行编程，通过递归缩减规模；

C． 龙贝格算法应该在做了的几个中最体现了“化复杂为简单的重复”的思想，多个循环通过变量的适当递增，和一个for循环语句来实现，循环主体只有一句话，但确是整个程序中的亮点和难点；

D． 三点高斯法直接通过一条简单的公式来编写程序，难度不大；

四. 实验体会

对实验过程进行分析总结，对比不同方法的精度，指出每种算法的设计要点及应注意的事项，以及自己通过实验所获得的对数值积分方法的理解。

Matlab中常用的数值计算方法

Matlab中常用的数值计算方法

Matlab中常用的数值计算方法

数值计算是现代科学和工程领域中的一个重要问题。Matlab是一种用于数值计算和科学计算的高级编程语言和环境，具有强大的数值计算功能。本文将介绍Matlab中常用的数值计算方法，包括数值积分、数值解微分方程、非线性方程求解和线性方程组求解等。

一、数值积分

数值积分是通过数值方法来近似计算函数的定积分。在Matlab中，常用的数值积分函数是'quad'和'quadl'。'quad'函数可以用于计算定积分，而'quadl'函数可以用于计算无穷积分。

下面是一个使用'quad'函数计算定积分的例子。假设我们想计算函数f(x) = x^2在区间[0, 1]上的定积分。我们可以使用如下的Matlab代码：

```

f = @(x) x^2;

integral = quad(f, 0, 1);

disp(integral);

```

运行这段代码后，我们可以得到定积分的近似值，即1/3。

二、数值解微分方程

微分方程是描述自然界各种变化规律的数学方程。在科学研究和工程应用中，常常需要求解微分方程的数值解。在Matlab中，可以使用'ode45'函数来求解常微分方程的数值解。'ode45'函数是采用基于Runge-Kutta方法的一种数值解法。 下面是一个使用'ode45'函数求解常微分方程的例子。假设我们想求解一阶常微分方程dy/dx = 2*x，初始条件为y(0) = 1。我们可以使用如下的Matlab代码：

```

fun = @(x, y) 2*x;

[x, y] = ode45(fun, [0, 1], 1);

plot(x, y);

```

运行这段代码后，我们可以得到微分方程的数值解，并绘制其图像。

三、非线性方程求解

非线性方程是指方程中包含非线性项的方程。在很多实际问题中，我们需要求解非线性方程的根。在Matlab中，可以使用'fsolve'函数来求解非线性方程的根。

matlab高斯数值积分

matlab高斯数值积分

matlab高斯数值积分

在MATLAB中，可以使用函数`integral`来计算高斯数值积分。`integral`函数可以用于计算定积分，也可以用于计算无穷积分。以下是使用`integral`函数计算高斯数值积分的示例代码：

```matlab

fun = @(x) exp(-x.^2); % 定义被积函数

a = -inf; % 积分下限（可以是负无穷）

b = inf; % 积分上限（可以是正无穷）

result = integral(fun, a, b);

在上面的示例中，我们定义了一个被积函数`fun`，然后使用`integral`函数计算了这个函数在负无穷到正无穷之间的高斯数值积分。计算结果会保存在变量`result`中。如果需要计算定积分，只需将积分上下限`a`和`b`设置为具体的数值即可：

```matlab

fun = @(x) exp(-x.^2); % 定义被积函数

a = 0; % 积分下限

b = 1; % 积分上限

result = integral(fun, a, b);

在实际使用中，还可通过设置`'RelTol'`和`'AbsTol'`选项来控制积分的精度。

matlab定义积分函数

matlab定义积分函数

matlab定义积分函数

MATLAB 是一种强大的科学计算软件，它有丰富的数学函数库，可以实现数值积分计算，便于工程和科学计算。在计算机编程中，数值积分计算是一个非常重要的部分。MATLAB 中定义积分的函数有：`integral`、`quad`、`quadl`、`quadgk` 等。本文将重点介绍 `integral` 函数的定义及使用。

一、 `integral` 函数的定义

`integral` 函数是 MATLAB 中求解一维定积分的函数，其函数调用格式如下：

```matlab

Q = integral(fun, a, b)

```

`fun` 是被积函数句柄 (函数句柄是指在 MATLAB 中定义的函数的名称)，变量 `a`

和 `b` 是积分区间的下限和上限，`Q` 是积分的结果。

二、 `integral` 函数的使用

1. 用户自定义函数

当我们需要求解用户自定义函数的积分时，可以采用以下方法：

1) 我们需要在 MATLAB 中定义被积函数。我们定义一个函数 $f(x)=\sqrt{1-x^2}$，用以下语句来实现：

```matlab

function y = myfun(x)

y = sqrt(1 - x.^2);

end

```

`myfun` 就是被积函数的名称。

2) 然后，我们就可以使用 `integral` 函数来计算积分的值。要求 $f(x)$ 在区间

$[-1, 1]$ 的积分值，我们可以用以下代码：

```matlab Q = integral(@myfun, -1, 1)

```

`@myfun` 表示把 `myfun` 函数的句柄传递给 `integral` 函数 ，`-1` 和 `1` 分别表示积分区间的下限和上限，`Q` 表示积分值。