mt2

clc
clear all
Lp=4.3e-3;  %Given value
Lp1=Lp*0.8;
Lp2=Lp*1.2;
Ra=13.51;   %Designed value
f=logspace(0,6,1000);  %
w=2*pi*f;
num=[1];

for i=0:4
    Lp=Lp*(0.8+i*0.1);
    den=[Lp/Ra 1];
    [mag,pha]=bode(num,den,w);
    loglog(f,mag,'LineWidth',2);
    hold on
    Lp=4.3e-3;
end

xlabel('f/Hz');
ylabel('Magnitude');
title('Dynamic Stiffness Frequency Response of Lp changing within 20%
from nominal' );
axis([10^0 10^6 10^-4 10^0.5])
grid on