空氣動力學(xué)試驗報告

1、實驗一邊界層流動測量實驗摘要 ：邊界層，又稱為流動邊界、附面層，它是流體流動過程中，緊貼壁面的粘性阻力不可忽略的一層薄薄的流體，它對主要流體運動的影響很大。自普朗特提出該概念起，邊界層研究就一直是流體力學(xué)研究中一個焦點和難點課題。本實驗通過熱線風(fēng)速儀測量距離凹口平板前緣不同位置點流體的速度分布情況，并對實驗數(shù)據(jù)加以分析處理，從而確定出在不同工況中的邊界層的厚度、位移厚度，以及避免粘性力等參數(shù)，最終分析邊界層的特性。關(guān)鍵詞：邊界層，熱線風(fēng)速儀，粘性力，雷諾數(shù)，擬合，標(biāo)定1. 實驗簡介此次實驗是在一個開口式風(fēng)洞中進行的，該風(fēng)洞試驗段截面尺寸為：500mm*500mm 。設(shè)置風(fēng)洞風(fēng)機的運行頻率為20

2、Hz 和 30Hz、 ，利用熱線風(fēng)速儀測量凹槽分離點20mm 的邊界層上的速度分布。然后用兩種不同的方法擬合熱線風(fēng)速儀實驗前后標(biāo)定曲線，得出標(biāo)定誤差值，從而分析比較這兩種擬合方法的優(yōu)缺點，并分析出實驗中熱線性能的穩(wěn)定性。2. 實驗步驟1)將皮托管固定在風(fēng)洞試驗段，軸線和來流速度方向平行。記錄皮托管標(biāo)定系數(shù)ko皮托管靜壓連接到壓力傳感器負壓接口，皮托管總壓連接到壓力傳感器通道1 ；2) 熱線風(fēng)速儀探頭安裝在二位坐標(biāo)架上，連接熱線探頭與恒溫控制器輸入、輸出。 此時熱線恒溫控制器切勿通電！將熱線探頭移至和皮托管同一高度；3)熱線輸出連接到數(shù)據(jù)采集卡AI0 ,皮托管輸出連接到數(shù)據(jù)采集卡AI1 ;4)

3、將熱線恒溫控制器通電，打開MATLAB 熱線風(fēng)速儀標(biāo)定程序“hw calibration.m ”，改變文件名運行程序；5) 將熱線移動至測量點(距離凹腔分離點X=20mm )上方自由來流中，調(diào)整風(fēng)洞風(fēng)速，風(fēng)機運行頻率f=30Hz, MATLAB運行熱線速度分布測量程序“hw measurements"改變文件存儲名稱。改變風(fēng)洞風(fēng)速，風(fēng)機運行頻率f=20Hz, 重復(fù)步驟4；6) 打開 MATLAB 熱線風(fēng)速儀標(biāo)定程序hw calibration.m ，改變標(biāo)定參數(shù)存儲文件名，重新運行標(biāo)定程序。7) 按要求處理實驗數(shù)據(jù)。3. 實驗設(shè)備1) 開口式風(fēng)洞，實驗段截面500mm*500mm 。

4、實驗段風(fēng)洞壁面裝有凹腔模型，測量模型分離點上游邊界層速度分布。2) 程控二位坐標(biāo)架：計算機軟件可精確控制探針在邊界層內(nèi)移動。3) 壓力傳感器：將壓力信號轉(zhuǎn)化為電信號，靈敏度為0.0245 V/mm.4)皮托管：測量來流靜壓Ps和總壓Pt,用于測量來流速度 U和標(biāo)定熱線風(fēng)速儀。5） 熱線風(fēng)速儀：測量氣體流場(尤其是湍流場)中速度大小，具有較高的頻率響應(yīng)( 10kHz) 。6） 數(shù)據(jù)采集卡：NI-USB 6212 BNC （16位，8差分通道）；7）計算機：控制二維坐標(biāo)架，控制和執(zhí)行數(shù)據(jù)采集；其中，熱線風(fēng)速儀的工作原理如下：熱線風(fēng)速儀前段的熱線探頭作為后端控制電路中電橋的一個橋端電阻。電橋工作時把

