music2.m

张贤达老师 通信信号处理4－9章的程序 这本书是张老师的经典之作。

clear
randn('state',0);
n = 0:1280;
x = sin(0.23*pi*n) + sin(0.2*pi*n) + 0.01*randn(size(n));
FONTSIZE=0.315;         % 单位为厘米;  
LINEWIDTH_AXES=1.2;%140/42;  % corresponding to 0.140mm 
LINEWIDTH_WAVE=0.5;%180/42;  % corresponding to 0.180mm
figNumber=figure('Name','      ', 'Visible','on');
set(gcf,'Color','white');
% plot(        ,'red','LineWidth',LINEWIDTH_WAVE);  
pmusic(x,4)

xlabel('Normalized Frequency  (\times\pi rad/sample)','FontUnits','centimeters','FontName','Times New Roman','FontSize',FONTSIZE);  
ylabel('dB','FontUnits','centimeters','FontName','Times New Roman','FontSize',FONTSIZE);
title('\bf MUSIC Pseudospectrum','Color','blue','FontUnits','centimeters','FontName','Times New Roman','FontSize',FONTSIZE*1.2);
set(gca,'XColor','blue','YColor','blue','FontUnits','centimeters','FontName','Times New Roman','FontSize',FONTSIZE,'LineWidth',LINEWIDTH_AXES);
set(figNumber,'PaperPosition',[0.634518 6.34518 20.3046*0.8 15.2284*0.8]);     %Shrink the figure size