Æû³µÀíÂÛmatlab×÷Òµ

plot(r,Law1,'r',r,Law2,'b',r,Law3,'g')

gtext('ez2') gtext('eFdG') gtext('efd') xlabel('r')

ylabel('efd/eFdG/ez2 dB')

title('fd,edFdG,ez2ËærµÄ±ä»¯') u=linspace(5,20,151) for i=1:151

f0=1.5;g=0.25;r=9;u1=10;fs=3;gs=0.25;v=20;Gqn0=2.56*10^(-8);n0=0.1; df=0.2;N=180; f=0.000001:0.2:36; ff0=f/f0; ffs1=f/fs; w=2*pi*f;

d=((1-ff0.^2).*(1+r-ff0.^2/u(i))-1).^2+4.*g^2.*ff0.^2.*(r-(1/u(i)+1).*ff0.^2).^2; z1q=r*sqrt(((1-ff0.^2).^2+4*g^2*ff0.^2)./d);

z2z1=sqrt((1+(2*g*ff0).^2)./((1-ff0.^2).^2+(2*g*ff0).^2)); pz2=sqrt((1+(2*gs*ffs1).^2)./((1-ffs1.^2).^2+(2*gs*ffs1).^2)); ez2(i)=sqrt(sum(z2z1.^2.*z1q.^2.*4.*pi^2.*Gqn0.*v.*n0^2*df)); FdG=(r.*w/9.8).^2.*((ff0.^2./(1+u(i))-1).^2+4*g^2.*ff0.^2)./d; eFdG(i)=sqrt(sum(FdG.*4.*pi^2.*Gqn0.*v.*n0^2*df)); fd=(r*ff0.^2./w).^2./d;

efd(i)=sqrt(sum(fd.*4.*pi^2.*Gqn0.*v.*n0^2*df)) ; end

for i=1:151

Law1(i)=20*log(ez2(i)/ez2(34)); Law2(i)=20*log(eFdG(i)/eFdG(34)); Law3(i)=20*log(efd(i)/efd(34)); end

plot(u,Law1,'r',u,Law2,'b',u,Law3,'g')

gtext('ez2') gtext('eFdG') gtext('efd') xlabel('u')

ylabel('efd/eFdG/ez2 dB') title('fd, dFdG,ez2ËæuµÄ±ä»¯')