1

%*********************%
%Program for Project 3
clc
close all
Lp=4.3e-3;  %Given value
Ra=13.51;   %Designed value
f=logspace(0,6,1000);  %
w=2*pi*f;
num=[1];
den=[Lp/Ra 1];
[mag,pha]=bode(num,den,w);
loglog(f,mag,'r','LineWidth',2);
xlabel('f/Hz');
ylabel('幅值');
title('Dynamic Stiffness Frequency Response of the Current Loop');


f=logspace(-2,6,1000);  %
w=2*pi*f;
Jp=0.015e-3;
br=5.84e-3;
Ksr=0.367;
Kisr=3.72;
num2=-[Jp br Ksr Kisr];
den2=[1 0];
figure(2);
[mag,pha]=bode(num2,den2,w);
%[mag,pha]=bode(den2,num2,w);
loglog(f,mag,'r','LineWidth',2);
xlabel('f/Hz');
ylabel('幅值');
title('Dynamic Stiffness Frequency Response');