Mt1

clc
clear all
a=20*pi;
b=100*pi;
c=500*pi;
Jp=0.015e-3;
f1=1/a;
f2=1/b;
f3=1/c;
K=Jp/(f1*f2+f1*f3+f2*f3);
bo=Jp^2/K*a*b*c;
Kso=K*bo/Jp*(f1+f2+f3);
Kiso=K*bo/Jp;
f=logspace(-1,4,1000);  
w=2*pi*f;
num=Jp*[Jp bo Kso Kiso];
den=[bo 0];
figure(1);
[mag,pha]=bode(num,den,w);
loglog(f,mag,'-b','LineWidth',1);
xlabel('f/Hz');
ylabel('Magnitude');
grid on