DSP实验Matlab程序实例.docVIP

2．离散时间信号（序列）的产生 利用MATLAB产生和绘制下列有限长序列： 单位脉冲序列?(n)，单位阶跃序列u(n)，矩形序列R8(n) 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)lt;0 y(i)=0; else y(i)=1; end end 在一幅图上绘出曲线 x?? 1(t)?Asin(5t?3)， x2(t)?Acos(t?? 3 ) subplot(3,1,2); %在子图画单位阶跃序列 stem(x,y,’*’); title(‘阶跃’); xlabel(‘n’); hold on; for i=1:21 if x(i)lt;-4 y(i)=0; else if x(i)gt;4 y(i)=0; else y(i)=1; end end end subplot(3,1,3); %在子图画矩形脉冲序列 stem(x,y,’*’); title(‘矩形’); xlabel(‘n’); hold on x3(t)?2sin(3t)cos(2t?) 3? x4(t)?Ae??tsin(2?ft)，A=2, α=0.5, f=2Hz A2=0.5, A3=0.2，x5(t)?A1sin(2?ft)?A2sin(2?2ft)?A3sin(2?3ft)，A1=1, 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．序列的运算 生成下列序列： x(n)??(n)?2?(n?1)?3?(n?2)?4?(n?3)?5?(n?4) h(n)??(n)?2?(n?1)??(n?2)?2?(n?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)对应的非零样