原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验1 常见离散信号的MATLAB产生和图形显示
1.单位抽样序列
x=zeros(1,50);?
x(1)=1;?
n=1:50;?
stem(n,x,.);
axis([0,50,0,1]);?
xlabel(n),ylabel(x)
2．单位阶跃序列
n=-25:25;?
x=ones(1,51);?
stem(n,x,.);?
xlabel(n);
ylabel(x(n));
axis([-25,25,0,1])
3．正弦序列
f=10;
Fs=250;
fai=0;
n=1:499;
N=500;
x=sin(2*pi*f*n/Fs+fai);
stem(n,x);
axis([0,100,-1.5,1.5]);
xlabel(n),ylabel(x)
4．复正弦序列
n=0:0.1:5;
w=2;
x=exp(j*w*n);
stem(n,x,.);?
xlabel(n);
ylabel(x(n)) ;
axis([0,5?-1,1])
5．指数序列
n=0:50;
a=2;
x=a.^n;
stem(n,x);
axis([0,50,0,10000]);
xlabel(n),ylabel(x)
实验2 离散系统的差分方程、单位脉冲响应和卷积分析
1.
A=[1,0.6,0.09];B=[1,-1];
x2n=ones(1,25);
x1n=[1,zeros(1,10)];
y1n=filter(B,A,x1n);
subplot(2,1,1);y=y1(n);
stem(y1n,g,.);
title(单位脉冲响应)
y2n=filter(B,A,x2n);
subplot(2,1,2);
y=y2(n);
stem(y2n,g,.);
title(阶跃响应)
2.
A=[1];B=[0,0.5,0.5,0.5,0.5,0.5];?
x2n=ones(1,15);?
x1n=[1,zeros(1,20)];?
y1n=filter(B,A,x1n);?
subplot(2,1,1);y=y1(n);?
stem(y1n,g,.);?
title(单位脉冲响应)?
y2n=filter(B,A,x2n);?
subplot(2,1,2);?
y=y2(n);?
stem(y2n,g,.);?
title(阶跃响应)?
实验3 离散系统的频率响应分析和零、极点分布
num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528];
den=[1 -1.8107 2.4947 -1.8801 0.9537 -0.2336];
[z,p,k]=tf2zp(num,den);
disp(零点);disp(z);
disp(极点);disp(p);
figure(1)
zplane(num,den)
figure(2)
freqz(num,den,128)
零点
-1.5870 + 1.4470i
-1.5870 - 1.4470i
0.8657 + 1.5779i
0.8657 - 1.5779i
-0.0669
极点
0.2788 + 0.8973i
0.2788 - 0.8973i
0.3811 + 0.6274i
0.3811 - 0.6274i
0.4910
实验4 离散信号的DTFT和DFT
1.
N=16;
n=1:16;
x=sin(5*pi*n/16);
X1=fft(x,16);
X11=abs(X1);
subplot(2,1,1);
stem(X11);
xlabel(频率);
ylabel(幅度);
title(16点);
X2=fft(x,32);
X22=abs(X2);
subplot(2,1,2);
stem(X22);
xlabel(频率);
ylabel(幅度);
title(32点);
2.
x= sin(5*pi*n/16);
X1=fft(x);
X11=abs(X1);
stem(X11);
xlabel(频率);
ylabel(幅度);
title( DTFT);
实验5 FFT算法的应用
1：
N=64;
n=0:2*N-1;
x=cos(2*pi*7*n/N)+1/2*cos(2*pi*19*n/N);
X=fft(x,128);
k=n;
stem(k,abs(X))
grid
xlabel(k);
ylabel(|X[k]|);
2.
k=0:63;
N=64;
X=1./(1-0.6*exp(-j*2*pi*k/N));
x
文档评论（0）