连续周期信号的分解与综合matlab,连续时间信号傅里叶级数分析及matlab实现

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

T=5;t=-2*T:0.01:2*T;

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])

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=subs(y,a1,a);

y=simple(y);

function x=fangbo_2

T=5;t=-2*T:0.01:2*T;duty=50;

x=square(t,duty);

源程序五：(典型周期脉冲的频谱分析——三角波)function [A_sym,B_sym]=sjb

syms t n k y