基于matlab的連續(xù)時(shí)間信號(hào)傅里葉級(jí)數(shù)分析及實(shí)現(xiàn)論文-資料下載頁(yè)
2024-11-07 21:53本頁(yè)面
【正文】 Ta。 e1=1/2+1/2.*sign(t+tao/2)。 e2=1/2+1/2.*sign(ttao/2)。 y=h.*(e1e2)。 function x=time_fun_x(t) %定義矩形脈沖函數(shù)（積分用） h=1。 x1=sym(39。Heaviside(t+)39。)*h。 x=x1sym(39。Heaviside()39。)*h。 源程序二：（吉布斯效應(yīng)） t=2::2。 N=input(39。N=39。)。 % N為輸入要達(dá)到的最高次諧波的次數(shù) c0=。 fN=c0*ones(1,length(t))。 for n=1:2:N fN=fN+cos(pi*n*t)*sinc(n/2)。 end figure plot(t,fN) axis([2 2 ]) 源程序三： （矩形脈沖的頻譜分析） function [A_sym,B_sym]=CTFS2 syms t n y if nargin3。Nf=input(39。pleas Input 所需展開(kāi)的最高諧波次數(shù) : Nf=39。)。end T=input(39。pleas Input 信號(hào)的周期 T=39。)。 if nargin5。Nn=32。end y=fun_in(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) 24 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=S2i*S4。 S6=fliplr(S5)。 N=Nf*2*pi/T。 k2=N:2*pi/T:N。 S7=[S6,S5(2:end)]。 x=fun_mc subplot(3,1,2) stem(k2,abs(S7))。 axis([150,150,0,]) subplot(3,1,3) stem(k2,angle(S7))。 axis([150,150,4,4]) end function y=fun_in(t) syms a a1 T=input(39。pleas Input 信號(hào)的周期 T=39。)。 M=input(39。周期與脈沖寬度之比 M=39。)。 A=1。tao=T/M。a=tao/2。 y1=sym(39。Heaviside(t+a1)39。)*A。 y=y1sym(39。Heaviside(ta1)39。)*A。 y=subs(y,a1,a)。 y=simple(y)。 function x=fun_mc T=5。tao=T/5。n=4。 t=n*T::n*T。 x=rectpuls(t,1)。 for i=1:n。 x=x+rectpuls(ti*T,1)+rectpuls(t+i*T,1)。 end subplot(3,1,1) plot(t,x) hold on axis([20,20,0,]) 源程序四： （典型周期脈沖的頻譜分析 —— 方波） function [A_sym,B_sym]=fangbo syms t n k y T=10。 if nargin4。Nf=input(39。pleas Input 所需展開(kāi)的最高諧波次數(shù) Nf= 39。)。end if nargin5。Nn=32。end 25 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=S2i*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::2*T。 plot(t,x) title(39。T=5，占空比為 50%的周期方波脈沖 39。) axis([10,10,]) subplot(2,1,2) stem(k2,abs(S7))。 title(39。連續(xù)時(shí)間函數(shù)周期方波脈沖的雙邊幅度譜 39。) axis([20,20,0,]) end function y=fangbo_1 syms a a1 T=5。a=T/2。 y1=sym(39。Heaviside(t)39。)*2sym(39。Heaviside(ta1)39。)。 y=y1sym(39。Heaviside(t+a1)39。)。 y=subs(y,a1,a)。 y=simple(y)。 function x=fangbo_2 T=5。t=2*T::2*T。duty=50。 x=square(t,duty)。 源程序五： （典型周期脈沖的頻譜分析 —— 三角波） function [A_sym,B_sym]=sjb syms t n k y 26 T=10。 if nargin4。Nf=input(39。pleas Input 所需展開(kāi)的最高諧波次數(shù) Nf= 39。)。end if nargin5。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=S2i*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=sjb_timefun T=10。t=2*T::2*T。 plot(t,x) title(39。T=10，脈寬為 1的周期三角波脈沖 39。) axis([20,20,1,]) subplot(2,1,2) stem(k2,abs(S7))。 title(39。連續(xù)時(shí)間函數(shù)周期三角脈沖的雙邊幅度譜 39。) axis([0,20,0,]) end function x=sjb_timefun T=10。t=2*T::2*T。 x=sawtooth(t2*T/3,1)。 function y=time_fun_s(t) syms a a1 T=10。a=T/2。 y=sym(39。Heaviside(t+a1)39。)*(2*t/a1)sym(39。Heaviside(ta1)39。)*(2*t/a1)。 %y=y1sym(39。Heaviside(t)39。)*(4*t/a1)。 y=subs(y,a1,a)。 y=simple(y)。
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1