1、 代码：t=0:10; f=t; plot（t,f）t=1:5 f=2-exp（-1.*t.）;80.1:2; f=exp（-1.*t.）*cos（10*pi*t）;（3） 三种奇异函数符号函数 代码: t=-5:0.05:5; f=sign（t）; plot（t,f）阶跃信号 t=-5: f=u（t）;单位冲激信号function chongji（t1,t2,t0）dt=0.01;t=t1:dt:t2;n=length（t）;x=zeros（1,n）;x（1,（-t0-t1）/dt+1）=1/dt;stairs（t,x）;axis（t1,t2,0,1.2/dt） title（t） （4）实验
2、三1234syms t f1=sym（-t+4）*（u（t）-u（t-4）; subplot（1,2,1）;ezplot（f1）;y1=subs（f1,t,-t）; f3=f1+y1;subplot（1,2,2）;ezplot（f3）;function f=u（t） f=（t0）;4、function f=u（t）f=（t0）syms tsubplot（1,3,1）;f2=sym（sin（2*pi*t）subplot（1,3,2）;ezplot（f2）;f6=f1.*f2;subplot（1,3,3）;ezplot（f6）;5、y6=subs（f6,t,t-2）;ezplot（y6）;f7=y
3、6+f2;ezplot（f7）;四、500;f=100.*abs（sin（2.*pi.*t./50）;plot（t,f,t,fD,t,fA）调用子程序：function fD=fDC（f）fD=mean（f）;function fA=fAC（f,fD）fA=f-fD; （5）求解信号的交直流分量 fD=mean（f）; fA=f-fD;f（t）=100|sin（2*PI*t/50）|;plot（t,fD,t,fA）实验二 线性系统时域分析（1）求解下面两个信号的卷积积分。要求：1） 在实验报告中推导出这两个信号卷积积分运算表达式；（手写） 2） 利用MATLAB 进行求解验证，附程序代码和波形
4、。（2）已知描述系统的微分方程和激励信号如下r（t ） + 3r（t ） +2r（t ） = e（t）+3 e（t） ，e（t）=u（t） 。1） 用解析法求系统的零状态响应r（t ）；2）利用MATLAB绘出系统零状态响应的时域仿真波形，并验证1）的结果是否正确，附程序代码和波形； 3）利用MATLAB绘出系统的冲激响应和阶跃响应波形，附程序代码和波形。在实验报告中推导出这两个信号卷积积分运算表达式；利用MATLAB 进行求解验证，附程序代码和波形。p=0.01;k1=-1/2:p:1;f1=1;k2=0:f2=0.5*k2;f,k=sconv（f1,f2,k1,k2,p）function
5、f,k=sconv（f1,f2,k1,k2,p）f=conv（f1,f2）;f=f*p;k0=k1（1）+k2（1）;k3=length（f1）+length（f2）-2;k=k0:（k3*p+k0）;subplot（2,2,1） plot（k1,f1） xlabel（t） ylabel（f1（t）subplot（2,2,2） plot（k2,f2） f2（t）subplot（2,2,3）plot（k,f）;h=get（gca,positionh（3）=2.5*h（3）;set（gca,h） title（f（t）=f1（t）*f2（t） xlabel（f（t） 用解析法求系统的零状态响应r（t
6、 ）；利用MATLAB绘出系统零状态响应的时域仿真波形，并验证的结果是否正确，附程序代码和波形； a=1,3,2; b=0,1,3; impulse（a,b）利用MATLAB绘出系统的冲激响应和阶跃响应波形，附程序代码和波形。b=0,1,3;step（b,a） 代码； p=0.01; t=0: x=exp（-t）; lsim（b,a,x,t）实验三：实验报告内容：对所给音频信号，进行时域压缩和扩展，画出时域波形与幅度谱，使其满足以下要求。（1）将music1.wav的音调变低a倍（0.8a1），及变高a倍（11.3），比较变换前后的10%带宽（最大值10%处的带宽）的变化情况。（2）将musi
7、c2.wav的音调变化到与儿童和男声相似的声音，分析变换前后的10%带宽的变化情况（变宽或变窄多少）。a=0.9 1.0 1.2;a1=30;b=11025;x0 b=wavread（F:信号与系统music1.wavfor m=1:3;a2=a1*a（m）;x=resample（x0,a1,a2） ;wavplay（x ,b）;1/b:（length（x）-1）/b;subplot（3,3,3*m-2）;plot（t,x）;axis（0 6 -inf inf）;title（a=,num2str（a（m）,的时域图）;N=（length（resample（x0,a1,a1*max（a）;f=（
8、-N/2:N/2-1）/N*b;X= fftshift（fft（x,N）;AMP=abs（X）;AMP0=AMP/max（AMP）;subplot（3,3,3*m-1）;plot（f,AMP0）;的频域图subplot（3,3,3*m）;axis（-inf inf 0 0.1）的10%带宽图enda=1 0.7 1.4;str=原 男 童;信号实验_2018music2.wav a2=a1*a（m）; x=resample（x0,a1,a2）; wavplay（x,b）; subplot（3,3,3*m-2）; plot（t,x）; axis（0 3 -inf inf） title（str（m
9、）,声的时域图 N=（length（resample（x0,a1,a1*max（a）; f=（-N/2: X=fftshift（fft（x,N）; AMP=abs（X）; AMP0=AMP/max（AMP）; subplot（3,3,3*m-1）; plot（f,AMP0）;声的频域图 subplot（3,3,3*m）; axis（-inf inf 0 0.1）声的10%带宽图 ）分析：扩展0.9倍，变换前后的10%带宽由2284.54Hz变为2056.16Hz，变化率：变窄9.997%；压缩1.2倍，变换前后的10%带宽由2284.54Hz变为2944.61Hz，变化率：变宽28.893%；
10、（2）男声：扩展0.7倍，变换前后的10%带宽由2394.9Hz变为1685.7Hz，变化率：变窄29.61%；童声：压缩1.4倍，变换前后的10%带宽由2394.9Hz变为3353.1Hz，变化率：变宽40.01%；实验七 连续系统零极点分析 实验报告内容：（1）根据例7-1、7-2、7-3的要求和提示完成。（2）完成第四部分实验内容中的（1）（3）。 a=8,2,3,1,5;b=1,3,2;q,p=sjdt（a,b） SJDT函数functionp,q=sjdt（A,B）p=roots（A）;q=roots（B）;p=p;q=qx=max（abs（p q）;x=x+0.1;y=x;clfh
11、old onaxis（-x x -y y）;axis（squareplot（-x x,0 0）plot（0 0,-y y）plot（real（p）,imag（p）,xplot（real（q）,imag（q）,otext（0.2,x-0.2,text（x-0.2,0.2,7-2（a） a=1 0; b=1; impulse（b,a）（b） a=1 2;（c）a=1 -2;b=1;impulse（b,a） （d） a=1 1 16.25;impulse（b,a,5）（e）a=1 0 16;（f） a=1 -1 16.25;7-3function splxy（f1,f2,k,p,q） q=q f=f
12、1:k:f2;w=f*（2*pi）;y=i*w; n=length（p）;m=length（q）;if n=0 yq=ones（m,1）*y;vq=yq-q*ones（1,length（w）;bj=abs（vq）;cosaij=angle（vq）./pi.*180;ai=1;thetai=0;elseif m=0 yp=ones（n,1）*y;vp=yp-p*ones（1,length（w）;ai=abs（vp）;thetai=angle（vp）./pi.*180;bj=1;cosaij=0;else end subplot（121）;Hw=prod（bj,1）./prod（ai,1）;plo
13、t（f,Hw）;连续系统幅频响应曲线） xlabel（频率 w（单位：赫兹） ylabel（F（jw） subplot（122）;Angw=sum（cosaij,1）-sum（thetai,1）;plot（f,Angw）;连续系统相频响应曲线Angle（jw）q=0;p=-100 -50;f1=0; f2=100;k=0.01;splxy（f1,f2,k,p,q）;a=3 5 4 -6;b=1 1 2;p,q=sjdt（a,b）不稳定q= -0.5000 - 1.3229i -0.5000 + 1.3229i;p= -1.1742 - 1.2472i -1.1742 + 1.2472i 0.6
14、816 + 0.0000i;f2=3;（3）a=1 0 20 0 64;b=3 0 -27 0;q=0 3 -3;p= -0.0000 - 4.0000i -0.0000 + 4.0000i 0.0000 - 2.0000i 0.0000 + 2.0000i;f2=1.5;实验八 离散系统z域分析（1）根据例1-例5的要求和提示完成。（2）完成第五部分的实验内容。8-1（1）A=1 -3 7 -5; B=3 -5 10 0; ljdt（A,B）调用的ljdt函数function ljdt（A,B）x=max（abs（p q 1）;w=0:pi/300:2*pi;t=exp（i*w）;plot（
15、t）text（0.1,x,jImztext（y,1/10,Rezpole-zero diagram for discrete systemhold off8-1（2） A=1 0.75 0.125;B=1 -0.5 0;ljdt（A,B）8-2例2：解：由例1绘出的零极点图可以看出两个系统的稳定性分别为：第（1）个系统不稳定；第（2）个系统稳定。8-3（1）a=1 -1;b=0 1; impz（b,a,10）8-3（2） a=1-2*cos（pi/8） 1; b=0 0 1;impz（b,a,50）8-2（3）a=1 -2 1;b=0 1 0;impz（b,a,10）8-2（4）a=1 -0.
16、8;8-2（5）a=1 -1 0.25;b=0 0 1;8-2（6）a=1 -1.2;例4：MATLAB 命令如下：B=1 -0.5;A =1 0;H,w=freqz（B,A,400,wholeHf=abs（H）;Hx=angle（H）;figure（1）plot（w,Hf）离散系统幅频特性曲线figure（2）plot（w,Hx）离散系统相频特性曲线例五：调用子程序;function dplxy（k,r,A,B）%The function to draw the frequency response of discrete systemljdt（A,B） r*pi/k:r*pi;y=exp（
17、i*w）;N=length（p）;M=length（q）;yp=ones（N,1）*y;yq=ones（M,1）*y;vp=yp-p*ones（1,k+1）;vq=yq-q*ones（1,k+1）;Ai=abs（vp）;Bj=abs（vq）;Ci=angle（vp）;Dj=angle（vq）;fai=sum（Dj,1）-sum（Ci,1）;H=prod（Bj,1）./prod（Ai,1）;plot（w,H）;角频率幅度figure（3）plot（w,fai）离散系统的相频特性曲线相位A=1 -1/4;B=5/4 -5/4;dplxy（500,2,A,B） 实验内容已知离散系统的系统函数分别为：试用MATLAB 分析：（1） 绘出系统的零极点图，根据零极点图判断系统的稳定性；（2） 如果系统稳定，绘出幅频特性和相频特性曲线。A=2 0 0 -1;B=0 1 -2 -1;稳定A=1 2 -4 1;B=0 1 0 2;
