mt3

clcclear all KT=4.3e-3;   KTe=KT;  Jp=0.015e-3; Jpe=Jp;      bo=0.0292;Kso=9.1787;Kiso=465.0942; f=logspace(0,6,1000);  w=2*pi*f;numtp=[Jp bo Kso Kiso];dentp=[Jpe bo Kso Kiso];figure(1)   for i=0:4    KTe=KT*(0.8+i*0.1);    num=[KTe/KT*Jp bo Kso Kiso]    den=[Jpe bo Kso Kiso];    [mag,pha]=bode(num,den,w);    loglog(f,mag,'LineWidth',2);    hold onend xlabel('f/Hz');ylabel('Magnitude');title('Dynamic Stiffness Frequency Response of KT changing within 20%from nominal');axis([10^0 10^6 0.6 1.4]);grid on  KTe=KT;      Jpe=Jp;      figure(2)    for i=0:4    Jpe=Jp*(0.8+i*0.1);    num=[KTe/KT*Jp bo Kso Kiso]    den=[Jpe bo Kso Kiso];    [mag,pha]=bode(num,den,w);    loglog(f,mag,'LineWidth',2);    hold onend xlabel('f/Hz');ylabel('Magnitude');title('Dynamic Stiffness Frequency Response of Jp changing within 20%from nominal');axis([10^0 10^6 0.6 1.4]);grid on