信号与系统实验matlab仿真!!.docVIP

信号与系统的亲们这可是实验利器哦！！ function f=Heaviside(t) f=(t0); function f=unitimpfun(t0,t) f=(t==t0); function f=unitimpseq(n0,n) if(n0==fix(n0)n0=10n0=-10) f=(n==n0); else printf(n0必须是-10到10整数); end function f=unitstepseq(n0,n) if(n0==fix(n0)n0=10n0=-10) f=(nn0); else printf(n0必须是-10到10的整数); end function f=rect(t) f=2*(Heaviside(t)-Heaviside(t-4)); function myfun1() t=-10:.001:10; n=-10:10; t0=input(单位连续冲击信号的平移参数 -10到10直接之间); n0=input(单位离散冲击信号和单位离散阶跃信号的平移参数 -10到10之间); figure(1) plot(t,Heaviside(t)); axis([-10 10 -0.1 1.1]); title(单位阶跃信号); figure(2) plot(t,rect(t)); axis([-10 10 -0.1 2.1]); title(矩形脉冲信号); figure(3) plot(t,unitimpfun(t0,t)); axis([-10 10 -0.1 1.1]); title(单位冲击信号); figure(4); stem(n,unitimpseq(n0,n)); axis([-10 10 -0.1 1.1]); title(单位冲击信号列 离散); figure(5); stem(n,unitstepseq(n0,n)); axis([-10 10 -0.1 1.1]); title(单位阶跃系列 离散); 实验一：单位阶跃信号 实验一：矩形脉冲信号 实验一：向右平移三个单位的单位冲击信号 实验一：向右平移两个单位的单位离散冲击信号 实验一：向右平移三个单位的离散阶跃信号 function f=f_signal2(r,n) f=r.^n.*exp(n*j*pi/4); function f=signal2(x,t) f=2*exp((-x+3*j*pi)*t); function myfun2() t=0:0.001:3.2; figure(1); subplot(221); plot(t,real(signal2(0,t))); ylabel(x为0实部); subplot(222); plot(t,imag(signal2(0,t))); ylabel(x为0虚部); subplot(223); plot(t,abs(signal2(0,t))); ylabel(x为0模); subplot(224); plot(t,angle(signal2(0,t))); ylabel(x为0角度); figure(2); subplot(221); plot(t,real(signal2(1,t))); ylabel(x=1实部); subplot(222); plot(t,imag(signal2(1,t))); ylabel(x=1虚部); subplot(223); plot(t,abs(signal2(1,t))); ylabel(x=1模); subplot(224); plot(t,angle(signal2(1,t))); ylabel(x=1角度); n=0:32; figure(3); subplot(211); stem(n,cos(pi/8*n)); ylabel(cos(nπ/8)); grid on; subplot(212); stem(n,cos(n/2)); ylabel(cos(n/2)); grid on; n=0:20; figure(4); subplot(311); stem(n,f_signal2(1.2,n)); ylabel(1.2^n*e^(jnπ/4)); subplot(312); stem(n,f_signal2(1,n)); ylabel(1^n*e^(jnπ/4)); subplot(313); stem(n,f_signal2(0.9,n)); ylabel(0.9^n*e^(jnπ/4)); 实验二：δ=0时候的所求图形 实验二：δ=1时候的所求图形 实验二：两个离散正弦信号