我想从多元taylor费用中滤除特定条件。我需要过滤掉受不同影响的术语。
例如:f(a,b) = a/b
过滤受以下因素影响的每个扩展
(a-a1)
(b-b1)
(a-a1)
和(b-b1)
现在,我可以将术语分开并“手动”选择它们。
这是用于方差分析。
syms a b a1 b1 test test2;
n = 150 % Taylor Expansion degree
for i = 1:n
test(i) = evalin(symengine,strcat('mtaylor(a/b , [a = a1, b = b1], ', num2str(i), ')'));
if i == 1
test2(i) = test(i);
else
test2(i) = test(i)-test(i-1);
end
end
filter = children(test2(2) + (a-a1) + (b-b1));
filter = [ a, -a1, b, -b1, (a - a1)/b1, -(a1*(b - b1))/b1^2]
filter(5) = (a - a1)/b1
filter(6) = -(a1*(b - b1))/b1^2
例如,我希望提取与受(a-1)
影响的标准有关的单个术语。