Matlab文件的編寫

Myfft.m是一個自編的matlab函數。其基本思路是，根據輸入的采樣頻率和采樣點數，建立數據庫連接后，取出指定表中的數據進行FFT分析，并繪制出時域和頻域波形進行對比分析。
function myfft(n,fs,sqlstring)  %建立myfft函數
conn=database('matlab_test','','');  %連接數據源matlab_test,用戶名為sa，密碼為空
ping(conn);
cursor1=exec(conn,sqlstring);  %獲得游標
cursor1=fetch(cursor1);                  %獲得數據 
A=cursor1.data;
a=cell2mat(A);                           %單元數組轉換為矩陣   
g=size(a,1);
t=1:1024;
subplot(2,1,1)
plot(t,a,'r');                      %繪制時域波形
title('--時域波形    幅值,時間(t)','fontsize',16,'fontweight','bold');
grid('on')
m=n/2;
y=fft(a,n);                             %調用malab的fft函數進行fft分析
Y=abs(y);
f=fs*(0:m)/n;                            %設置x軸坐標
subplot(2,1,2)
plot(f,Y(1:(m+1)),'r');                   %繪制頻域波形
xlabel('--頻域波形圖  幅值 , 頻率(Hz)','fontsize',16,'fontweight','bold');
grid('on');