如何在matlab裡面向訊號裡面加入白噪聲

2021-03-27 05:29:07 字數 6203 閱讀 9660

1樓:go風影

y=awgn(x,snr),加性高斯白du噪函式,x原訊號,zhisnr是訊雜比,不知道你dao這個2%是能量

專還是幅屬值,如果是能量的話,snr=10*log(100/2);如果是幅值,snr=20*log(100/2);

**如下:

clc; clear ;

t=0:0.01:1;

y1=(pi/2)*cos(2*pi*t)+(pi/64)*cos(pi*t/8+pi/3);

snr=10*log(100/2); % 2%是能量比x=awgn(y1,snr);

subplot(211), plot(t,y1);

title('原訊號');

% 加入噪聲後有毛刺,但2%的噪聲有點小,毛刺不明顯。

subplot(212), plot(t,x);

title('加高斯白噪訊號');

2樓:匿名使用者

在matlab中無論是wgn還是baiawgn函式,實質du都是由randn函式產生

zhi的噪聲。即,wgn函式中調

dao用了randn函式,而版awgn函式中呼叫了wgn函式。

權根據awgn的實現**可以知道「向已知訊號新增某個訊雜比(snr)的高斯白噪聲」,

即:awgn(x,snr,』measured』,'linear』),命令的作用是對原訊號x新增訊雜比(比值)為snr的噪聲,在新增之前先估計訊號x的強度。

直接對原始訊號新增噪聲:

y=x+rand(length(x),1)y=x+randn(length(x),1))

3樓:匿名使用者

t=0:0.001:1;

y1=(pi/2)*cos(2*pi*t)+(pi/64)*cos(pi*t/8+pi/3);

temp=y1+0.02*y1.*rand(size(y1));

plot(t,y1,t,temp,'r')

4樓:日向淳正

l=length(y1);

noise=rand(1,l); %生成跟你的訊號一樣長度的白噪聲.

y=y1+noise; %加入白噪聲

matlab中如何給我的一段訊號加高斯白噪聲

5樓:匿名使用者

你這個是訊號波形資料嗎? 加噪聲一般是對模擬訊號加噪聲matlab有一個函式很方便加入高斯白噪聲x是訊號序列, y = awgn(x,snr) 在訊號x中加入高斯白噪聲,訊雜比snr以db為單位。x的強度假定為0dbw。

比如:x = sqrt(2)*sin(0:pi/8:6*pi);

y = awgn(x,10);

plot(x,'b');hold on; plot(y,'r')

matlab中給訊號加入白噪聲,如何利用傅立葉變換去噪

6樓:我是舞三

總體思路:設定篩選條件,將陣列中符合條件的元素篩選出來(可以篩選出元素值也可以篩選出元素對應的下標),然後將元素值賦值給另外一個陣列的元素

參考**:

clcclear all

close all

x = randi(10, 1, 5); % 原始資料

disp('x=');

disp(x);

y1 = -1 * ones(1, length(x)); % 儲存篩選出的結果,並初始化成不可能出現的結果

%-- 篩選符合條件的元素下標,並將對應元素存到y1 --%

ix = find(x>=3); % 例如,篩選出原始資料中大於3的元素的下標

if(isempty(ix))

disp('x的元素均不符合條件');

else

y1(1:length(ix)) = x(ix); % 將符合條件的元素存到y1中

y1( find(y1==-1) ) = ; % 剩餘元素刪除

disp('y1=');

disp(y1);

y2 = -1 * ones(1, length(x)); % 儲存篩選出的結果,並初始化成不可能出現的結果

%-- 篩選符合條件的元素值,並存入y2 --%

for i=1:length(x)

if(x(i)>=3) % 逐個檢測x的元素是否符合條件,把符合條件的元素值賦值給y2

y2(i) = x(i);

endend

y2( find(y2==-1) ) = ; % 剩餘元素刪除

disp('y2=');

disp(y2);

%----------------------------------%

