%附录 1 不同母线结构的反射波幅频特性
%电压反射系数幅频特性
num1=[-1 0]; %n=1 时，阶跃电压的输入幅频特性
den1=[1 1e6];
[mag1,phase1,w1]=bode(num1,den1);
den2=[1 5e5];
[mag2,phase2,w2]=bode(num1,den2);
den3=[1 333333];
[mag3,phase3,w3]=bode(num1,den3);
f1=w1/(2*pi);f2=w2/(2*pi);f3=w3/(2*pi);
subplot(1,2,1);
plot(f1,mag1,f2,mag2,f3,mag3);axis([0 2e6 0 1]);
title('n=1');
xlabel('f(HZ)');
ylabel('Gain(DB)');
text(1e6,0.7, '\leftarrowc=5000pF');
text(1e6,0.8,'\leftarrowc=10000pF');
text(0,1e6,0.99 ,'\leftarrowc=15000pF');
subplot(1,2,2);
plot(f1,phase1,f2,phase2,f3,phase3); axis([0 2e6 0 1]);
title('n=1');
xlabel('f(HZ)');
ylabel('phase(deg)');
text(4e6,-160,'\leftarrowc=5000pF');
text(3.5e6,-170,'\leftarrowc=10000pF');
text(2e6,-175,'c =15000pF\rightarrow');
%n=3 时，阶跃电压的输入幅频特性
a1=[-1 -1e6];a2=[-1 -0.5e6];a3=[-1 -0.33e6];
b1=[1 2e6];b2=[1 1e6];b3=[1 0.667e6];[mag1,phase1,w1]=bode(a1,b1);
[mag2,phase2,w2]=bode(a2,b2);
[mag3,phase3,w3]=bode(a3,b3);
f1=w1/(2*pi);f2=w2/(2*pi);f3=w3/(2*pi);
subplot(1,2,1);
plot(f1,mag1,f2,mag2,f3,mag3); axis([0 2e6 0.4 1]);
title ('n=3');xlabel('f(HZ)');ylabel('Gain(db)');
text(0.5e6,0.87,'\leftarrowc=5000pF');
text(0.22e6,0.9, '\leftarrowc=10000pF');
text(0.3e6,0.95,'\leftarrowc=15000pF');
subplot(1,2,2);plot(f1,phase1,f2,phase2,f3,phase3); axis([0 2e6 0.4 1]);
title('n=3');xlabel('f(HZ)');ylabel('Gain(db)');
text(1e6,-171,'\leftarrowc=5000pF');
text(0.75e6,-174,'\leftarrowc=10000pF');
text(0.2e6,-177,'c =l5000pF\rightarrow');
