MATLAB多项式

MATLAB将多项式作为行向量执行，包括按幂次降序排列的系数。例如，方程P(X)= X 4 + 7×3 – 5×+ 9可以被表示为-

p = [1 7 0 -5 9];

多项式函数

多元拟合

给定两个向量x和y，命令a = polyfit(x，y，n)通过数据点(x i ，y i )拟合n阶多项式，并返回(n + 1)x in的幂系数行向量系数以x的幂的降序排列，即a = [a n a n-1 …a 1 a 0 ]。

多重

给定数据向量x和行向量中多项式的系数，命令y = polyval(a，x)评估数据点x i处的多项式并生成值y i ，使得

y i = a(1)x n i + a(2)x i (n-1) +…+ a(n)x + a(n + 1)。

因此，向量a的长度为n + 1，因此，求出的多项式的阶次为n。因此，如果a是5个元素长，则要求值的多项式将自动确定为四阶。

如果需要误差估计，则polyfit和polyval都使用最佳参数。

示例：直线(线性)拟合

以下数据是从旨在测量给定弹簧的弹簧常数的实验中收集的。弹簧悬挂有不同的质量m，并测量了弹簧从其未拉伸状态产生的相应偏移δ。

根据物理学，我们有F =kδ，这里F = mg。因此，我们可以从关系k = mg /δ中找到k。

在这里，我们将通过绘制实验数据，通过数据拟合最佳直线(我们知道δ和F之间的关系是线性的)，然后测量最佳拟合线的斜率来找到k。

m(g)

5.00         10.00         20.00         50.00         100.00

δ(mm)

15.5         33.07         53.39         140.24         301.03

通过数据拟合一条直线意味着我们想要找到多项式系数a 1和0 (一阶多项式)，使得1 x i + a 0给出y i的“最佳”估计。在步骤中，我们需要执行以下操作：

步骤1：找出系数a k 's：

a = polyfit(x，y，1)

步骤2：使用拟合多项式以更精细(更紧密分布)的x j评估y：

y_fitted = polyval(a，x_fine)

步骤3：绘制并查看。将给定输入绘制为点，将拟合数据绘制为线：

情节(x，y，'o'，x_fine，y_fitted);

以下脚本文件显示了通过弹簧实验的给定数据进行直线拟合并找到弹簧常数所包含的所有步骤。

并绘制下图：

最小二乘曲线拟合

最小二乘曲线拟合的过程可以简单地在MATLAB中实现，因为该技术产生了一组需要求解的线性方程。

大多数曲线拟合是多项式曲线拟合或指数曲线拟合(包括幂律，例如y = ax b )。

两个最常用的功能是：

y = ae bx

y = cx d

1.ln(y)= ln(a)+ bx或y = a 0 + a 1 x，其中y = ln(y)，a 1 = b和0 = ln？(a)。

2. ln(y)= ln(c)+ d ln(x)或y = a 0 + a 1 x ，其中y = ln(y)，a 1 = d和0 = ln？(c)。

现在我们可以在两种方法中仅使用一阶多项式来使用polyfit来确定未知常数。

涉及的步骤如下：

第1步：开发新输入：通过获取原始输入的对数，开发分配的新输入向量y和x。例如，要拟合y = ae bx类型的曲线，请创建

ybar = log(y)并保留x并拟合y = cx d类型的曲线，创建ybar = log(y)和xbar = log(x)。

步骤2：进行线性拟合：使用polyfit找出线性曲线拟合的系数a 0和1。

步骤3：绘制曲线：根据曲线拟合系数，计算原始常数的值(例如a，b)。根据获得的关系重新计算给定x处的y值，并绘制曲线和原始数据。

下表显示了真空泵的时间与压力变化读数。我们将通过数据拟合一条曲线P(t)= P 0 e -t /τ ，并确定未知常数P 0和τ。

t

0             0.5         1.0         5.0         10.0         20.0

P

760         625         528         85         14         0.16

通过记录双方关系的日志，我们可以

其中P = ln(P)， 1 =

-并且0 ＝ ln？(P 0 )。因此，一旦我们有1和0 ，我们就可以很容易地计算出P 0和τ。

示例：创建脚本文件

并绘制下图：

0

相关文章:MATLAB-循环类型 MATLAB-循环类型 在某 […]...

讨论MATLAB 讨论MATLAB MATLAB是MathWorks开& […]...

MATLAB简介 什么是MATLAB？ MATLAB是用于高性能Math 计算，可视化和编程环境的软件包。它提供了具有数百种内置 […]...

MATLAB中的工作区，变量和函数 MATLAB中的工作区，变量和函数 工作空间 工作区包含我们在MATLAB中工作时创建的所有变量。 每当我们为 […]...

MATLAB while循环 MATLAB while循环 while循环重复执行语句，而指定的语句为true。 句法：...

MATLAB 2D Bar() MATLAB Bar() 条形图是其中每个点都由垂直条或水平条表示的图。 句法...

MATLAB 2D Silhouette() MATLAB轮廓线() 等高线图显示了矩阵Z的等值线。 句法...

MATLAB特征值和特征向量 特征值和特征向量 方阵A的特征值和特征向量为标量λ和满足条件的非零向量v Av =λv 在这个等式中， A是一 […]...