数字信号处理实验报告和自主设计实验.docxVIP

北京科技大学《信号系统与信号处理综合实验》实 验 报 告专业班级：学生姓名：学 号：指导教师：实验成绩： 年 月 日计算机与通信工程学院目录一、参考实验21、实验122、实验273、实验名称（实验三： 用FFT对信号作频谱分析）124、结论及思考18二、自主设计实验18三、总结与收获18一、参考实验1、实验1实验目的（1）掌握求系统响应的方法。（2）掌握时域离散系统的时域特性。（3）分析、观察及检验系统的稳定性。实验内容实验报告要求1）实验1-（1）程序清单图1-1 2）实验1-（2）程序清单3）实验1-（3）程序清单4）实验1-（4）程序清单图1-4 输入u(n)和x(n)的输出波形2、实验2实验目的实验内容实验报告要求1）时域采样理论的验证程序清单A=444.128;a=50*sqrt(2)*pi;w0=50*sqrt(2)*pi; Tp=50/1000;F1=1000;F2=300;F3=200;T1=1/F1;T2=1/F2;T3=1/F3;n1=0:Tp*F1-1;n2=0:Tp*F2-1;n3=0:Tp*F3-1; x1=A*exp(-a*n1*T1).*sin(w0*n1*T1);x2=A*exp(-a*n2*T2).*sin(w0*n2*T2);x3=A*exp(-a*n3*T3).*sin(w0*n3*T3);f1=fft(x1,length(n1));f2=fft(x2,length(n2));f3=fft(x3,length(n3));k1=0:length(f1)-1; fk1=k1/Tp; k2=0:length(f2)-1; fk2=k2/Tp; k3=0:length(f3)-1; fk3=k3/Tp; subplot(3,2,1) stem(n1,x1,.) title((a)Fs=1000Hz); xlabel(n);ylabel(x1(n)); subplot(3,2,3) stem(n2,x2,.) title((b)Fs=300Hz); xlabel(n);ylabel(x2(n)); subplot(3,2,5) stem(n3,x3,.) title((c)Fs=200Hz); xlabel(n);ylabel(x3(n)); subplot(3,2,2) plot(fk1,abs(f1)) title((a) FT[xa(nT)],Fs=1000Hz); xlabel(f(Hz));ylabel(幅度) subplot(3,2,4) plot(fk2,abs(f2)) title((b) FT[xa(nT)],Fs=300Hz); xlabel(f(Hz));ylabel(幅度) subplot(3,2,6) plot(fk3,abs(f3)) title((c) FT[xa(nT)],Fs=200Hz); xlabel(f(Hz));ylabel(幅度)图2-2 时域采样理论的验证结果（2）频域采样理论的验证M=27;N=32;n=0:M; xa=0:floor(M/2); xb= ceil(M/2)-1:-1:0; xn=[xa,xb]; Xk=fft(xn,1024); X32k=fft(xn,32) ; x32n=ifft(X32k); X16k=X32k(1:2:N); x16n=ifft(X16k,N/2); subplot(3,2,2);stem(n,xn,.);box on title((b) 三角波序列x(n));xlabel(n);ylabel(x(n));axis([0,32,0,20]) k=0:1023;wk=2*k/1024; % subplot(3,2,1);plot(wk,abs(Xk));title((a)FT[x(n)]); xlabel(\omega/\pi);ylabel(|X(e^j^\omega)|);axis([0,1,0,200]) k=0:N/2-1; subplot(3,2,3);stem(k,abs(X16k),.);box on title((c) 16点频域采样);xlabel(k);ylabel(|X_1_6(k)|);axis([0,8,0,200]) n1=0:N/2-1; subplot(3,2,4);stem(n1,x16n,.);box on title((d) 16点IDFT[X_1_6(k)]);xlabel(n);ylabel(x_1_6(n));axis([0,32,0,20]) k=0:N-1; subplot(3,2,5);stem(k,abs(X32k),.);box on title((e) 32点频域采样);xlabel(k);ylabe