周期信号频域解析及MATLAB实现-全文可读.pptxVIP

;;7.1 周期信号的傅里叶级数与信号的频谱;7.1 周期信号的傅里叶级数与信号的频谱;7.1 周期信号的傅里叶级数与信号的频谱;7.1 周期信号的傅里叶级数与信号的频谱;7.1 周期信号的傅里叶级数与信号的频谱;7.2 Matlab分析周期信号的频谱;7.2.2 单边频谱分析;例题7.1 周期矩形脉冲信号幅度为1，脉宽为1，周期为4，求该信号的三角形式的其各阶傅里叶系数，并绘出各次谐波叠加的傅里叶综合波形图。 1.编写函数文件 function [y,tao]=SqureWaves(t,T) %该函数用于产生对称的、脉宽可调的矩形脉冲； %它由符号变量和表达式写成。 % M——为脉冲占空比； % tao——为脉冲宽度； % a——为脉冲半宽度； % y——周期函数在一个周期的符号函数表达式。 M=input(duty cycle脉冲占空比M=); A=1;tao=T*M; y=A*heaviside(t+tao/2)-A*heaviside(t-tao/2); y=simple(y);;绘制函数图 figure(3) Syms t; T=4; [y,tao]=SqureWaves(t,T); yn=subs(y,t,t-T)+y+subs(y,t,t+T); %定义有原函数3个周期的函数 ezplot(yn,[-1.5*T 1.5*T]); 调用m文件SquarePulse_CTFS. m求傅里叶;A0=0.25;;调用SquarePulse计算傅里叶系数;原函数的各阶截断傅里叶级数;7.2.3 双边频谱 周期信号可以分解成一系列虚指数信号之和，并可以求得相应的傅里叶系数;思路：利用CTFShchsym.m先求出周期信号的傅里叶系数，再利用系数关系求An，Fn，φn。;绘制周期函数;编写函数文件 子函数y=time_fun_s(t)，表示出信号在一个周期内的符号表达式 傅里叶复指数系数及绘制频谱图 调用函数time_fun_s(t)，获取周期信号的符号表达式; 求出信号的三角级数形式的傅里叶系数an和bn 。 求出信号的复指数形式的傅里叶级数展开系数Fn 。 绘制Fn的振幅频谱图。 利用y的符号表达式绘制信号3个周期的波形图;7.2.4 用FFT实现周期信号的频谱分析;2. 用FFT 法对x(n)傅里叶分析的计算公式;三、用FFT 实现周期信号的频谱分析;例7.2 用fft分析淹没在噪声中的50Hz和150H在信号;7.2 Matlab分析周期信号的频谱;例7.3 用FFT 周期矩形脉冲函数的频谱， 幅度为 1，脉冲宽度为tao = 1 function x=fun_ex0706(t,T) % 该函数是CTFSfft.m的子函数。 %t是时间数组 %T是周期 A=1;tao=1;T=5; x=zeros(size(t)); ii=find(t=-tao /2 t=tao/2); x(ii)=ones(size(ii)).*A;x(t==0)=1; return 调用CTFSfft;7.2 Matlab分析周期信号的频谱;7.3 用MATLAB分析典型周期信号的频谱;例题7.4 周期方波的理论推算;调用Ex07_04.m;7.3.2周期三角波脉冲的频谱 例题7.5 周期方波的频谱分析；T=5 s，幅度=±1 function y=fun_ex0705(t,T) syms a a1 % T=input(please Input 信号的周期T=); 调用Ex07_05.m global T M=input(duty cycle脉冲占空比M=); A=1;tao=T*M;a=tao/2; y1=sym(Heaviside(t+a1))*A; y=y1-sym(Heaviside(t-a1))*A; y=subs(y,a1,a); y=simple(y); return;7.3 用MATLAB分析典型周期信号的频谱