数字信号处理(00002).docx

数字信号处理实验报告 实验一 1.分析并绘出常用函数（a）锯齿波；（b）三角波; （c）方波；（d）抽样函数 的时域特性波形 程序： clear alia O t=0:0.0001:0.1; x1=sawtooth(2*pi*50*t); subplot(221) plot(t,x1) title( '锯齿波 110900623') xlabel( 't') x2=sawtooth(2*pi*50*t,0.5); subplot(222) plot(t,x2) xlabel( 't') title( '三角波 110900623') x3=square(2*pi*50*t); subplot(223) plot(t,x3) xlabel( 't') title( '方波 110900623') axis([0,0.1,-1.2,1.2]) t=-4:0.1:4; x4=s in c(t); subplot(224) plot(t,x4) xlabel( 't') axis([-5,5,-0.5,0.5]) title( '抽样信号 110900623') % 产生五个周期锯齿波 %产生五个周期三角波 % 产生十个周期方波 %产生抽样信号 0.1 t 0.05 0 t 2■分析并绘出常用窗函数时域特性波形 clear all y1=boxcar(80); % 矩形窗 plot(y1,'k') axis([-1,81,-0.2,1.2]) gtxst(矩形窗' hold on y2=triang(80); % 三角窗 plot(y2,'m.') hold on y3=hanning(80); % 汉宁窗 plot(y3,'y*') gtxst(汉宁窗' hold on y4=hamming(80); % 海明窗 plot(y4,'r-') gtxst(海明窗' hold on y5=blackman(80); % 布莱克曼窗 plot(y5,'g:') gtxst(布莱克曼窗' hold on y6=kaiser(80,7.865);% 凯塞-贝尔窗 plot(y6,'b-.') gtxst(凯塞-贝尔窗' title('常用窗函数 110900623') 常用窗函数110900623 矩形窗 0.8 0.6 0.4 布莱克曼窗 凯塞-贝尔窗 0.2 海明窗??三角窗 海明窗 汉■宁窗 -0.2」0『r -0.2」 0 10 20 30 40 50 60 70 80 实验二内容： 1?计算序列 x(n)=[1,2,3,4,5],与序列 h(n)=[2,-2,3,5] 的线性卷积和6点、8点和10点圆周卷积. xn=[123,4,5]; hn=[2,-2,3,5]; yln=con v(x n,hn); yen 1=circ on v2(x n,hn ,6); yen 2=circ on v2(x n,hn ,8); yen 3=circ on v2(x n,hn ,10); ny0=[0:1:le ngth(yl n)-1]; ny 1= [0:1:le ngth(yc n1)-1]; ny2=[0:1:le ngth(yc n2)-1]; ny3=[0:1:le ngth(yc n3)-1]; subplot(2,2,1); stem( ny 0,yl n) ylabel('线性卷积') subplot(2,2,2); stem( ny1,yc n1); ylabel('圆周卷积6') subplot(2,2,3); stem( ny 2,yc n2) ylabel('圆周卷积8') subplot(2,2,4)； stem( ny 3,yc n3) ylabel('圆周卷积10') axis([0,10,0,40]); fun cti on yc=circo nv2(x1,x2,N) % 子程序 1 if len gth(x1)>N error( 'N must not be less tha n len gth of x1' ); end if len gth(x2)>N error( 'N must not be less tha n len gth of x2' ); end x1=[x1,zeros(1,N-le ngth(x1))]; x2=[x2,zeros(1,N-le ngth(x2))]; n=[0:1:N-1]; x2=x2(mod(-n,N)+1); H=zeros(N,N); for n=1:1:N H( n,:)=cirshiftd(x2, n-1,N); end yc=x1*H'; function y=cirshiftd(x,m,N) % 子程序 2 if len gth(x)>N error(