大連理工大學(xué)隨機(jī)信號實(shí)驗(yàn)報(bào)告材料完整

1、大連理工大學(xué)實(shí)驗(yàn)預(yù)習(xí)報(bào)告學(xué)院（系）：信息與通信工程學(xué)院專業(yè)：電子信息工程班級：姓名：學(xué)號：組：_實(shí)驗(yàn)時(shí)間：2015.12.14實(shí)驗(yàn)室：C221實(shí)驗(yàn)臺(tái)：指導(dǎo)教師：實(shí)驗(yàn) I ：隨機(jī)信號的產(chǎn)生、相關(guān)分析及其應(yīng)用實(shí)驗(yàn)實(shí)驗(yàn) 1均勻分布隨機(jī)數(shù)的產(chǎn)生，統(tǒng)計(jì)特性分析及計(jì)算機(jī)仿真一、實(shí)驗(yàn)?zāi)康暮鸵笳莆站鶆蚍植茧S機(jī)信號的基本產(chǎn)生方法二、實(shí)驗(yàn)原理和內(nèi)容較簡單的偽隨機(jī)序列產(chǎn)生方法是采用數(shù)論中基于數(shù)環(huán)理論的線性同余法（乘同余法、混合同余法），其迭代公式的一般形式為f(x) = (r ·x+bM，其離散形式為sn) Mod(+1)=rs nbMsn為n時(shí)刻的隨機(jī)數(shù)種子，r為擴(kuò)展因子，·( ) + Mo

2、d。其中，( )b 為固定擾動(dòng)項(xiàng)， M 為循環(huán)模， ModM 表示對 M 取模。為保證s( n) 的周期為 M，r的取值應(yīng)滿足 r= 4k+ 1 ，Mp， k 與 p 的選取應(yīng)滿足： r<M， r M31-1。s2， r， bM( -1)+1<2通常公式中參數(shù)常用取值為(0) =12357= 2045= 1，。=1048576三、實(shí)驗(yàn)步驟1. 編程實(shí)現(xiàn)產(chǎn)生 10000個(gè)在 (0, 1) 區(qū)間均勻分布隨機(jī)數(shù)。2. 計(jì)算生成隨機(jī)數(shù)的 14 階矩，最大值，最小值，頻度直方圖。實(shí)驗(yàn) 2高斯分布隨機(jī)數(shù)的產(chǎn)生，統(tǒng)計(jì)特性分析及計(jì)算機(jī)仿真一、實(shí)驗(yàn)?zāi)康暮鸵笳莆崭咚拱自肼暤幕井a(chǎn)生方法二、實(shí)驗(yàn)原理和

3、內(nèi)容1. 變換法2. 較簡單的高斯白噪聲產(chǎn)生方法是基于概率論中的中心極限定理。即無窮多個(gè)同分布隨機(jī)變量之和構(gòu)成隨機(jī)變量服從高斯分布。方便起見，可以使用多個(gè)均勻分布隨機(jī)變量之和近似高斯分布隨機(jī)變量。三、實(shí)驗(yàn)步驟1. 編程實(shí)現(xiàn)產(chǎn)生10000個(gè)N(3,4)高斯隨機(jī)數(shù)。2. 計(jì)算生成隨機(jī)數(shù)的 14 階矩，最大值，最小值，頻度直方圖。實(shí)驗(yàn) 3隨機(jī)信號相關(guān)函數(shù)計(jì)算、相關(guān)分析及計(jì)算機(jī)仿真一、實(shí)驗(yàn)?zāi)康暮鸵笳莆针S機(jī)信號相關(guān)函數(shù)計(jì)算、相關(guān)分析及實(shí)現(xiàn)二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)自相關(guān)和互相關(guān)的定義， 自相關(guān)，互相關(guān)計(jì)算隨機(jī)信號的自相關(guān)和互相關(guān)。三、實(shí)驗(yàn)步驟1. 產(chǎn)生高斯隨機(jī)信號。2. 計(jì)算其自相關(guān)函數(shù)。3. 計(jì)算兩個(gè)高