5、探頭上的電阻絲加熱，當(dāng)流體以速度U流過探頭時帶走一部分熱量Q電橋設(shè)計是為了保持探頭熱絲溫度恒定，探頭熱絲損失熱量Q后，兩端電壓E發(fā)生變化從而提供相應(yīng)的熱量以彌補損失掉的那部分熱量，最終保證探頭熱絲的溫度恒定。通過標(biāo)準(zhǔn)速度源可以將速度U與電橋電壓E之間的關(guān)系擬合出來。這樣在測量時，只要測量到電橋電壓 E,就可以知道當(dāng)?shù)?風(fēng)速U。熱線風(fēng)速儀的標(biāo)定：用風(fēng)洞自由來流速度作為標(biāo)準(zhǔn)速度，自由來流速度U由皮托管測量得到。擬合關(guān)系式用：1） King，s Law： E = A + BU n,n = 0.45;2342） U = A + BE +CE +DE + FE 。4.實驗數(shù)據(jù)處理1）分別用兩種不同方法擬

6、合熱線風(fēng)速儀實驗前后標(biāo)定曲線，給出標(biāo)定誤差值，分析擬合方法優(yōu)缺點。將實驗前后熱線風(fēng)速儀標(biāo)定曲線畫在同一張圖上，將實驗后熱線標(biāo)定加上溫度矯正和前面兩組曲線畫在同一張圖上，分析實驗中熱線性能是否穩(wěn)定及溫度矯正的作用。表4-1實驗前標(biāo)定時采集的數(shù)據(jù)流速/U電壓/E溫度/T3.08361.696221.2363.99461.742521.2944.94061.782421.3925.86891.817121.5216.7911.847521.6547.7121.876121.7678.62041.901521.8729.55781.925621.98910.4681.947322.173表4-2實驗后

7、標(biāo)定時采集的數(shù)據(jù)流速/U電壓/E溫度/T3.05831.692422.0293.99231.738822.0094.93641.778922.0335.84981.813322.0916.78111.844422.1167.69911.872622.1918.61191.89822.2829.54051.921622.402按照 King' s Law E=A + BU n,n = 0.45關(guān)系擬合程序：clc;clear;標(biāo)定實驗 hw calibration-xxxxxx(U_E)_0527.txt');% 載入數(shù)據(jù) v=data (: 1);w=v.人0.45;e=dat

8、a (:, 2);f=inline ('a(1)+a(2)*w','a','w');b,c=lsqcurvefit(f,1.35,0.205,w,e);%b 是(A,B)向量組合，c是目標(biāo)函數(shù)值； b,c for i=1:xv(i)=data(i,1)h(i戶b*1;v(i)A0.45;endfigure;plot(v,h,'-Ar')Xlabel('流體流速U');Ylabel('電橋電壓E');title('電壓E與流速U的關(guān)系');按照U = A+BE+CE2+DE3 + FE

9、4關(guān)系擬合程序：clc;clear;標(biāo)定實驗 hw calibration-xxxxxx(U_E)_0527.txt');% 載入數(shù)據(jù) v=data(:,1);e=data(:,2);f=inline('a(1)+a(2)*e+a(3)*e.A2+a(4)*e.A3+a(5)*e.A4','a','e');b,c=lsqcurvefit(f,1 1 1 1 1,e,v);%b 是(A,B)向量組合，c是目標(biāo)函數(shù)值； b,c for i=1:xe(i)=data(i,2);h(i)=b*1 e(i) e(i).A2 e(i3 e(i)A4&

10、#39;end figure;plot(e,h,'-.g')Xlabel('電橋電壓E');Ylabel('流體流速U');title('電壓E與流速U的關(guān)系');擬合曲線為：5 -Hu 5- -Jf 11.7-1澧巾- ih.J圖1 King;s Law法擬合曲線由圖1可以看出，未擬合前實驗前后的曲線間距明顯，這就表明，在實驗結(jié)束后熱線探頭發(fā)生了一定距離的偏移。而且使用 King' s Law以合后，擬合前后的曲線還是存在一定的1 N 誤差。根據(jù)=一£ y(ui(data) -ui(mod el),實驗前曲線的

11、標(biāo)定誤差為2.618%,實驗N y ,后曲線的標(biāo)定誤差為 2.433%。1 ?I«111 H電麗電J5E圖2擬合方法二標(biāo)定曲線使用同樣的誤差公式計算得到：對實驗前的曲線的標(biāo)定誤差為：0.598%,對實驗后的曲線的標(biāo)定誤差為：0.803%。由上圖可以看出，使用擬合方法二擬合后的曲線可以與未擬合的曲線比較好地重合，說明擬合誤差很小。通過比較上面兩圖，可以得出，第二種擬合方法明顯比第一種方法得到的擬合效果好，這就說，第二種擬合方法的標(biāo)定誤差較小。同時，從實驗前后的兩條曲線可以看出，實驗前后熱線探頭發(fā)生了一定距離的偏移，這證明該實驗使用的熱線性能還不夠穩(wěn)定。2)利用文獻方法，預(yù)測測量點距離壁

