资源描述：
《dsp实验matlab程序实例》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、2．离散时间信号（序列）的产生利用MATLAB产生和绘制下列有限长序列：单位脉冲序列，单位阶跃序列，矩形序列x=-10:10;y=[];fori=1:21ifx(i)==0y(i)=1;elsey(i)=0;endendsubplot(3,1,1);%在子图画单位冲激序列stem(x,y,'*');title('冲激');xlabel('n');holdon;fori=1:21ifx(i)<0y(i)=0;elsey(i)=1;endendsubplot(3,1,2);%在子图画单位阶跃序列stem(x,y,'*');title('阶跃');xlabel('n');hold
2、on;fori=1:21ifx(i)<-4y(i)=0;elseifx(i)>4y(i)=0;elsey(i)=1;endendendsubplot(3,1,3);%在子图画矩形脉冲序列stem(x,y,'*');title('矩形');xlabel('n');holdon在一幅图上绘出曲线，，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,
3、1);%在子图中画出一系列正弦余弦曲线plot(t,x1,'-red',t,x2,'-gr',t,x3,'-bl');title('正弦余弦曲线');xlabel('t');holdon;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');holdon;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
4、);subplot(3,1,3);%在子图中画出谐波信号plot(t,x5)title('谐波信号');xlabel('t');holdon;3．序列的运算生成下列序列：(1)利用MATLAB编程完成上述两序列的卷积，并绘制运算后序列的波形。x=[12345];h=[1212];n1=[01234];n2=[2345];y=conv(x,h);%计算序列x与h的卷积和序列yn0=n1(1)+n2(1);%计算卷积和序列y的起点位置n3=length(x)+length(h)-2+n0;%计算卷积和序列y的终止位置n=n0:1:n3;%确定卷积和y非零样值的时间向量subpl
5、ot(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]=
6、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;fori=1:n4%求卷积和yifi>length(x)x(i)=0;endforj=1:iifj>length(h)h(j)=0;endm=m+x(i+1-j)*h(j
7、);endy(i)=m;m=0;endn1=n1(1):1:n1(1)+n4-1;n1:%x(n)对应的新的时间向量n2=n2(1):1:n2(1)+n4-1;%h(n)对应的新的时间向量n=n3:1:n3+n4-1;%y(n)对应的时间向量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
