Mt1

来源:互联网 发布:义乌淘宝美工培训学校 编辑:程序博客网 时间:2024/06/05 20:25
clcclear alla=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