12、面精確距離，測量點摩擦速度，計算壁面粘性應(yīng)力，邊界層厚度，邊界層位移厚度，邊界層動量厚度，基于動量厚度的邊界層雷諾數(shù)等參數(shù)。根據(jù)文獻指導(dǎo)并參考徐進同學(xué)的程序，然后進行必要的校正完善，得到以下程序及圖形：clc;clear;邊界層試驗 x20f20 profile.txt'data=dlmread(file);junk=sortrows(data);yi=junk(:,1)./1000;u=junk(:,2);umean=mean(data(1:20,4);u99=umean*0.99;%自由來流風(fēng)速；yd=-9.9:0.01:-4;for j=1:length(yd) %截止距離循環(huán)迭

13、代；y=yi+yd(j)./1000;ut=0.20:0.01:1;lp=1.5e-5./ut;k=0.41;B=5.0;for i=1:length(ut) %壁面摩擦應(yīng)力循環(huán)迭代 ；yp=y./lp(i);up=u./ut(i);ye=up+exp(-k*B)*(exp(k*up)-1-k*up-(k*up).A2/2-(k*up).A3/6);err(i,j戶sum(abs(yp(1:35)-ye(1:35)./yp(1:35)/35;% 殘余函數(shù)計算; endendm indexi=min(err);% 行數(shù)；n indexj尸min(m);% 列數(shù)；yoff=yd(indexj)%

14、截止距離；utf=ut(indexi(indexj)% 壁面摩擦應(yīng)力； lpf=1.5e-5./utf;yp=(yi+yoff./1000)./lpf；upp=1：0.5：25；up=junk(：,2)./utf;ue=(1./k).*log(yp)+B;ye=upp+exp(-k*B)*(exp(k*upp)-1-k*upp-(k*upp).A2/2-(k*upp).A3/6);figure。)；semilogx(yp,up,'-ob','linewidth',1); hold on;semilogx(ye,upp,'-or','lin

15、ewidth',1);hold on;semilogx(yp,ue,'-k','linewidth',1);xlabel('yA+')；ylabel('UA+')；legend('Experimental data','Spalding profile','Log-law');axis(1 1000 0 26);yabs=yi+yoff./1000;ul=1-u/umean;yabsp=trapz(yabs,ul).*1000%邊界層位移厚度；uul=(u/umean).*(

16、1-u./umean);theta=trapz(yabs,uul).*1000% 邊界層動量厚度； Re=1.2*umean*theta/1000./(1.5e-5)% 基于動量厚度的邊界層雷諾數(shù)；figure(2);Yreal=yi.*1000+yabsp;Ureal=junk(:,2);plot(Ureal,Yreal,'-oc')Xlabel('Ureal'),Ylabel('Yreal');axis(0 8 0 120);取f=20Hz時，其風(fēng)速近似為 6.26m/s,經(jīng)循環(huán)迭代得到，測點距壁面的誤差距離為 9.34mm；測量點摩擦速度為

17、0.34m/s；計算出壁面的剪切應(yīng)力為 0.149N/m2；邊界層厚度從上圖可以近似得到，為12mm；邊界層位移厚度為 2.3755mm；邊界層動量厚度為 2.0819mm；基于動量厚度的雷諾數(shù)為1051;取f=30Hz時，其風(fēng)速近似為9.3m/s,經(jīng)循環(huán)迭代得到，測點距壁面的誤差距離為7.99mm；測量點摩擦速度為 0.46m/s；計算出壁面的剪切應(yīng)力為0.273N/m2;邊界層厚度從上圖可以)；,邊界層速度rms分。（圖例為20Hz）仃Q>壁工引生工鳳華生HU近似得到，為11mm；邊界層位移厚度為 2.0217mm；邊界層動重厚度為 1.5098mm；基于動 量厚度的雷諾數(shù)為1122

18、。3）畫出邊界層速度分布曲線（基于邊界層厚度和自由來流無量綱化）布曲線。畫出邊界層速度分布曲線（基于邊界層內(nèi)層參數(shù)無量綱化）clc;clear;邊界層試驗 x20f20 profile.txt'data=dlmread(file);junk=sortrows(data);yi=junk(:,1);ui=junk(:,2);y=yi./12;u=ui./6.26;plot(u,y,'-ob');Xlabel（'基于自由來流風(fēng)速無量綱化）,Ylabel（'基于邊界層厚度無量綱化 axis(0.45 1.1 0 10) clc;clear;邊界層試驗 x20f

19、20 profile.txt'data=dlmread(file);junk=sortrows(data);yi=junk(:,1);rms=junk(:,3);y=yi-9.34;plot(rms,y,'-ob');Xlabel('rms'),Ylabel('邊界層');axis(0 0.8 0 100)clc;clear;邊界層試驗 x20f20 profile.txt'data=dlmread(file);junk=sortrows(data);yi=junk(:,1);ui=junk(:,2);y=(yi-9.34)./1

20、2;u=ui./0.34;plot(u,y,'-ob');Xlabel('基于邊界層摩擦速度無量綱化),Ylabel('基于邊界層厚度無量綱化');axis(8 20 0 10)c 81Q1*14,-桑理-& +以上圖形均為在20Hz下的數(shù)據(jù)得到的，30Hz的相應(yīng)圖形用同樣的方法即可得到，在 此不再贅述。4）對比分析討論不同風(fēng)洞運行速度下邊界層速度分布曲線特性，邊界層基本參數(shù)。«* 立 4fl- 4 4 _2 4 首瞿日雷項目20Hz30Hz自由來流速度m/s6.269.3測量點摩擦速度 m/s0.340.46邊界層厚度/mm1211邊

21、界層位移厚度/mm2.37552.0217邊界層動量厚度/mm2.08191.5098邊界層雷諾數(shù)10511122通過mat lab編程將兩種頻率下的三類圖形放在一起進行比較，由數(shù)據(jù)文件獲知，在風(fēng)洞運行頻率f=20Hz時，自由來流速度約為6.26m/s,而在風(fēng)洞運行頻率f=30Hz時，自由來流速度為9.3m/s,從圖上可以看到，無量綱化后的圖形幾乎一致，邊界層速度分布曲線其基本走勢均為距離壁面較近的區(qū)域速度逐漸增加，到一定程度后就不再增加，這段距離即為邊界層厚度。5)對比分析討論不同風(fēng)洞運行速度下剪切層速度大小，rms分布曲線特性。clc;clear;邊界層試驗 x20f20 profile.

22、txt'邊界層試驗 x20f30 profile.txt'data1=dlmread(file1);data2=dlmread(file2);junk1=sortrows(data1);junk2=sortrows(data2);yi1=junk1(:,1);ui1=junk1(:,2);yi2=junk2(:,1);ui2=junk2(:,2);y1=yi1-9.34;u1=ui1;y2=yi2-7.99;u2=ui2;figure。)；plot(u1,y1,'-og',u2,y2,'-*c');legend('20Hz',&#

23、39;30Hz');Xlabel('V'),Ylabel('Y');title('不同頻率下剪切層速度大小比較);hold onyi1=junk1(:,1);rms1=junk1(:,3);yi2=junk2(:,1);rms2=junk2(:,3);y1=yi1-9.34;y2=yi2-7.99;figure(2);plot(rms1,y1,'-og',rms2,y2,'-*c');legend('20Hz','30Hz');Xlabel('rms'),Ylabel

24、('Y');title('不同頻率下rms大小比較);用他工I的的層鹿底大.1，也必20Hz下的剪切層速度小于 30Hz下的剪切層速度rms圖上可以看到，基本趨勢兩種風(fēng)速下是一致的，但是距壁面相同距離時，30Hz下的rms值要大于20Hz時的值，而且 30Hz時的rms最大值為1.18,而20Hz時的最大值為 0.704。6）在Urms最大值處，速度的能量譜圖，分析能量譜特征。從圖上可以看到峰值出現(xiàn)的地方，對應(yīng)的頻率即為該信號主要能量來源，從低到高頻率階次依次為100Hz,160Hz,220Hz,280Hz ,而且由圖形趨勢看到其主要能量集中于低頻處，頻 率越高處能量越

25、低。實驗二數(shù)據(jù)采集與處理1 .實驗儀器1） 信號發(fā)生器及示波器:DSOX2012A ;2） 數(shù)據(jù)采集卡： NI-USB 6212 BNC （16 位，8 差分通道，400 kS/s）;3）計算機：控制和執(zhí)行數(shù)據(jù)采集；2 .實驗步驟1）連接信號發(fā)生器，數(shù)據(jù)采集卡，電腦等設(shè)備；2）設(shè)置信號發(fā)生器產(chǎn)生正弦波信號，信號頻率為f=200Hz,信號幅值為A=0.05v;3）信號發(fā)生器產(chǎn)生信號連接到數(shù)據(jù)采集卡AI0 ;4) 打開MATLAB數(shù)據(jù)采集程序 'dataacquisition.m ,'設(shè)置采樣頻率Fs=100Hz, 通道輸入范圍(InputRange')為-0.2 0.2。

26、更改數(shù)據(jù)存儲文件名，運行程序；5) 檢查步驟4中數(shù)據(jù)是否成功保持。重復(fù)步驟4改變采樣頻率Fs=200Hz。更改數(shù)據(jù)存儲文件名，運行程序；6) 檢查步驟5中數(shù)據(jù)是否成功保持。重復(fù)步驟4改變采樣頻率Fs=400Hz。更改數(shù)據(jù)存儲文件名，運行程序；7) 檢查步驟6中數(shù)據(jù)是否成功保持。重復(fù)步驟4改變采樣頻率 Fs=1000Hz。更改數(shù)據(jù)存儲文件名，運行程序；8)檢查步驟7中數(shù)據(jù)是否成功保持。重復(fù)步驟7,改變通道輸入范圍('InputRange')為-1010。更改數(shù)據(jù)存儲文件名，運行程序。3.實驗數(shù)據(jù)處理1)畫出數(shù)據(jù)隨時間變化曲線(帶點線) x(t);不同參數(shù)下所得到的曲線畫到同一張圖

