bpf.m

我做了很久的

a=[1 150 5000];b=[1 0] % H(s)=s/(s^2+150s+5000)
ljd(a,b);
q=[0]
p=[-50 -100]
f1=0
f2=100
k=0.1
splxy(f1,f2,k,p,q)