南邮通达数字信号实验报告.docVIP

实验一：熟悉MATLAB环境 实验要求 )熟悉MATLAB的主要操作命令。 学会简单的矩阵输入和数据读写。 掌握简单的绘图命令。 用MATLAB编程并学会创建函数。 观察离散系统的频率响应。 实验内容 用MATLAB实现以下序列。 x(n)=e(0.2+3j)n 0≤n≤15 实验程序： n=0:15;x=exp((0.2+3*j)*n); stem(n,x,fill); xlabel (时间序列n );ylabel(x(n)=exp((0.2+3*j)*n)); 实验结果： x(n)=3cos(0.125πn+0.2π)+2sin(0.25πn+0.1π) 0≤n≤15 实验程序： n=0:1:15; x=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi); stem(n,x); xlabel(n); ylabel(幅度); 实验结果： 将c)中的x(n)扩展为以16为周期的函数x16(n)=x(n+16),绘出四个周期 实验程序： n=0:1:63; x=3*cos(0.125*pi*rem(n,16)+0.2*pi)+2*sin(0.25*pi*rem(n,16)+0.1*pi); stem(n,x,fill); xlabel (时间序列n ); ylabel(x16(n)); 实验结果： n=0:1:63; x=3*cos(0.125*pi*rem(n,16)+0.2*pi)+2*sin(0.25*pi*rem(n,16)+0.1*pi); stem(n,x,fill); xlabel (时间序列n );ylabel(x16(n)); 绘出时间函数的图形，对x轴、y轴图形上方均须加上适当的标注。 x(t)=sin(2πt) 0≤t≤10s 实验程序： clc; t1=0:0.001:10;t2=0:0.01:4; xa=sin(2*pi*t1);xb=cos(100*pi*t2).*sin(pi*t2); subplot(2,1,1); plot(t1,xa);xlabel (t);ylabel(x(t));title(x(t)=sin(2*pi*t)); subplot(2,1,2); plot(t2,xb);xlabel (t);ylabel(x(t));title(x(t)=cos(100*pi*t2).*sin(pi*t2)); 实验结果： 给一定因果系统求出并绘制H(z)的幅频响应与相频响应。 实验程序： a=[1 -0.67 0.9]; b=[1 sqrt(2) 1]; [h w]=freqz(b,a); fp=20*log(abs(h)); subplot(2,1,1); plot(w,fp);xlabel(时间序列t);ylabel(幅频特性); xp=angle(h); subplot(2,1,2); plot(w,xp);xlabel(时间序列t);ylabel(相频特性); 实验结果： 求以下差分方程所描述系统的单位脉冲响应 实验程序： N=20; a=[1 -2]; b=[1 0.1 -0.06]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(a,b,x); stem(k,y); title(单位脉冲响应); xlabel(n);ylabel(幅度); plot() 实验结果： 实验二： 实验目的 在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉FFT子程序。 熟悉应用FFT对典型信号进行频谱分析的方法。 了解应用FFT进行信号频谱分析过程中可能出现的问题以便在实际中正确应用FFT。 熟悉应用FFT实现两个序列的线性卷积的方法。 对DCT变换用作数据压缩有初步的认识。 实验内容 观察高斯序列的时域和幅频特性。 实验程序： clc; n=0:1:15; % p=8,q=2 p=8;q=2; x1=exp(-(n-p).^2./q); fp1=fft(x1);fp1=abs(fp1); subplot(5,2,1);plot(x1);xlabel(n);ylabel(时域特性);title(p=8,q=2); subplot(5,2,2);plot(fp1);xlabel(n);ylabel(幅频特性);title(p=8,q=2); % p=8,q=4 p=8;q=4; x2=exp(-(n-p).^2./q); fp2=fft(x2);fp2=abs(fp2); subplot(5,2,3);plot(x2);xlabel(n);ylabel(时域特性);title(p=8,q=4); subplot(5,2,4);plot(fp2);xlabel(n);ylabel(幅频特性);t