How can I solve this expression for Qi_s/thetha_s = ....?
1 view (last 30 days)
Show older comments
clc,clear
syms theta_s s C Qi_s R theta_a
Transfer=solve(Qi_s/C-theta_s/(R*C)+theta_a/(R*C)==s*theta_s,...)
Hi, I was wondering if you can help me with this exercise, I must get my transfer function, I need to obtain "Qi_s/thetha_s" and the expression, can anyone help me please. I think solve can't do this so I am looking for a new tool that could help me.
2 Comments
Answers (1)
Sulaymon Eshkabilov
on 2 Feb 2023
Edited: Torsten
on 2 Feb 2023
To get a transfer function, you should divide both sides of the equation by theta_s, and you will get this:
% Qi_s/C-theta_s/(R*C)+theta_a/(R*C)==s*theta_s
% Qi_s/(C*theta_s)-1/(R*C)+theta_a/(R*C) = s
% Qi_s/(C*theta_s) = s+(1/R*C)-theta_a/(R*C)
% Qi_s/theta_s = s/C + (1-theta_a)/(R*C^2)
% Qi_s/theta_s = (s*R*C+1-theta_a)/(R*C^2)
% TF = (s*R*C+1-theta_a)/(R*C^2)
% Analytical solution of the found TF and its reverse using the inverse and
% Laplace transforms
% Some numerical values assigned to the unknowns:
R = 100;
C = 5;
theta_a = pi;
syms s
Solution1 = ilaplace((s*R*C+1-theta_a)/(R*C^2))
fplot(Solution1, [-10, 10])
% Note that this is an improper TF. Presuming it to be reverse:
% TF = (R*C^2)/(s*R*C+1-theta_a)
% Some numerical values
R = 100;
C = 5;
theta_a = pi;
s = tf('s');
TF = (R*C^2)/(s*R*C+1-theta_a);
% Step response
step(TF, 10)
% analytical solution of the found TF and its reverse using the inverse and
% Laplace transforms
syms s
Solution2 = ilaplace((R*C^2)/(s*R*C+1-theta_a))
fplot(Solution2, [-10, 10])
1 Comment
Torsten
on 2 Feb 2023
Edited: Torsten
on 2 Feb 2023
@Jhoao Venegas comment moved here:
I think in the second line there should be: %Qi_s/(C*theta_s)-1/(R*C)+theta_a/(theta_s*R*C) = s
Although, thanks.
See Also
Categories
Find more on Calculus in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!