数字信号处理功率谱估计Matlab源代码

数字信号处理 丁玉美著 第四章答案 希望大家会有用处

close all
clear
clc
%%%%%%%%%%%%%%%%%%%    输入信号   %%%%%%%%%%%%%%%%%%%
N = 512;fs=100;IP=80;;
f1 = 0.2*fs;f2 = 0.3*fs;
fan1 = 0; fan2 = 0;
t = 0:1./fs:（N-1）/fs;
x = sin（2*pi*f1*t + fan1） + sin（2*pi*f2*t + fan2）;
y = awgn（x10）;
plot（tx）;title（‘无噪声原始信号‘）;
figure;plot（ty）;title（‘加噪输入信号‘）;
%%%%%%%%%%%%%%%%%%%   周期图法求功率谱 %%%%%%%%%%%%%%%%
y_z = fftshift（abs（fft（yN））.^2/N）;
df=fs./N;
f=（0:df:df*（N-1））-fs./2;
figure;plot（f10*log10（y_z））;title（‘周期图法（dB）谱估计‘）;

%%%%%%%%%%%%%%%%%%%%%%% p=0     %%%%%%%%%%%%%%%%%%%%%%
ef=zeros（IPN）;
eb=zeros（IPN）;
ef（1:）=y;      % ef（0:）
eb（1:）=y;      % eb（0:）
row（1）=y*y‘/N;  %ρ=row
k（1）=0;         %k（0）
%%%%%%%%%%%%%计算a（pi） power（p） p=1 2 3 ...IP%%%%%%%%%% 
a=zeros（IP+1IP+1）;
power（1）=row（1）;
for p=2:IP+1
    numerator=0;
    denominator=0;
    for n=p:N
    numerator= numerator+ （-2）*ef（p-1n）*eb（p-1n-1）;
    denominator=denominator+（ef（p-1n））^2+（eb（p-1n-1））^2;
    end
    k（p）=numerator./（denominator+0.0001）;  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1732  2006-11-29 13:57  功率普估计\burg.m

     文件       1332  2006-11-10 09:21  功率普估计\levenson.m

     文件        361  2006-11-10 00:05  功率普估计\zhouqitu.m

     目录          0  2009-06-07 19:36  功率普估计

----------- ---------  ---------- -----  ----

                 3425                    4