DSP作业讲解.pptVIP

* * 连续时间周期信号的matlab实现 其中 其中 幅频： 相频： 傅里叶级数三角函数形式： 傅里叶级数指数函数形式： 二者之间的关系： 由此可知，只要我们求出了傅里叶级数的系数an、bn ，进而我们就可以求出信号的幅度谱和相位谱。 由于matlab中没有现成的求傅里叶级数的函数，需要自己编写，这里我们用到了matlab中的符号工具箱。另外，由于周期信号的频谱是离散的，故在绘制频谱图时，我们采用的是stem命令，而不是plot命令。 连续周期信号matlab实现流程： 1.用符号表达式表示出周期信号在第一个周期内的符号表达式，并赋值给返回符号变量y； 2.求出信号的三角级数形式的傅里叶级数展开系数an、bn ； 3.求出信号的复指数形式的傅里叶级数展开系数Fn 4.绘制信号的频谱图 连续时间周期信号的matlab实现 参考资料：《信号与线性系统分析：基于matlab的方法与实现》 梁虹等编著 例题：已知周期信号f(t)如下，画出其频谱图。 其中 解：将f(t)整理为标准形式 function [A,B,F]=fouriers(f,t,T,a,b,k) w=2*pi/T; A=1/T*int(f,t,a,b);%计算A0 B=[]; F=A; if k==0 syms k integer ak=2/T*int(f*cos(k*w*t),t,a,b); bk=2/T*int(f*sin(k*w*t),t,a,b); A=[A,ak]; B=[B,bk]; F=[]; else for i=1:k ak=2/T*int(f*cos(k*w*t),t,a,b); bk=2/T*int(f*sin(k*w*t),t,a,b); A=[A,ak]; B=[B,bk]; F=F+ak*cos(i*w*t)+bk*sin(i*w*t); end end end 根据前面的公式，可以自己编写求解傅里叶级数系数的函数，然后调用绘制频谱图。 %本函数实现周期方波信号的频谱图绘制，采用傅里叶级数方法 function [A_sym,B_sym]=fangbo syms t n k y; T=10; if nargin4;Nf=input(‘pleas Input 所需展开的最高谐波系数Nf= );end if nargin5;Nn=32;end y=fangbo_1; A0=2*int(y,t,0,T)/T; 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)); 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 if nargout==0 S1=fliplr(A_sym); S1(1,k+1)=A_sym(1); S2=fliplr(1/2*S1); S3=fliplr(1/2*B_sym); S3(1,k+1)=0; S4=fliplr(S3); S5=S2-i*S4; S6=fliplr(S5); N=Nf*2*pi/T; k2=-N:2*pi/T:N; S7=[S6,S5(2:end)]; subplot(2,1,1) x=fangbo_2; plot(t,x) title(‘T=5,占空比为50%的周期方波信号) axis([-10,10,-1.2,1.2]) subplot(2,1,2),stem(k2,abs(S7)); title(‘连续时间函数周期方波信号的双边幅度谱) axis([-20,20,0,0.6]) figure,stem(angle(S7)); end end function y=fangbo_1 syms a a1 T=5;a=T/2; y1=sym(Heaviside(t))*2-sym(Heaviside(t-a1)); y=y1-sym(Heaviside(t+a1)); y=