# NumericalAnalysisCode

**Repository Path**: iamttp/NumericalAnalysisCode

## Basic Information

- **Project Name**: NumericalAnalysisCode
- **Description**: 数值分析第二版 Timothy Sauer Matlab/Octave测试代码
- **Primary Language**: Matlab
- **License**: Not specified
- **Default Branch**: master
- **Homepage**: None
- **GVP Project**: No

## Statistics

- **Stars**: 1
- **Forks**: 0
- **Created**: 2020-03-27
- **Last Updated**: 2021-09-25

## Categories & Tags

**Categories**: Uncategorized

**Tags**: None

## README

# 数值分析 Matlab/Octave代码

> 验证平台：`https://octave-online.net/`
>
> 编写平台：Ubuntu VSCode
>
> 参考书籍：Numerical Analysis(Second Edition) Timothy Sauer

## 第0章 基础知识（nest.m）

### 程序0.1 嵌套乘法：

* 使用霍纳方法计算

* 调用示例：

```matlab
    nest(4,[-1 5 -3 3 2],1/2)
```

* 最终结果：

```matlab
    ans =  1.2500
```

## 第1章 求解方程  (bisect.m  fpi.m  f.m)

### 程序1.1 二分法

* 调用示例：

```matlab
    f = @(x) x^3+x-1
    xc = bisect(f,0,1,0.00005)
```

* 最终结果：

```matlab
    xc =  0.68234
```

### 程序1.2 不动点迭代

* 调用示例：

```matlab
    g = @(x) (1+2.*x^3)/(1+3.*x^2)
    xc = fpi(g,0.5,10)
```

* 最终结果：

```matlab
    xc =  0.68233
```

### 事实验证1 Stewart平台运动学

* Stewart平台运动学

> 参数意义：
>
> p1\p2\p3为臂长
>
> v\L1\L2\L3定义三角形的形状(等腰直角）
>
> x1/x2/y2定义支杆位置
>

* 使用方法：

```matlab
    theta = -pi:0.05:pi;
    fx = f(theta);
    plot(theta,fx)
    grid on
```

* 结果：

```matlab
    解为pi/4和-pi/4
```

![avatar](picture/f.png)
![avatar](picture/f3.png)

## 第2章 方程组 (Euler_Bernoulli.m)

## 第3章 插值 (newtdd.m  clickinterp.m)

### 程序3.1 牛顿差商插值方法

* 调用示例

```matlab
    x0 = [0 2 3]
    y0 = [1 2 4]
    c = newtdd(x0,y0,3);
    x = 0:0.01:4
    y = nest(2,c,x,x0)
    plot(x0,y0,'o',x,y)
```

![avatar](picture/newtdd.png)

## 第4章 最小二乘  (test.m)

### 程序4. 抛物线拟合

* 调用示例

```matlab
    直接运行test.m
```

![avatar](picture/test.png)

## 第5章 数值微分和积分

## 第6章 常微分方程  (euler.m)

### 程序6.1 求解初值问题的欧拉方法

* 调用示例

```matlab
euler([0,1],1,10);
```

![avatar](picture/euler.png)