tt2

最新推荐文章于 2018-08-20 16:22:00 发布

转载最新推荐文章于 2018-08-20 16:22:00 发布 · 606 阅读

clc
clear all
Lp=4.3e-3; 
Ra=13.51;   
f=logspace(1,5,200);  
w=2*pi*f;
num=[1];
den=[Lp/Ra 1];
[mag,pha]=bode(num,den,w);
 
figure(2);
num=[1];
for i=0:8
    Lp=Lp*(0.8+i*0.05);
    den=[Lp/Ra 1];
    [mag,pha]=bode(num,den,w);
    loglog(f,mag,'LineWidth',1);
    hold on
    Lp=4.3e-3;
end