splxy.m

我做了很久的

function splxy(f1,f2,k,p,q)
p=p'
q=q'
f=f1:k:f2
w=f*(2*pi)
y=i*w
n=length(p)
m=length(q)
if n==0
    yq=ones(m,1)*y
    vq=yq-q*ones(1,length(w))
    bj=abs(vq)
    ai=1
else if m==0
    yp=ones(n,1)*y
    vp=yp-p*ones(1,length(w))
    ai=abs(vp)
    bj=1
else
    yp=ones(n,1)*y
    yq=ones(m,1)*y
    vp=yp-p*ones(1,length(w))
    vq=yq-q*ones(1,length(w))
    ai=abs(vp)
    bj=abs(vq)
end
Hw=prod(bj,1)./prod(ai,1)
plot(f,Hw)
title('连续系统幅频响应曲线')
xlabel('频率w(单位:HZ)')
ylabel('F(jw)')
end