【文章內(nèi)容簡介】 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。%產(chǎn)生16點DFT對應(yīng)的采樣點頻率（以零頻率為中心）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。%產(chǎn)生16點DFT對應(yīng)的采樣點頻率（以零頻率為中心）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èi)容N=82018161412|X(f)|10864200123Hz456圖 實驗內(nèi)容1采樣點數(shù)N=8時的頻譜N=612108|X(f)|64200123Hz456圖 實驗內(nèi)容1采樣點數(shù)N=6時的頻譜N=3280706050|X(f)|40302010001234Hz5678圖 實驗內(nèi)容1采樣點數(shù)N=32時的頻譜N=21151050200幅度譜150100500HzN=1*** 實驗內(nèi)容2頻譜圖（1）實驗內(nèi)容1，根據(jù)給定模擬f1=2Hz，f=4Hz,所以滿足采樣定理2fs179。2fh，采取樣頻率fs=16Hz,F=42=2Hz,根據(jù)N179。`2f0/F,所以N179。8，當(dāng)N=6時，不滿足采樣點數(shù)，所以頻譜在時域混疊，而當(dāng)N=32時，滿足采樣點數(shù),頻譜中能夠看到有兩個頻率成分，在f1=2Hz,f2=4Hz處有譜線，。(2)實驗內(nèi)容2,f1=120Hz,f2=140Hz,取采樣頻率fs=400Hz,F=140120=20Hz,所以N=20，顯然，當(dāng)N=21時，能把f1和f2兩個譜峰分辨出來，當(dāng)N=11時，兩個譜峰分辨不出來。16點|DFT[x(nT)]|10幅度503020100f(Hz)32點|DFT[x(nT)]|102030幅度1003020100f(Hz)64點|DFT[x(nT)]|102030幅度2003020100f(Hz) 實驗內(nèi)容3頻譜圖（3）實驗內(nèi)容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倍，這種結(jié)果正好驗證了用DFT對中期序列譜分析的理論。（1）完成各個實驗任務(wù)和要求，附上程序清單和有關(guān)曲線。（2）簡要回答思考問題第二篇：MATLAB仿真技術(shù)實驗教案《MATLAB仿真技術(shù)》 實驗教案實驗一實驗名稱：熟悉Matlab交互工作界面一、實驗?zāi)康氖煜atlab各種工作界面的操作要旨掌握Matlab的基本操作命令二、實驗步驟命令窗口（1）體驗命令窗口的菜單及各項功能（2）嘗試命令窗口編輯特殊功能鍵和設(shè)置工作空間窗口與當(dāng)前路徑窗口（1）在工作空間窗口查看及修改變量（2）添加新的路徑為Matlab路徑圖形窗口和文本編輯窗口（1）練習(xí)圖形窗口中修改圖形的方法（2）在文本編輯窗口調(diào)試程序體會Matlab的基本操作命令三、實驗儀器PC機 MATLAB軟件四、實驗結(jié)果五、結(jié)論實驗二實驗名稱：Matlab在符號計算方面的應(yīng)用一、實驗?zāi)康恼莆諛?biāo)識符的生成和使用掌握矩陣及變量的賦值熟悉三類運算符及其功能二、實驗內(nèi)容標(biāo)識符的生成和使用1）、計算y=x+()/(x+)5(x+1/x)，當(dāng)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矩陣魔方陣和相應(yīng)的隨機矩陣，將兩個矩陣并接起來，然后提取任意兩個列向量。 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）設(shè)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軟件四、實驗結(jié)果五、結(jié)論實驗三實驗名稱：一維二維數(shù)組創(chuàng)建與尋訪一、實驗?zāi)康恼莆諗?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ù)組 保存并調(diào)用變量a,b比較數(shù)組運算與矩陣運算操作的區(qū)別 ⑴分別對變量a,b 進行矩陣乘和數(shù)組乘 ⑵分別對變量a,b 進行矩陣除和數(shù)組除三、實驗結(jié)果直接輸入法創(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ù)組 保存并調(diào)用變量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軟件四、實驗結(jié)果五、結(jié)論實驗四實驗名稱：Matlab在數(shù)值計算方面的應(yīng)用一、實驗?zāi)康膶W(xué)會使用MATLAB關(guān)系運算符學(xué)會使用MATLAB邏輯運算符掌握多項式的表示方法及運算熟練掌握精度可控的數(shù)值積分二、實驗步驟關(guān)系運算符生成一個包含10個 元素的行向量A，要求A中的元素服從均值為0，方差為1的正態(tài)分布，求sin(A）./A的極限。注：當(dāng)A中的元素為0時MATLAB不發(fā)出警告邏輯運算符分別生成兩個包含10個 元素的行向量A、B，要求A、B中的元素在（0，1）之間服從均勻分布，，并將結(jié)果分別賦值給變量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)關(guān)于x1求從0到1的積分，要求精度為1e6三、實驗結(jié)果關(guān)系運算符生成一個包含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）之間服從均勻分布，，并將結(jié)果分別賦值給變量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)關(guān)于x1求從0到1的積分，要求精度為1e6 MATAB源程序： y=quad(39。x1.*log(1+x1)39。,0,1,39。1e639。)y =三、實驗儀器PC機 MATLAB軟件四、實驗結(jié)果五、結(jié)論實驗五實驗名稱：數(shù)據(jù)和函數(shù)的可視化一、實驗?zāi)康恼莆斩S圖形的繪制熟悉二維圖形的修飾掌握三維圖形的繪制能從不同角度觀察試圖二、實驗步驟二維圖形的繪制Y=sinx2ex分別在四個子圖中給出棒狀圖（紅）、階梯圖（藍(lán)）、條形圖（黃）、填充圖形（綠），并加網(wǎng)絡(luò)線。二維圖形的修飾已知y1=sin(x)，y2=cos(x)，在一個圖形繪制y1，y2的曲線。⑴給出圖形題目（“正弦曲線和余弦曲線”）⑵給出橫坐標(biāo)（“弧度值”），縱坐標(biāo)（“函數(shù)值”）⑶分別加曲線標(biāo)注（“正弦曲線”），（“余弦曲線”）⑷將x，y坐標(biāo)軸的單位設(shè)置為相等 ⑸將當(dāng)前圖形設(shè)置為正方形三維曲線的繪制f(x,y)=2sinx2+y2x+y22分別繪制其網(wǎng)線圖及曲面圖三維曲線的修飾及視角 對上題的曲面圖作⑴用截面顏色分布方式填充 ⑵用插列式顏色分布方式填充 ⑶用平面式顏色分布方式填充 ⑷用300方位角，－ ⑸用1800方位角，00仰角觀察三、實驗結(jié)果二維圖形的繪制Y=sinx2ex分別在四個子圖中給出棒狀圖（紅）、階梯圖（藍(lán)）、條形圖（黃）、填充圖形（綠），并加網(wǎng)絡(luò)線。MATAB源程序： x=0::2*pi。y=sin(x.^2).*exp(x)。figure(1)subplot(2,2,1)stem(x,y,39。r39。)gri