Reduce block diagrams symbolically

18 次查看（过去 30 天）
Aleem Andrew
Aleem Andrew 2021-2-6
回答： Paul 2021-2-6
Can the code below be modified to find the transfer function in terms of symbolic functions g1 to g5 instead of numeric functions? I tried passing symbolic arrays as arguments to tf but I get an error message stating that the values of num and den must be numeric arrays. Is there a way to write code similar to the code below but with g1, g2, g3, g4, and g5 being passed as symbolic arguments?
g1 = tf(1,[1 0 0]);
num2 = 50; den2 = [1 1];
g2 = tf(num2,den2);
g3 = tf(2,[1,0]); g4 = tf([1 0],1);
g5 = tf(2,[1]);
sys1 = parallel(g4,-g5);
sys2 = feedback(g2,g3,-1);
sys4 = series(sys2,sys1);
sys5 = series(sys4,g1);
sys6 = feedback(sys5,[1])
answer = minreal(sys6)

请先登录，再回答此问题。

回答（1 个）

Paul
Paul 2021-2-6
Assuming you want final result in terms of the gi(s):
syms g1(s) g2(s) g3(s) g4(s) g5(s)
symsys1 = g4(s)-g5(s);
symsys2 = g2(s)/(1 + g3(s)*g2(s));
symsys4 = symsys2*symsys1;
symsys5 = symsys4*g1(s);
symsys6 = symsys5/(1 + symsys5);
symsys6 = simplify(symsys6,50);
[num,den] = numden(symsys6);
symsys6 = num/den;
symsys6
symsys6 =
(g1(s)*g2(s)*g4(s) - g1(s)*g2(s)*g5(s))/(g2(s)*g3(s) + g1(s)*g2(s)*g4(s) - g1(s)*g2(s)*g5(s) + 1)

类别

在 Help Center 和 File Exchange 中查找有关 Dynamic System Models 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by