%??先從示波器數據中取樣，因示波器采樣點數太多，每3個數據取1個
clear;

%?load?PHA_i.tmp;
%?N=length（PHA_i）;
%?x1=zeros（1N）;
%?x1=PHA_i（:2）‘;
%?x=x1（200:800）;

%?
%?load?ia1.txt;?
%?b=ia1（250:3:55002）;?????%取整周期

%?load?ia2.txt;?
%?b=ia2（625:3:56752）;?????%取整周期

%?
load?u1.txt;?
b=u1（222:3:53442）;?????%取整周期

%?load?emf.txt;?
%?b=emf（730:1:61502）;?????%取整周期


%?
%?load?u2.txt;?
%?b=u2（450:3:58102）;?????%取整周期

%?x=v0（2:100001）;

M=length（b）;????????????%M為2pi區間內的采樣點數
M_end=50*M;
?for?k=0:M_end/M-1;
?????for?i=1:M;
????????x（k*M+i）=b（i）;?????%原始信號
?????end
?end
%----------------------------------------------??
n=length（x）;?%N=1024;?%總采樣點
Ts=9.6e-5*2/3;???%Ts=tp/N;?%采樣周期
tp=n*Ts;?????%tp=0.2;?%采樣總時間
detf=1/tp;???%detf=1/tp;?%頻域分辨率
t=0:Ts:Ts*（n-1）;??????%t=0:Ts:Ts*（N-1）;?%時域采樣序列
f=0:detf:（n-1）*detf;??%f=0:detf:（N-1）*detf;?%頻域采樣序列
fftx=fft（x）/（n/2）;????%FFT運算
y=abs（fftx）;

figure（2）
axes（‘FontSize‘[18]）;
plot（fy‘LineWidth‘2）;????????%半對數圖顯示
ylabel（‘Vol?line?to?line（V）‘‘fontsize‘[25]）;
%?ylabel（‘ia（A）‘‘fontsize‘[25]）
xlabel（‘harmonic（Hz）‘‘fontsize‘[25]）;
grid;
%?print?-djpeg?‘uasSpetrum.jpeg‘;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1245??2019-01-17?16:35??實用MATLABfft程序\Ifft_fromScope.m

?????目錄??????????0??2019-02-25?14:52??實用MATLABfft程序

-----------?---------??----------?-----??----

?????????????????1245????????????????????2