窄帶隨機過程的產生

2021-10-01 15:25:14 字數 2153 閱讀 9677

1、窄帶隨機過程的產生過程

隨機過程的產生過程如圖1所示,在利用萊斯表示式生成窄帶隨機過程時,需要通過高斯白雜訊生成兩個低頻的隨機過程,即首先使高斯白雜訊通過低通濾波器(本實驗中的使用的低通濾波器為巴特沃斯低通濾波器)得到兩個低頻隨機過程,接著讓兩個隨機過程分別通過乘法器,最後再通過累加器即可得到所需要的窄帶隨機過程。

2、**程式

(1)主程式

clear all

clcfs=10000;

t=1/fs;

t=0:t:1;

n=length(t);

%-----------------------------------------生成低頻隨機過程

fa=0.005; %低頻過程a(t)截止頻率

wa=2*pi*fa;

fb=0.0055; %低頻過程b(t)截止頻率

wb=2*pi*fb;

f0=2000; %中心頻率

at=lowfrequency(n,wa); %低頻隨機過程

figure

subplot(2,1,1),plot(t,at)

title('低頻過程a'),xlabel('t'),ylabel('b(t)')

bt=lowfrequency(n,wb);

subplot(2,1,2),plot(t,bt)

title('低頻過程b'),xlabel('t'),ylabel('b(t)')

%-----------------------------------------窄帶隨機過程及性質

x=at.*cos(2*pi*f0*t)-bt.*sin(2*pi*f0*t); %窄帶隨機過程

figure,plot(t,x),title('窄帶隨機過程')

rtau=xcorr(x); %自相關函式

tt=-n+1:n-1;

figure,plot(tt,rtau),title('自相關函式r_x(\tau)')

sx=fft(rtau); %功率譜密度

len=length(sx);

k=0:len-1;

w=2*pi*(k/len-1/2)*fs;

figure,plot(w/2/pi,abs(fftshift(sx)));xlim([-2300,2300])

title('功率譜密度s_x(\omega)'),xlabel('f/hz')

(2)產生低頻隨機過程的函式

function [m_t] = lowfrequency(n,omega)

%lowfrequency 產生低頻高斯隨機訊號

% t-時間,w-頻率,m_t -低頻雜訊

x=wgn(1,n,5); %產生高斯白雜訊

[b,a]=butter(10,omega,'low');

m_t=filter(b,a,x);

end

實驗所得窄帶隨機過程如圖3所示,窄帶隨機過程的性質如圖4、圖5所示,在本實驗中窄帶隨機過程的中心頻率為

由圖5可知窄帶隨機過程的功率譜密度中心頻率

邊界頻率分別為

滿足

隨機產生漢字

region 產生n個隨機的漢字 產生n個隨機的漢字 產生的隨機漢字個數 public static string getrandomchinese intstrlength 此函式在漢字編碼範圍內隨機建立含兩個元素的十六進製制位元組陣列,每個位元組陣列代表乙個漢字,並將 四個位元組陣列儲存在obj...

隨機數的產生

方法 標頭檔案中的rand 函式可以產生乙個閉區間 0,rand max 的均勻隨機數,其中rand max在中被定義為常數32767,嚴格地說,這裡的隨機數是 偽隨機數 因為它是由數學公式計算出來的,不過在演算法領域,多數情況下可以吧她當作真正地隨機數。示例 產生n個區間 1,m 上的隨機數。in...

隨機數的產生

1 隨機數的使用 1 arc4random 比較精確不需要生成隨即種子 使用方法 通過arc4random 獲取0到x 1之間的整數的 如下 int value arc4random x 獲取1到x之間的整數的 如下 int value arc4random x 1 2 ccrandom 0 1 c...