DSP实验Matlab程序实例.doc

2．离散时间信号（序列）的产生 利用MATLAB产生和绘制下列有限长序列： 单位脉冲序列，单位阶跃序列，矩形序列 x=-10:10; y=[]; for i=1:21 if x(i)==0 y(i)=1; else y(i)=0; end end subplot(3,1,1); %在子图画单位冲激序列 stem(x,y,*); title(冲激); xlabel(n); hold on; for i=1:21 if x(i)0 y(i)=0; else y(i)=1; end end subplot(3,1,2); %在子图画单位阶跃序列 stem(x,y,*); title(阶跃); xlabel(n); hold on; for i=1:21 if x(i)-4 y(i)=0; else if x(i)4 y(i)=0; else y(i)=1; end end end subplot(3,1,3); %在子图画矩形脉冲序列 stem(x,y,*); title(矩形); xlabel(n); hold on  在一幅图上绘出曲线 ， ，A=2, α=0.5, f=2Hz ，A1=1, A2=0.5, A3=0.2，f=2Hz。 t=[0:0.5:720]*pi/180; x1=5*sin(pi*t/5+pi/3); x2=5*cos(t-pi/3); x3=2*sin(3*t).*cos(2*t-pi/3); subplot(3,1,1); %在子图中画出一系列正弦余弦曲线 plot(t,x1,-red,t,x2,-gr,t,x3,-bl); title(正弦余弦曲线); xlabel(t); hold on; A=2; a=0.5; f=2; x4=A*exp(-a*t).*sin(2*pi*f*t); subplot(3,1,2);%在子图中画出正弦衰减信号 plot(t,x4) title(正弦衰减信号); xlabel(t); hold on; A1=1; A2=0.5; A3=0.2; f=2; x5=A1*sin(2*pi*f*t)+A2*sin(2*pi*2*f*t)+A3*sin(2*pi*3*f*t); subplot(3,1,3);%在子图中画出谐波信号 plot(t,x5) title(谐波信号); xlabel(t); hold on; 3．序列的运算 生成下列序列： (1) 利用MATLAB编程完成上述两序列的卷积，并绘制运算后序列的波形。 x=[1 2 3 4 5]; h=[1 2 1 2]; n1=[0 1 2 3 4]; n2=[2 3 4 5]; y=conv(x,h); %计算序列x与h的卷积和序列y n0=n1(1)+n2(1); %计算卷积和序列y的起点位置 n3=length(x)+length(h)-2+n0; %计算卷积和序列y的终止位置 n=n0:1:n3; %确定卷积和y非零样值的时间向量 subplot(3,1,1); stem(n1,x,.); %在子图绘x(n) title(x(n)); xlabel(n); ylabel(x(n)); subplot(3,1,2); stem(n2,h,.); %在子图绘h(n) title(h(n)); xlabel(n); ylabel(h(n)); subplot(3,1,3); stem(n,y,.); %在子图绘卷积和y(n) title(y(n)=x(n)*h(n)); xlabel(n); ylabel(y(n)); （2）编写一个完成卷积的函数，输入是两个序列，输出是结果。 function [y,n]=jsjuanji(x,n1,h,n2) %计算序列卷积和y(n)=x(n)*h(n) % y: 卷积和y(n)对应的非零样值向量 % n：卷积和y(n)对应的时间向量 % x: x(n)对应的非零样值向量 % n1: x(n)对应的时间向量 % h: h(n)对应的非零样值向量 % n2: h(n)对应的时间向量 n3=n1(1)+n2(1); %卷积和序列y的起始位置 n4=length(x)+length(h)-1;%卷积和序列y的长度 m=0; for i=1:n4 %求卷积和y if ilength(x) x(i)=0; end for j=1:i if jlength(h) h(j)=0; end m=m+x