【文章內容簡介】 pi*n*T)。%對x(t)32點采樣 Xk32=fft(xnT)。%計算xnT的32點DFT Xk32=fftshift(Xk32)。%將零頻率移到頻譜中心 Tp=N*T。F=1/Tp。%頻率分辨率F k=N/2:N/21。fk=k*F。%產生16點DFT對應的采樣點頻率（以零頻率為中心）subplot(3,1,2)。stem(fk,abs(Xk32),39。.39。)。box on %繪制8點DFT的幅頻特性圖 title(39。32點|DFT[x(nT)]|39。)。xlabel(39。f(Hz)39。)。ylabel(39。幅度39。)。axis([N*F/21,N*F/21,0,*max(abs(Xk32))])N=64。n=0:N1。%FFT的變換區(qū)間N=16 xnT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T)。%對x(t)64點采樣 Xk64=fft(xnT)。%計算xnT的64點DFT Xk64=fftshift(Xk64)。%將零頻率移到頻譜中心 Tp=N*T。F=1/Tp。%頻率分辨率F k=N/2:N/21。fk=k*F。%產生16點DFT對應的采樣點頻率（以零頻率為中心）subplot(3,1,3)。stem(fk,abs(Xk64),39。.39。)。box on%繪制8點DFT的幅頻特性圖 title(39。 64點|DFT[x(nT)]|39。)。xlabel(39。f(Hz)39。)。ylabel(39。幅度39。)。axis([N*F/21,N*F/21,0,*max(abs(Xk64))])內容N=82018161412|X(f)|10864200123Hz456圖 實驗內容1采樣點數(shù)N=8時的頻譜N=612108|X(f)|64200123Hz456圖 實驗內容1采樣點數(shù)N=6時的頻譜N=3280706050|X(f)|40302010001234Hz5678圖 實驗內容1采樣點數(shù)N=32時的頻譜N=21151050200幅度譜150100500HzN=1*** 實驗內容2頻譜圖（1）實驗內容1，根據(jù)給定模擬f1=2Hz，f=4Hz,所以滿足采樣定理2fs179。2fh，采取樣頻率fs=16Hz,F=42=2Hz,根據(jù)N179。`2f0/F,所以N179。8，當N=6時，不滿足采樣點數(shù)，所以頻譜在時域混疊，而當N=32時，滿足采樣點數(shù),頻譜中能夠看到有兩個頻率成分，在f1=2Hz,f2=4Hz處有譜線，。(2)實驗內容2,f1=120Hz,f2=140Hz,取采樣頻率fs=400Hz,F=140120=20Hz,所以N=20，顯然，當N=21時，能把f1和f2兩個譜峰分辨出來，當N=11時，兩個譜峰分辨不出來。16點|DFT[x(nT)]|10幅度503020100f(Hz)32點|DFT[x(nT)]|102030幅度1003020100f(Hz)64點|DFT[x(nT)]|102030幅度2003020100f(Hz) 實驗內容3頻譜圖（3）實驗內容3，對莫逆周期信號進行譜分析x(t)=cos8Pt+cos16Pt+cos20Pt x(t)有3個頻率成分，f1樣頻率=4Hz,f2=8Hz，f3=10Hz,所以x(t)，采fs=64Hz。變換區(qū)間N=16時，觀察時間Tp=16T=,不是x(t)的整數(shù)倍周期，所以所得頻譜不正確，變換區(qū)間N=32，64時，觀察時間Tp=,1s,是x(t)的整數(shù)周期，所以所得頻譜正確，圖中3根譜線正好位于4Hz,8Hz,10Hz處。變換區(qū)間N=64時，頻譜幅度是變換區(qū)間N=32時的2倍，這種結果正好驗證了用DFT對中期序列譜分析的理論。（1）完成各個實驗任務和要求，附上程序清單和有關曲線。（2）簡要回答思考問題第二篇：MATLAB仿真技術實驗教案《MATLAB仿真技術》 實驗教案實驗一實驗名稱：熟悉Matlab交互工作界面一、實驗目的熟悉Matlab各種工作界面的操作要旨掌握Matlab的基本操作命令二、實驗步驟命令窗口（1）體驗命令窗口的菜單及各項功能（2）嘗試命令窗口編輯特殊功能鍵和設置工作空間窗口與當前路徑窗口（1）在工作空間窗口查看及修改變量（2）添加新的路徑為Matlab路徑圖形窗口和文本編輯窗口（1）練習圖形窗口中修改圖形的方法（2）在文本編輯窗口調試程序體會Matlab的基本操作命令三、實驗儀器PC機 MATLAB軟件四、實驗結果五、結論實驗二實驗名稱：Matlab在符號計算方面的應用一、實驗目的掌握標識符的生成和使用掌握矩陣及變量的賦值熟悉三類運算符及其功能二、實驗內容標識符的生成和使用1）、計算y=x+()/(x+)5(x+1/x)，當x=2和x=4時的值。x=[2 4]。y=x.^3+().^2./(x+).^35*(x+1./x)。y y = 32）、計算cos6092。o323 y=cos(pi/3)(9sqrt(2))^(1/3)y =矩陣及變量的賦值21)、已知a=3,A=4，b=a,B=b21,c=a+A2B,C=a+2B+c,求C a=3。A=4。b=a^2。B=b^21。c=a+A2*B。C=a+2*B+c。C C = 2）、創(chuàng)建34矩陣魔方陣和相應的隨機矩陣，將兩個矩陣并接起來，然后提取任意兩個列向量。 A=magic(4)。A(4,:)=[]。B=rand(3,4)。C=[A B]。D=C(:,3)。E=C(:,4)。D,E D = 3 10 6E =8 12 3）、創(chuàng)建一個55隨機陣并求其逆。 A=rand(5,5)。B=inv(A)B = 4）利用上題的矩陣，計算矩陣的5次方。 A=rand(5,5)。C=A^5 C = 三類運算符及其功能1）創(chuàng)建一個44隨機矩陣，提取對角線以上部分。 A=rand(5,5)。D=triu(A,1)D =0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2）創(chuàng)建一個45隨機矩陣。 A=rand(4,5)。B=A([1 2],:)。C=B.*A([1 2],:)C = 0 4813249。432249。233。1233。5234。，B=234。6238，求C=A*B,D=A.*B。593）設A=36234。234。234。234。235。27128235。1397 A=[1 4 8 13。3 659。27128]。B=[5 4 32。62 38。1 39 7]。C=A*B39。 D=A.*B C =1982 30 12 27 338 54 29D =2426181215 72221 10856三、實驗儀器PC機 MATLAB軟件四、實驗結果五、結論實驗三實驗名稱：一維二維數(shù)組創(chuàng)建與尋訪一、實驗目的掌握數(shù)組的創(chuàng)建和尋訪掌握數(shù)組的生成函數(shù)和操作運算掌握數(shù)組操作運算與矩陣操作運算二、實驗步驟直接輸入法創(chuàng)建數(shù)組創(chuàng)建一個數(shù)組并賦值給變量a 233。1234249。234。5678 使a=a=234。234。9101112234。13141516235。使用MATLAB函數(shù)創(chuàng)建數(shù)組創(chuàng)建一個數(shù)組并賦值給變量b，使b中元素服從均值為0，方差為1的正態(tài)分布（b為44階）利用外部數(shù)據(jù)文件（MAT）保存和裝載數(shù)組 保存并調用變量a,b比較數(shù)組運算與矩陣運算操作的區(qū)別 ⑴分別對變量a,b 進行矩陣乘和數(shù)組乘 ⑵分別對變量a,b 進行矩陣除和數(shù)組除三、實驗結果直接輸入法創(chuàng)建數(shù)組創(chuàng)建一個數(shù)組并賦值給變量a 233。1234249。234。5678 使a=a=234。234。9101112234。13141516235。MATAB源程序： a=[1 2 3 4。5 6 7 8。9 10 11 12。13 14 15 16] a =13使用MATLAB函數(shù)創(chuàng)建數(shù)組創(chuàng)建一個數(shù)組并賦值給變量b，使b中元素服從均值為0，方差為1的正態(tài)分布（b為44階）MATAB源程序： b=randn(4,4)b = 利用外部數(shù)據(jù)文件（MAT）保存和裝載數(shù)組 保存并調用變量a,b MATAB源程序： save D:55 a b clear load D:55 a b比較數(shù)組運算與矩陣運算操作的區(qū)別 ⑴分別對變量a,b 進行矩陣乘和數(shù)組乘 MATAB源程序： a*bans = a.*bans = ⑵分別對變量a,b 進行矩陣除和數(shù)組除 MATAB源程序： a/bans = ab Warning: Matrix is close to singular or badly may be = =+016 * a./bans = ans =三、實驗儀器PC機 MATLAB軟件四、實驗結果五、結論實驗四實驗名稱：Matlab在數(shù)值計算方面的應用一、實驗目的學會使用MATLAB關系運算符學會使用MATLAB邏輯運算符掌握多項式的表示方法及運算熟練掌握精度可控的數(shù)值積分二、實驗步驟關系運算符生成一個包含10個 元素的行向量A，要求A中的元素服從均值為0，方差為1的正態(tài)分布，求sin(A）./A的極限。注：當A中的元素為0時MATLAB不發(fā)出警告邏輯運算符分別生成兩個包含10個 元素的行向量A、B，要求A、B中的元素在（0，1）之間服從均勻分布，，并將結果分別賦值給變量a、b，求變量a、b的同或（只利用與、或、非三個基本邏輯運算符）。多項式的表示方法及運算令m(x)=x6+6x5+20x4+50x3+75x2+84x+64；n(x)=x412x3＋25x+16 求m(x)與 n(x)的卷積。精度可控的數(shù)值積分已知y=x1log(1＋x1)關于x1求從0到1的積分，要求精度為1e6三、實驗結果關系運算符生成一個包含10個 元素的行向量A，要求A中的元素服從均值為0，方差為1的正態(tài)分布，求sin(A）./A的極限。MATAB源程序： A=randn(1,10)。A=A+(A==0)*eps。sin(A)./Aans =邏輯運算符分別生成兩個包含10個 元素的行向量A、B，要求A、B中的元素在（0，1）之間服從均勻分布，，并將結果分別賦值給變量a、b，求變量a、b的同或（只利用與、或、非三個基本邏輯運算符）。MATAB源程序： A=rand(1,10)B=rand(1,10)c=A。a=c.*A d=BB = a =00 0b =0000 0ans =0000000多項式的表示方法及運算令m(x)=x6+6x5+20x4+50x3+75x2+84x+64；n(x)=x412x3＋25x+16 求m(x)與 n(x)的卷積。MATAB源程序： m=[1 6 20 50 75 84 64]。n=[112 0 25 16]。w=conv(m,n)w =Columns 1 through 1052165626190733002944035900220Column 111024精度可控的數(shù)值積分已知y=x1log(1＋x1)關于x1求從0到1的積分，要求精度為1e6 MATAB源程序： y=quad(39。x1.*log(1+x1)39。,0,1,39。1e639。)y =三、實驗儀器PC機 MATLAB軟件四、實驗結果五、結論實驗五實驗名稱：數(shù)據(jù)和函數(shù)的可視化一、實驗目的掌握二維圖形的繪制熟悉二維圖形的修飾掌握三維圖形的繪制能從不同角度觀察試圖二、實驗步驟二維圖形的繪制Y=sinx2ex分別在四個子圖中給出棒狀圖（紅）、階梯圖（藍）、條形圖（黃）、填充圖形（綠），并加網絡線。二維圖形的修飾已知y1=sin(x)，y2=cos(x)，在一個圖形繪制y1，y2的曲線。⑴給出圖形題目（“正弦曲線和余弦曲線”）⑵給出橫坐標（“弧度值”），縱坐標（“函數(shù)值”）⑶分別加曲線標注（“正弦曲線”），（“余弦曲線”）⑷將x，y坐標軸的單位設置為相等 ⑸將當前圖形設置為正方形三維曲線的繪制f(x,y)=2sinx2+y2x+y22分別繪制其網線圖及曲面圖三維曲線的修飾及視角 對上題的曲面圖作⑴用截面顏色分布方式填充 ⑵用插列式顏色分布方式填充 ⑶用平面式顏色分布方式填充 ⑷用300方位角，－ ⑸用1800方位角，00仰角觀察三、實驗結果二維圖形的繪制Y=sinx2ex分別在四個子圖中給出棒狀圖（紅）、階梯圖（藍）、條形圖（黃）、填充圖形（綠），并加網絡線。MATAB源程序： x=0::2*pi。y=sin(x.^2).*exp(x)。figure(1)subplot(2,2,1)stem(x,y,39。r39。)gri