matlab的傅里葉變換小實驗

1、三項電壓波形程序如下：A=xlsread('yl.xls'); %讀取yl.xls中的數(shù)據(jù)x=A(:,1);y1=A(:,3);y2=A(:,4);y3=A(:,5); %x為矩陣A的第一列；y1為矩陣A的第三列；y2為矩陣A的第%四列；y3為矩陣A的第五列figure();plot(x,y1,'r'); %畫出A項波形圖形hold on;

2、 %保持圖形plot(x,y2,'g'); %畫出B項波形圖形hold on; %保持圖形plot(x,y3,'b'); %畫出C項波形圖形hold on; %保持圖形xlabel('時間/t'); %x軸標(biāo)注為時間/tylabel('電壓/Kv');

3、 %y軸標(biāo)注為電壓/Kv三項電壓波形如下圖所示：圖一 三項電壓波形圖四 C 項電壓波形將上述三項電壓波形中 A,B,C 項分別進行傅立葉變換，程序如下：A=xlsread('yl.xls'); %讀取yl.xls中的數(shù)據(jù)x=A(:,1);y1=A(:,3);y2=A(:,4);y3=A(:,5); %x為矩陣A的第一列；y1為矩陣A的第三列；y2為矩陣A的第四列；y3為矩陣A的第五

4、列dt=x(2)-x(1); %用diff（t）看了下，采樣時間間隔一樣，就只算了一個差值fs=1/dt; %傅里葉變換后的坐標(biāo)軸的最大值就是采樣頻率fs，與dt是倒數(shù)關(guān)系ay1=abs(fft(y1));ay2=abs(fft(y2));ay3=abs(fft(y3)); %fft是matlab中傅里葉變換的命令，abs是求模f=(0:length(x)-1)*fs/len