%求信號x的功率譜
%時域信號--（FFT變換）-->幅度譜--（平方）-->功率譜??直接法
%時域信號--->相關(guān)函數(shù)--（FFT變換）-->功率譜--（除以頻率分辨率）-->功率譜密度??間接法
clear?all
%a函數(shù)及其頻譜
m=-20:0.2:20;
index1=find（m<0）;
a?=?10*sin（m（index1））./（m（index1））;
a=[a?10];
index2?=?find（m?>?0）;
a?=?[a?10*sin（m（index2））./（m（index2））];
figure（1）;subplot（411）;stem（ma）;title（‘原函數(shù)a‘）;
M=length（m）;
p=0:M-1;q=0:M-1;
A=a*exp（-1i*2*pi/M）.^（p‘*q）;
subplot（412）;plot（fftshift（abs（A）））;title（‘原函數(shù)a的幅度譜‘）;
%調(diào)制函數(shù)b
t=0:1/100:2;
b=100*sin（2*pi*20*t）;
f=linspace（0100128）;
subplot（413）;plot（tb）;title（‘其中的一個調(diào)制函數(shù)b=50sin（2pi20t）‘）;
B=fft（b128）;
subplot（414）;plot（fabs（B））;title（‘其中的一個調(diào)制函數(shù)b的幅度譜‘）;
%調(diào)制函數(shù)c
c=80*sin（2*pi*40*t）;
ba=a.*b;ca=a.*c;
x=ba+ca+randn（size（ba））;
X=abs（fft（x128））;?
nfft=512;
index=0:nfft-1;
k=index*100/nfft;%100采樣頻率
[r?tao]=xcorr（x‘unbiased‘）;
subplot（412）;plot（taor）;title（‘x的自相關(guān)函數(shù)‘）;
Pxx=abs（fft（rnfft））;
Sf=10*log10（Pxx（index+1））;??%功率譜密度函數(shù)
n0=length（X）;
n1=length（Sf）;
%??多尺度
m=-4:-1;
for?i=1:4
????delta（i）=2.^m（i）;
end
N=16;?%?濾波器長度（需要調(diào)整，必須是偶數(shù)）
A=-1/sqrt（2*pi）;??%??幅度
%??構(gòu)造高斯函數(shù)
for?index_x=1:N;
????x=index_x-（N