輸出結果

x=7 1 9 10 7

y1=7 9 10 7

y2=7 9 10 7

7樓:我想呆會兒

小波的閾值去噪效果比較好。可以呼叫wden函式,比較簡單

本題中前面改一下比較直觀,可以把原始訊號輸出來subplot(2,2,1)後面去噪訊號subplot(2,2,2)頻譜圖subplot(2,2,3)後面可加程式

xd=wden(x,'rigrsure','s','sln',5,'sym8');

subplot(2,2,4);plot(xd);

title('去噪訊號');

err=x-xd;

函式wden的呼叫格式有以下兩種:

(1)[xd,cxd,lxd]=wden(x,tptr,sorh,scal,n,'wname')

(2)[xd,cxd,lxd]=wden(c,l,tptr,sorh,scal,n,'wname')

函式wden用於一維訊號的自動消噪。x為原始訊號,[c,l]為訊號的小波分解,n為小波分解的層數。

thr為閾值選擇規則:

*tptr='rigrsure',自適應閾值選擇使用stein的無偏風險估計原理。

*tptr='heursure',使用啟發式閾值選擇。

*tptr='sqtwolog',閾值等於sqrt(2*log(length(x))).

*tptr='minimaxi',用極大極小原理選擇閾值。

sorh是軟閾值或硬閾值的選擇(分別對應's'和'h')。

scal指所使用的閾值是否需要重新調整,包含下面三種:

*scal='one' 不調整;

*scal='sln' 根據第一層的係數進行噪聲層的估計來調整閾值。

*scal='mln' 根據不同的噪聲估計來調整閾值。

xd為消噪後的訊號,[cxd,lxd]為消噪後訊號的小波分解結構。

8樓:匿名使用者

山科同跟張老師上小波分析的路過

怎樣用matlab產生白噪聲訊號?

9樓:匿名使用者

1. wgn:產生高斯白噪聲

y = wgn(m,n,p) 產生一個m行n列的高斯白噪聲的矩陣,p以dbw為單位指定輸出噪聲的強度。

y = wgn(m,n,p,imp) 以歐姆(ohm)為單位指定負載阻抗。

y = wgn(m,n,p,imp,state) 重置randn的狀態。

在數值變數後還可附加一些標誌性引數:

y = wgn(…,powertype) 指定p的單位。powertype可以是'dbw', 'dbm'或'linear'。線性強度(linear power)以瓦特(watt)為單位。

y = wgn(…,outputtype) 指定輸出型別。outputtype可以是'real'或'***plex'。

2. awgn:在某一訊號中加入高斯白噪聲

y = awgn(x,snr) 在訊號x中加入高斯白噪聲。訊雜比snr以db為單位。x的強度假定為0dbw。如果x是複數,就加入復噪聲。

y = awgn(x,snr,sigpower) 如果sigpower是數值,則其代表以dbw為單位的訊號強度;如果sigpower為'measured',則函式將在加入噪聲之前測定訊號強度。

y = awgn(x,snr,sigpower,state) 重置randn的狀態。

y = awgn(…,powertype) 指定snr和sigpower的單位。powertype可以是'db'或'linear'。如果powertype是'db',那麼snr以db為單位,而sigpower以dbw為單位。

如果powertype是'linear',那麼snr作為比值來度量,而sigpower以瓦特為單位

10樓:匿名使用者

xi=sqrt(0.1)*randn(100,1);

matlab 給訊號加高斯白噪聲

11樓:匿名使用者

一般是這樣的

如果一個訊號為x (向量),產生一個均值為0,方差為a的高斯白噪聲,noise = a*randn(size(x));

則噪聲訊號為

y = x + noise;

12樓:匿名使用者

給訊號x加白噪聲:

y=awgn(x,snr);函式

請問如何在matlab裡在一個訊號曲線上加上1到3個高斯白噪聲?謝謝

13樓:匿名使用者

matlab中產生高斯白

噪聲的兩個函式

