第二部分 周期信号傅里叶级数MATLAB实现.docVIP

第二部分 周期信号的傅里叶级数的MATLAB实现 【实验目的】 了解傅里叶变换的基本原理； 掌握周期信号傅里叶级数的编程思想。 【原理】 设有周期信号，它的周期是T，角频率，它可以分解为 式中： 傅里叶系数 【示例】 已知周期方波脉冲信号如图所示，其幅度为1，脉冲宽度（占空比）duty=1/2，周期T=5。试用MATLAB编程绘出该周期信号的频谱。(红色部分为说明) 解： function [A_sym,B_sym]=CTFsdbfb 采用符号计算求[0,T]内时间函数的三角级数展开系数。 函数的输入输出都是数值量 Nf 谐波的阶数 Nn 输出数据的准确位数 A_sym 第1元素是直流项，其后元素依次是1,2,3...次谐波cos项展开系数 B_sym 第2,3,4,...元素依次是1,2,3...次谐波sin项展开系数 syms t n k y 创建符号函数，相当于y=f(x)中的x T=5; 周期T=5 Nf=60; 谐波的阶数 Nn=32; 输出数据的准确位数，就是有效位数 y=time_fun_s; 调用time_fun_s函数，和C语言相当 A0=2*int(y,t,0,T)/T; 计算，int为积分函数，格式为int（需要积分的函数，积分变量，下线，上线） As=int(2*y*cos(2*pi*n*t/T)/T,t,0,T); 计算 Bs=int(2*y*sin(2*pi*n*t/T)/T,t,0,T); 计算 A_sym(1)=double(vpa(A0,Nn)); vpa:计算变量的精度，就是有效位数 double：强制转换为双精度浮点型 for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn)); B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn)); end for循环，计算，的精度值 subs是将符号函数转换为数值计算，调用格式subs(转换函数，转换的变量，将要计算变量) k=Nf; S1=fliplr(A_sym); 对A_sym阵左右对称交换，fliplr为左右对称交换函数，相当于y=[1 2 3 4],fliplr(y)=[-4 -3 -2 -1] S1(1,k+1)=A_sym(1); 上面两句语句相当于将A_sym的1*k阵扩展为1*(k+1)阵 S2=fliplr(1/2*S1); S3=fliplr(1/2*B_sym); 乘以1/2为公式需要 S3(1,k+1)=0; 这句是为了让S3和S2的向量中的个数相等，使后面相减时不出错 S4=fliplr(S3); S5=S2-i*S4; N=Nf*2*pi/T; 计算 k2=0:2*pi/T:N; 计算横坐标的值 subplot(2,1,1) T=5;t=-2*T:0.01:2*T; x=squ_timefun(t); 调用squ_timefun函数，计算周期方波脉冲 plot(t,x); title(周期方波脉冲); axis([-10,10,-1,1.2]); line([-10,10],[0,0]); 相当于在y=0时画条辅助线 subplot(2,1,2) stem(k2,abs(S5)); title(周期方波脉冲的单边频谱); axis([0,60,0,0.6]); %--------------------------------------------------------- function y=time_fun_s(t) syms a a1 T=5;a=T/2; y1=heaviside(t)*2-heaviside(t-a1); y=y1-heaviside(t+a1); heaviside为MATLAB自带的单位跃阶信号函数，这条语句是用单位跃阶信号形成单个方波脉冲信号 例： t=-5:0.001:5; T=5;a1=T/2; y=heaviside(t+a1).*((2*t+a1)/a1);%相当于直线在a1点上乘以它的斜率tan(60°) plot(t,y);axis([-3,3,-2,2]); ps:与作业有关，大家好好好想想，不会可以来提问 如图： y=subs(y,a1,a); y=simple(y); 简化符号函数y,是成为多项式的表示形式，可以help simple试试 %---------------------------------------------------------- function x=squ_timefun(t,T) duty=50; x=square(t,duty); square为MATLA