自己编写的一个MATLAB FFT程序

clear all clc % 读取从示波器保存下来:的文档数据，通过下面的for循环来改善有些文档数据量过大 x1=csvread('G:\ArcData\lunwenshuju\shiyan\20151225\tek0005.csv',21,1); %读取通道1:通道2：通道3：通道4： % 通道1：LEM LA25-p霍尔电流传感器 采样电阻阻值99.5欧姆 % 通道2：YHDC 电流互感器 采样电阻阻值99.40欧姆 % 通道3：PA3655nl 电流互感器 采样电阻阻值196.35欧姆 变比为200： % 通道4：电流钳 TCP0030 % 实验采样率采用500k/s 5M个点数据采集 [afrow1,afcol1]=size(x1); %afrow1是采样点数；afcol1是采样通道数 for i=1:afcol1 for j=1:afrow1 if (x1(j,i)==inf)|(x1(j,i)==-inf) x1(j,i)=0; end end end freq=1000000; %采样频率 t1=1/freq; %采样周期 t=0:t1:(afrow1/freq-t1); %0~采样时间窗口 % %特殊时段处理 qishidian=1; jieshudian=1000000; x2=x1(qishidian:jieshudian,1); t2=0:t1:t1*(jieshudian-qishidian); % % %********************原始数据转换*************************** x11=x2(:,1); % x22=x2(:,2); % x33=x2(:,3); % x44=x2(:,4); % x11=(x11*1000)/20; % x22=x22; % x33=x33/196.35; % x44-x44; %**********************显示******************************** figure; % subplot(4,1,1) plot(t2,x11,'b-'); grid on; % figure;% subplot(4,1,2) % plot(t2,x44,'g'); % grid on; % figure;% subplot(4,1,3) % plot(t2,x22,'r'); % grid on; % figure;% subplot(4,1,4) % plot(t2,x33,'k'); % grid on; % ylabel('电流互感器输出电流/A'); % xlabel('时间/s'); % title('直流串联型电弧'); % [w1,YY11]=hua_fft_return(x11,1e6,1); % % [w2,YY22]=hua_fft_return(x22,5e5,1); % % [w3,YY33]=hua_fft_return(x33,5e5,1); % % [w4,YY44]=hua_fft_return(x44,5e5,1); % % y11=20*log(YY11); % % y22=20*log(YY22); % % y33=20*log(YY33); % % y44=20*log(YY44); figure; % % % subplot(3,1,1) semilogx(w1,YY11,'b'); % % set(gca,'xscale','log'); % % % axis([0 200000]); grid on; % % figure;% subplot(3,1,2) % % plot(w4,YY44,'r'); % % % set(gca,'xscale','log'); % % % % axis([0 200000]); % % grid on; % % figure;% subplot(3,1,3) % % plot(w2,YY22,'g'); % % % % plot(w2,y22,'b',w3,y33,'g'); % % % set(gca,'xscale','log'); % % grid on; % % figure;% subplot(3,1,3) % % plot(w3,YY33,'k'); % % % % plot(w2,y22,'b',w3,y33,'g'); % % % set(gca,'xscale','log'); % % grid on;