《数字信号处理实验讲义.docVIP

王实验一 离散信号的matlab实现 实验目的 熟悉matlab软件，学会matlab语言的编写 使用matlab软件产生一些常见的离散信号 掌握用matlab软件作信号的相关分析 实验环境 计算机操作系统、matlab软件 实验内容 1、用matlab程序产生下列离散信号或连续信号，并画出其波形。 a单位抽样序列的产生 参考程序：N=100; x=zeros(1,N); 产生一个1行N列值全为0的矩阵，如看成数组x（1）-x（100）都为0 x(1)=1; n=0:N-1; stem(n,x); 产生序列 参考程序：N=100; x=zeros(1,N); k=20; x(k+1)=1; xn=0:N-1; stem(xn,x); b单位阶跃序列的产生 参考程序：N=32; x=ones(1,N);产生一个1行N列值全为1的矩阵 n=0:N-1; stem(n,x); 产生序列 参考程序：N=32; k=20; x1=zeros(1,k); x2=ones(1,N-k); x=[x1,x2]; xn=0:N-1; stem(xn,x); c模拟信号，以t=0.01n（n=0：N-1）进行采样后的离散信号。 参考程序： N=128；n=[0:N-1]; t=0.01*n; x=2*sin(4*pi*t)+5*cos(8*pi*t); figure(1); subplot(211); stem(t,x); subplot(212); stem(n,x); d产生一个sinc(t)=sint/t抽样函数 参考程序：n=200; step=4*pi/n; t=-2*pi:step:2*pi; y=sinc(t); plot(t,y,t,zeros(size(t)));%同时画出y(t）和横轴 grid on; plot(t,y,t,zeros(size(t)),zeros(size(y)),y);%同时画出y(t）和横轴、纵轴 e方波信号square(t) square(t,duty) 产生周期是2pi，幅度为正负1的方波，duty占空比,高电平跟整个周期的比值 参考程序：t=0:0.01:2*pi; y=square(t,3); plot(t,y); 试产生一个周期为1，高低电平分别为半个周期的方波信号 2、相关分析去除噪声 x(n)=sin(2*pi*n)+u(n) 噪声为高斯分布白噪声，使用相关分析去除噪声，噪声1功率为1，噪声2功率为0.1 %rxy=xcorr(x,y); %rx=xcorr(x,Mlag,flag) Mlag表示rx的单边长度，总长度为2Mlag+1, flag---biased rx(m)/N --unbiased rx(m)/(N-abs(m)) 参考程序：N=500; p1=1; p2=0.1; f=1/8; Mlag=60; u=randn(1,N); u2=u*sqrt(p2); n=[0:N-1]; s=sin(2*pi*f*n); x1=u(1:N)+s; rx1=xcorr(x1,Mlag,biased); subplot(211); plot(-Mlag:Mlag,rx1); x2=u2(1:N)+s; rx2=xcorr(x2,Mlag,biased); subplot(212); plot(-Mlag:Mlag,rx2); 3.书后习题1.15 实验二、离散信号的傅里叶变换 实验目的 进一步熟悉matlab软件的使用，熟悉matlab的编程语言 用matlab语言编写程序进行离散信号的傅里叶分析 实验原理 设离散序列，长度为N，其DTFT定义为： 在实际计算中无法取到无限长序列，通常通过无限长序列加窗作有限长序列的DTFT。 实验内容 试求序列，其中n=0，1，……，N-1，N=12的DTFT。并画出其幅频特性曲线和相频特性曲线。 若N=24,36,120,其幅频特性曲线和相频特性曲线如何变化，为什么？ 参考程序：N=12; n=[0:1:N-1]; xn=cos(n*pi/6); w=0:0.01:2*pi; M=length(w); xejw=zeros(1,M); for i=0:1:N-1 xejw=xejw+cos(i*pi/6)*exp(-j*w*i); end; xr=abs(xejw); xphase=angle(xejw); subplot(211); plot(