信号与系统实验3周期信号频域分析.docVIP

实验三 周期信号频域分析 一、目的 （1）掌握周期信号傅立叶级数分解与合成的计算公式 （2）掌握利用MATLAB实现周期信号傅立叶级数分解与综合方法 （3）理解并掌握周期信号频谱特点 二、周期信号傅立叶级数 周期信号是定义在区间内，按一定时间间隔（周期）不断重复的信号。可表示为 式中为任意整数，为周期，周期的倒数成为该信号频率。 实验内容 仿照例程，实现下述周期信号的傅立叶级数分解与合成： 要求： （a）首先，推导出求解，，的公式，计算出前10次系数； （b）利用MATLAB求解，，的值，其中，求解前10次系数，并给出利用这些系数合成的信号波形。 （a）设周期信号的周期为，角频率，且满足狄里赫利条件，则该周期信号可以展开成傅立叶级数。 （1）三角形式傅立叶级数 （2）指数形式傅立叶级数 （b）求解，，及合成信号波形所用程序： function [A_sym,B_sym]=CTFShchsym % 采用符号计算求一个周期内连续时间函数f的三角级数展开系数,再用这些 % 展开系数合成连续时间函数f.傅立叶级数 % 函数的输入输出都是数值量 % Nf=6 谐波的阶数 % Nn 输出数据的准确位数 % A_sym 第1元素是直流项，其后元素依次是1,2,3...次谐波cos项展开系数 % B_sym 第2,3,4,...元素依次是1,2,3...次谐波sin项展开系数 % tao=1 tao/T=0.2 syms t n k x T=4; tao=T/4; a=-1.5; if nargin4 Nf=10; end if nargin5 Nn=32; end x=time_fun_x(t); A0=int(x,t,a,T+a)/T; %求出三角函数展开系数A0 As=2/T*int(x*cos(2*pi*n*t/T),t,a,T+a); %求出三角函数展开系数As Bs=2/T*int(x*sin(2*pi*n*t/T),t,a,T+a); %求出三角函数展开系数Bs A_sym(1)=double(vpa(A0,Nn)); %获取串数组A0所对应的ASC2码数值数组 for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn)); %获取串数组A所对应的ASC2码数值数组 B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn)); %获取串数组B所对应的ASC2码数值数组 end ; if nargout==0 c=A_sym; disp(c); %输出c为三角级数展开系数:第1元素是直流项，其后元素依次是1,2,3...次谐波cos项展开系数 d=B_sym; disp(d); %输出d为三角级数展开系数: 第2,3,4,...元素依次是1,2,3...次谐波sin项展开系数 t=-3*T:0.01:3*T; f0=c(1); %直流 f1=c(2).*cos(2*pi*1*t/T)+d(2).*sin(2*pi*1*t/T); % 基波 f2=c(3).*cos(2*pi*2*t/T)+d(3).*sin(2*pi*2*t/T); % 2次谐波 f3=c(4).*cos(2*pi*3*t/T)+d(4).*sin(2*pi*3*t/T); % 3次谐波 f4=c(5).*cos(2*pi*4*t/T)+d(5).*sin(2*pi*4*t/T); % 4次谐波 f5=c(6).*cos(2*pi*5*t/T)+d(6).*sin(2*pi*5*t/T); % 5次谐波 f6=c(7).*cos(2*pi*6*t/T)+d(7).*sin(2*pi*6*t/T); %