使用Chebfun求解Blasius方程(一)

Chebfun的特点:

1. 基于Chebyshev展开，展开项数由机器精度自适应控制;

2. 将符号计算和数值计算结合，以处理数值的速度处理函数；

3. 在Matlab上实现，将Matlab处理向量和矩阵的命令重载，以处理函数和算子；

4. 基于Newton迭代法求解非线性微分方程;

5. 使用自动微分技术计算Frechet导数；

6. Chebop的实现利用了谱方法和惰性求值的思想

7. 能表示具有可去奇点的函数

 

Chebfun仍然在持续开发中，后期对较复杂的求解过程进行了封装，使得用户将更多的精力放在自己的问题上。下面从三个层面使用chebfun系统求解Blasius方程，“麻雀虽小五脏俱全”，以期对Chebfun求解非线性边值问题有较深入的理解。

 

Blasius方程模拟了一个半无限长平板上的二维粘性层流流动，控制方程是

边界条件是