周期信号傅里叶分析

《周期信号傅里叶分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验三周期信号傅里叶分析1、将fl(t)按照最高谐波次数为7,21,41进行分解，然后叠加恢复原信号fl(t)的波形变化。n=7:symstnf=0;forn=l:2:7;f二f+(4/pi)*(l/n)・*sin(n*2*pi*t);endezplot(t,f,[0,10]);n=21:symstnf=0;forn=l:2:21;f二f+(4/pi)*(l/n).*sin(n*2*pi*t);enclezplot(t,f,[0,10]);n=41:symstnf二0;forn=l:2:41;f二f+(4/pi)*(l/n)・*sin(n*2*pi*t);e

2、ndezplot(t,f,[0,10]);2、周期信号f2(t)的单边频谱图:closeal1;clear;symsnN二25;fw=zcros(1,N);w二pi/2;forn二1:1:N;holdon;f=(j*(l/w))*(l/n).*(exp(-j*n*w)T);fw(l,n)=abs(f);endplot(fw(l,:));bar(fw(l,:),0.01);将前25项合成原信号f2(t)的波形：symstnf=2;forn=l:1:25;f二f+(2*j/pi)*(1/n)•*(exp(-j*n*(pi/2))T).*(exp(j*n*(pi/

3、2)*t));endezplot(t,f,[0,10]);closeallclear八symsnNgfwHzerospN)八WHPiaforrlHl二Nholdon-f"l、8Ysinc((n*pi)、8)-fwprlTabs(f)八endplot(fw(lmbar(fw(l二)o・01)-(1)Fnu(t、16)关Sa(n*pi*t、16)(T=16);4:closeallclearsymsnNgfwuzcrospN)八WHPi、?0forrlHl二Nholdon-f"l、4Ysinc((n*pi)、4)-fwprlTabs(f)八endploc(fw

4、(l二))八bar(fw(l二)o・01)八T=1：closeall;clear;symsnN二25;fw=zcros(1,N);w=pi/2;forn=l:1:N;holdon;f=(1/16)*sinc((n*pi)/16);fw(l,n)=abs(f);endplot(fw(l,:));bar(fw(1,:),0.01);closeallclear八symsnNgfwHzerospN)八WH2关pi、5-fornul二Nholdon.f"l、(2^5)¥sinc((I^pi)、(2^5))-fwpnvabs(f)-endploc+(fw(l二))八b

5、ar(fw(l二)0.01)-(2)Fnu(l、(2*T))*sa(n*e、4)(；0・5)TH2“closeall.clear-symsnNgfwnzcrosPN)八WHPiifornul二Nholdon八f"l、2Ysinc((n*pi)、4)-fwprlTabs(f)八endploc(fw(l二))八bar(fw(l二)o・01)八T=9：closeall;clear;symsnN=25;fw=zeros(1,N);w二2*pi/9;forn=l:1:N;holdon;f=(l/(2*9))*sinc((n*pi)/(2*9));fw(l,n)=abs

6、(f);endplot(fw(l,:));bar(fw(l,:),0.01);4、3=pi,C0=0,Cn=Bn=(4/(n*pi厂2)(2*sin(n*pi/2)-sin(n*pi)),(

7、>n=-pi/2,Fn=-j/2Bnf4(t)的幅度谱：closeall;clear;symsnN二15;fw=zeros(1,N);w二pi/2;forn=l:1:N;holdon;f二(-j/2)*(4/((n*pi)2)).*(2*sin(n*pi/2)-sin(n*pi));fw(l,n)=abs(f);endplot(fw(l,:))；bar(fw(l,:),

8、0.01);0.45f4(t)的相位谱：closeall;clear;symsnN二15;fw=zeros(1,N);w二pi/2;forn=l:1:N;holdon;f=(-j/2)*(4/((n*pi厂2))・*(2*sin(n*pi/2)-sin(n*pi));fw(l,n)二angle(f);endbar(fw(l,:),0.01);25>3=2*pi,a0=2/pi,An二(T厂(n+l)*4/(pi*(4*rT2-l)),(

9、>n=0,Fn=An/2f5(t)在正频率域的幅度谱：closeall;clear;symsnN二15;fw=zeros(

10、1,N);w二2*pi;forn=l:1:N;hol