来做道matlab编程题,与信号与系统有关
已知系统函数H(s)=(s^2+s+2)/(3s^3+5s^2+4s-6),试用MATLAB绘制系统的零极点图,并根据零极点图判断系统的稳定性,写出完整程序,最好带点小说...
已知系统函数H(s)=(s^2+s+2)/(3s^3+5s^2+4s-6),试用MATLAB绘制系统的零极点图,并根据零极点图判断系统的稳定性,写出完整程序,最好带点小说明解释。
展开
1个回答
展开全部
这个式计算N阶的B(s)的分母的多项式;
N=11;
for p=1:N
s(p)=exp(i*pi*(0.5+(2*p-1)/2/N));
end
a=[1];
for p=1:N
c=[1,-1*s(p)];
a=conv(a,c);
end
a=abs(a);
自己改N即可;
1)
把for p=1:N改为for p=1:2*N即可
最后加一句画所有极点 pzmap([1],a)
2)令N=1
N=1;wc=10*pi;
for p=1:N
s(p)=wc*exp(i*pi*(0.5+(2*p-1)/2/N));
end
a=[1];
for p=1:N
c=[1,-1*s(p)];
a=conv(a,c);
end
a=abs(a);
b=wc.^N;
w=linspace(0,1000);
freqs(b,a,w)
4)在1)中令N=3
5)在4)中已经实现 分母的计算;没有指定wc;
加一句bode(b,a);
时间关系,不写了
N=11;
for p=1:N
s(p)=exp(i*pi*(0.5+(2*p-1)/2/N));
end
a=[1];
for p=1:N
c=[1,-1*s(p)];
a=conv(a,c);
end
a=abs(a);
自己改N即可;
1)
把for p=1:N改为for p=1:2*N即可
最后加一句画所有极点 pzmap([1],a)
2)令N=1
N=1;wc=10*pi;
for p=1:N
s(p)=wc*exp(i*pi*(0.5+(2*p-1)/2/N));
end
a=[1];
for p=1:N
c=[1,-1*s(p)];
a=conv(a,c);
end
a=abs(a);
b=wc.^N;
w=linspace(0,1000);
freqs(b,a,w)
4)在1)中令N=3
5)在4)中已经实现 分母的计算;没有指定wc;
加一句bode(b,a);
时间关系,不写了
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询