现代信号处理技术(上机实验报告)
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一系统响应与系统稳定性一实验程序：
1.调用filter函数接差分方程：
>> A=[1,-0.9];B=[0.05,0.05];
>> x1n=[1 1 1 1 1 1 1 1 zeros(1,50)];
>> x2n=ones(1,128);
>> hn=impz(B,A,58);
>> subplot(1,1,1);y='h(n)';tstem(hn,y);
>> title(' (a)系统单位脉冲响应h（n）')
>> y1n=filter(B,A,x1n);
>> title('(b)系统对R8(n)的响应y1n');
>> y2n=filter(B,A,x2n);
>> subplot(1,1,1);y='y2(n)';tstem(y2n,y);
>> title('(c) 系统对u(n)的响应y2(n)');
2.调用conv函数计算卷积：
>> x1n=[1 1 1 1 1 1 ];
>> h1n=[ones(1,10) zeros(1,10)];
>> h2n=[1 2.5 2.5 1 zeros(1,10)];
>> y21n=conv(h1n,x1n);
>> y22n=conv(h2n,x1n);
>> subplot(1,1,1);y='h1(n)';tstem(h1n,y);
>> title('(1) 系统单位脉冲响应h1(n)')
>> subplot(1,1,1);y='y21n';tsttem(y21n,y);
>> subplot(1,1,1);y='y21n';tstem(y21n,y);
>> title('(2) h1n与R8(n)的卷积y21(n)')
>> subplot(1,1,1);y='h2(n)';tstem(h2n,y);
>> title('(3) 系统单位脉冲响应h2(n)')
>> subplot(1,1,1);y='y22(n)';tstem(y22n,y);
>> title('(4) h2(n)与R8(n)的卷积y22(n)')
3.谐振器分析：
>> un=ones(1,256);
>> n=0:255;
>> xsin=sin(0.014*n)+sin(0.4*n);
>> A=[1,-1.8237,0.9801];B=[1/100.49,0,-1/100.49]; >> y31n=filter(B,A,un);
>> y32n=filter(B,A,xsin);
>> subplot(1,1,1);y='y31(n)';tstem(y31n,y);
>> title('(1) 谐振器对u(n)的响应y31(n)')
>> subplot(1,1,1);y='y32(n)';tstem(y32n,y);
>> title('(2) 谐振器对正弦信号的响应y32(n)') 二试验程序运行结果:
实验二用FFT对信号作频谱分析一实验程序：
1.调用函数fft计算序列x(n)的DFT
>>x1n=[ones(1,4)];
>> M=8;xa=1:(M/2);xb=(M/2):(-1):1;x2n=[xa,xb];
>> M=8;xa=1:(M/2);xb=(M/2):-1:1;x2n=[xa,xb];
>> x3n=[xb,xa];
>> X1k8=fft(x1n,8);
>> X1K16=fft(x1n,16);
>> X1k16=fft(x1n,16);
>> X2k8=fft(x2n,8);
>> X2k16=fft(x2n,16);
>> X3k8=fft(x3n,8);
>> X3k16=fft(x3n,16);
%以下为绘图程序
>> subplot(1,1,1);mstem(X1k8)
>> title('(a)8点DFT[x_1(n)]');
>> axis([0,2,0,1.2*max(abs(X1k8))])
>> subplot(1,1,1);mstem(X1k16);
>> title('(b)16点DFT[x_1(n)]');
>> axis([0,2,0,1.2*max(abs(X1k16))])
>> subplot(1,1,1);mstem(X2k8)
>> title('(c) 8点DFT[X_2(n)]');
>> axis([0,2,0,1.2*max(abs(X2k8))])
>> subplot(1,1,1);mstem(X2k16)
>> title('(d) 16点DFT[x_2(n)]');
>> axis([0,2,0,1.2*max(abs(X2k16))])
>> subplot(1,1,1);mstem(X3k8);
>> title('(e) 8点DFT[x_3(n)]');
>> axis([0,2,0,1.2*max(abs(X3k8))])
>> subplot(1,1,1);mstem(X3k16);
>> title('(f) 16点DFT[x_3(n)]');
>> axis([0,2,0,1.2*max(abs(X3k16))])
2.周期序列谱分析：
>> N=8;n=0:N-1;
>> x4n=cos(pi*n/4);
>> x5n=cos(pi*n/4)+cos(pi*n/8);
>> X4k8=fft(x4n);
>> x5k8=fft(x5n);
>> N=16;n=0:N-1;