Matlab讲义实验报告连续时间信号的分析.docx

1、Matlab讲义实验报告连续时间信号的分析连续时间信号的分析一、实验目的1.学习使用MATLAB产生基本的连续信号、绘制信号波形。2.实现信号的基本运算，为信号分析和系统设计奠定基础。二、实验原理1、基本信号的产生时间间隔代替连续信号。连续指数信号的产生连续矩形脉冲信号(门信号)的产生。连续周期矩形波信号的产生。2、信号的基本运算相加、相减、相乘、平移、反折、尺度变换。三、实验内容1. 用MATLAB编程产生正弦信号，并画图。 代码如下：clcclearf0=5;w0=2*pi*f0;t=0:0.001:1;x=2*sin(w0*t+pi/3);plot(t,x) title(正弦信号)2.

2、用MATLAB编程产生信号，画出波形。代码如下：clcclearf0=2;t=0:0.0001:2.5;y=square(w0*t,50);plot(t,y);axis(0 2.5 -1.5 1.5)title(周期方波);图形如下：3. 分别画出2中移位3个单位的信号、反折后的信号、尺度变换后的信号。代码如下：clccleart=-10:0.001:10;subplot(3,1,1)plot(t,f(t-3)axis(-7 7 -2 2)xlabel(t)ylabel(f(t-3)title(移位)grid onsubplot(3,1,2)plot(t,f(-t)axis(-7 7 -2 2

3、)xlabel(t)ylabel(f(-t)title(反折)grid onsubplot(3,1,3)plot(t,f(3*t)axis(-7 7 -2 2)xlabel(t)ylabel(f(3t)title(尺度变换)grid on图形如下：4. 用MATLAB编程画出下图描述的函数。分别画出的图形。代码如下：clccleart=-10:0.001:10;subplot(5,1,1)plot(t,fc(t)axis(-2 4 -1 3)xlabel(t)ylabel(f(t)title(f(t)grid onsubplot(5,1,2)plot(t,fc(1-t)axis(-4 3 -1

4、 3)xlabel(t)ylabel(f(1-t)grid onsubplot(5,1,3)plot(t,fc(3*t+2)axis(-2 2 -1 3)xlabel(t)ylabel(f(3t+2)grid onsubplot(5,1,4)plot(t,fc(2-t/3)axis(-6 10 -2 3)xlabel(t)ylabel(f(2-t/3)grid onF=(fc(t)+f(2-t).*heaviside(1-t)subplot(5,1,5)plot(t,fc(2-t/3)axis(-8 10 -1 2)xlabel(t)grid on图形如下：四、实验思考通过实验我熟悉了各个函数

5、功能并掌握各种变换有助于进一步了解信号的各种变化的方式和意义。利用Matlab可以求连续、有限时域的零状态响应，求连续、有限时域的零状态响应、冲激响应和阶跃响应。Matlab可以大大提高效率，对于某些信号或响应的给出形象的物理解释，给自己一个直观的感受。实验二 连续时间系统的时域分析一、实验目的1. 加深对线性时不变系统中零状态响应概念的理解，掌握其MATLAB的求解方法。2. 加强对卷积运算的理解，验证卷积的一些重要性质。3. 会求给定连续系统的单位冲击响应和阶跃响应。二、实验原理1. 卷积积分定义 (1) (2)MATLAB conv()函数：计算两个离散时间序列的卷积。性质：交换律 结合

6、律 分配律C=conv(a,b)。连续卷积MATLAB求解functiony,k=sconv(f1,f2,k1,k2,T)y=conv(f1,f2);y=y*T;k_start=k1(1)+k2(1);k_end=length(f1)+length(f2)-2;k=k_start:T:(k_start+k_end*T);给定两个函数 和 求。clcclearT=0.001;k1=-1:T:5;f1=2*(k10);k2=-1:T:3;f2=(k20&k22);y,k=sconv(f1,f2,k1,k2,T);subplot(3,1,1)plot(k1,f1);axis(-1,5,0,2.2);

7、title(f1);subplot(3,1,2)plot(k2,f2);axis(-1,3,0,1.2);title(f2);subplot(3,1,3)plot(k,y);axis(min(k),max(k),min(y),max(y)+0.2);title(y=f1*f2);2. 零状态响应、单位冲击响应、单位阶跃响应概念：输出响应=零输入响应(初始状态)+零状态响应(输入信号)LTI系统：。阶跃响应：单位阶跃信号输入+零状态 (4) MATLAB语句：传递函数建立 sys=tf(b,a) b=b3,b2,b1,b0 a=a3, a2, a1, a0;脉冲响应：h=impulse(sys,

8、 t); t-采用时间任意输入零状态响应：h=lsim(sys,f,t);f-任意输入函数单位阶跃响应：h=step(sys,t);求微分方程的单位冲击响应。比较理论值和计算值。clear clct=0:0.01:5;b=1,3;a=1,3,2;sys=tf(b,a);y=impulse(sys,t);subplot(2,1,1);plot(t,y)title( MATLAB计算解)y1=2*exp(-t)-exp(-2*t);subplot(2,1,2)plot(t,y1);title(理论解)figure(2)plot(t,abs(y-y1);title(误差)讲解 plot语句可以画横向

9、量也可以画纵向量。向量、矩阵转置 “”。abs 求绝对值。三、实验内容1.已知 和 求，并且画出、以及波形。代码如下：clcclearT=0.001;k1=-3:T:3;f1=1*(k10);k2=-3:T:3;f2=2*(k10);y,k=sconv(f1,f2,k1,k2,T);subplot(3,1,1)plot(k1,f1);axis(-3,3,0,2);subplot(3,1,2)plot(k2,f2);axis(-3,3,0,3);subplot(3,1,3)plot(k,y);axis(-3,3,0,5)图形如下：2.已知系统的微分方程为，求该系统的单位冲激响应并画出相应的波形。

10、代码如下：clearclct=0:0.01:5;b=0,1;a=1,6,9;sys=tf(b,a);y=impulse(sys,t);subplot(2,1,1);plot(t,y)title( MATLAB 计算解)y1=2*exp(-t)-exp(-2*t);subplot(2,1,2)plot(t,y1);title(理论解)图形如下：3.已知系统的微分方程为，当外加激励信号为时，系统的零状态响应并画出相关的信号波形。代码如下：T=0.001;t1=0:T:5;b=3;a=1,5,6;sys=tf(b,a);x1=1*impulse(sys,t1);t2=-1:T:11x2=1*exp(

11、-t2).*(t2=0&t2=10);y,k=sconv(x1,x2,t1,t2,T);subplot(3,1,1)plot(t2,x2);axis(min(t2),max(t2),min(x2),max(x2)+0.2);subplot(3,1,2)plot(t1,x1);axis(min(t1),max(t1),min(x1),max(x1)+0.2);subplot(3,1,3)plot(k,y);图形如下：4.已知和如图3.1所示，求，并且画出波形。代码如下：T=0.001t1=0:T:3;t2=0:T:4;f1=1*(t10);f2=1*(t20).*(t2/4);y,k=sconv

12、(f1,f2,t1,t2,T)t3=-1:T:1;f3=2*(t30);plot(t3,f3)图形如下：四、实验思考通过实验我熟悉了各个函数功能并掌握各种变换有助于进一步了解信号的各种变化的方式和意义。加深了对线性时不变系统中零状态响应概念的理解，掌握其MATLAB的求解方法。加强了对卷积运算的理解，验证卷积的一些重要性质。会求给定连续系统的单位冲击响应和阶跃响应实验三 连续时间系统的频域分析一、 实验目的4. 理解周期信号的傅里叶分解，掌握傅里叶系数的计算方法。5. 深刻理解和掌握非周期信号的傅里叶变换和其计算方法。6. 熟悉傅里叶变换的性质，并能应用其性质实现信号的幅度调制。7. 理解连续

13、时间系统的频域分析原理和方法，掌握连续系统的频域响应求解方法，并画出相应的幅频和相频响应曲线。二、 实验原理三、周期信号的傅里叶分解三角函数形式 (1)式中傅里叶系数，指数函数形式 (2)式中傅里叶系数 (3)MATLAB实现：数值积分quad()函数 quadl()函数 y=quad(func,a,b)/func为字符串，表示被积函数.m文件的文件名y=quad(myfun,a,b)/表示函数句柄。符号函数int()函数四、周期信号的频谱 幅度谱和相位谱傅里叶系数直接获得。4-1给定周期为4，脉冲宽度为2，幅值为0.5的矩形信号，计算傅里叶级数，绘制幅度谱和相位谱。将其带入，求N=10实的合

14、成波形。clcclearT=4;width=2;A=0.5;t1=-T/2:0.001:T/2;ft1=0.5*abs(t1)width/2;t2=t1-2*T t1-T t1 t1+T t1+2*T;ft=repmat(ft1,1,5);subplot(4,1,1);plot(t2,ft);xlabel(t);ylabel(时域波形)grid on;w0=2*pi/T;N=10;K=-N:N;for k=-N:Nfactor=exp(-j*t*, num2str(w0),*,num2str(k),);f_t=num2str(A),*rectpuls(t,2) ;Fn(k+1+N)=quad(

15、f_t,.*,factor,-T/2,T/2)/T;endsubplot(4,1,2)stem(K*w0,abs(Fn);xlabel(nw0);ylabel(幅度谱);grid on;ph=angle(Fn);subplot(4,1,3)stem(K*w0,ph);xlabel(nw0);ylabel(相位谱);grid on;t=-2*T:0.01:2*T;K=-N:N;Ft=Fn*exp(j*w0*K*t);subplot(4,1,4)plot(t,Ft)xlabel(t);ylabel(合成波形)grid on;讲解1. B = repmat(A,M,N) creates a larg

16、e matrix B consisting of an M-by-N tiling of copies of A. The size of B is size(A,1)*M, size(A,2)*N.The statement repmat(A,N) creates an N-by-N tiling. 复制矩阵函数2. T = Num2str(X) converts the matrix X into a string representation T with about 4 digits and an exponent if required. This is useful forlabe

17、ling plots with the TITLE, XLABEL, YLABEL, and TEXT commands. 数值到字符转换。五、非周期信号的傅里叶变换和性质傅里叶变换 频谱密度函数 (4)傅里叶逆变换 (5)Matlab实现：符号函数syms 定义的变量fourier()函数, ifourier()函数。4-2求单边函数的傅里叶变换，画出其幅频特性和相频特性图。clcclearsyms t w f;f=exp(-2*t)*sym(Heaviside(t);F=fourier(f);subplot(2,1,1);ezplot(f,0,2,0,1.2);subplot(2,1,2)

18、;ezplot(abs(F),-10:10);讲解：1.ezplot(FUN,A,B) plots FUN(X) over A X B.2. ezplot(FUN2,XMIN,XMAX,YMIN,YMAX) plots FUN2(X,Y) = 0 overXMIN X XMAX and YMIN Y YMAX.3. Heaviside Step function. Heaviside (X) is 0 for X 0, and NaN for X = 0.4-3求的傅里叶逆变换。clcclearsyms t wF=1/(1+w2);f=ifourier(F,w,t);ezplot(f);讲解：

19、1. ezplot(FUN) plots the function FUN(X) over the default domain -2*PI X 2*PI.2. ezplot()函数无法画出函数。调制 (6) 7)4. 连续系统的频域分析和频域响应 (8)， (9)MATLAB实现：freqs()函数H= freqs(B,A,W) returns the complex frequency response vector H of the filter B/A: nb-1 nb-2 B(s) b(1)s + b(2)s + . + b(nb) H(s) = - = - na-1 na-2 A(

20、s) a(1)s + a(2)s + . + a(na)given the numerator and denominator coefficients in vectors B and A. The frequency response is evaluated at the points specified in vector W (in rad/s). The magnitude and phase can be graphed by calling FREQS(B,A,W) with no output arguments. H,W = FREQS(B,A) automatically

21、 picks a set of 200 frequencies W on which the frequency response is computed. FREQS(B,A,N) picks N frequencies.三、实验内容1. 求下列信号的傅里叶变换表达式并画图。代码如下：clcclearsyms t w f;f=exp(-2*abs(t);F=fourier(f);subplot(2,1,1);ezplot(f,0,2);subplot(2,1,2);ezplot(abs(F),-10:10)图形如下：2.求的傅里叶逆变换式并且画图。代码如下：clcclearsyms t w

22、F=2*sin(w)/w;f=ifourier(F,w,t);ezplot(f);图形如下：4.画出N=30时，例4-1的傅里叶级数，绘制幅度谱和相位谱。将其带入，求N=30实的合成波形。代码如下：clcclearT=4;width=2;A=0.5;t1=-T/2:0.001:T/2;ft1=0.5*abs(t1)width/2;t2=t1-2*T t1-T t1 t1+T t1+2*T;ft=repmat(ft1,1,5);subplot(4,1,1);plot(t2,ft);xlabel(t);ylabel(时域波形)grid on;w0=2*pi/T;N=30;K=-N:N;for k=

23、-N:Nfactor=exp(-j*t*, num2str(w0),*,num2str(k),);f_t=num2str(A),*rectpuls(t,2) ;Fn(k+1+N)=quad(f_t,.*,factor,-T/2,T/2)/T;endsubplot(4,1,2)stem(K*w0,abs(Fn);xlabel(nw0);ylabel(幅度谱);grid on;ph=angle(Fn);subplot(4,1,3)stem(K*w0,ph);xlabel(nw0);ylabel(相位谱);grid on;t=-2*T:0.01:2*T;K=-N:N;Ft=Fn*exp(j*w0*K

24、*t);subplot(4,1,4)plot(t,Ft)xlabel(t);ylabel(合成波形)grid on;图形如下：四、实验思考通过实验我熟悉了各个函数功能并掌握各种变换有助于进一步了解信号的各种变化的方式和意义。理解周期信号的傅里叶分解，掌握傅里叶系数的计算方法。深刻理解和掌握了非周期信号的傅里叶变换和其计算方法。熟悉了傅里叶变换的性质，并能应用其性质实现信号的幅度调制。理解了连续时间系统的频域分析原理和方法，掌握了连续系统的频域响应求解方法。实验四 连续时间系统的复频域分析一、 实验目的1.深刻理解和掌握拉普拉斯变换的运算方法和其性质。2.熟练掌握利用部分分式展开的方法求解拉普拉

25、斯逆变换，并能用MATLAB实现。3.掌握复频域系统函数的意义，并能够熟练画出其频谱。4.掌握利用复频域系统函数的零、极点分布对连续时间系统进行复频域分析的原理和方法。二、 实验原理1. 拉普拉斯变换 (1)MATLAB实现：F=laplace(f) f=ilaplace(F)F和f都是符号函数。5-1 用laplace和ilaplace求：(1) 的拉普拉斯变换 (2) 的拉普拉斯逆变换(1)clcclearsyms a t;F=laplace(exp(-2*t)*cos(a*t)F=(s+2)/(s2+4*s+4+a2)(2)clc clearsyms sF=1/(s+1)*(s+2);f

26、=ilaplace(F);f=exp(-t)-exp(-2*t)2. 部分分式展开法求拉普拉斯逆变换 (2)特征多项式，特征方程，极点的概念。MATLAB实现： r p k=residue(num,den) r系数向量p极点k常数项B(s) R(1) R(2) R(n) - = - + - + . + - + K(s) A(s) s - P(1) s - P(2) s- P(n)。重根情况R(j) R(j+1) R(j+m-1) - + - + . + - s - P(j) (s - P(j)2 (s - P(j)m3. 系统传递函数和频域响应函数 (3) (4)4. 系统传递函数零极点与系统

27、稳定性稳定极点在左边平面。MATLAB实现： roots()函数 p=roots(a)，p为多项式系数向量。pzmap函数 pzmap(sys) sys系统模型 syst=tf(b,a) b, a分别是H(s)分子分母多项式系数函数。5-5求传递函数的零极点分布图。clcclearnum=1 -1;den=1 2 2;zs=roots(num);ps=roots(den);figure(1);plot(real(zs),imag(zs),o,real(ps),imag(ps),kx, markersize,12);axis(-2 2 -2 2);grid onsys=tf(num,den);figure(2)pzmap(sys);axis(-2 2 -2 2);讲解1. 函数real()和函数imag()计算一个复数实部和虚部。2. figure画另一个图形3. 再分析plot函数 x = -pi:pi/10:pi; y = tan(sin(x) - sin(tan(x); plot(x,y,-rs,LineWidth,2,. MarkerEdgeColor,k,. MarkerFaceColor,g,. MarkerSize,10