4、斯隨機(jī)信號的互相關(guān)函數(shù)。大連理工大學(xué)實(shí)驗(yàn)報(bào)告學(xué)院（系）：信息與通信工程學(xué)院專業(yè)：電子信息工程班級：電子 1303姓名：實(shí)驗(yàn)時(shí)間：指導(dǎo)教師：李彤學(xué)號 ：201383081組：_2015.12.14實(shí)驗(yàn)室：C221實(shí)驗(yàn)臺(tái)：實(shí)驗(yàn) I ：隨機(jī)信號的產(chǎn)生、相關(guān)分析及其應(yīng)用實(shí)驗(yàn)實(shí)驗(yàn) 1均勻分布白噪聲的生成一、實(shí)驗(yàn)?zāi)康暮鸵蠡诰鶆蚍植紓坞S機(jī)數(shù)，掌握均勻分布白噪聲典型生成方法。二、實(shí)驗(yàn)原理和內(nèi)容較簡單的偽隨機(jī)序列產(chǎn)生方法是采用數(shù)論中基于數(shù)環(huán)理論的線性同余法（乘同余法、混合同余法），其迭代公式的一般形式為f(x) = (r ·x+bModM，其離散形式為sn)(+1)=rs nbMsn為n時(shí)刻的隨機(jī)

5、數(shù)種子，r為擴(kuò)展因子，·( ) + Mod。其中，( )b 為固定擾動(dòng)項(xiàng)， M 為循環(huán)模， ModM 表示對 M 取模。為保證s( n) 的周期為 M，r的取值應(yīng)滿足r= 4k+ 1，M2p， k 與 p 的選取應(yīng)滿足：r<M， rM31-1。，r， b， M(-1)+1<2本實(shí)驗(yàn)中參數(shù)取值為s(1) =12357= 2025= 1。=1048576三、主要儀器設(shè)備微型計(jì)算機(jī)、 Matlab 開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法1. 編程實(shí)現(xiàn)產(chǎn)生 10000 個(gè)在 (0, 1) 區(qū)間均勻分布隨機(jī)數(shù)。2. 計(jì)算生成隨機(jī)數(shù)的 14 階矩，最大值，最小值，頻度直方圖五、實(shí)驗(yàn)數(shù)據(jù)記錄和

