大連理工大學(xué)第二次隨機(jī)信號(hào)實(shí)驗(yàn)報(bào)告

大連理工大學(xué)實(shí)驗(yàn)報(bào)告學(xué)院（系）：電子信息與電氣工程專業(yè)：電子信息工程班級(jí)：姓名：學(xué)號(hào)：組：實(shí)驗(yàn)時(shí)間：實(shí)驗(yàn)室：實(shí)驗(yàn)臺(tái)：指導(dǎo)教師簽字：成績(jī)：實(shí)驗(yàn)二7／8／9／10實(shí)驗(yàn)7一、實(shí)驗(yàn)?zāi)康暮鸵笳莆罩苯臃ü烙?jì)隨機(jī)信號(hào)功率譜。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)隨機(jī)信號(hào)的Fourier變換結(jié)果，計(jì)算信號(hào)的功率譜。功率譜與頻譜關(guān)系：三、主要儀器設(shè)備微型計(jì)算機(jī)，matlab開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法(1)生成高斯白噪聲的隨機(jī)信號(hào)。(2)通過Fourier變換計(jì)算高斯白噪聲的功率譜，并繪圖。代碼： 五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理clcclearx1=rand(1,1000);x2=rand(1,1000);y=sqrt((-2)*log(x1)).*cos(2*pi*x2);N=length(y);n=0:(N-1);k=0:(N-1);z=y*exp(-1j*2*pi/N).^(n'*k);G=z.*conj(z)/(2*N);plot(G);六、實(shí)驗(yàn)結(jié)果與分析見實(shí)驗(yàn)數(shù)據(jù)和處理。實(shí)驗(yàn)8一、實(shí)驗(yàn)?zāi)康暮鸵笳莆臻g接法估計(jì)隨機(jī)信號(hào)功率譜。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)維納-辛欽定理，計(jì)算信號(hào)的功率譜。即：三、主要儀器設(shè)備微型計(jì)算機(jī)，matlab開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法(1)生成高斯白噪聲的隨機(jī)信號(hào)。(2)計(jì)算高斯白噪聲的自相關(guān)函數(shù)。(3)通過計(jì)算高斯白噪聲自相關(guān)函數(shù)的Fourier變換，得到噪聲功率譜并繪圖。代碼：clcclearx1=rand(1,1000);x2=rand(1,1000);r=zeros(1,1000);y=sqrt((-2)*log(x1)).*cos(2*pi*x2);figurehist(y);form=1:500fork=1:1000-mr(m)=r(m)+y(k)*y(k+m-1);endendr=r./500;form=1:500r(500+m)=r(m);endform=1:499r(501-m)=r(501+m);endr(1)=0;g=fft(r);G=abs(g);figureplot(G);實(shí)驗(yàn)數(shù)據(jù)記錄和處理六、實(shí)驗(yàn)結(jié)果與分析見實(shí)驗(yàn)數(shù)據(jù)和處理。實(shí)驗(yàn)9實(shí)驗(yàn)?zāi)康暮鸵笳莆諘r(shí)域上系統(tǒng)對(duì)隨機(jī)信號(hào)響應(yīng)的統(tǒng)計(jì)特性分析及仿真實(shí)現(xiàn)。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)系統(tǒng)卷積性質(zhì)，計(jì)算系統(tǒng)輸出信號(hào)的統(tǒng)計(jì)特性。有如下性質(zhì)：三、主要儀器設(shè)備微型計(jì)算機(jī)，matlab開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法(1)生成均勻分布的隨機(jī)信號(hào)。(2)計(jì)算均勻分布隨機(jī)信號(hào)通過平均低通濾波器。(3)計(jì)算輸出信號(hào)均值、方差等統(tǒng)計(jì)特性。代碼： clcclearx1=rand(1,1000);x2=rand(1,1000);y=sqrt((-2)*log(x1)).*cos(2*pi*x2);t=0:0.001:0.999;h=0.9.^t;Y=conv(y,h);plot(Y);my=0;fori=1:1000my=my+y(i);endmy=my/1000;mY=0;fori=1:1000mY=mY+h(i);endmY=my*mYRx=zeros(1,1000);forj=0:999fori=1:1000-abs(j-500)Rx(j+1)=Rx(j+1)+y(i)*y(i+abs(j-500));endRx(j+1)=Rx(j+1)/1000;endfigure(2);plot(Rx);Ry=zeros(1,1000);forj=1:1000fori=j:1000fork=i:1000Ry(j)=Ry(j)+h(k)*h(i)*Rx(j+k-i);endendendRy=Ry/1000;figure(3);plot(Ry);Gy=Ry(1)-mY^2五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理x自相關(guān)函數(shù)y自相關(guān)函數(shù)六、實(shí)驗(yàn)結(jié)果與分析見實(shí)驗(yàn)數(shù)據(jù)和處理。實(shí)驗(yàn)10實(shí)驗(yàn)?zāi)康暮鸵笳莆疹l域上系統(tǒng)對(duì)隨機(jī)信號(hào)響應(yīng)的統(tǒng)計(jì)特性分析及仿真實(shí)現(xiàn)。二、實(shí)驗(yàn)原理和內(nèi)容根據(jù)卷積定理，計(jì)算系統(tǒng)輸出信號(hào)的統(tǒng)計(jì)特性。即：三、主要儀器設(shè)備微型計(jì)算機(jī)，matlab開發(fā)環(huán)境四、實(shí)驗(yàn)步驟與操作方法(1)生成均勻分布的隨機(jī)信號(hào)。(2)頻域上計(jì)算均勻分布隨機(jī)信號(hào)通過平均低通濾波器。(3)計(jì)算輸出信號(hào)均值、方差等統(tǒng)計(jì)特性。代碼： clcclearx1=rand(1,1000);x2=rand(1,1000);y=sqrt((-2)*log(x1)).*cos(2*pi*x2);G=zeros(1,1000);Y=fft(y);fork=1:1000G(k)=1/1000*(abs(Y(k)))^2;endh=zeros(1,1000);forn=1:1000h(n)=0.9^n;endHH=fft(h);H=abs(HH).^2;Gy=G.*H;figure;plot(Gy);disp(['方差=',num2str(G)]);五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理功率譜密度六、實(shí)驗(yàn)結(jié)果與分