The step function only delivers unit amplitudes.
One option would be to cascade an ‘amplifier’ with your system.
Example —
a = [-0.5572,-0.7814;0.7814,0];
b = [1,-1;0,2];
c = [1.9691,6.4493];
sys = ss(a,b,c,0);
[y1, t1] = step(sys);
a = [-0.5572,-0.7814;0.7814,0];
b = [1,-1;0,2] * 91;
c = [1.9691,6.4493];
sys = ss(a,b,c,0);
[y2, t2] = step(sys);
figure(1)
subplot(2,1,1)
plot(t1, squeeze(y1))
title('Original')
grid
subplot(2,1,2)
plot(t2, squeeze(y2))
title('Amplified')
grid