6、處理程序如下：M=1048576;b=1;r=2025;n(2)=n(2)+(s(i)-avr)2;%求 二階s=zeros(1,10000);矩s(1)=12357;n(3)=n(3)+(s(i)-avr)3;%求三階s=zeros(1,10000);n(4)=n(4)+(s(i)-avr)4;%求四階for i=2:10000ends(i)=mod(s(i-1)*r+b,M);m=m/10000;endn=n/10000;s=s/M;% 均勻分布隨機(jī)生成10000 數(shù)據(jù)disp('均值 =',num2str(m(1);figure,plot(s)%全部數(shù)據(jù)畫線disp(&#

7、39;二階原點(diǎn)矩title('全部數(shù)據(jù)連線 ')=',num2str(m(2);figure,plot(s,'.')%全部數(shù)據(jù)畫點(diǎn)disp('三階原點(diǎn)矩title('全部數(shù)據(jù)畫點(diǎn) ')=',num2str(m(3);%畫直方圖disp('四階原點(diǎn)矩hist(s,40)=',num2str(m(4);title('40個(gè)區(qū)間 ')disp('均值 =',num2str(n(1);sum=0;disp('二階中心矩for i=1:10000=',num2str(n

8、(2);sum=sum+s(i);% 求所有數(shù)的總和disp('三階中心矩end;=',num2str(n(3);avr=sum/10000;% 求所有數(shù)的平均數(shù)disp('四階中心矩m=zeros(1,4);=',num2str(n(4);for i=1:10000disp('方差 =',num2str(var(s);m(1)=m(1)+s(i); %求均值c=0;m(2)=m(2)+s(i)2;%求二階矩d=1;m(3)=m(3)+s(i)3;%求三階for i=1:10000%求最大值m(4)=m(4)+s(i)4;%求四階if c<

9、s(i)endc=s(i);n=zeros(1,4);end;for i=1:10000end;n(1)=n(1)+(s(i)-avr); %求均值for i=1:10000%求最小值if d>s(i)end;d=s(i);cend;d六、實(shí)驗(yàn)結(jié)果與分析運(yùn)行程序， command窗口中顯示的結(jié)果如下：均值 =0.49699二階原點(diǎn)矩 =0.32916三階原點(diǎn)矩 =0.24551四階原點(diǎn)矩 =0.19546均值 =-8.3666e-017二階中心矩 =0.08217三階中心矩 =0.00024392四階中心矩 =0.012193方差 =0.082179c =1.0000d =0實(shí)驗(yàn)得到圖表

10、如下：圖 1.1圖 1.2七、討論、建議、質(zhì)疑本實(shí)驗(yàn)中編寫了生成隨機(jī)序列的程序，通過設(shè)定不同的參數(shù)值可以得到不同的隨機(jī)序列，通過計(jì)算多階原點(diǎn)矩、中心矩和繪制直方圖可以更加清楚的看到產(chǎn)生的隨機(jī)序列的特點(diǎn)，對隨機(jī)序列的理解更加深刻。實(shí)驗(yàn) 2 高斯分布白噪聲的生成一、實(shí)驗(yàn)?zāi)康暮鸵蠡诰鶆蚍植紓坞S機(jī)數(shù)，掌握高斯分布白噪聲典型生成方法。二、實(shí)驗(yàn)原理和內(nèi)容1. 變換法2. 較簡單的高斯白噪聲產(chǎn)生方法是基于概率論中的中心極限定理。即無窮多個(gè)同分布隨機(jī)變量之和構(gòu)成隨機(jī)變量服從高斯分布。方便起見，可以使用多個(gè)均勻分布隨機(jī)變量之和近似高斯分布隨機(jī)變量。三、主要儀器設(shè)備微型計(jì)算機(jī)、 Matlab 開發(fā)環(huán)境四、實(shí)驗(yàn)

11、步驟與操作方法1. 編程實(shí)現(xiàn)產(chǎn)生10000個(gè)N(3,4)高斯隨機(jī)數(shù)。2. 計(jì)算生成隨機(jī)數(shù)的 14 階矩，最大值，最小值，頻度直方圖。五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理實(shí)驗(yàn)程序如下：m=3;n(3) = n(3) + (s(i)-avr)3;%求三n=4;階n1=sqrt(n);n(4) = n(4) + (s(i)-avr)4;%求四pi=3.1416;階s=zeros(1,10000);endfor i=1:10000m=m/10000;a=sqrt(-2*log(rand);n=n/10000;b=2*pi*rand;disp('均值 = ',num2str(m(1) );s(i)=n

12、1*a*cos(b)+m;%生成10000 個(gè)disp('二階原點(diǎn)矩=N(3, 4) 高斯隨機(jī)數(shù)',num2str(m(2) );enddisp('三階原點(diǎn)矩=figure',num2str(m(3) );plot(s)disp('四階原點(diǎn)矩=sum=0;',num2str(m(4) );for i=1:10000disp('一階中心矩 = ',num2str(n(1) );sum=sum+s(i);% 求所有數(shù)總數(shù)disp('二階中心矩=end;',num2str(n(2) );avr=sum/10000;% 求

13、平均數(shù)disp('三階中心矩=m = zeros(1,4);',num2str(n(3) );for i = 1 : 10000disp('四階中心矩=m(1) = m(1) + s(i);%求均值',num2str(n(4) );m(2) = m(2) + s(i)2;%求二階矩hist(s,100)m(3) = m(3) + s(i)3;%求三階title('100個(gè)區(qū)間 ')% 顯示頻率m(4) = m(4) + s(i)4;%求四階c=0;endd=1;n = zeros(1,4);for i=1:10000%求最大for i = 1 :

14、 10000if c<s(i)n(1) = n(1) + (s(i)-avr);%求一c=s(i);階矩end;n(2) = n(2) + (s(i)-avr)2;%求end;二階矩for i=1:10000%求最小if d>s(i)end;d=s(i);cend;d六、實(shí)驗(yàn)結(jié)果與分析運(yùn)行程序， command窗口中顯示的結(jié)果如下：均值 = 2.993二階原點(diǎn)矩 = 12.9888三階原點(diǎn)矩 = 62.522四階原點(diǎn)矩 = 341.0162一階中心矩 =-5.8249e-015二階中心矩 = 4.0306三階中心矩 = -0.48106四階中心矩 = 49.8857c =10.21

15、22d =-4.5628實(shí)驗(yàn)得到圖像如下：圖 1.3七、討論、建議、質(zhì)疑本實(shí)驗(yàn)和上一個(gè)實(shí)驗(yàn)類似， 但需要用到高斯函數(shù)的編寫方法，利用循環(huán)結(jié)構(gòu)可以求得若干個(gè)隨機(jī)數(shù)，通過直方圖可以對高斯分布有更直觀實(shí)際的了解。實(shí)驗(yàn) 3隨機(jī)信號相關(guān)函數(shù)估計(jì)一、實(shí)驗(yàn)?zāi)康暮鸵笳莆针S機(jī)信號相關(guān)函數(shù)計(jì)算、相關(guān)分析及實(shí)現(xiàn)二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)自相關(guān)和互相關(guān)的定義， 自相關(guān)，互相關(guān)計(jì)算隨機(jī)信號的自相關(guān)和互相關(guān)。三、主要儀器設(shè)備微型計(jì)算機(jī)、 Matlab 開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法4. 產(chǎn)生高斯隨機(jī)信號。5. 計(jì)算其自相關(guān)函數(shù)。6. 計(jì)算兩個(gè)高斯隨機(jī)信號的互相關(guān)函數(shù)。五、 實(shí)驗(yàn)數(shù)據(jù)記錄和處理實(shí)驗(yàn)程序如下：m=3;N(3,

16、4)高斯隨機(jī)數(shù) sn=4;endn1=sqrt(n);q=zeros(1,10000);pi=3.1416;for i=1:10000Fs=1000;c=sqrt(-2*log(rand);s=zeros(1,10000);d=2*pi*rand;for i=1:10000q(i)=n*c*cos(d)+m;%生成10000 個(gè)a=sqrt(-2*log(rand);N(3, 4)高斯隨機(jī)數(shù) qb=2*pi*rand;ends(i)=n1*a*cos(b)+m;%生成 10000 個(gè)figure(1);subplot(211),plot(s),title('s(n)');%co

17、r2 lag2=xcorr(q);畫出 s 的圖像 sfigure(3);plot(lag2/Fs,cor2),title('q的自相關(guān)subplot(212),plot(q),title('q(n)');%函數(shù) ');% 畫出 q 的自相關(guān)函數(shù)畫出 q 的圖像cor3 lag3=xcorr(s,q);cor1 lag1=xcorr(s);figure(4);figure(2);plot(lag3/Fs,cor3),title('s與 q 的互相關(guān)函數(shù) ');% 畫出 s 與 q 的互相關(guān)函數(shù)plot(lag1/Fs,cor1),title(&#

18、39;s的自相關(guān)函數(shù) ');% 畫出 s 的自相關(guān)函數(shù)六、實(shí)驗(yàn)結(jié)果與分析運(yùn)行程序，得到各圖像如下：圖 1.4 s和 q 的數(shù)據(jù)圖像圖 1.5 s的自相關(guān)函數(shù)圖 1.6 q的自相關(guān)函數(shù)圖 1.7 s與 q 的互相關(guān)函數(shù)七、討論、建議、質(zhì)疑本實(shí)驗(yàn)與實(shí)驗(yàn)2 相似，也要產(chǎn)生高斯隨機(jī)數(shù)，關(guān)鍵在于自相關(guān)程序以及互相關(guān)程序的編寫。通過對自相關(guān)以及互相關(guān)圖像的觀察，可以讓我們對于自相關(guān)以及互相關(guān)有更全面的了解。大連理工大學(xué)實(shí)驗(yàn)預(yù)習(xí)報(bào)告學(xué)院（系）：信息與通信工程學(xué)院專業(yè)：電子信息工程班級：電子 1303姓名：實(shí)驗(yàn)時(shí)間：指導(dǎo)教師：李彤學(xué)號 ：201383081組：_2015.12.15實(shí)驗(yàn)室：C221實(shí)驗(yàn)

19、臺(tái)：實(shí)驗(yàn) II ：系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析、功率譜分析及應(yīng)用實(shí)驗(yàn)實(shí)驗(yàn) 4 隨機(jī)信號的功率譜分析1一、實(shí)驗(yàn)?zāi)康暮鸵笳莆罩苯臃ü烙?jì)隨機(jī)信號功率譜。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)隨機(jī)信號的 Fourier變換結(jié)果，計(jì)算信號的功率譜。功率譜與頻譜關(guān)系：三、實(shí)驗(yàn)步驟1. 生成高斯白噪聲的隨機(jī)信號。2. 通過 Fourier 變換計(jì)算高斯白噪聲的功率譜，并繪圖。實(shí)驗(yàn) 5 隨機(jī)信號的功率譜分析方法2一、實(shí)驗(yàn)?zāi)康暮鸵笳莆臻g接法估計(jì)隨機(jī)信號功率譜。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)維納 - 辛欽定理，計(jì)算信號的功率譜。即三、實(shí)驗(yàn)步驟1. 生成高斯白噪聲的隨機(jī)信號。2. 計(jì)算高斯白噪聲的自相關(guān)函數(shù)。3. 通過計(jì)算高斯白噪

20、聲自相關(guān)函數(shù)的 Fourier 變換，得到噪聲功率譜并繪圖。實(shí)驗(yàn) 6 系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析及仿真1一、實(shí)驗(yàn)?zāi)康暮鸵笳莆諘r(shí)域上系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析及仿真實(shí)現(xiàn)。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)系統(tǒng)卷積性質(zhì)，計(jì)算系統(tǒng)輸出信號的統(tǒng)計(jì)特性。有如下性質(zhì)：三、實(shí)驗(yàn)步驟1. 生成均勻分布的隨機(jī)信號。2. 計(jì)算均勻分布隨機(jī)信號通過平均低通濾波器。3. 計(jì)算輸出信號均值、方差等統(tǒng)計(jì)特性。實(shí)驗(yàn) 7 系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析及仿真2一、實(shí)驗(yàn)?zāi)康暮鸵笳莆疹l域上系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析及仿真實(shí)現(xiàn)。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)卷積定理，計(jì)算系統(tǒng)輸出信號的統(tǒng)計(jì)特性。即：三、實(shí)驗(yàn)步驟1. 生成均勻分

21、布的隨機(jī)信號。2. 頻域上計(jì)算均勻分布隨機(jī)信號通過平均低通濾波器。3. 計(jì)算輸出信號均值、方差等統(tǒng)計(jì)特性。大連理工大學(xué)實(shí)驗(yàn)報(bào)告學(xué)院（系）：信息與通信工程學(xué)院專業(yè)：電子信息工程班級：電子 1303姓名：實(shí)驗(yàn)時(shí)間：指導(dǎo)教師：李彤學(xué)號 ：201383081組：_2015.12.15實(shí)驗(yàn)室：C221實(shí)驗(yàn)臺(tái)：實(shí)驗(yàn) II ：系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析、功率譜分析及應(yīng)用實(shí)驗(yàn)實(shí)驗(yàn) 4 隨機(jī)信號的功率譜分析1一、實(shí)驗(yàn)?zāi)康暮鸵笳莆罩苯臃ü烙?jì)隨機(jī)信號功率譜。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)維納 - 辛欽定理，計(jì)算信號的功率譜。即三、主要儀器設(shè)備微型計(jì)算機(jī)、 Matlab 開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法1. 生成高斯白

22、噪聲的隨機(jī)信號。2. 通過 Fourier 變換計(jì)算高斯白噪聲的功率譜，并繪圖。五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理程序如下：n=65536;+a(2)+randn(1,n);fs=2000;pi=3.14159;b1=fft(xn);t=(0:n-1)/fs;%b2=cos(2*pi/n)a=random('unif',0,1,1,2)*2*pi;b2=abs(b1).2/n;xn=cos(2*pi*30*t+a(1)+3*cos(2*pi*100*tb3=log(b2(1:n/2);f=(0:n/2-1)*fs/n;title('功率譜密度 ');figure,plot(

23、f,10*b3);六、實(shí)驗(yàn)結(jié)果與分析實(shí)驗(yàn)得到圖表如下：圖 2.1七、討論、建議、質(zhì)疑通過本實(shí)驗(yàn)，我掌握了直接估計(jì)隨機(jī)信號功率譜的方法。實(shí)驗(yàn) 5 隨機(jī)信號的功率譜分析2一、實(shí)驗(yàn)?zāi)康暮鸵笳莆臻g接法估計(jì)隨機(jī)信號功率譜。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)隨機(jī)信號的 Fourier變換結(jié)果，計(jì)算信號的功率譜。功率譜與頻譜關(guān)系：三、主要儀器設(shè)備微型計(jì)算機(jī)、 Matlab 開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法1. 生成高斯白噪聲的隨機(jī)信號。2. 計(jì)算高斯白噪聲的自相關(guān)函數(shù)。3. 通過計(jì)算高斯白噪聲自相關(guān)函數(shù)的 Fourier 變換，得到噪聲功率譜并繪圖。五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理程序如下：n=65536;b1=xcorr(xn

24、,'biased');fs=2000;b2=fft(b1);t=(0:n-1)/fs;b3=abs(b2);a=random('unif',0,1,1,2)*2*pi;f=(0:n-1)*fs/n/2;xn=cos(2*pi*30*t+a(1)+3*cos(2*pi*1figure,plot(f,10*log(b3(1:n);00*t+a(2)+randn(1,n);title('功率譜密度 ');六、實(shí)驗(yàn)結(jié)果與分析實(shí)驗(yàn)得到圖表如下：圖 2.2七、討論、建議、質(zhì)疑本實(shí)驗(yàn)與上一個(gè)實(shí)驗(yàn)結(jié)果相同，但是用的方法卻不同。通過本實(shí)驗(yàn)，我掌握了間接法估計(jì)隨機(jī)信

25、號功率譜。實(shí)驗(yàn) 6 系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析及仿真1一、實(shí)驗(yàn)?zāi)康暮鸵笳莆諘r(shí)域上系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析及仿真實(shí)現(xiàn)。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)系統(tǒng)卷積性質(zhì)，計(jì)算系統(tǒng)輸出信號的統(tǒng)計(jì)特性。有如下性質(zhì)：三、主要儀器設(shè)備微型計(jì)算機(jī)、 Matlab 開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法1. 生成均勻分布的隨機(jī)信號。2. 計(jì)算均勻分布隨機(jī)信號通過平均低通濾波器。3. 計(jì)算輸出信號均值、方差等統(tǒng)計(jì)特性。五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理程序如下：m=1048576;b=1;r=2045;x=zeros(1,10000); s(1)=12357;for i=2:10000s(i)=mod(s(i-1)*r+b,m)

26、;endmy=ms*sum(s); for i=1:10000h(i)=0.9i;endn1=1024;rs=zeros(1,1024);for m=0:(n1-1)for k=0:(n1-1)s=s/m;hist(s)r=zeros(1,1001);t=zeros(1,1001);ms=sum(s)/length(s);rs(m+1)=rx(m+1)+s(k+1)*s(k+m+1);endendp=0;rs=rs/1000;endfigureplot(rs(1:1000)ry(i+1)=ry(i+1)+rx(a+1)*h(i+1)*h(k+ry=zeros(1,100);1);for i=0

27、:99endfor j=0:99endfor k=0:99enda=i+j-k;ry=ry/200;if(a<0)figurea=-a;plot(ry(1:100)均值=2.4725e+03方差 =-6.1134e+06六、實(shí)驗(yàn)結(jié)果與分析實(shí)驗(yàn)得到圖表如下：圖 2.3圖 2.4圖 2.5七、討論、建議、質(zhì)疑通過本實(shí)驗(yàn)我對時(shí)域上系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析的了解更加深刻，并掌握了對隨機(jī)信號響應(yīng)的仿真。實(shí)驗(yàn) 7 系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析及仿真2一、實(shí)驗(yàn)?zāi)康暮鸵笳莆疹l域上系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析及仿真實(shí)現(xiàn)。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)卷積定理，計(jì)算系統(tǒng)輸出信號的統(tǒng)計(jì)特性。即：三、

28、主要儀器設(shè)備微型計(jì)算機(jī)、 Matlab 開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法1. 生成均勻分布的隨機(jī)信號。2. 頻域上計(jì)算均勻分布隨機(jī)信號通過平均低通濾波器。3. 計(jì)算輸出信號均值、方差等統(tǒng)計(jì)特性。五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理程序如下：n=500;w=(1:n)/n;xt=random('norm',0,1,1,n);t=(-n:n-1)/n*(n/20000);ht=fir1(500,0.3 0.4);subplot(4,1,1);plot(w,abs(sxx(1hw=fft(ht,2*n);:n);rxx=xcorr(xt,'biased');subplot(4,1,2

29、);plot(w,abs(hw2(1sxx=abs(fft(xt,2*n).2)/(2*n);:n);hw2=abs(hw).2;subplot(4,1,3);plot(w,abs(syy(1syy=sxx.*hw2;:n);ryy=fftshift(ifft(syy);subplot(4,1,4);plot(t,ryy);六、實(shí)驗(yàn)結(jié)果與分析實(shí)驗(yàn)得到圖表如下：圖 2.6七、討論、建議、質(zhì)疑通過本實(shí)驗(yàn)我對頻域上系統(tǒng)對隨機(jī)信號響應(yīng)的統(tǒng)計(jì)特性分析的了解更加深刻，并掌握了對隨機(jī)信號響應(yīng)的仿真。大連理工大學(xué)實(shí)驗(yàn)預(yù)習(xí)報(bào)告學(xué)院（系）：信息與通信工程學(xué)院專業(yè)：電子信息工程班級：電子 1303姓名：實(shí)驗(yàn)時(shí)間：指

30、導(dǎo)教師：李彤學(xué)號 ：201383081組：_2015.12.17實(shí)驗(yàn)室：C219實(shí)驗(yàn)臺(tái)：實(shí)驗(yàn) III ：通信信號統(tǒng)計(jì)特性分析、仿真實(shí)現(xiàn)及應(yīng)用實(shí)驗(yàn) 8 Hilbert變換算法實(shí)現(xiàn)1一、實(shí)驗(yàn)?zāi)康暮鸵笳莆誋ilbert變換的時(shí)域?qū)崿F(xiàn)。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)信號的卷積性質(zhì)，進(jìn)行Hilbert變換。Hilbert變換等效系統(tǒng)單位響應(yīng)為三、實(shí)驗(yàn)步驟生成正弦波信號，卷積得到Hilbert變換信號。生成余弦波信號，卷積得到Hilbert變換信號。實(shí)驗(yàn) 9 Hilbert變換算法實(shí)現(xiàn)2一、實(shí)驗(yàn)?zāi)康暮鸵笳莆誋ilbert變換的頻域?qū)崿F(xiàn)。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)信號的卷積定理，在頻域進(jìn)行Hilbert變換。 Hi

31、lbert變換等效系統(tǒng)函數(shù)為三、實(shí)驗(yàn)步驟生成正弦波信號，在頻域相乘得到Hilbert變換信號。生成余弦波信號，在頻域計(jì)算得到Hilbert變換信號。實(shí)驗(yàn)10窄帶平穩(wěn)隨機(jī)信號統(tǒng)計(jì)特性分析一、實(shí)驗(yàn)?zāi)康暮鸵笳莆諘r(shí)域窄帶平穩(wěn)隨機(jī)信號統(tǒng)計(jì)特性分析方法。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)窄帶平穩(wěn)隨機(jī)信號的特點(diǎn)，進(jìn)行信號分析。窄帶過程可以表示為 Ac(n)cos( 0n)-As(n)sin( 0n) ?？梢酝ㄟ^提取載波，相干解調(diào)得到x(n)=Ac(n) 和 As(n) ，進(jìn)一步求得包絡(luò)和相位三、實(shí)驗(yàn)步驟1. 按照實(shí)驗(yàn) 2中方法得到高斯白噪聲 xc(n) 和 xs(n) ，設(shè)定這兩個(gè)信號采樣頻率為10MHz。2. 用M

32、ATLAB的fdatool 函數(shù)生成采樣頻率為 10MHz，通帶頻率為 1KHz、阻帶頻率為10KHz、通帶波動(dòng)為 1dB、阻帶衰減為 60dB 的低通濾波器 h(n) 。3. 用低通濾波器 h1(n) 分別對 xc(n) 和 xs(n) 濾波得到 Ac(n) 和As(n) 。4. 用載波頻率為 0.2MHz的余弦、正弦波 ( 等效數(shù)字角頻率 0 = 0.04 ) 分別對 Ac(n) 和As(n) 調(diào)制，并組合為 x(n) 。5. 用載波頻率為 0.2MHz 的余弦、正弦波對 x(n) 進(jìn)行解調(diào)，并分別通過 h(n) 得到 Ac(n)和 As(n) 的估計(jì)值 Ac?(n) As?(n)。6.根

33、據(jù)，計(jì)算包絡(luò)、相位。7. 計(jì)算包絡(luò)分布和相位概率分布。大連理工大學(xué)實(shí)驗(yàn)報(bào)告學(xué)院（系）：信息與通信工程學(xué)院專業(yè)：電子信息工程班級：電子 1303姓名：實(shí)驗(yàn)時(shí)間：指導(dǎo)教師：李彤學(xué)號 ：201383081組：_2015.12.17實(shí)驗(yàn)室：C219實(shí)驗(yàn)臺(tái)：實(shí)驗(yàn) III ：通信信號統(tǒng)計(jì)特性分析、仿真實(shí)現(xiàn)及應(yīng)用實(shí)驗(yàn) 8 Hilbert變換算法實(shí)現(xiàn)1六、實(shí)驗(yàn)?zāi)康暮鸵笳莆誋ilbert變換的時(shí)域?qū)崿F(xiàn)。七、實(shí)驗(yàn)原理和內(nèi)容根據(jù)信號的卷積性質(zhì)，進(jìn)行Hilbert變換。Hilbert變換等效系統(tǒng)單位響應(yīng)為八、主要儀器設(shè)備微型計(jì)算機(jī)、 Matlab開發(fā)環(huán)境九、實(shí)驗(yàn)步驟與操作方法生成正弦波信號，卷積得到Hilbert變

34、換信號。生成余弦波信號，卷積得到Hilbert變換信號。十、實(shí)驗(yàn)數(shù)據(jù)記錄和處理程序如下：h=zeros(1,15);x=cos(pi/10*n);g=zeros(1,15);for i=1:15;n=1:400;h(i)=1/(i*pi)-(-1)*i*1/(i*pi);endg(j)=1/(j*pi)-(-1)*j*1/(j*pi);y=conv(x,h);endsubplot(2,1,1)b=conv(a,g);plot(x);subplot(2,1,1)title('余弦信號 ');figuresubplot(2,1,2)plot(a);plot(y);title(

35、9;正弦信號 ');title('余弦 Hilbert變換 ');subplot(2,1,2)a=sin(pi/10*n);plot(b);for j=1:15;title('正弦 Hilbert變換 ');十一、實(shí)驗(yàn)結(jié)果與分析生成的余弦信號及其Hilbert變換信號如下圖：圖 3.1余弦信號及其Hilbert變換生成的正弦信號及其Hilbert變換信號如下圖：圖 3.2 正弦信號及其Hilbert變換七、討論、建議、質(zhì)疑本次實(shí)驗(yàn)加深了我對希爾伯特變換的理解，并且編程實(shí)現(xiàn)了正弦和余弦信號希爾伯特變換的時(shí)域?qū)崿F(xiàn)。實(shí)驗(yàn) 9 Hilbert變換算法實(shí)現(xiàn)2一、實(shí)

36、驗(yàn)?zāi)康暮鸵笳莆誋ilbert變換的頻域?qū)崿F(xiàn)。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)信號的卷積定理，在頻域進(jìn)行Hilbert變換。 Hilbert變換等效系統(tǒng)函數(shù)為三、主要儀器設(shè)備微型計(jì)算機(jī)、 Matlab開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法生成正弦波信號，在頻域相乘得到Hilbert變換信號。生成余弦波信號，在頻域計(jì)算得到Hilbert變換信號。五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理程序如下：Y=zeros(1,400);H=zeros(1,400);n=1:400;a=sin(pi/10*n);x=cos(pi/10*n);subplot(2,1,1)subplot(2,1,1)figure,plot(a);plot(x);t

37、itle('正弦信號 ');title('正弦 Hilbert變換 ');G=fft(a,400);X=fft(x,400);for k=2:200for k=2:200H(k)=-1i*G(k);Y(k)=-1i*X(k);endendfor k=202:400for k=202:400H(k)=1i*G(k);Y(k)=1i*X(k);endendH(1)=G(1);Y(1)=X(1);H(201)=G(201);Y(201)=X(201);b=ifft(H);y=ifft(Y);subplot(2,1,2)subplot(2,1,2)plot(b);plot(y);title('正弦 Hilbert變換 ');title('余弦 Hilbert變換 ');六、實(shí)驗(yàn)結(jié)果與分析生成的余弦信號及其Hilbert變換信號如下圖：圖 3.3余弦信號及其Hilbert變換生成的正弦信號及其Hilbert變換信號如下圖：圖 3.4正弦信號及其Hilbert變換七、討論、建議、質(zhì)疑本實(shí)驗(yàn)與上一個(gè)實(shí)驗(yàn)類似，也是實(shí)現(xiàn)正弦和余弦信號的