Halo

A magic place for coding

0%

介绍

  在日常生活的很多测试中,我们得到的数据都是由不同频率不同振幅的波形叠加起来的数据。在计算的时候需要逼近这种波形,选取具有周期性的三角函数作为基函数是合适的。对于计算傅里叶逼近系数问题,都可以统一地归结为:
$$
c_j = \sum^{N-1}_{k=0}x_k\omega^{kj}_N,j=0,1,\dots,N-1
$$,这就是N点DFT
其中${x_k}^{N-1}_0$为已知的输入数据点(采样点),${c_j}^{N-1}_0为输出数据$。

Read more »

介绍

  我们在以前可能会在直线拟合的时候使用过最小二乘法,常用于处理实现数据,为了获得一条直线或者曲线。而在这里我们对更一般的情况做讨论。
  对于给定的超定线性方程组$Ax = b$,直接通过矩阵求解是困难的。这里就用到了逼近的思想。假设上述方程存在一个噪声$e$,则可以写成$b-Ax = e$,接下来要做的就是求解出能够极小化这个噪声的$x$。这就是递推最小二乘法需要解决的问题。

Read more »

介绍

  在前面我们提到了牛顿法求解非线性方程,但是牛顿法有一个明显的缺点,就是在每一步迭代中都要计算函数值和一阶导数,这个计算量是巨大而且困难的,因此我们希望减少这个计算,对牛顿法进行了改进–简化牛顿法

Read more »

介绍

  当我们对一个非线性方程进行求根的时候,一般采用的是迭代法,通过迭代法的收敛性质,最终求得一个逼近精确解的近似解。这里介绍的是其中一种方法–二分法

Read more »

介绍

  在实际生活中,我们常常需要通过观测来获取某些数据,这些数据有时候是连续的,有时候却是离散的,对于这些数据,我们希望找到一个方便运算的函数$P(x)$能接近真实的函数$f(x)$,求得这个$P(x)$的过程就叫插值。插值的方法有很多种,这里重点介绍拉格朗日插值

Read more »

Introduction

  在之前我实现了求解低阶矩阵的一些算法。考虑线性方程组Ax=b,其中A为非奇异矩阵,当A是低阶稠密矩阵时,可以使用之前提到的中的高斯消元法以及列主
元消元法(详情请参考高斯消元法之讲解与代码实现)来进行求解。但是实际生活中工程技术产生的矩阵都是大型的稀疏矩阵(阶数很大,但零元素很多),当A是高阶稀疏矩阵时,就可以利用迭代法来进行求解。迭代法利用了A中零元素较多的特点,这对于计算机的存储和运算是很有利的。

Read more »