%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;??


am_n1=zeros（1041）;??????????????????????????????%定義曲線條數
am_n2=zeros（1041）;??????????????????????????????%定義曲線條數
N=10;
while（N）
for?snr=-10:30?
?????y=awgn（usnr）;???????????????????????????????%加信噪比為snr的噪聲
?????Yi=fftshift（fft（y））;
?????Ei=（abs（Yi））.^2;?????????????????????????????%求能量
?????E=sum（Ei）;???????????????????????????????????%總能量
?????Pi=Ei/E;
?????H1（snr+11）=-sum（Pi.*log10（Pi））;
?????H2（snr+11）=sum（Pi.*exp（1-Pi））;
?????SNR_N（snr+11）=snr;????????
end
am_n1（N:）=H1;
am_n2（N:）=H2;
N=N-1;
end
H1=mean（real（am_n1）1）;
H2=mean（real（am_n2）1）;

figure（22）
plot（SNR_NH1‘-b.‘）;title（‘功率譜香農熵‘）;
xlabel（‘SNR/dB‘）;ylabel（‘H（X）‘）;
hold?on
figure（33）
plot（SNR_NH2‘-b.‘）;title（‘功率譜指數熵‘）;
xlabel（‘SNR/dB‘）;ylabel（‘H（X）‘）;
hold?on

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1390??2014-05-25?14:44??奇異譜熵\qiyi_entropy.m

?????目錄??????????0??2014-05-25?14:37??奇異譜熵

-----------?---------??----------?-----??----

?????????????????1390????????????????????2