fft.rar_FFT谐波提取_MATLAB 谐波提取_傅里叶 谐波_提取各次谐波_谐波提取

%原始信号图形输出 %clear clc Fs=8652; t=(1:8652)/Fs; % load Cigre_06.out % load Cigre_07.out % load Cigre_16.out % load Cigre_2.out % load Cigre_04.out % iy0=Cigre_06(:,11); % is0=Cigre_07(:,11); % iy=Cigre_16(:,11); % is=Cigre_2(:,11); % e=Cigre_04(:,11); % dy=iy-iy0; % ds=is-is0; %filename='F:\xsw\新建文件夹\pinpu.emt\noname_01.out'; %读入数据,只能读入.out； %y=textread(filename); %delt_a1=sin(2*pi*200*t); % delt_a=y(:,2); % delt_t=y(:,1); % subplot(2,1,1); % plot(delt_t,delt_a); % title('A相电流分量信号'); % set(gcf,'color','white'); % xlabel('Time(seconds)'); ylabel('Current waveform'); % subplot(2,1,2); % plot(delt_t,delt_a1); % title('对侧相电流分量信号'); % set(gcf,'color','white'); % xlabel('Time(seconds)'); ylabel('Voltage waveform'); % subplot; id1=ia444(:,2); N=173; %N=length(id1); for i=2500:8479; %figure; Ya=fft(id1(i:(i+N)),N);%做FFT变换 %Yc=fft(delt_a1,N); for t=1:N; Mag_a(i-2499,t)=(abs(Ya(t))); end p(i-2499)=Mag_a(i-2499,3)./Mag_a(i-2499,2); end plot(p); %Mag_c=(abs(Yc)); % subplot(2,1,1); % plot(Mag_a(1:N)); % title('A相电流幅频图'); % subplot(2,1,2); % plot(Mag_a1(1:N)); % title('对策a相电流幅频图'); % subplot; %figure; % Mag_a=Mag_a/(N/2);%换算成实际的幅度 % %Mag_c=Mag_a1/(N/2); % Mag_a(1)=Mag_a(1)/2; %直流分量的换算 % %Mag_c(1)=Mag_a1(1)/2; % F=([1:Fs]-1)*Fs/N;%换算成实际的频率值 % %subplot(2,1,1); % stem(F(1:N/2),Mag_a(1:N/2));%只计算半波即可 % title('换算后的A相电流幅频图'); % subplot(2,1,2); % plot(F(1:N/2),Mag_a1(1:N/2)); % title('换算后的C相电流幅频图'); % subplot; % figure; % subplot(2,1,1); % Ph_a=[1:N/2]; % for i=1:N/2; % Ph_a(i)=phase(Ya(i)); %计算相位 % Ph_a(i)=Ph_a(i)*180/pi; %换算为角度 % end; % plot(F(1:N/2),Ph_a(1:N/2)); %显示相位图 % title('A相相位-频率曲线图'); % subplot(2,1,2); % Ph_a1=[1:N/2]; % for i=1:N/2; % Ph_a1(i)=phase(Yc(i)); %计算相位 % Ph_a1(i)=Ph_c(i)*180/pi; %换算为角度 % end; % plot(F(1:N/2),Ph_a1(1:N/2)); %显示相位图 % title('C相相位-频率曲线图');