mt3

clc
clear 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 on
end
 
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 on
end
 
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