matlab的傅里葉變換小實(shí)驗(yàn)

1、三項(xiàng)電壓波形程序如下：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'); %畫(huà)出A項(xiàng)波形圖形hold on;

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

3、 %y軸標(biāo)注為電壓/Kv三項(xiàng)電壓波形如下圖所示：圖一 三項(xiàng)電壓波形圖四 C 項(xiàng)電壓波形將上述三項(xiàng)電壓波形中 A,B,C 項(xiàng)分別進(jìn)行傅立葉變換，程序如下：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）看了下，采樣時(shí)間間隔一樣，就只算了一個(gè)差值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