3/2/2015
AE_157_Homework_2_Problem_3_c
Contents
deletesallthevariablespace,closeallfigures,clearthescreenandopenanewpage.
MATRIX
Parta:Ifwewouldliketohaveboth(x1,x2)statesasoutputs,outputmatrixCshouldbe
Partb:Ifwedonotwantanytypeofinputeffectontheoutput,theshape(size)andcontentofthestatetransmission
matrix(D)shouldbe
Partc:statespacerepresentationofthesystem
Partd:givenmatrixrepresentations,obtain4transferfunctionswhichdefineinputoutputdynamics,fromeachofthe
inputstoeachoftheoutputs.(NamelyY1/U1,Y1/U2,...etc.)
Parte:stepcommand,evaluatethestepresponseofthesystem.Commentonyourresults,whatdoyouobserve?
WhataretheMpandsteadystatecharacteristics.
Partf:impulseresponse,evaluatetheimpulseresponseofthesystem.
END
deletesallthevariablespace,closeallfigures,clearthescreenandopenanewpage.
clearall,closeall,clc;
MATRIX
A=[0.81.5;5.70]
B=[1.11.1;1.10]
A=
0.80001.5000
5.70000
B=
1.10001.1000
1.10000
Parta:Ifwewouldliketohaveboth(x1,x2)statesasoutputs,outputmatrixCshouldbe
C=[10;01]%eye(2)
C=
10
01
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
1/11
3/2/2015
AE_157_Homework_2_Problem_3_c
Partb:Ifwedonotwantanytypeofinputeffectontheoutput,theshape(size)andcontentof
thestatetransmissionmatrix(D)shouldbe
D=[00;00]%zeros(2)
D=
00
00
Partc:statespacerepresentationofthesystem
StateSpace=ss(A,B,C,D)
StateSpace=
a=
x1x2
x10.81.5
x25.70
b=
u1u2
x11.11.1
x21.10
c=
x1x2
y110
y201
d=
u1u2
y100
y200
Continuoustimestatespacemodel.
Partd:givenmatrixrepresentations,obtain4transferfunctionswhichdefineinputoutput
dynamics,fromeachoftheinputstoeachoftheoutputs.(NamelyY1/U1,Y1/U2,...etc.)
[numerator1,denominator1]=ss2tf(A,B,C,D,1)
[numerator2,denominator2]=ss2tf(A,B,C,D,2)
Y1U1=tf(numerator1(1,:),denominator1)
Y2U1=tf(numerator1(2,:),denominator1)
Y1U2=tf(numerator2(1,:),denominator2)
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
2/11
3/2/2015
AE_157_Homework_2_Problem_3_c
Y2U2=tf(numerator2(2,:),denominator2)
numerator1=
01.10001.6500
01.10007.1500
denominator1=
1.00000.80008.5500
numerator2=
01.10000
006.2700
denominator2=
1.00000.80008.5500
Y1U1=
1.1s1.65
s^2+0.8s+8.55
Continuoustimetransferfunction.
Y2U1=
1.1s+7.15
s^2+0.8s+8.55
Continuoustimetransferfunction.
Y1U2=
1.1s
s^2+0.8s+8.55
Continuoustimetransferfunction.
Y2U2=
6.27
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
3/11
3/2/2015
AE_157_Homework_2_Problem_3_c
s^2+0.8s+8.55
Continuoustimetransferfunction.
Parte:stepcommand,evaluatethestepresponseofthesystem.Commentonyourresults,
whatdoyouobserve?WhataretheMpandsteadystatecharacteristics.
figure
step(Y1U1)
Info=stepinfo(Y1U1)
title('StepResponseforY1/U1')
Info=
RiseTime:0.1688
SettlingTime:10.1213
SettlingMin:0.4439
SettlingMax:0.0309
Overshoot:130.0284
Undershoot:101.0928
Peak:0.4439
PeakTime:1.3967
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
4/11
3/2/2015
AE_157_Homework_2_Problem_3_c
figure
step(Y2U1)
Info=stepinfo(Y2U1)
title('StepResponseforY2/U1')
Info=
RiseTime:0.3389
SettlingTime:9.7631
SettlingMin:0.4494
SettlingMax:1.4316
Overshoot:71.1921
Undershoot:0
Peak:1.4316
PeakTime:0.9670
figure
step(Y1U2)
Info=stepinfo(Y1U2)
title('StepResponseforY1/U2')
Info=
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
5/11
3/2/2015
AE_157_Homework_2_Problem_3_c
RiseTime:0
SettlingTime:10.3152
SettlingMin:0.1991
SettlingMax:0.3063
Overshoot:Inf
Undershoot:Inf
Peak:0.3063
PeakTime:0.5372
figure
step(Y2U2)
Info=stepinfo(Y2U2)
title('StepResponseforY2/U2')
Info=
RiseTime:0.3975
SettlingTime:9.7883
SettlingMin:0.4259
SettlingMax:1.2083
Overshoot:64.7728
Undershoot:0
Peak:1.2083
PeakTime:1.0744
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
6/11
3/2/2015
AE_157_Homework_2_Problem_3_c
Partf:impulseresponse,evaluatetheimpulseresponseofthesystem.
figure
impulse(Y1U1)
title('ImpulseResponseforY1/U1')
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
7/11
3/2/2015
AE_157_Homework_2_Problem_3_c
figure
impulse(Y2U1)
title('ImpulseResponseforY2/U1')
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
8/11
3/2/2015
AE_157_Homework_2_Problem_3_c
figure
impulse(Y1U2)
title('ImpulseResponseforY1/U2')
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
9/11
3/2/2015
AE_157_Homework_2_Problem_3_c
figure
impulse(Y2U2)
title('ImpulseResponseforY2/U2')
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
10/11
3/2/2015
AE_157_Homework_2_Problem_3_c
END
PublishedwithMATLABR2012b
file:///C:/Users/komal/Documents/MATLAB/html/AE_157_Homework_2_Problem_3_c.html
11/11
