谐波分析C语言程序

谐波分析C语言程序谐波分析C语言程序谐波分析C语言程序谐波分析C语言程序编制仅供参考审核批准生效日期地址：电话：传真:邮编:谐波分析1）作业要求编程求出图一所示波形的最大值、平均值、有效值、求出1~8次谐波的幅值和频率。图一被分析波形2）解题思路、拟定方法、使用工具软件描述解题思路：上图为周期20ms的正弦波，该波形为奇函数。（1）n=0,1,2，…（2）n=1,2,…（3）因为f（t）为奇函数，故=0。（4）故要求得的1~8次谐波可分为两个步骤：其一，求的值；其二，求。拟定方法：先构造上述函数，取一个周期分成1000等份。每个等份对应一个数值，在对这些数值进行一定的运算后便可求得最大值、平均值、有效值和1~8次谐波的幅值和频率。使用工具软件描述：用c语言软件求得最大值、平均值、有效值和1~8次谐波的幅值和频率，并显示出来。3）子函数描述f(t)函数：f（t）为所要构建的周期性方波信号波形，在第一个周期内（）当时f（t）=100,当时f（t）=-100，其余时刻f（t）=0。再通过一个整除函数进行周期延拓。f(t){ doublei;if(t>=&&t<= i=;elseif(t>=&&t<= i=;else i=;return(i);}jf(a,b)函数：jf(a,b)函数为求f（t）有效值,设a，b为一个周期的起点和终点。有效值=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)结束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);}//*/}//构造f（t）函数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;