如何在符号数学 MATLAB 中提取系数。
例如,我有像这样的传递函数
H(S) = (a*S^2 + b*S^ + a*c*S + d*a) / (S^3 + b*c*S^2 + d*S^1 + a*d)
我想获取 H(s) 中“S”项的系数作为向量数组,如
num = [a b a*c d*a]
den = [1 b*c d a*d]
最佳答案
可以申请numden
提取分子和分母多项式:
[numexpr, denexpr] = numden(sym(H)) %// 'sym' makes sure that H is symbolic
然后您可以使用 coeffs
提取 S
的符号系数命令(记住首先对每个表达式应用 expand
以获得它们的多项式形式)。
但是,请注意 coeffs
仅返回非零系数。为了克服这个问题,我提出以下建议:
%// Extract numerator coefficients
[numcoef, numpow] = coeffs(expand(numexpr), S);
num = rot90(sym(sym2poly(sum(numpow))), 2);
num(num ~= 0) = coeffs(expand(numexpr), S);
%// Extract denominator coefficients
[dencoef, denpow] = coeffs(expand(denexpr), S);
den = rot90(sym(sym2poly(sum(denpow))), 2);
den(den ~= 0) = coeffs(expand(denexpr), S);
P.S: 你也可以申请sym2polys
tool来自 numexpr
和 denexpr
上的 MATLAB Exchange。
另请注意,系数向量中的最后一个元素与 S
的最高次幂相关联更为常见,因此此解决方案的结果将与您的顺序相反在你的问题中描述。
例子
%// Create symbolic function
syms a b S
H = b * S / (a + S^2)
[numexpr, denexpr] = numden(sym(H));
%// Extract numerator coefficients
[numcoef, numpow] = coeffs(expand(numexpr), S);
num = rot90(sym(sym2poly(sum(numpow))), 2);
num(num ~= 0) = coeffs(expand(numexpr), S);
%// Extract denominator coefficients
[dencoef, denpow] = coeffs(expand(denexpr), S);
den = rot90(sym(sym2poly(sum(denpow))), 2);
den(den ~= 0) = coeffs(expand(denexpr), S);
结果是:
num =
[0, b]
den =
[a, 0, 1]
关于matlab - 如何在MATLAB中提取传递函数的符号系数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/15498684/