meizaosheng.m

这是一些自己几个月的劳动成果

clear,clc,clf;
t=0:10^(-9):0.04*10^(-4);
N=length(t);
lamde=0.03;
Pj=50;
Lj=7;
Lr=10;
Br=2*10^6;
Bj=10*10^6;
mae=1/3;
U0=sqrt(2*Pj/(mae^2+1));
Gj=13;
Ga=sinc(0.3-0.15)+sinc(0.3+0.15);
Gb=sinc(0.3-0.15)-sinc(0.3+0.15);
eta=0.7;
gmax=178;
Gaa=eta*gmax*Ga;
Gbb=eta*gmax*Gb;
R=30*10^3;
Pjc=(lamde^2*Br)/((4*pi*R)^2*10^(Lj/10)*10^(Lr/10)*Bj);
U01=-U0*sqrt(Pjc)*10^(Gj/20)*Gaa;
U02=-U0*sqrt(Pjc)*10^(Gj/20)*Gbb;
M=square(2*pi*1e3*t,0.1);
x=M.*cos(2*pi*1e8*t+pi/3);
s1=x.*U01;
s2=x.*U02;
local=cos(2*pi*1.6*10^8*t);
ylocal=fft(local);
n=0:N-1;
MUT1=s1.*local;
MUT2=s2.*local;
hhlen=length(MUT1);
r1=fir1(hhlen,[0.118,0.122],'bandpass');
r2=fir1(hhlen,[0.118,0.122],'bandpass');
h1=conv(r1,MUT1);
h1=h1(hhlen-(hhlen-1)/2:hhlen+(hhlen-1)/2);
h2=conv(r2,MUT2);  
h2=h2(hhlen-(hhlen-1)/2:hhlen+(hhlen-1)/2);
j=fir1(hhlen,0.002,'low');
J1=h1.*h2;
J2=h1.*h1;
J11=conv(j,J1);
J11=J11(hhlen-(hhlen-1)/2:hhlen+(hhlen-1)/2);
J21=conv(j,J2);
J21=J21(hhlen-(hhlen-1)/2:hhlen+(hhlen-1)/2);
Phase1=h2./h1;
Phase2=J11./J21;
r3=fir1(hhlen,[0.118,0.122],'bandpass');
freqz(r3);
r4=fir1(hhlen,[0.118,0.123],'bandpass');
h3=conv(r3,MUT1);
h3=h3(hhlen-(hhlen-1)/2:hhlen+(hhlen-1)/2);
h4=conv(r4,MUT2);  
h4=h4(hhlen-(hhlen-1)/2:hhlen+(hhlen-1)/2);
j=fir1(hhlen,0.002,'low');
J3=h3.*h4;
J4=h3.*h3;
J31=conv(j,J3);
J31=J31(hhlen-(hhlen-1)/2:hhlen+(hhlen-1)/2);
J41=conv(j,J4);
J41=J41(hhlen-(hhlen-1)/2:hhlen+(hhlen-1)/2);
Phase3=h4./h3;
Phase4=J31./J41;
subplot(221);plot(t,Phase1);grid;
xlabel('(a)');title('第一类鉴相处理','color','b');
subplot(222);plot(t,Phase2);grid;
xlabel('(b)');
title('第二类鉴相处理','color','b');
subplot(223);plot(t,Phase3);grid;
xlabel('(c)');
subplot(224);plot(t,Phase4);grid;
xlabel('(d)');