evalit.mv3 (Jul 2013)

评估一个数学表达式和它相应的误差传播理论。

354下载

更新2013年7月22日

查看许可协议

你可曾想过有多少位数应该使用计算器的结果?

这个函数将为你做这些。

例如,让我们来计算地球的引力:

> > %符号输入:
> > sym G M R
> > g = g * M / R ^ 2;
> >
> > %数值输入:
顾> > = 6.67384 e-11;eGu = 0.00080 e-11;%牛顿引力常数及其误差。
> >我= 5.97 e24;高速= 3;% 3有效数字地球的质量。
> > Re = 6.371 e6;之前= 4;% 4有效数字地球的半径。
> >
> > %和计算:
> > evalit (g, {g M R},[顾我再保险],[eGu高速eRe))

所以我们得到:

EVALIT: FUNC (G、M R) = (G * M) / R ^ 2
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
值+ / -错误(参数的误差贡献)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
9.8160 + / - 0.0085 (0.0012 + 0.0082 + 0.0015)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

那就是:
g = (9.8160±0.0085) N /公斤
而不是
g = 9.816008178047202……N /公斤

正如您可以看到的,你也会得到每个因素的误差的贡献。
在这种情况下,地球的质量(2)结果与最大的错误:0.0082 N /公斤。
在G贡献最小的错误。

但最重要的是,所有错误从每个因子(G、M和R)有相同的数量级。

除了……
——这与矩阵输入。
——非function_handles(@)而不是符号表达式。
——它给数值和/或打印结果。
——它详尽检查输入它有助于输入错误。
——它使用误差最小的双精度浮点(因为v2)。
——输入数组(几乎):任何细胞,矩阵,向量(v3)。
——空,零,正或南错误更改为双精度浮点错误(v3)。
——它使用一阶误差约。但检查它们是否满足小于二阶的(v3)。
——它使用高阶误差约。当第一个是零(v3)。

享受它!
卡洛斯。

引用作为

卡洛斯阿德里安·巴尔加斯·阿奎莱拉(2023)。evalit.mv3 (Jul 2013)(//www.tianjin-qmedu.com/matlabcentral/fileexchange/42307 - evalit - m - v3 - 7 - 2013), MATLAB中央文件交换。检索

MATLAB版本兼容性
创建R2008b
兼容任何释放
平台的兼容性
窗户 macOS Linux

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!
版本 发表 发布说明
1.5.0.0

现在处理复杂的瓦尔斯。固定的错误当var dos不配合函数的参数。var和PDS输入数组可以是细胞或矩阵。其他一些小的变化。

1.2.0.0

现在处理双精度浮点数。其他一些小的变化。

1.0.0.0