諧波分析C語(yǔ)言程序

諧波分析C語(yǔ)言程序諧波分析C語(yǔ)言程序諧波分析C語(yǔ)言程序諧波分析C語(yǔ)言程序編制僅供參考審核批準(zhǔn)生效日期地址：電話：傳真:郵編:諧波分析1）作業(yè)要求編程求出圖一所示波形的最大值、平均值、有效值、求出1~8次諧波的幅值和頻率。圖一被分析波形2）解題思路、擬定方法、使用工具軟件描述解題思路：上圖為周期20ms的正弦波，該波形為奇函數(shù)。（1）n=0,1,2，…（2）n=1,2,…（3）因?yàn)閒（t）為奇函數(shù)，故=0。（4）故要求得的1~8次諧波可分為兩個(gè)步驟：其一，求的值；其二，求。擬定方法：先構(gòu)造上述函數(shù)，取一個(gè)周期分成1000等份。每個(gè)等份對(duì)應(yīng)一個(gè)數(shù)值，在對(duì)這些數(shù)值進(jìn)行一定的運(yùn)算后便可求得最大值、平均值、有效值和1~8次諧波的幅值和頻率。使用工具軟件描述：用c語(yǔ)言軟件求得最大值、平均值、有效值和1~8次諧波的幅值和頻率，并顯示出來(lái)。3）子函數(shù)描述f(t)函數(shù)：f（t）為所要構(gòu)建的周期性方波信號(hào)波形，在第一個(gè)周期內(nèi)（）當(dāng)時(shí)f（t）=100,當(dāng)時(shí)f（t）=-100，其余時(shí)刻f（t）=0。再通過(guò)一個(gè)整除函數(shù)進(jìn)行周期延拓。f(t){ doublei;if(t>=&&t<= i=;elseif(t>=&&t<= i=;else i=;return(i);}jf(a,b)函數(shù)：jf(a,b)函數(shù)為求f（t）有效值,設(shè)a，b為一個(gè)周期的起點(diǎn)和終點(diǎn)。有效值=jf(a,b){intj,n=1000;doublei;doubleh;h=(b-a)/1000;i=(f(a)*f(a)+f(b)*f(b))*h/2;for(j=1;j<n;j++)i+=f(a+j*h)*f(a+j*h)*h;i=sqrt(i/(b-a))；returni;}max=f(a)max=f(a),i=0i++≠1000f(a+i×h)≥maxmax=f(a+i×h)結(jié)束NNYYf\n",max);f\n",ave);f\n",value);f，頻率為%d\n",k,d,k*50);}///*求理論值for(j=1;j<8;){m=100*(cos(j**-cos(j**-cos(j**+cos(j**)/j/;j=j+2;printf("%f/n",m);}//*/}//構(gòu)造f（t）函數(shù)doublef(doublet){ doublei;if(t>=&&t<= i=;elseif(t>=&&t<= i=;else i=;return(i);}doublejf(doublea,doubleb){intj,n=1000;doublei=0,k=0;doubleh;h=(b-a)/(double)n;i=(f(a)*f(a)+f(b)*f(b))*h/;for(j=1;j<n;++j)i+=f(a+j*h)*f(a+j*h)*h;i=sqrt(i/(b-a));returni;}//求傅立葉變換doublejf1(doublea,doubleb,intk){intj,n=1000;doublei,h;h=(b-a)/;i=(f(a)*sin(100**k*a)+f(b)*f(100**k*b))*h/2;for(j=1;