信号与系统部分实验matlab代码-供参考分析.doc

信号与系统第三版112学期部分实验代码 目录 实验一 常用信号的实现和时域变换部分实验答案 1 实验二 LTI时间系统的时域分析 3 实验四 连续信号与系统的频域分析 4 实验七 连续时间信号的采样与重构 15 实验一 常用信号的实现和时域变换部分实验答案 （1） …………………………………………………………………………………………… t=0:0.1:30; x=exp(-0.1*t).*sin(2/3*t); plot(t,x); axis([0 30 -1 1]); grid; xlabel(Time(sec)); ylabel(f(t)); …………………………………………………………………………………………… 对上面的f(t)进行抽样f(n)=f(t)|t=nT=x(nT)就可以得到抽样信号，将上述命令里的时间增量改为1s，plot(t,x)命令改为 可以得到下图。 4．已知离散序列x(n)如图所示，试用Matlab编程绘出以下的离散序列波形。 （1）x(-n+3) clc; close all; clear all; n=-3:2 x=[3 3 3 3 2 1]; [x1,n1]=xlpy(x,n,-3); [x2,n2]=xlfz(x1,n1); subplot(311) stem(n,x,filled); axis([min(n)-1,max(n)+1,min(x)-0.5,max(x)+0.5]) subplot(312) stem(n1,x1,filled); axis([min(n1)-1,max(n1)+1,min(x1)-0.5,max(x1)+0.5]) subplot(313) stem(n2,x2,filled); axis([min(n2)-1,max(n2)+1,min(x2)-0.5,max(x2)+0.5]) (2) x(n-3)x(n+2) clc; close all; clear all; n=-3:2 x=[ 3 3 3 3 2 1 ]; [x1,n1]=xlpy(x,n,3); [x2,n2]=xlpy(x,n,-2); nn=-5:5; x3=[0 0 0 0 0 x1]; x4=[x2 zeros(1,5)]; x=x3.*x4; subplot(311) stem(nn,x3,filled); subplot(312) stem(nn,x4,filled); subplot(313) stem(nn,x,filled); 实验二 LTI时间系统的时域分析 6.对于因果和稳定的LTI系统，对于下列二阶微分方程确定其单位冲激响应是否是欠阻尼、过阻尼或临界阻尼，画出系统的h(t)和频率响应模的波特图。 (1) (2) (3) (4) …………………………………………………………………………………………… clc; close all; clear all; num=[1]; den=[1 4 4]; H=tf(num,den); bode(H); hold on; num=[7]; den=[5 4 5]; H=tf(num,den); bode(H); hold on; num=[1/3 7]; den=[5 4 5]; H=tf(num,den); bode(H); …………………………………………………………………………………………… 实验四 连续信号与系统的频域分析 1.利用那个fourier函数求下列信号的傅里叶变换F(jw)，并用ezplot函数绘出其幅度谱和相位谱. （1） close all; clear all; syms t v w phase im re; % 定义变量t,v,w,phase,im re f=exp(-3*abs(t))*sin(2*t); % Fw=fourier(f); subplot(311); ezplot(f); % 画-2*pi到2*pi内函数 axis([-0.01 2 0 0.5]); subplot(312); ezplot(abs(Fw)); im=imag(Fw); re=real(Fw); phase=atan(im/re); subplot(313); ezplot(phase); axis([-6 6 -0.5 0.5]); （2） close all; clear all; syms t v w phase im re; % 定义变量t,v,w,phase,im re f=t*exp(-2*t)*sin(4*t)*sym(Heaviside(