matlab中產生高斯白噪聲非常方便,可以直接應用兩個函式,一個是wgn,另一個是awgn。wgn用於產生高斯白噪聲,awgn則用於在某一訊號中加入高斯白噪聲。

1. wgn:產生高斯白噪聲

y = wgn(m,n,p) 產生一個m行n列的高斯白噪聲的矩陣,p以dbw為單位指定輸出噪聲的強度。

y = wgn(m,n,p,imp) 以歐姆(ohm)為單位指定負載阻抗。

y = wgn(m,n,p,imp,state) 重置randn的狀態。

在數值變數後還可附加一些標誌性引數:

y = wgn(…,powertype) 指定p的單位。powertype可以是'dbw', 'dbm'或'linear'。線性強度(linear power)以瓦特(watt)為單位。

y = wgn(…,outputtype) 指定輸出型別。outputtype可以是'real'或'***plex'。

2. awgn:在某一訊號中加入高斯白噪聲

y = awgn(x,snr) 在訊號x中加入高斯白噪聲。訊雜比snr以db為單位。x的強度假定為0dbw。如果x是複數,就加入復噪聲。

y = awgn(x,snr,sigpower) 如果sigpower是數值,則其代表以dbw為單位的訊號強度;如果sigpower為'measured',則函式將在加入噪聲之前測定訊號強度。

y = awgn(x,snr,sigpower,state) 重置randn的狀態。

y = awgn(…,powertype) 指定snr和sigpower的單位。powertype可以是'db'或'linear'。如果powertype是'db',那麼snr以db為單位,而sigpower以dbw為單位。

如果powertype是'linear',那麼snr作為比值來度量,而sigpower以瓦特為單位。

註釋 1. 分貝(decibel, db):分貝(db)是表示相對功率或幅度電平的標準單位,換句話說,就是我們用來表示兩個能量之間的差別的一種表示單位,它不是一個絕對單位。

例如,電子系統中將電壓、電流、功率等物理量的強弱通稱為電平,電平的單位通常就以分貝表示,即事先取一個電壓或電流作為參考值(0db),用待表示的量與參考值之比取對數,再乘以20作為電平的分貝數(功率的電平值改乘10)。

2. 分貝瓦(dbw, db watt):指以1w的輸出功率為基準時,用分貝來測量的功率放大器的功率值。

3. dbm (db-milliwatt):即與1milliwatt(毫瓦)作比較得出的數字。

0 dbm = 1 mw

10 dbm = 10 mw

20 dbm = 100 mw

也可直接用randn函式產生高斯分佈序列,例如:

y=randn(1,2500);

y=y/std(y);

y=y-mean(y);

a=0.0128;

b=sqrt(0.9596);

y=a+b*y;

就得到了 n ( 0.0128, 0.9596 ) 的高斯分佈序列。

如何在matlab中呼叫matlab中已有的函式

呼叫函式檔案 myfile.m clear clcfor t 1 10 y mylfg t end 自定義函式檔案 mylfg.m function y mylfg x 注意 函式名 mylfg 必須與檔名 mylfg.m 一致 y x 1 3 注 這種方法要求自定義函式必須單獨寫一個m檔案,不能與...

如何在matlab中的labelylabellegend和tet函

1 首先設定畫圖的引數。close all 語句清空之前的圖,set 語句將背景設定為白色,定義一個fontsize變數方便以後使用。2 畫圖之前先用figure語句產生一個空圖,這樣每次生成一幅圖都會自動順序編號。否則用figure n 指定生成圖的編號或名字也可以。3 用plot函式在同一圖中畫...

如何在matlab中運用ezplot顯示多個函式的影象

plot可以用一個函式copy繪製多個影象,ezplot可以使用hold on 畫完一個影象後用hold on 然後畫另一個。sym x ezplot sin x hold on ezplot cos x hold on命令不行嗎?在matlab中如何輸出一個函式的影象 如果要畫y x 3 x 2 ...