27、上， 用不同顏色標(biāo)示。討論區(qū)別，及采樣參數(shù)影響。clc;clear;數(shù)據(jù)采集試驗 g1y100r02.bin','rt');數(shù)據(jù)采集試驗 g1y200r02.bin','rt');數(shù)據(jù)采集試驗 g1y400r02.bin','rt');數(shù)據(jù)采集試驗 g1y1000r02.bin','rt');數(shù)據(jù)采集試驗 g1y3000r02.bin','rt');數(shù)據(jù)采集試驗 g1y3000r10.bin','rt');A1,NUM1=fread(fip1,inf,

28、'double')A2,NUM2=fread(fip2,inf,'double')A3,NUM3=fread(fip3,inf,'double')A4,NUM4=fread(fip4,inf,'double')A5,NUM5=fread(fip5,inf,'double')A6,NUM6=fread(fip6,inf,'double')t1=(0:1:9)*10A-2;figure。)；plot(t1,A1(1:10),'- Or');% 第一 hold ont2=(0:1:19)/

29、2*10A-2;plot(t2,A2(1:20),'-+m');hold ont3=(0:1:39)/4*10A-2;plot(t3,A3(1:40),'-*g');hold ont4=(0:1:99)*10A-3;plot(t4,A4(1:100),'->c'); hold ont5=(0:1:299)/3*10A-3;plot(t5,A5(1:300),'-<b');Xlabel('時間 t/s');Ylabel('幅值 V/v');title('不同采樣頻率下的曲線比較)f

30、igure(2);t6=(0:1:299)/3*10A-3;plot(t6,A5(1:300),'-<b');hold ont7=(0:1:299)/3*10A-3;plot(t7,A6(1:300),'-+k');Xlabel('時間 t/s');Ylabel('幅值 V/v');title('不同采樣范圍下的曲線比較)Sis惜尚M從圖上可以清楚的看到：采樣頻率取100Hz和200Hz時完全失真，采到的為兩條橫線，當(dāng)取400Hz （原始信號頻率為 200Hz）時，采樣信號有一點變化，但還離原始信號差很遠， 因為不滿

31、足香農(nóng)采樣定理（采樣頻率>400Hz）,當(dāng)把頻率取到1000Hz時，采樣信號有了原始信號的輪廓，但還是失真比較大，當(dāng)把頻率取為 3000Hz時，采樣信號已經(jīng)很大程度的逼近了原始信號，然后將采樣幅值范圍取為-10, 10,在峰值處有了平臺信號失真，所以幅值取為-0.2,0.2非常接近真實信號。因此在采集數(shù)據(jù)時，必須滿足采樣定理，同時信號的靈敏 度要求幅值與真實信號幅值接近。2）計算采集到數(shù)據(jù)的平均值 xm,均方差xrmso比較不同采樣參數(shù)下得到值的區(qū)別。討論 區(qū)別，及采樣參數(shù)影響。clc;clear;數(shù)據(jù)采集試驗 g1y100r02.bin','rt');A1,NU

32、M1=fread(fip1,inf,'double');sum1=0;for i=1:NUM1sum1=A1(i)+sum1;endx1=sum1/NUM1sum2=0;for i=1:NUM1cha=A1(i)-x1;fang=cha.A2;sum2=sum2+fang;endxrms1=sqrt(sum2/NUM1)依次算出6種情況下的均值和均方差(即標(biāo)準(zhǔn)差)。表2-1不同情況下采集數(shù)據(jù)的均值和均方差列表比較項目100Hz0.2200Hz0.2400Hz0.21000Hz0.23000Hz0.23000Hz10平均0.0110-0.0168-3.7972e-4-3.5217

33、e-4-3.5032e-4-0.0010值均方0.00209.5665e-40.00320.01770.01770.0177差編程計算求得各自的均值和均方差。從上表可以看出，100Hz和200Hz時均方差較小，而均值較大，這與實際情況一致，因為其圖形近乎為一直線； 400Hz時嚴(yán)重失真，圖形呈鋸 齒波形，均值較小近乎為0;而后面三種高頻采集的信號均值都很小，而且均方差值接近，較好地反映了原始信號，相比之下低的采樣頻率采集的信號均方差的誤差較大。3) 用mat lab編程計算所得數(shù)據(jù)的傅里葉變化X( k),驗證Parsevals Theorem。畫出|X(k)|隨頻率變化曲線。不同采樣參數(shù)得到的

34、曲線畫到同一張圖上，用不同顏色標(biāo)示。討論區(qū)別，及采樣參數(shù)影響。clc;clear;數(shù)據(jù)采集試驗 g1y100r02bin','rt');數(shù)據(jù)采集試驗 g1y200r02bin','rt');數(shù)據(jù)采集試驗 g1y400r02bin','rt');數(shù)據(jù)采集試驗 g1y1000r02bin','rt');數(shù)據(jù)采集試驗 g1y10000r02bin','rt');數(shù)據(jù)采集試驗 g1y10000r10.bin','rt');data1=fread(fip1,51

35、2,'double')data2=fread(fip2,512,'double')data3=fread(fip3,512,'double') data4=fread(fip4,512,'double') data5=fread(fip5,512,'double') data6=fread(fip6,512,'double') N=512; %512個采樣點 figure。)；Fs=100; %采樣頻率t=(0:N-1)/Fs; % 時域Y=fft(data1,N); %快速傅里葉變換 f=(0:N

36、-1)*Fs/N; % 頻域的點數(shù) plot(f,abs(Y),'-sr') % 頻譜圖 hold onFs=200;t=(0:N-1)/Fs;Y=fft(data2,N);f=(0:N-1)*Fs/N;plot(f,abs(Y),'-dm') hold onFs=400;t=(0:N-1)/Fs;Y=fft(data3,N);f=(0:N-1)*Fs/N; plot(f,abs(Y),'-hg') hold onFs=1000;t=(0:N-1)/Fs;Y=fft(data4,N);f=(0:N-1)*Fs/N;plot(f,abs(Y),&#

37、39;-*c') hold onFs=10000;t=(0:N-1)/Fs;Y=fft(data5,N);f=(0:N-1)*Fs/N;plot(f,abs(Y),'-oy') hold onFs=10000;t=(0:N-1)/Fs;Y=fft(data6,N);f=(0:N-1)*Fs/N;plot(f,abs(Y),'-pk')43210420SO10D501001502D010844204 口。4noF400060008D00TQQQ 004I 口位10000從上圖6種采樣信號單獨的傅立葉變換后的圖形，可以很明顯地看到，100Hz和200Hz2000-000200600080D010000200-3兩種情況下，只有 0Hz處有幅值，屬于直流分量，而其余部分近乎為零，沒有能量，沒有采集到信號，或者說采集到的信號完全失真。一般來說，滿足采樣定理的信號頻譜圖是對稱的，后面的幾種情況就是這樣，明顯看到在200Hz處出現(xiàn)峰值，表示能量集中在此頻率下， 說明原始信號的頻率為 200Hz,這與實際情況是相符的。在后兩