数字信号处理实验全部程序MATLAB.docVIP

单位采样序列 function [x,n]=impseq(n1,n2,n0) n=[n1:n2]; x=[(n-n0)==0]; [x,n]=impseq(-2,8,2); stem(n,x); title(电信1201) n0=-2; n=[-10:10]; nc=length(n); x=zeros(1,nc); for i=1:nc if n(i)==n0 x(i)=1 end end stem(n,x) ; title(电信1201 采样序列第二种方法) 单位阶跃序列 function [x,n]=stepseq(n1,n2,n0) n=[n1:n2]; x=[(n-n0)=0]; [x,n]=stepseq(-2,8,2); stem(n,x); title(电信1201) x=0.9.^n; stem(n,x); title(电信1201) 复数指数序列 n=[-10:10]; alpha=-0.1+0.3*j; x=exp(alpha*n); real_x=real(x); image_x=imag(x); mag_x=abs(x); phase_x=angle(x); subplot(2,2,1); stem(n,real_x) ;title(电信1201) subplot(2,2,2); stem(n,image_x);title(电信1201) subplot(2,2,3); stem(n,mag_x);title(电信1201) subplot(2,2,4); stem(n,phase_x);title(电信1201) 正余弦 n=[0:10]; x=3*cos(0.1*pi*n+pi/3); stem(n,x); title(电信1201) 例： 求出下列波形 x1(n)=2x(n-5)-3x(n+4) function [y,n]=sigadd(x1,n1,x2,n2) n=min(min(n1),min(n2)) : max(max(n1),max(n2)); y1=zeros(1,length(n)); y2=y1; y1(find((n=min(n1))(n=max(n1))==1))=x1; y2(find((n=min(n2))(n=max(n2))==1))=x2; y=y1+y2; function [y,n]=sigshift(x,m,n0) n=m+n0; y=x; n=[-2:10]; x=[1:7,6:-1:1]; [x11,n11]=sigshift(x,n,5); [x12,n12]=sigshift(x,n,-4); [x1,n1]=sigadd(2*x11,n11,-3*x12,n12); stem(n1,x1); title(电信1201) 已知某一系统方程为：y[n]-y[n-1]+0.9y[n-2]=x[n]计算并画出脉冲响应h(n),n=(-20,100) n=(-20:100); num=(1); den=[1 -1 0.9]; x=impseq(-20,100,0); h=filter(num,den,x); stem(n,h) xlabel(时间序号N); ylabel(脉冲响应h); title(脉冲响应 电信1201); 实现两个序列a,b的卷积 x=[3,11,7,0,-1,4,2]; h=[2,3,0,-5,2,1]; c=conv(x,h); stem(c); title(电信1201) 自己给的数 x=[0.5,1,1.5,0,0]; h=[1,1,1]; c=conv(x,h); stem(c); title(电信1201) 试求卷积C(t)=f1(t)*f2(t),并绘制出f1、f2、及卷积以后的波形。 function [y,ny]= conv_m(x,nx,h,nh,p) % 信号处理的改进卷积程序 nyb=nx(1)+nh(1); nyc=nx(length(x))+nh(length(h)); ny=(nyb :p: nyc); y=conv(x , h); p=0.1; t1= [0:p:1]; f1=t1.*(t10); t2= [-1:0.1:2]; f2=t2.*exp(t2).*(t2=0)+exp(t2).*(t20); [y