信号与系统实验

《信号与系统实验》由会员分享，可在线阅读，更多相关《信号与系统实验（5页珍藏版）》请在装配图网上搜索。

1、实验六 利用MATLAB实现周期信号的分解与合成一、实验目的1复习巩固信号与系统课程中关于周期信号展开成傅里叶级数等有关理论知识； 2验证、设计并掌握利用MATLAB实现周期信号的分解与合成的方法。二、实验内容1问题一周期方波信号如图所示，求出该信号的三角函数形式的傅里叶级数，并用MATLAB编程实现其各次谐波叠加情况的观察与分析。三角函数形式的傅里叶级数展开：该信号的周期T=2，基频为因为该方波信号为奇函数，所以 而 所以 程序代码：%观察周期方波信号的分解与合成%m：傅里叶级数展开的项数display(Please input the value of m（傅里叶级数展开的项数）); %在

2、命令窗口中显示提示m=input(m=); %键盘输入傅里叶级数展开的项数t=-2*pi:0.01:2*pi %时域波形的时间范围-22，采样间隔0.01n=round(length(t)/4); %根据周期方波信号的周期，计算1/2周期的数据点数f=ones(n,1);-1*ones(n,1);ones(n,1);-1*ones(n+1,1); %构造周期方波信号y=zeros(m+1,max(size(t);y(m+1,:)=f;figure(1);plot(t/pi,y(m+1,:); %绘制方波信号grid on;axis(-2 2 -1.5 1.5);title(周期方波);xlab

3、el( 单位： pi,Fontsize,8);x=zeros(size(t);kk=1;for k=1:2:2*m-1 %循环显示谐波叠加图形 pause; x=x+sin(k*t)/k; y(k+1)/2,:)=4/pi*x; %计算各次谐波叠加和 plot(t/pi,y(m+1,:); hold on; plot(t/pi,y(k+1)/2,:); %绘制谐波叠加信号 hold off; grid; axis(-2 2 -1.5 1.5); title(strcat(第,kk,次谐波叠加); xlabel(单位： pi,Fontsize,8); kk=strcat(kk,、,num2str

4、(k+2);endpause;plot(t/pi,y(1:m+1,:)grid;axis(-2 2 -1.5 1.5);title(各次谐波叠加波形);xlabel( 单位： pi,Fontsize,8);%end执行结果：观察实验结果，理解：“当所包含的谐波分量越多时，合成波形越接近于原来的周期信号”。结合实验过程，理解上述程序中的命令和相应显示结果的对应关系。2问题二求下图所示矩形脉冲信号f(t)的傅里叶级数表示式，并用MATLAB方法求出N=7和N=21时的合成图。傅里叶级数展开：信号f(t)的离散频谱函数（复傅里叶系数） （参考教材97页式48） 得傅里叶级数表达式 即f(t)的前N项

5、合成表达式为 t=-2:0.001:2;N=input(N=);F0=0.5;fN=F0*ones(1,length(t);for n=1:2:N fN=fN+cos(pi*n*t)*sinc(n/2);endplot(t,fN);title(N= num2str(N)axis(-2 2 -0.2 1.2);程序代码：执行结果：3设计信号f(t)如图所示，设T=2s、A=1，试求其复指数形式的傅里叶级数，设计一个MATLAB程序，通过运行程序观察此信号的分解与合成。三、作业 1进一步理解周期信号分解为傅里叶级数的理论；2完成问题二中f(t)的前N项合成表达式的详细推导。附加说明 关于dsolv

6、e函数调用格式 yp=dsolve(方程1，方程2，条件1，条件2，变量)简要说明（1）函数返回微分方程1（或方程1、方程2、构成的微分方程组）的解；（2）各个微分方程中的一阶导数符号用D表示，D2、D3、Dn分别表示二阶、三阶、n阶导数；（3）各个微分方程中的自变量默认为t，否则，可用格式中的 变量 重新指明；（4）各个条件项可用于指明初始条件。若不指明或指明不全，则函数的返回值中可能出现C1、C2、等待定常数。 关于tf函数 tf函数根据描述系统的微分方程的系数生成系统函数对象（TF对象）。调用格式 sys=tf(b,a)简要说明输入参量a 微分方程左边多项式系数构成的行向量；输入参量b

7、微分方程右边多项式系数构成的行向量；输出参量sys 系统函数对象对一般形式的系统微分方程: 或a和b应分别定义成：、。a=1 3 2;b=-1 2;sys=tf(b,a) %调用tf函数生成系统函数对象sys例如，对由微分方程描述的系统，由tf函数生成其系统对象sys的命令为Transfer function:-s+2-s2+3s+2上述命令运行结果为 调用tf函数生成系统函数对象sys,并用向量f和t定义了系统激励信号后，即可调用lsim函数对连续系统的响应进行仿真。 关于lsim函数调用格式 lsim函数的一种调用格式为： lsim(sys,f,t)简要说明 这种调用格式对向量t定义的时间

8、范围内的系统响应进行仿真，即绘制LTI连续系统的响应的时域波形，同时还绘出系统的激励信号对应的时域波形。 在这种调用格式中，输入参量f和t是两个表示输入信号的行向量，其中t为表示输入信号时间范围的向量，f则是输入信号在向量t定义的时间点上的采样值。 例如，命令 t=0:0.01:5; f=cos(2*t); 就定义了05秒时间范围内的余弦输入信号cos(2*t)（采样时间间隔为0.01秒）。 输入参量sys是由tf函数根据描述系统的微分方程的系数生成的系统函数对象。 关于impulse函数 Impulse函数主要用于求连续系统冲激响应并绘制其时域波形。调用格式 （1）impulse(b,a)

9、（2）impulse(b,a,t) （3）impulse(b,a,t1:ts:t2) （4）y=impulse(b,a,t1:ts:t2)简要说明 （1）该调用格式以默认的方式绘制由向量a和b定义的LTI连续系统的冲激响应时域波形，并提供交互式功能查看任意时刻冲激响应的信号样值。其中输入参量a和b分别是描述系统的微分方程左边和右边系数构成的行向量。 例如，对由微分方程描述的系统，则定义该系统的向量a和b应使用命令： a=1 0 2 4; b=-1 6; 。 （2）该调用格式绘制由向量a和b定义的连续系统在时间参数t所指定的0t时间范围内系统冲激响应的时域波形。（3）该调用格式绘制由向量a和b定义的连续系统在时间范围t1t2内，且以时间间隔ts均匀采样的系统冲激响应的时域波形。（4）该调用格式并不绘制系统冲激响应的波形，而是求出由向量a和b定义的连续系统在时间范围t1t2内，且以时间间隔ts均匀采样的系统冲激响应的数值解。 关于step函数调用格式 （1）step(b,a) （2）step (b,a,t) （3）step (b,a,t1:ts:t2) （4）y= step (b,a,t1:ts:t2) 上述4种调用格式的使用方法和调用过程与impulse函数完全相同，只是该函数绘制的是系统阶跃响应的时域波形而不是冲激响应波形。