dsp32.m

现代信号处理里的一些内容

%DSP32.m
N=8096;
NR1=1024;
NR2=8096;
p=100;
pr=[10,100];
f1=100;
f2=110;
Fs=1000;
Ph0=0;
sn1=-20;
sn2=-30;
t=[0:N-1]./Fs;
FLAG=2;%tt{FLAG}通过选择FLAG的值显示使用的估计方法；
tt={'自相关法(N=','Burg法(N=','协方差法(N=','修正协方差法(N='};
A1=sqrt(2*10.^(sn1/10));
A2=sqrt(2*10.^(sn2/10));
%w=normrnd(0,1,1,N);%产生干扰为均值为0,方差为1的高斯噪声;
x=w+A1.*cos(2*pi*f1.*t+Ph0)+A2.*cos(2*pi*f2.*t+Ph0);
%plot(t,x);
%title('信号时域波形');
%xlabel('t/s');
%ylabel('x(t)');
%axis([min(t),max(t),min(x),max(x)]);
%grid on;
figure;
subplot(2,2,1);

[Pxx,f]=pburg(x(1:NR1),pr(1),[],Fs);%11111
plot(f,10*log10(Pxx));
title([tt{FLAG},num2str(NR1),',p=',num2str(pr(1)),')']);
xlabel('f/Hz');
ylabel('P_{xx}/dB');
grid on;
subplot(2,2,2);

[Pxx,f]=pburg(x(1:NR2),pr(1),[],Fs);%11111
plot(f,10*log10(Pxx));
title([tt{FLAG},num2str(NR2),',p=',num2str(pr(1)),')']);
xlabel('f/Hz');
ylabel('P_{xx}/dB');
grid on;
subplot(2,2,3);

[Pxx,f]=pburg(x(1:NR1),pr(2),[],Fs);%111111
plot(f,10*log10(Pxx));
title([tt{FLAG},num2str(NR1),',p=',num2str(pr(2)),')']);
xlabel('f/Hz');
ylabel('P_{xx}/dB');
grid on;
subplot(2,2,4);

[Pxx,f]=pburg(x(1:NR2),pr(2),[],Fs);%111111
plot(f,10*log10(Pxx));
title([tt{FLAG},num2str(NR2),',p=',num2str(pr(2)),')']);
xlabel('f/Hz');
ylabel('P_{xx}/dB');
grid on;