实验三、四程序.docVIP

实验三、信号基本运算 1、已知，绘制和的图形。 MATLAB程序代码： t=linspace(-5,5,100); y1=exp(-0.5*t); y2=sin(10*t); y3=y1+y2; y4=y1.*y2; plot(t,y3,r:) hold on; plot(t,y4) 2、利用impz函数，计算以下差分方程表示的离散系统的单位脉冲响应，并绘制前30点的波形。 a=[1,0.7,-0.45,-0.6]; b=[0.8,-0.44,0.36,0.02]; k=0:30; x=impz(b,a,k); stem(x); grid on; 单位脉冲响应的波形如下： 3、已知输入信号为，系统的单位冲激响应为，利用卷积和求系统的响应y[k],并绘制、、的波形 程序如下： n=[1:100]; x=[1,1,1,1,1,1,1,1,1,1,zeros(1,90)];%输入信号 h=0.9.^n;%系统冲激响应 y=conv(x,h);%系统输出 k=[1:199]; subplot(3,1,1);stem(n,x); subplot(3,1,2);stem(n,h); subplot(3,1,3);stem(k,y); 输入信号、冲激响应、输出信号的波形分别如下： 4 、已知某连续时间系统的微分方程为，输入信号，初始状态为，，试求系统的零输入响应、零状态响应和完全响应，并绘制响应的波形。 解 由微分方程可知，系统的极点为-1和-3，再结合初始条件，可以求解系统的零输入响应 则零输入响应为： 程序如下： sys=tf([2 1],[1 4 3])%系统函数对于的系数向量 t=1:0.05:10;kt=length(t); f=ones(1,kt) %输入信号 yf=lsim(sys,f,t);%系统的零状态响应 subplot(3,1,1),plot(t,yf);title(零状态响应); yx=2.5*exp(-t)-1.5*exp(-3*t);%系统的零输入响应 subplot(3,1,2),plot(t,yx);title(零输入响应); y=yx+yf;%lsim计算得到的零状态响应为列向量，零输入相应为行向量，故需要转置 subplot(3,1,3),plot(t,y);title(完全响应); 系统的零输入响应、零状态响应、完全响应如下： 5、已知周期三角信号，请分解为傅里叶级数，并用前面7次谐波合成。(幅值在正负1之间交替，周期4) MATLAB程序代码： %由已知： %周期三角函数为偶函数，bn=0; %an=1/n^2*(1-cos(2n)) %f=-2+1/n^2*(1-cos(2*n))*cos(n*t) %f=-2+1/1^2*(1-cos(2*1))*cos(t)+1/2^2*(1-cos(2*2))*cos(2t)+1/3^2*(1-cos(2*3))* %cos(3t)+1/4^2*(1-cos(2*4))*cos(4t)+1/5^2*(1-cos(2*5))*cos(5t)+... %程序代码如下： %绘制三角波型 t1=-2:0.001:2; tt=-1:0.001:1; ft=tripuls(t1,4,0); y=2.*ft-1; %计算傅里叶级数 t=-4:0.001:4; c1=1./1^2.*(1-cos(2*1)).*cos(1*t);%一次谐波 c3=1./3^2*(1-cos(2*3))*cos(3*t);%三次谐波 c5=1./5^2*(1-cos(2*5))*cos(5*t);%五次谐波 c7=1./7^2*(1-cos(2*7))*cos(7*t);%七次谐波 subplot(2,2,1) %绘第1张图 hold on; plot(t,0,b:); plot(0,tt,b:); plot(t1,y,r); plot(t,c1);title(基波+1次谐波) subplot(2,2,2) %绘第2张图 hold on plot(t,0,b:); plot(0,tt,b:); plot(t1,y,r); plot(t,c1+c3);title(基波+1次谐波+3次谐波) subplot(2,2,3) %绘第3张图 hold on plot(t,0,b:); plot(0,tt,b:); plot(t1,y,r); plot(t,c1+c3+c5);title(基波+1次谐波+3次谐波+5次谐波) subplot(2,2,4) %绘第4张图 hold on plot(t,0,b:); plot(0,tt,b:); plot(t1,y,r); plot(t,c1+c3+c5+c7);title(基波+1次谐波+3次谐波+5次谐波+7次谐波) MAT