www.pudn.com > work_2.rar > figure_S.m, change:2005-06-12,size:1355b

a=polyfit(SNR1,10*log10(S1),3); 
mymin=min(SNR1); 
mymax=max(SNR1); 
xx=linspace(mymin,mymax,100); 
yy=polyval(a,xx); 
plot(SNR1,10*log10(S1),'o',xx,yy); 
hold on 
 
a=polyfit(SNR2,10*log10(S2),3); 
mymin=min(SNR2); 
mymax=max(SNR2); 
xx=linspace(mymin,mymax,100); 
yy=polyval(a,xx); 
plot(SNR2,10*log10(S2),'o',xx,yy); 
hold on 
 
a=polyfit(SNR3,10*log10(S3),3); 
mymin=min(SNR3); 
mymax=max(SNR3); 
xx=linspace(mymin,mymax,100); 
yy=polyval(a,xx); 
plot(SNR3,10*log10(S3),'o',xx,yy); 
hold on 
 
a=polyfit(SNR4,10*log10(S4),3); 
mymin=min(SNR4); 
mymax=max(SNR4); 
xx=linspace(mymin,mymax,100); 
yy=polyval(a,xx); 
plot(SNR4,10*log10(S4),'o',xx,yy); 
hold on 
 
a=polyfit(SNR5,10*log10(S5),3); 
mymin=min(SNR5); 
mymax=max(SNR5); 
xx=linspace(mymin,mymax,100); 
yy=polyval(a,xx); 
plot(SNR5,10*log10(S5),'o',xx,yy); 
hold on 
 
a=polyfit(SNR6,10*log10(S6),3); 
mymin=min(SNR6); 
mymax=max(SNR6); 
xx=linspace(mymin,mymax,100); 
yy=polyval(a,xx); 
plot(SNR6,10*log10(S6),'o',xx,yy); 
hold on 
 
a=polyfit(SNR7,10*log10(S7),3); 
mymin=min(SNR7); 
mymax=max(SNR7); 
xx=linspace(mymin,mymax,100); 
yy=polyval(a,xx); 
plot(SNR7,10*log10(S7),'o',xx,yy); 
hold on 
 
 
a=polyfit(SNR,10*log10(Shannon),3); 
mymin=min(SNR); 
mymax=max(SNR); 
xx=linspace(mymin,mymax,100); 
yy=polyval(a,xx); 
plot(SNR,10*log10(Shannon),'o',xx,yy); 
hold on 
%axis([2,8,-120,-50])