及起始状态r(k)(0-) （k=0,1,2,…,n-1）的解；它是齐次解中的一部分：
rzi(t)??Azikeakt，由于没有外界激励，系统状态不会发生变化，即：r(k)(0-)=r(k)(0+)
k?1n零状态响应：不考虑起始时刻（t=0-）系统储能的作用，仅由系统的外加激励信号产生响应，满足：
dndn?1dC0nrzs(t)?C1n?1rzs(t)???Cn?1rzs(t)?Cnrzs(t)dtdtdt
dmdm?1d?E0me(t)?E1m?1e(t)???Em?1e(t)?Eme(t)dtdtdt及起始状态 r(k)(0-)=0（k=0,1,…,n-1）的解；它是由特解和一部分齐次解构成：
rzs(t)??Azskeakt?B(t)。
k?1n零状态响应也可由卷积的方法得出：rzs(t)?e(t)?h(t)
五、试验报告：
试验报告应包含以下内容
试验目的、内容、步骤、原理、程序、图形（可打印）等。 附录
%试验二,内容1 %产生单位冲激信号
clear all;t0=0;t1=1;dt=0.05; t=0:0.05:5;st=length(t); n1=floor((t1-t0)/dt); x1=zeros(1,st); x1(n1)=1/dt;
subplot(4,1,1),stairs(t-0.95,x1),grid; axis([-1,4,-1,22]) %产生单位阶跃信号
x2=[zeros(1,n1-1),ones(1,st-n1+1)]; subplot(4,1,2),stairs(t-0.95,x2),grid; axis([-1,4,-0.1,1.1]) %产生周期的三角波 clear all;
t=-4*pi:0.01:4*pi; a=sawtooth(t,0.5);
subplot(4,1,3),plot(t,a),axis([-4*pi,4*pi,-1.1,1.1]),grid; %产生周期的方波
5
b=square(t,50);
subplot(4,1,4),plot(t,b),axis([-4*pi,4*pi,-1.1,1.1]),grid; %实验二，内容2
a=input('输入分母系数向量a＝[a1,a2,...]='); n=length(a)-1;
y0=input('输入初始条件向量 y0＝[y0,dy0,...]='); p=roots(a);
v=rot90(vander(p)); c=v/y0; dt=0.2; tf=8; t=0:dt:tf;
y=zeros(1,length(t));
for k=1:n y=y+c(k)*exp(p(k)*t);end; plot(t,y),grid;
%实验二，内容3
a=input('输入多项式分母系数向量 a＝[a1,a2,...]='); b=input('输入多项式分子系数向量 b＝[b1,b2,...]='); u=input('输入信号 u(t)='); dt=0.2;tf=8; t=0:dt:tf;
[r,p,k]=residue(b,a); h=zeros(1,length(t)); for i=1:length(a)-1;
h=h+r(i)*exp(p(i)*t);end
subplot(2,1,1),plot(t,h),grid;title('冲激响应'); y=conv(u,h);
subplot(2,1,2),plot(t,y(1:length(t))),grid;title('零状态响应');
6
试验三 连续系统的频域分析
一、试验目的：
分析傅里叶级数以及利用傅氏变换研究信号频谱。
二、试验内容及要求：
编制Matlab程序，实现以下内容
1 傅氏级数：将周期为T=2?，幅度为?的对称方波分解为多次正弦波之和，用Matlab演示谐波合成情况。（取到7次谐波） 2 利用Matlab分析矩形脉冲信号的频谱； 求f(t)=u(t+0.5)-u(t-0.5)的傅立叶变换
三、实验步骤：
打开电脑，双击桌面Matlab程序图标，打开文件菜单，点击新建M文件，在M文件编辑器中编写程序，编写完毕后按F5键运行调试。
四、实验原理：
傅式级数所选项数越多，相加后波形越逼进原信号，二者方均误差越小；当信号是脉冲信号时，其高频分量主要影响脉冲的跳变沿，低频分量主要影响脉冲的顶部，即信号波形变化越剧烈，包含的高频成分越丰富；变化越缓慢，包含的低频成分越丰富；当信号中任一频谱分量的幅度或相位发生相对变化时，输出波形一般要发生失真。
傅立叶变换对：
F(?)??f(t)e?j?tdt——傅立叶正变换；
???1?j?t??f?t??F?ed?——傅立叶逆变换，简写f?t??F??? ???2? 7
傅立叶变换的物理意义在于信号可分解为无穷多个幅度为无穷小的连续指数信
?1???F?d??之和，整个频域-∞→∞，由于幅度无穷小，所以频谱不能再用幅号??2??度表示，而改用密度函数表示。
五、试验报告要求：
试验报告应包含以下内容
试验目的、内容、步骤、原理、程序、图形（可打印）等。 附录
%实验三，内容1 function [h]=my301(a); t=0:0.01:2*pi; nn=a; y1=0; for i=1:nn; xb=i*2-1;
subplot(2,1,1), y=sin(xb*t)/xb;
plot(t,y),grid,hold on;title('对称方波傅式级数'); subplot(2,1,2) y1=y1+y; end plot(t,y1),grid h=0;%实验三，内容2 clear,tf=1; N=256;
dt=1/N;t=[-N+1:N]*dt;
f=[zeros(1,N/2),ones(1,N),zeros(1,N/2)]; wf=50; nf=64;
w1=linspace(0,wf,nf);dw=wf/(nf-1);
8
f1=f*exp(-j*t'*w1)*dt; w=[-fliplr(w1),w1(2:nf)]; f2=[-fliplr(f1),f1(2:nf)];
subplot(1,2,1),plot(t,f,'linewidth',1.5),axis([-1,1,-0.1,1.1]),grid;title('时域方波'); subplot(1,2,2),plot(w,abs(f2),'linewidth',1.5),grid;title('方波频谱');
9