内容发布更新时间 : 2024/11/17 4:24:50星期一 下面是文章的全部内容请认真阅读。
plot(t,y,'r:') title('基波+2 次谐波+...+5次谐波+6次谐波 ') axis([-4,4.5,-0.5,1.5]) subplot(2,2,4) plot(t,f33), hold on y=time_fun_e(t) plot(t,y,'r:') title(' 基波+2 次谐波+...+9次谐波+10次谐波 ') axis([-4,4.5,-0.5,1.5]) end function x=time_fun_x(t) % 该函数是 CTFShchsym.m的子函数。它由符号变量和表达式写成。 h=1; x1=sym('heaviside(t+0.5)')*h; x=x1-sym('heaviside(t-0.5)')*h; function y=time_fun_e(t) % 该函数是 CTFShchsym.m的子函它由符号函数和表达式写成 a=0.5; T=5; h=1; tao=0.2*T; t=-8*a:0.01:T-a; e1=1/2+1/2.*sign(t+tao/2); e2=1/2+1/2.*sign(t-tao/2); y=h.*(e1-e2); %连续时间函数-周期矩形脉冲 2、 实验程序： %testFynction2 WaveNum=10; Nn=5; a=2; T=4; h=1; t0=-2:0.01:4-a; x=sym('abs(2-ts)-1'); A0=double(int(x,ts,0,2*a))*h/T; An=int(2*x*cos(2*pi*n*ts/T)/T,ts,0,2*a); Bn=int(2*x*sin(2*pi*n*ts/T)/T,ts,0,2*a); clear A_List B_List; for i=1:WaveNum 实验结果： A_List(i)=double(vpa(subs(An,n,i),Nn))*h; B_List(i)=double(vpa(subs(Bn,n,i),Nn))*h; end t=mod(t0,T); f0=zeros(1,length(t)); f=zeros(WaveNum,length(t)); for i=1:WaveNum f(i,:)=A_List(i).*cos(2*pi*i*t/T)+B_List(i).*sin(2*pi*i*t/T); if i==1 f(i,:)=f(i,:)+A0; end f0=f0+f(i,:); end y=testFunction2(t); figure(2)