scalar=4;
%--------------------------------------------------------------------------
%AM信號
%clear
%clc
t0=0.000008;??????????????????????????????????????%定義t0信號的持續時間
tz=0.000000001;???????????????????????????????????%定義抽樣時間
fc=100000000;?????????????????????????????????????%定義載波頻率
a=0.8;????????????????????????????????????????????%定義調制系數
fz=1/tz;????????????????????????????????????????????????????????
t=[0:tz:t0];??????????????????????????????????????%定義出抽樣點數據
m=sin（2*pi*10000000*t）;???????????????????????????%定義信號m??
c=cos（2*pi*fc.*t）;????????????????????????????????%載波信號??
u=（1+a*m）.*c;??
U=fftshift（fft（u））;
scalar=4;%小波變換的尺度

snr=20;
y=awgn（usnr）;
s=cwt（yscalar‘haar‘）;
s=abs（s）;????????????????????????????%取小波系數的絕對值
%--------------------------------------------------------------------------
%中值濾波函數。x是需要濾波的信號n是模板大?。?×n）
n=5;??????????????????????????????????%奇數，便于計算

L=length（s）;??????????????????????????%輸入圖像是1×p的且p>n
s2=s;
for?i=1:L-n+1
??????????c=s（i:i+（n-1））;?????????????%取出x1中從（i）開始的n個元素?
??????????[caindex?]=sort（c）;????????%升序排列C=fliplr（B）;反轉操作，即倒序
??????????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信號的持續時間
tz=0.000000001;?????????????????????????????????????????%定義抽樣時間
t=0:tz:t0;?????????????????????????????%定義橫軸范圍
fz=1/tz;
y=sin（2*pi*10000000.*t）;?????????????????????????????????????%定義信號m??
x=[ones（11000）zeros（11000）ones（11000）ones（11000）zeros（11000）zeros（11000）ones（11000）zeros（11001）];?%定義一個與二元序列對應的時間序列10110010
u=x.*y;??????????????????????????????????%幅頻鍵控
U=fftshift（fft（u））;?

y=awgn（usnr）;
s=cwt（yscalar‘haar‘）;
s=abs（s）;
%--------------------------------------------------------------------------
%中值濾波函數。x是需要濾波的信號n是模板大?。?×n）

L=length（s）;??????????????????????????%輸入圖像是1×p的且p>n
s2=s;
for?i=1:L-n+1
??????????c=s（i:i+（n-1））;?????????????%取出x1中從（i）開始的n個元素?
??????????[caindex?]=sort（c）;????????%升序排列C=fliplr（B）;反轉操作，即倒序
??????????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信號的持續時間
tz=0.000000001;?????????????????????????????????????????%定義抽樣時間
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）;?????????????????????????????????????%定義信號m??
U=fftshift（fft（u））;?

y=awgn（usnr）;
s=cwt（yscalar‘haar‘）;
s=abs（s）;
%---------------------------------------------------------------------