ctfssjbshbd.m

我做了很久的

function [A_sym,B_sym]=CTFSsjbshbd
syms t n  k y 
T=5;
if nargin<4;Nf=input('please input 所需展开的最高谐波的次数:');end
if nargin<5;Nn=32;end
y=time_fun_s(t);
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(3,3,3)
x=sjb_timefun
T=5;t=-2*T:0.01:2*T;
plot(t,x)
title('连续时间函数周期三角波脉冲')
axis([-10,10,-1,1.2])
line([-10,10],[0,0])
subplot(3,1,3)
stem(k2,abs(S7))
     title('续时间函数周期三角波脉冲的双边幅度谱')
     axis([-80,80,0,0.25])
 end