原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
[45];b=[555];[zpk]=tf2zp(a,b)z=-1.2500p=-0.500=buttord(wp,ws,rp,rs,s);[B,A]=butter(N,wc,s)Rp=20*log10((1+dev1)/(1-dev1));Rs=40;[N,wc]=e列中含有两种频率成分，
[45];b=[555];[zpk]=tf2zp(a,b)z=-1.2500p=-0.500
=buttord(wp,ws,rp,rs,s);[B,A]=butter(N,wc,s)
Rp=20*log10((1+dev1)/(1-dev1));Rs=40;[N,wc]=e
列中含有两种频率成分，f1=2hz，f2=2.05hz，采样频率为fs=10Hz，即x(n)=sin
信号处理与系统分析课程设计内容
一、※三个系统的系统函数分别为h1(
h2(s)=(4s+5)/[s(5s+1)];h3(s)=(4s+5)/[5s2+5s+5)]用simulink来仿真三个系统
的阶跃响应，并分析系统的稳定性。
1、a=[5];
b=[51];
[zpk]=tf2zp(a,b)
z=Emptymatrix:0-by-1
p=-0.2000
k=1
极点p=-0.2，在s域的左半开平面，所以系统稳定！
2、a=[45];
b=[510];
(wp,ws,rp,rs,s);[b,a]=butter(n,ws,s);[bz,az]8;[N,wp]=cheb1ord(wp/pi,ws/pi,rp,rs);[Bz,Az]=chebyTs=1/fs;Tp=20;N=fs*Tp;n=[0:N-1];xn=sin(2*pi*f1,
(wp,ws,rp,rs,s);[b,a]=butter(n,ws,s);[bz,az]
8;[N,wp]=cheb1ord(wp/pi,ws/pi,rp,rs);[Bz,Az]=cheby
Ts=1/fs;Tp=20;N=fs*Tp;n=[0:N-1];xn=sin(2*pi*f1
,阻带截止频率0.3*pi,通带波动1dB；在频率0.3π到π之间的阻带衰减大于10dB。t=1
[zpk]=tf2zp(a,b)
z=-1.2500
p=0
-.2000
k=0.8000
极点p=-0.2在s域的左半开平面，但p=0不在左半开平面，所以系统不稳定！
3、a=[45];
b=[555];
[zpk]=tf2zp(a,b)
z=-1.2500
p=-0.5000+0.8660i
-0.5000-0.8660i
k=0.8000
由极点p可知，-0.5000+0.8660i和0.5000-0.8660i并不都在s域的左半开平面，所以该系统不稳定！
信号：x(t)=s(t)+u(t),其中u(t)是干扰信号，
信号：x(t)=s(t)+u(t),其中u(t)是干扰信号，s(t)+u(t)的频谱不混叠，其幅度谱
pi);ylabel(相位(rad));Wp=3/8;dev1=0.02;Ws=0.5;
,w]=freqz(b,a);subplot(2,1,1)plot(w/pi,20*log10(ab
=buttord(wp,ws,rp,rs,s);[B,A]=butter(N,wc,s)
1s2s
三、频谱分析
设一序列中含有两种频率成分，f1=2hz，f2=2.05hz，采样频率为fs=10Hz，即x(n)=sin(2
πfn/f)+sin(2πfn/f)，分析其频谱。
f1=2;f2=2.05;
fs=10;Ts=1/fs;Tp=20;N=fs*Tp;
n=[0:N-1];
xn=sin(2*pi*f1*Ts*n)+sin(2*pi*f2*Ts*n);
Xk=fft(xn,N);
stem(n,abs(Xk),*)
xlabel(k);ylabel(|X|)
Ts=1/fs;Tp=20;N=fs*Tp;n=[0:N-1];xn=sin(2*pi*f11,2)plot(w,angRp=20*log10((1+dev1)/(1-dev1));Rs=40;[N,wc]=ek=0.8000
Ts=1/fs;Tp=20;N=fs*Tp;n=[0:N-1];xn=sin(2*pi*f1
1,2)plot(w,ang
Rp=20*log10((1+dev1)/(1-dev1));Rs=40;[N,wc]=e
k=0.8000极点p=-0.2在s域的左半开平面，但p=0不在左半开平面，
文档评论（0）