帕德逼近matlab算法,帕德逼近算法.doc

《MATLAB程序设计实践》课程作业

一、用MATLAB编程实现“帕德逼近”的科学计算算法，及举例应用。

1)帕德逼近算法说明如下：

帕德逼近是一种有理分式逼近，逼近公式如下：

大量实验表明，当L+M为常数时，取L=M，帕德逼近精确度最好，而且速度最快。此时，分子与分母中的系数可通过以下方式求解。

首先，求解线形方程Aq=b，得到(…)的值，其中

，，

然后，通过下式求出的值。

注意，函数的帕德逼近不一定存在。

在MATLAB中编程实现的帕德逼近法函数为：Pade。

功能：用帕德形式的有理分式逼近已知函数。

调用格式：f=Pade(y,n)或f=Pade(y,n,x0)。

其中，y为已知函数；

n为帕德有理分式的分母多项式的最高次数；

x0为逼近点的x坐标；

f为求得的帕德有理分式或在x0处的逼近值。

2)程序源代码如下：

①在m文件中编写实现函数的Pade逼近的代码如下：

function f=Pade(y,n,x0)

%用帕德形式的有理分式逼近已知函数

%已知函数：y

%帕德有理分式的分母多项式的最高次数：n

%逼近点的坐标：x0

%求得的帕德有理分式或在x0处的逼近值：f

syms t;

A=ze