关于CZT对谐波的检测

cckkcclove

　　电力系统谐波信号进行检测，采用频率域插值小波变换法，谐波信号为x=0.3*cos(2*pi*f1*t+70)+0.7*cos(2*pi*f2*t+80)+1.0*cos(2*pi*f3*t)+0.5*cos(2*pi*f4*t+90)+0.4*cos(2*pi*f5*t+40);
　　由采样定理，信号中最高频率为150，fs=100
　　程序：
　　fs=500;
　　N=1024;
　　f1=25,f2=35.85,f3=50,f4=86.6,f5=150;
　　m=0:N-1;
　　t=m/fs;
　　x=0.3*cos(2*pi*f1*t+70)+0.7*cos(2*pi*f2*t+80)+1.0*cos(2*pi*f3*t)+0.5*cos(2*pi*f4*t+90)+0.4*cos(2*pi*f5*t+40);
　　w=exp(-j*2*pi*200/(N*fs));
　　a=1;
　　z=czt(x,N,w,a);
　　fz=(0:length(z)-1)'*200/length(z);
　　plot(fz,abs(z)*2/N);
　　axis([0 200 0 1.2]); grid;




　　请简单描述一下该程序该图
　　答辩用 急用 谢谢