Let omega0 10 Plot the gain function Gs for c 05 1 2 3 4 5

Let omega_0 = 10. Plot the gain function G(s) for c = 0.5, 1, 2, 3, 4, 5 on the interval s elementof [0, 2]. G(s) = 1/omega^2_0 squareroot (1 - s^2) + D^2s^2 where s = omega_omega_0 and D = 2c/omega_0. Plot the different curves on the same plot with different color. Label your axes.

Solution

syms s;
w0 = 10;

c = 0.5;
D = (2*c)/w0;
G05 = 1/(100*sqrt((1-s^2)^2 + D^2*s^2));

c = 1;
D = (2*c)/w0;
G1 = 1/(100*sqrt((1-s^2)^2 + D^2*s^2));

c = 2;
D = (2*c)/w0;
G2 = 1/(100*sqrt((1-s^2)^2 + D^2*s^2));

c = 3;
D = (2*c)/w0;
G3 = 1/(100*sqrt((1-s^2)^2 + D^2*s^2));

c = 4;
D = (2*c)/w0;
G4 = 1/(100*sqrt((1-s^2)^2 + D^2*s^2));

c = 5;
D = (2*c)/w0;
G5 = 1/(100*sqrt((1-s^2)^2 + D^2*s^2));

plotgraph05 = ezplot(G05,[0,1]);
hold on

plotgraph1 = ezplot(G1,[0,1]);
hold on

plotgraph2 = ezplot(G2,[0,1]);
hold on

plotgraph3 = ezplot(G3,[0,1]);
hold on

plotgraph4 = ezplot(G4,[0,1]);
hold on

plotgraph5 = ezplot(G5,[0,1]);
hold on

set(plotgraph05,\'color\',[1 0 0])
set(plotgraph1,\'color\',[0 0 1])
set(plotgraph2,\'color\',[0 1 0])
set(plotgraph3,\'color\',[1 1 0])
set(plotgraph4,\'color\',[0 1 1])
set(plotgraph5,\'color\',[1 1 1])