Pade近似值
找到Padé的近似值sin (x)
.默认情况下,pade
返回一个三阶Padé近似。
信谊x pade (sin (x))
Ans = -(x*(7*x^2 - 60))/(3*(x^2 + 20))
如果不指定扩展变量,symvar
选择它。找到Padé的近似值sin (x) + cos (y)
.的symvar
功能选择x
作为扩展变量。
Syms x y pade(sin(x) + cosy)
ans = (- 7 * x ^ 3 + 3 * cos (y) * x ^ 2 + 60 * 60 * x + cos (y)) / (3 * (x ^ 2 + 20))
将扩展变量指定为y
.的pade
函数返回关于的Padé近似y
.
pade (sin (x) + cos (y), y)
ans = (12 * sin (x) + y ^ 2 * sin (x) - 5 * y ^ 2 + 12) / (y ^ 2 + 12)
找出…的价值谭(3 *π/ 4)
.使用pade
找到Padé的近似值谭(x)
然后代入潜艇
找到谭(3 *π/ 4)
.
Syms x f = tan(x);P = pade (f);y =潜艇(P, x, 3 *π/ 4)
Y =(9*²/16 - 15)/(4*(9*²)/8 - 5)
使用vpa
转换y
转换成一个数值。
vpa (y)
ans = -1.2158518789569086447244881326842
您可以通过增加订单来提高Padé近似的准确性。如果膨胀点是极点或零点,也可以通过设置来提高精度OrderMode
来相对
.的OrderMode
如果膨胀点不是极点或零点,选项没有影响。
找到Padé的近似值谭(x)
使用pade
扩展点为0
和订单
的[1]
.找出…的价值谭(1/5)
用Padé近似代入潜艇
,并使用vpa
转换1/5
转换成一个数值。
信谊x侯= pade (tan (x), x, 0,“秩序”,[1])侯=潜艇(侯,x, vpa (1/5))
P11 = x P11 = 0.2
用减法求近似误差赛
的实际值谭(1/5)
.
y = tan (vpa (1/5));错误= y - p11
错误= 0.0027100355086724833213582716475345
通过增加使用的顺序来增加Padé近似的准确性订单
.集订单
来(2 - 2)
,找出错误。
第22位= pade (tan (x), x, 0,“秩序”,2[2])第22位=潜艇(第22位,x, vpa (1/5));错误= y - p22
P22 = -(3*x)/(x^2 - 3) error = 0.0000073328059697806186555689448317799
精度随阶数的增加而增加。
如果膨胀点为极点或零,则Padé近似的精度会降低。设置OrderMode
选项相对
补偿精度的下降。有关详细信息,请参见Pade近似值.因为棕褐色
函数的at为零0
,设置OrderMode
来相对
提高准确性。如果膨胀点不是极点或零,这个选项没有效果。
p22Rel = pade(tan(x),x,0,'Order',[2 2],'OrderMode','relative');错误= y - p22Rel
p22Rel = (x*(x^2 - 15))/(3*(2*x^2 - 5)) error = 0.0000000084084014806113311713765317725998 . p22Rel = (x*(x^2 - 15))/(3*(2*x^2 - 5)
如果膨胀点是极点或零点,则精度增加OrderMode
被设置为相对
.
画出exp (x)
以及它的Padé近似值[1]
通过(4 - 4)
.使用轴
专注于感兴趣的领域。从图中可以看出,准确率随着Padé近似的阶数的增加而增加。
信谊xexpr = exp (x);持有在网格在为I = 1:4 fplot(expr - pade(expr,“秩序”,我)结束轴([-4 4 -4 4])的订单(1),的订单(2,2),“订单(3 3)”,“订单(4,4)”,...“位置”,“最佳”)标题(“exp(x)与Pade近似的区别”) ylabel (“错误”)
如果你同时使用第三个参数一个
和ExpansionPoint
要指定扩展点,通过ExpansionPoint
生活的全部。
的参数一个1、……bn的级数展开,使得Padé近似的级数展开与f可能的最大阶。
扩展点±∞和±我∞是不允许的。
当pade
如果找不到Padé近似函数,则返回函数调用。
为pade
的泰勒或劳伦级数展开f必须存在于膨胀点。