[工学]信号与线性系统 答案.doc

实验一 信号的MATLAB表示 实验内容： 1. 用MATLAB表示连续信号：，，。 t=0:001:10; A=1; a=-0.4; ft=A*exp(a*t); plot(t,ft) t=0:0.1:10; A=1; a=1; b=pi/4; ft=A*sin(a*t+b); plot(t,ft) t=0:0.1:10; A=1; a=1; b=pi/4; ft=A*cos(a*t+b); plot(t,ft) 2. 用信号处理工具箱提供的函数表示抽样信号、矩形脉冲信号及三角脉冲信号。 y=sinc(t) y=sinc(t); plot(t,y) y=rectpuls(t, width) t=0:0.01:4; T=1; y=rectpuls(t-2*T, 2*T); plot(t,y) y=tripuls(t, width, skew) t=-5:0.01:5; width=2;skew=0.6; y=tripuls(t, width, skew); plot(t,y) 3. 编写如图所示的MATLAB函数，并画出，的图形。 t=-2:0.01:3; ft=rectpuls(t+0.5, 1)+(1-t).*rectpuls(t-0.5,1)-rectpuls(t-1.5, 1); plot(t,ft) function ft=f(t) ft=rectpuls(t+0.5, 1)+(1-t).*rectpuls(t-0.5,1)-rectpuls(t-1.5, 1); plot(t,ft) t=-5:0.01:5; y=f(0.5*t); plot(t,y) t=-2:0.01:8; y=f(2-0.5*t); plot(t,y) 4. 用MATLAB表示离散信号：，。 k=0:0.2:10; a=2; y=a.^k; stem(k,y) k=0:0.2:10; A=2; y=A*sin(k); stem(k,y) 实验二 连续时间系统的时域分析 三、实验内容： 1. 分别用函数lsim和卷积积分两种方法求如图所示系统的零状态响应。其中L=1，R=2，，。 用函数lsim L=1;R=2; sys=tf([1],[L R]); t=0:0.1:10; f=exp(-t); y=lsim(sys,f,t); plot(t,y); xlabel(Time(sec)) ylabel(y(t)) 用卷积积分 L=1;R=2; sys=tf([1],[L R]); t=0:0.001:10; f=exp(-t); y=impulse(sys,t); z=conv(y,f); plot(z) 2. 求上述系统的冲激响应与阶跃响应。 冲激响应 L=1;R=2; sys=tf([1],[L R]); t=0:0.1:10; f=exp(-t); y=impulse(sys,t); plot(t,y) xlabel(Time(sec)) ylabel(y(t)) 阶跃响应 L=1;R=2; sys=tf([1],[L R]); t=0:0.1:10; f=exp(-t); y=step(sys,t); plot(t,y) xlabel(Time(sec)) ylabel(y(t)) 实验三 信号的Fourier分析 三、实验内容： 1. 求如图所示周期矩形脉冲信号的Fourier级数表达式，画出频谱图，并用前N次谐波合成的信号近似。 画出周期矩形脉冲信号的频谱图 N=8; n1=-N:-1; c1=-2*j*sin(n1*pi/2)/pi./n1; c0=0; n2=1:N; c2=-2*j*sin(n2*pi/2)/pi./n2; cn=[c1 c0 c2]; n=-N:N; subplot(1,2,1); stem(n,abs(cn)); ylabel(Cn的幅度); xlabel(\Omega); subplot(1,2,2); stem(n,angle(cn)); ylabel(Cn的相位); xlabel(\Omega); 用前N次谐波合成的信号近似 建立M函数f.m function s=f(m) s=0;p=0; t=0:0.01:10; for n=1:2:m p=p+1/n*(sin(n*t)); end s=(4/pi)*p; plot(t,s) f(20); xlabel(Time(sec)) ylabel(f(t)) 2. 求信号的幅度谱。 w=-2*pi:0.2:2*pi; c=(2+j*w).^(-2); plot(w,abs(c)); ylabel(C的幅度); xlabel(\Omega); 实验四 连续系统的频域分析 三、实验内容： 如图所示系统： (1) 对不同的RC值，用freqs函