信号与系统MATLAB仿真-2012.pptVIP

《信号与系统的MATLAB仿真》 东南大学成贤学院? 2011级 范梦恒 2012.7. [4-1] 单位冲激信号的产生 t=-5:0.01:5; y=(t==0); subplot(121); plot(t,y,r); n=-5:5; x=(n==0); subplot(122); stem(n,x); [4-2] 比较不同的实指数信号 clear; n=0:1:10; a1=0.5;a2=-0.5;a3=1.2;a4=-1.2; x1=a1.^n;x2=a2.^n; x3=a3.^n;x4=a4.^n; subplot(2,2,1),stem(n,x1); title(‘实指数序列(0a11)); xlabel(n);ylabel(x1(n)); subplot(2,2,2),stem(n,x2); title(实指数序列(-1a20)); line([0,10],[0,0]); xlabel(n);ylabel(x2(n)); subplot(2,2,3),stem(n,x3); title(实指数序列(1a3)); xlabel(n);ylabel(x3(n)); subplot(2,2,4),stem(n,x4); line([0,10],[0,0]); title(实指数序列(a4-1)); xlabel(n);ylabel(x4(n)); [4-3] 生成负指数信号 clear; n=1:10; delta=-0.2;w0=0.7; x=exp((delta+j*w0)*n); subplot(211),stem(n,real(x)); line([0,10],[0,0]); title(‘负指数序列); ylabel(‘负指数序列实部); subplot(212),stem(n,imag(x)); line([0,10],[0,0]); ylabel(‘负指数序列虚部); xlabel(n); [4-4] 生成正弦信号 clear; n=0:10; w0=pi/5;w1=pi/4; x=sin(w0*n+w1); stem(n,x); title(‘正弦型序列); line([0,10],[0,0]); xlabel(n);ylabel(x(n)); [4-5] 产生周期为0.2的三角波，width取值分别为0,1,0.5。 td=1/100000; t=0:td:1; x1=sawtooth(2*pi*5*t,0); x2=sawtooth(2*pi*5*t,1); x3=sawtooth(2*pi*5*t,0.5); subplot(311);plot(t,x1); subplot(312);plot(t,x2); subplot(313);plot(t,x3); [4-6] 仔细观察由下面代码产生的3个三角波信号之间的区别 t=-3:0.001:3; x1=tripuls(t,4,0); subplot(131); plot(t,x1); axis([-4 4 0 1]); grid t=-6:0.001:6; x2=tripuls(t,4,0.5); subplot(132); plot(t,x2); axis([-4 4 0 1]); grid x3=tripuls(t+2,4,0.5); subplot(133); plot(t,x3); axis([-4 4 0 1]); grid [4-7] 产生频率为40HZ，占空比分别为25%，50%，75%的周期性方波 clear; td=1/100000; t=0:td:1; x1=square(2*pi*40*t,25); x2=square(2*pi*40*t,50); x3=square(2*pi*40*t,75); subplot(311);plot(t,x1); title(‘占空比25%);axis([0 0.2 -1.5 1.5]); subplot(312);plot(t,x2); title(占空比50%);axis([0 0.2 -1.5 1.5]); subplot(313);plot(t,x3); title(‘占空比75%);axis([0 0.2 -1.5 1.5]); [4-8] 生成幅度为2，宽度T=4，中心在t=0的矩形波x(t)以及x(t-T/2). t=-4:0.0001:4; T=4; x1=2*rectpuls(t,T); subplot(121);plot(t,x1); title(x(t));axis([-4 6 0 2.2]); x2=2*rectpuls(t-T/2,T); subplot(122);plot(t,x2); title(t-T/2);axis([-4 6 0