%%%%%頻差噪聲的影響
close
Tc=0.8;
delt=600;
N=128;
E=-19.5:0.5:0
D=1+N*10^27*N*E*(sin(pi*N*delt*Tc)./(pi*N*delt*Tc)).^2;%%(128,600,0.8)

delt1=600;
N1=128.1;
E1=-19.5:0.5:0
D1=1+N1*10^27*N1*E1*(sin(pi*N1*delt1*Tc)./(pi*N1*delt1*Tc)).^2;%%(256,600,0.8)

delt2=600;
N2=128.2;
E2=-19.5:0.5:0
D2=1+N2*10^27*N2*E2*(sin(pi*N2*delt2*Tc)./(pi*N2*delt2*Tc)).^2;%%(512,600,0.8)

%delt3=100;
%N3=124.5;
%E3=-19.5:0.5:0
%D3=1+N3*10^27*N3*E3*(sin(pi*N3*delt3*Tc)./(pi*N3*delt3*Tc)).^2;%%(128,10,0.8

delt4=10;
N4=43.5;
E4=-19.5:0.5:0
D4=1+N4*10^27*N4*E4*(sin(pi*N4*delt4*Tc)./(pi*N4*delt4*Tc)).^2;%%(256,1,0.8

delt5=10;
N5=65;
E5=-19.5:0.5:0
D5=1+N5*10^27*N5*E5*(sin(pi*N5*delt5*Tc)./(pi*N5*delt5*Tc)).^2;%%(512,1,0.8

delt6=5000;
N6=64;
E6=-19.5:0.5:0
D6=1+N6*10^27*N6*E6*(sin(pi*N6*delt6*Tc)./(pi*N6*delt6*Tc)).^2;%%(256,5000,0.8

delt7=10000;
N7=64.6;
E7=-19.5:0.5:0
D7=1+N7*10^27*N7*E7*(sin(pi*N7*delt7*Tc)./(pi*N7*delt7*Tc)).^2;%%(256,10000,0.8

%delt8=10000;
%N8=64.5;
%E8=-19.5:0.5:0
%D8=1+N8*10^27*N8*E8*(sin(pi*N8*delt8*Tc)./(pi*N8*delt8*Tc)).^2;%%(256,10000,0.8

figure(1);
plot(E,D,'-.b',E1,D1,':b')%,E2,D2,'-r')
%hold on
%plot(E3,D3,'-.r','LineWidth',2)
hold on
plot(E4,D4,'--r','LineWidth',1)
hold on
plot(E5,D5,':r','LineWidth',0.8)
hold on
plot(E6,D6,'--b','LineWidth',1)
hold on
plot(E7,D7,'-.r','LineWidth',1)
%hold on
%plot(E8,D8,'b','LineWidth',1.5)


set(gca,'Xscale','log')

axis([-20  -0.5  0.15  1])
set(gca,'XTick',[0.001,0.1])
grid on
text(-20,0.11,'-20')
text(-8.4,0.11,'-15')
text(-3.5,0.11,'-10')
text(-1.4,0.11,'-5')
text(-0.55,0.11,'0')