-
大小: 12KB文件類型: .m金幣: 1下載: 0 次發(fā)布日期: 2021-06-09
- 語(yǔ)言: Matlab
- 標(biāo)簽: 小波??調(diào)制識(shí)別??
資源簡(jiǎn)介
小波方差用于調(diào)制識(shí)別,matlab程序,已經(jīng)調(diào)試通過(guò)。
代碼片段和文件信息
scalar=4;
%--------------------------------------------------------------------------
%AM信號(hào)
%clear
%clc
t0=0.000008;??????????????????????????????????????%定義t0信號(hào)的持續(xù)時(shí)間
tz=0.000000001;???????????????????????????????????%定義抽樣時(shí)間
fc=100000000;?????????????????????????????????????%定義載波頻率
a=0.8;????????????????????????????????????????????%定義調(diào)制系數(shù)
fz=1/tz;????????????????????????????????????????????????????????
t=[0:tz:t0];??????????????????????????????????????%定義出抽樣點(diǎn)數(shù)據(jù)
m=sin(2*pi*10000000*t);???????????????????????????%定義信號(hào)m??
c=cos(2*pi*fc.*t);????????????????????????????????%載波信號(hào)??
u=(1+a*m).*c;??
U=fftshift(fft(u));
scalar=4;%小波變換的尺度
snr=20;
y=awgn(usnr);
s=cwt(yscalar‘haar‘);
s=abs(s);????????????????????????????%取小波系數(shù)的絕對(duì)值
%--------------------------------------------------------------------------
%中值濾波函數(shù)。x是需要濾波的信號(hào)n是模板大小(即1×n)
n=5;??????????????????????????????????%奇數(shù),便于計(jì)算
L=length(s);??????????????????????????%輸入圖像是1×p的且p>n
s2=s;
for?i=1:L-n+1
??????????c=s(i:i+(n-1));?????????????%取出x1中從(i)開(kāi)始的n個(gè)元素?
??????????[caindex?]=sort(c);????????%升序排列C=fliplr(B);反轉(zhuǎn)操作,即倒序
??????????mm=median(ca);??????????????%mm是中值
??????????s2(i+(n-1)/2)=mm;???????????%將中值賦給模板中心位置的元素
end?
d=uint8(s2);?????????????????????????%未被賦值的元素取原值
d=double(d);
v1=var(d);
%--------------------------------------------------------------------------
%2ASK
t0=0.000008;??????????????????????????????????????%定義t0信號(hào)的持續(xù)時(shí)間
tz=0.000000001;?????????????????????????????????????????%定義抽樣時(shí)間
t=0:tz:t0;?????????????????????????????%定義橫軸范圍
fz=1/tz;
y=sin(2*pi*10000000.*t);?????????????????????????????????????%定義信號(hào)m??
x=[ones(11000)zeros(11000)ones(11000)ones(11000)zeros(11000)zeros(11000)ones(11000)zeros(11001)];?%定義一個(gè)與二元序列對(duì)應(yīng)的時(shí)間序列10110010
u=x.*y;??????????????????????????????????%幅頻鍵控
U=fftshift(fft(u));?
y=awgn(usnr);
s=cwt(yscalar‘haar‘);
s=abs(s);
%--------------------------------------------------------------------------
%中值濾波函數(shù)。x是需要濾波的信號(hào)n是模板大小(即1×n)
L=length(s);??????????????????????????%輸入圖像是1×p的且p>n
s2=s;
for?i=1:L-n+1
??????????c=s(i:i+(n-1));?????????????%取出x1中從(i)開(kāi)始的n個(gè)元素?
??????????[caindex?]=sort(c);????????%升序排列C=fliplr(B);反轉(zhuǎn)操作,即倒序
??????????mm=median(ca);??????????????%mm是中值
??????????s2(i+(n-1)/2)=mm;???????????%將中值賦給模板中心位置的元素
end?
d=uint8(s2);?????????????????????????%未被賦值的元素取原值
d=double(d);
v2=var(d);
%--------------------------------------------------------------------------
%2FSK:
t0=0.000008;??????????????????????????????????????%定義t0信號(hào)的持續(xù)時(shí)間
tz=0.000000001;?????????????????????????????????????????%定義抽樣時(shí)間
t=0:tz:t0;
fz=1/tz;
?f=[ones(11000)5*ones(11000)ones(11000)ones(11000)5*ones(11000)5*ones(1?1000)ones(11000)5*ones(11001)];%定義序列10110010
u=sin(2*pi*10000000*f.*t);?????????????????????????????????????%定義信號(hào)m??
U=fftshift(fft(u));?
y=awgn(usnr);
s=cwt(yscalar‘haar‘);
s=abs(s);
%---------------------------------------------------------------------
評(píng)論
共有 條評(píng)論