MATLAB第一次作业..docVIP

班级：2011211123 姓名：黄昭龙 学号：2011211512 实验一：数字信号的产生和基本运算 a) δ(n) imqsep.m： function [x,n] = impseq(n0,n1,n2) % Generates x(n) = delta(n-n0); n1 = n,n0 =n 2; % [x,n] = impseq (n0,n1,n2) if ((n0n1)||(n0n2)||(n1n2)) error(Arguments must satisfy n1 = n0 = n2) end n = [n1:n2]; x = [(n-n0) == 0]; 主程序： n0 = 0; n1 = -9; n2 = 9; [x,n] = impseq(n0,n1,n2); stem(n,x) title(δ(n)) xlabel(n (sample number)) ylabel(δ(n)) axis([-10 10 -1 2]) impseq.fig: b) 单位阶跃序列 2 u(n-5) stepseq.m: function [x,n] = stepseq(n0,n1,n2) % Generates x(n) = u(n-n0); n1 = n,n0 =n 2; % [x,n] = impseq (n0,n1,n2) if ((n0n1)||(n0n2)||(n1n2)) error(Arguments must satisfy n1 = n0 = n2) end n = [n1:n2]; x = [(n-n0) = 0]; 主程序： n0=5; n1=-9; n2=9; [x,n] = stepseq(n0,n1,n2);x=2*x; stem(n,x) title(2u(n-5)) xlabel(n (sample number)) ylabel(2u(n-5)) axis([-10 10 -1 3]) stepseq.fig: c) 矩形序列 R(n) recseq.m: function [x,n] = recseq(n1,n2,n3,n4) % R(n) will begin with n3, and end with n4; if (n3n4||n1n2||n1n3||n4n2) error(Arguments must satisfy n1 = n3 = n4 = n2) end n=[n1:n2]; x=[stepseq(n3,n1,n2)-stepseq(n4,n1,n2)]; 主程序： n1=-9; n2=9; n3=-3; n4=4; [x,n]=recseq(n1,n2,n3,n4); stem(n,x) title(R(n)) xlabel(n (sample number)) ylabel(R(n)) axis([-10 10 -1 2]) recseq.fig: d) y(n)＝2sin(0.3πn)+ 0.5cos2(0.6πn) yseq.m: function [x,n] = yseq(n1,n2) % R(n) will begin with n3, and end with n4; if (n1n2) error(Arguments must satisfy n1 = n3 = n4 = n2) end n=[n1:n2]; x=2*sin(0.3*pi*n)+0.5*cos(0.6*pi*n).^2; 主程序： n1=-9; n2=9; [x,n]=yseq(n1,n2); stem(n,x) title(2sin(0.3πn)+ 0.5cos2(0.6πn)) xlabel(n (sample number)) ylabel(y(n)) axis([-10 10 -2 3]) yseq.fig: (2) 数字信号的基本运算： 加、减、尺度（乘除）和移位是数字信号处理中最基本的算术运算，将上述基本序列进行这些基本运算，得到多个序列构成的组合序列。 n1=-9; n2=9; n=[n1:n2]; y1=impseq(0,n1,n2); y2=2*stepseq(5,n1,n2); y3=recseq(n1,n2,-3,4); y4=yseq(n1,n2); y5=y1+y2; subplot(2,2,1);stem(n,y5); title(δ(n)+2 u(n-5)) xlabel(n (sample number)) ylabel(y(n)) axis([-10 10 -1 3]) y6=y4-y2;