error using ss2tf function

Question

Tenzing Thiley 2022-2-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1649430-error-using-ss2tf-function

评论： Paul 2022-2-14

clear all;
clc;
%% Declaration with parameter
I = 205113.07; % kg-m^2
m= 15118.35;   % kg
c = 3.511;     % m
rho = 1.225;   % kg/m^3
Tm = 49817.6;  % N
S = 37.16;     % m^2
%% Equation of motion(Longitudinal modes)
syms alpha etta delta theta q V 
f1 = ((etta*Tm*cos(alpha))/m) - ((0.5*rho*(V^2)*S*((0.0013*alpha^2)-(0.00438*alpha)+0.1423))./m)-(9.81*sin(theta-alpha));
f2 = q-(etta*Tm*sin(alpha)/(m*V))-(0.5*rho*V*S*((0.0751*alpha)+(0.0144*delta)+0.732)/m)+(9.81*cos(theta-alpha)/V);
f3 = (0.5*rho*V^2*S*c*(-0.00437*alpha-0.0196*delta-0.123*q-0.1885)/I);
f4 = q;
%% Linearization
A = jacobian([f1,f2,f3,f4],[V alpha q theta]);
B = jacobian([f1,f2,f3,f4],[etta delta]);
%% Equillibrium
V = 94.5077;  %m/s
etta = 0.58;
delta = -0.1678; %rad
theta = 0;
alpha = 0;
q = 0;
%% Evaluate System and Control Matrix
format shortG; 
fprintf("State Matrix = ");
A = eval(A)
fprintf("Control Matrix = ");
B = eval(B)

I want to convert State Space to Single input Single Output Transfer function Transfer function but I get error when using ss2tf function.

can any one help me out?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Paul 2022-2-14

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1649430-error-using-ss2tf-function#answer_895655

在 MATLAB Online 中打开

Define the nonlinear output equation y = h(x,u), take the jacobians wrt x and u, and evaluate them at the equilibirium point. These results will be C and D. Same exact approach to derive A and B from xdot = f(x,u). Once you have C and D then if you really need the transfer function using ss2tf()

[n,q]=ss2tf(A,B,C,D)

or better

hzpk = zpk(ss(A,B,C,D))

Or just leave the model in state space form

hss = ss(A,B,C,D);

unless you have to have the explicit transfer function for some reason.

5 个评论
显示 3更早的评论隐藏 3更早的评论

Tenzing Thiley 2022-2-14

transfer function needed was in this form

v(s)/etta(s) =

alpha(s)/etta(s) =

q(s)/etta(s) =

theta(s)/etta(s) =

v(s)/delta(s) =

alpha(s)/delta(s) =

q(s)/delta(s) =

theta(s)/delta(s) =

how can I achieve is using function ss2tf

Paul 2022-2-14

It looks like the output vector y = [v; alpha; q; theta] is exactly the same as the state variable vector, x. In this case, what should C and D be in this equation, where u = [etta; delta]?

y = C*x + D*u

Once you have C and D you can get all eight transfer functions at once using the second or third method in my comment above.

If sst2tf must be used, you can use

[num,den] = sst2tf(A,B,C,D,1)

to get the first four transfer functions and

[num,den] = sst2tf(A,B,C,D,2)

to get the second four transfer functions.

请先登录，再进行评论。