一十种概率密度函数

一十种概率密度函数 function zhifangtu(x,m) %画数据的直方图,x表示要画的随机数,m表示所要画的条数 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a=min(x); b=max(x); l=length(x); h=(b-a)/m; %量化x x=x/h; x=ceil(x); w=zeros(1,m); for i=1:l for j=1:m if (x(i)==j) %x(i)落在j的区间上，则w(j)加1 w(j)=w(j)+1; else continue end end end w=w/(h*l); z=a:h:(b-h); bar(z,w); title('直方图') function y=junyun(n) %0-1的均匀分布，n代表数据量，一般要大于1024 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% y=ones(1,n); x=ones(1,n); m=100000; x0=mod(ceil(m*rand(1,1)),m); x0=floor(x0/2); x0=2*x0+1; u=11; x(1)=x0; for i=1:n-1 x(i+1)=u*x(i)+0; x(i+1)=mod(x(i+1),m); x(i)=x(i)/m; end %x(n)单位化 x(n)=x(n)/m; y=x; function y=zhishu(m,n) %指数分布,m表示指数分布的参数,m不能为0.n表示数据量,n一般要大于1024 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% x=junyun(n); for i=1;n if (x(i)==0) x(i)=0.0001; else continue; end end u=log(x); y=-(1/m)*u; function y=ruili(m,n) %瑞利分布,m是瑞利分布的参数,n代表数据量,n一般要大于1024 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% x=junyun(n); for i=1:n if (x(i)==0) x(i)=0.0001; else continue; end end u=(-2)*log(x); y=m*sqrt(u); function y=weibuer(a,b,n) %韦布尔分布,a,b表示参数,b不能为0.n表示数据量,一般要大于1024 %a=1时,是指数分布 %a=2时,是瑞利分布 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% x=junyun(n); for i=1:n if (x(i)==0) x(i)=0.0001; else continue; end end u=-log(x); y=b*u.^(1/a); function y=swerling(n) %swelingII分布 %%%%%%%%%%%%%%%%%%%%%% r=ones(1,n); u=junyun(n); v=junyun(n); for i=1:n if (u(i)==0) u(i)=0.0001; else continue end end for i=1:n if (u(i)==v(i)) u(i)=u(i)+0.0001 else continue end end t=-2*log(u); h=2*pi*v; x=sqrt(t).*cos(h); z=sqrt(t).*sin(h); y=(r/2).*(x.^2+z.^2); function y=bernoulli(p,n) %产生数据量为n的贝努利分布,其中p属于(0-1)之间。 %----------------------- % u=junyun(n); y=zeros(1,n); for i=1:n if(u(i)<=p) y(i)=1; else y(i)=0; end end function y=duishuzhengtai(a,b,n) %产生对数正态分布，a,b为随机分布的参数，n为数据量 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% x=gaussian(n); u=sqrt(b)*x+a; y=exp(u); function y=kaifeng(m,n) %产生开丰分布，其中m代表开丰分布的自由度，n表示产生的点数量 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% y=zeros(1,n); if(floor(m/2)==m/2) for i=1:m/2 [x1,x2]=gaussian(n); for j=1:n y(j)=x1(j)^2+x2(j)^2+y(j); end end else for i=1:floor(m/2) [x1,x2]=gaussian(n); for j=1:n y(j)=x1(j)^2+x2(j)^2+y(j); end end x=gaussian(n); for j=1:n y(j)=y(j)+x(j)^2; end end function y=dajiama(a,b,n) %产生伽马随机分布的数据，a、b为随机分布的参数，数据量为n %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % k=1; if(a<1) while(k<=n) x1=junyun(1); x2=junyun(1); y2=(exp(1)+a)/exp(1)*x2; if(y2>1) p=-log(((exp(1)+a)/exp(1)-y2)/a); if(x1=1) while(k<=n) x1=junyun(1); x2=junyun(1); v=(2*a-1)^(-0.5)*log(x1/(1-x2)); x=a*exp(v); z=x1^2*x2; w=a-log(4)+(a+sqrt(2*a-1))*v-x; if(w>=log(z)) y(k)=x; k=k+1; else continue; end end end y=b*y; function y=beitafenbu(a1,a2,n) %产生贝他分布的随机数，其中a1、a2是贝他分布的参数，n代表数据量 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% x1=dajiama(a1,1,n); x2=dajiama(a2,1,n); y=x1./(x1+x2); function [y1,y2]=gaussian(n) %产生数据量为n的两个相互独立高斯分布y1、y2 %--------------------------------------- % k=1; y1=zeros(1,n); y2=zeros(1,n); while(k<=n) u1=junyun(1); u2=junyun(1); v1=2*u1-1; v2=2*u2-1; s=v1^2+v2^2; if(s>=1) continue; elseif(s==0) k=k+1; else y1(k)=v1*sqrt(-2*log(s)/s); y2(k)=v2*sqrt(-2*log(s)/s); k=k+1; end end function y=canshu(x); y=ones(1,2); n=length(x); y(1)=sum(x)/n; z=x-y(1); z=z.^2; y(2)=sum(z)/(n-1); function y=correlation(x) %计算x的自相关函数 %%%%%%%%%%%%%%%%%%%%%%%%% n=length(x); for i=1:n x1(i)=x(n+1-i); end y=conv(x,x1); 二(三种相关杂波 function y=gaussianpu(x) %由数据量为n的高斯白噪声产生向量为n，功率谱为高斯型的高斯随机向量 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% N=0:20; f=20; T=1/256; c=2*f*T*sqrt(pi)*exp(-4*f^2*pi^2*T^2*N.^2); n=length(x); y=zeros(1,n); for k=1:n for i=20:-1:0 if ((k-i)<=0) continue; else y(k)=y(k)+c(21-i)*x(k-i); end end for i=20:40 if ((k-i)<=0) continue; else y(k)=y(k)+c(i-19)*x(k-i); end end end y=0.5*y; function y=weibuerpu(a,b,n) %由数据量为n的高斯白噪声产生向量为n，功率谱为高斯型的韦布尔分布的随机向量 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [z1,z2]=gaussian(n); z1=5*z1; z2=5*z2; y1=sqrt(b^a/2)*z1; y2=sqrt(b^a/2)*z2; x1=gaussianpu(y1); x2=gaussianpu(y2); x1=sqrt(b^a/2)*x1; x2=sqrt(b^a/2)*x2; y=x1.^2+x2.^2; b=canshu(y); y=y-b(1); function y=duishuzhengtaipu(a,b,n) %由数据量为n的高斯白噪声产生向量为n，功率谱为高斯型的对数正态随机向量 %a表示

标准

excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载

方差,b表示均值 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% z1=gaussian(n); x=gaussianpu(z1); y=a*x; y=exp(y); y=b*y; b=canshu(y); y=y-b(1); %去掉直流分量 function y=swerling2pu(n) %由数据量为n的高斯白噪声产生向量为n，功率谱为高斯型的斯维凌II型随机向量 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% r=6; [z1,z2]=gaussian(n); x1=gaussianpu(z1); x2=gaussianpu(z2); y=x1.^2+x2.^2; y=r*y; b=canshu(y); y=y-b(1); %去掉直流分量 function y=kexipu(m,n) %由数据量为n的高斯白噪声产生向量为n，功率谱为柯西谱的高斯随机向量 wc=2*pi*256; T0=1/(256*m); x=gaussian(n); y=zeros(1,n); y(1)=wc*T0*x(1); for i=2:n y(i)=wc*T0*x(i)+exp(-wc*T0)*y(i-1); end b=canshu(y); %y=y-b(1); %去掉直流分量 y=conv(y,y); y=fft(y); y=abs(y); i=1:2*n-1; plot(i,y) function plotpu(x) %绘出随机数的功率谱密度函数频域的图形。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% w=fft(x); w=abs(w); v=2*pi/length(w); i=0:v:(2*pi-v); plot(i,w); 三.雷达系统仿真 function [t,s,g,f0,fs,f1]=huibo %产生目标回波信号x,系统噪声y，地物杂波z以及回波p%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% f0=3*10^7; %发射信号频率 w=0; %发射信号初始相位 c=3*10^8; %光速 l=c/f0; %雷达信号波长(载波波长) R=40000; %目标范围 Vd=200; %雷达与目标之间的径向速度 fd=2*Vd/l; %多普勒频率 Tr=600/f0; %脉冲重复周期 N=10; %雷达脉冲串长度 f1=f0/5; %调频带宽是发射信号频率的1/5 k=10*f1/Tr; fs=3*f0; %仿真采样频率 Ts=1/fs; %Tt=2*R/c; Btar=4*pi*R/l; % M=floor(Tr*fs); %一个脉冲重复周期内的采样点数M=600 mt1=floor(2*Tr*fs/5);%mt1=720 mt2=floor(3*Tr*fs/10);%mt2=540 mt3=floor(7*Tr*fs/10);%mt3=1260 mt4=floor(3*Tr*fs/5);%mt4=1080 mt5=mt1-mt2;%mt5=180 Vgain=6; s=zeros(M,N); %回波幅度起伏 for m=1:M for n=1:N v(m,n)=(u(mt1-m)-u(mt2-m))*Vgain*cos(2*pi*f0*(m-mt2)*Ts+2*pi*k*Ts*(m-mt2)/2*(m-mt2) *Ts+2*pi*fd*n*Tr); g(m,n)=(u(mt3-m)-u(mt4-m))*Vgain*cos(2*pi*(f0*(m-mt4)*Ts+k*Ts*(m-mt4)/2*(m-mt4)*Ts) ); %u(t)是发射信号包络 end end for i=1:M/10 t(i)=Vgain*cos(2*pi*(k*Ts*(M/10-i)/2*(M/10-i)*Ts)); end for n=1:N y(1:M,n)=gaussian(M)'; %系统噪声服从高斯分布 end for n=1:N z(1:M,n)=swerling2pu(M)'; %地物杂波服从swerling2型分布 end s=v+g+y+z; %s=v+g; [m,n]=size(s); x=zeros(1,m*n); q=zeros(1,m*n); x=s(:)'; q=g(:)'; p=x; i=0:length(p)-1; subplot(2,1,1); plot(i,x),title('目标回波信号'); %目标回波信号x subplot(2,1,2); l=canshu(x); b=x-l(1); plotpu(b),title('目标回波频谱'); function y=gaofang(s,f0,fs,f1) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%高频放大器 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [M,N]=size(s); Vgain=3; %高放增益w1=(f0-4*f1/3)*2*pi/fs; %w1=26pi/45 w3=(f0-f1/3)*2*pi/fs; %w3=29pi/45 w2=(f0+7*f1/3)*2*pi/fs; %w2=37pi/45 w4=(f0+4*f1/3)*2*pi/fs; %w4=34pi/45 th=min((w3-w1),(w2-w4));%w3-w1=pi/15,w2-w4=pi/15 M1=ceil(6.6*pi/th)+1; %M1=99 %n=[0:M1-1]; w5=(w1+w3)/2;w6=(w4+w2)/2;%w5=11pi/18,w6=71pi/90 h=wide(w6,M1)-wide(w5,M1); w=(hamming(M1))'; h=h.*w; L=length(h); s=Vgain*s; z=zeros(M+L-1,N); y=zeros(M,N); K=ceil(L/2); for i=1:N for j=1:10 z(1+(j-1)*M/10:j*M/10+L-1,i)=conv(h,s(1+(j-1)*M/10:j*M/10,i)')'; y(1+(j-1)*M/10:j*M/10,i)=z(1+(j-1)*M/10+K:j*M/10+K,i); end end w=y(:)'; subplot(3,1,1); i=0:length(w)-1; plot(i,w),title('高放'); subplot(3,1,2); plotpu(w),title('频谱'); subplot(3,1,3); plotpu(h),title('幅频特性'); function y=hunpin(s,f0,fs,f1) %%%%%%%%%%%%%%%%%%%%%%%%%%进行混频，输出为中频信号 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [M,N]=size(s); i=0:M-1; z1=cos(2*pi*2*f0*i/3/fs); y1=zeros(M,N); for i=1:N y1(1:M,i)=(z1.*s(1:M,i)')'; end w1=(f0-5*f1)*2*pi/(3*fs);%w1=pi/9 w3=(f0-2*f1)*2*pi/(3*fs);%w3=1.6pi/9 w2=(f0+8*f1)*2*pi/(3*fs);%w2=0.4pi w4=(f0+5*f1)*2*pi/(3*fs);%w4=pi/3 th=min((w3-w1),(w2-w4)); %w3-w1=0.2pi/3,w2-w4=0.2pi/3 M1=ceil(6.6*pi/th)+1; %M1=99 %n=[0:M1-1]; w5=(w1+w3)/2;w6=(w4+w2)/2;%w5=1.3pi/9,w6=1.1pi/3 h=wide(w6,M1)-wide(w5,M1); wth=(hamming(M1))'; h=h.*wth; L=length(h); z=zeros(M+L-1,N); y=zeros(M,N); K=ceil(L/2); for i=1:N for j=1:10 z(1+(j-1)*M/10:j*M/10+L-1,i)=conv(h,s(1+(j-1)*M/10:j*M/10,i)')'; y(1+(j-1)*M/10:j*M/10,i)=z(1+(j-1)*M/10+K:j*M/10+K,i); end end w=y(:)'; v=y1(:)'; subplot(4,1,1); i=0:length(w)-1; plot(i,w),title('混频'); subplot(4,1,2); plotpu(v),title('频谱'); subplot(4,1,3); plotpu(w); subplot(4,1,4); plotpu(h),title('幅频特性'); function y=zhongfang(s,f0,fs,f1) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%中频放大器 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [M,N]=size(s); Vgain=5; %中放增益 w1=(f0-5*f1)*2*pi/(3*fs);%w1=pi/9 w3=(f0-2*f1)*2*pi/(3*fs);%w3=1.6pi/9 w2=(f0+8*f1)*2*pi/(3*fs);%w2=0.4pi w4=(f0+5*f1)*2*pi/(3*fs);%w4=pi/3 th=min((w3-w1),(w2-w4)); %w3-w1=0.2pi/3,w2-w4=0.2pi/3 M1=ceil(6.6*pi/th)+1; %M1=99 %n=[0:M1-1]; w5=(w1+w3)/2;w6=(w4+w2)/2;%w5=1.3pi/9,w6=1.1pi/3 h=wide(w6,M1)-wide(w5,M1); w=(hamming(M1))'; h=h.*w; L=length(h); s=Vgain*s; z=zeros(M+L-1,N); y=zeros(M,N); K=ceil(L/2); for i=1:N for j=1:10 z(1+(j-1)*M/10:j*M/10+L-1,i)=conv(h,s(1+(j-1)*M/10:j*M/10,i)')'; y(1+(j-1)*M/10:j*M/10,i)=z(1+(j-1)*M/10+K:j*M/10+K,i); end end w=y(:)'; subplot(3,1,1); i=0:length(w)-1; plot(i,w),title('中放'); subplot(3,1,2); plotpu(w),title('频谱'); subplot(3,1,3); plotpu(h),title('幅频特性'); function [I,Q]=xiangganjianbo(s,fs,f0,f1) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%相位相干检波 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [M,N]=size(s); i=0:M-1; z1=cos(2*pi*f0*i/(3*fs)); z2=sin(2*pi*f0*i/(3*fs)); y1=zeros(M,N); y2=zeros(M,N); for i=1:N y1(1:M,i)=(z1.*s(1:M,i)')'; y2(1:M,i)=(z2.*s(1:M,i)')'; end wt1=2*pi*(4*f1)/(3*fs);%wt1=4pi/45 wt2=2*pi*(7*f1)/(3*fs);%wt2=7pi/45 th=wt2-wt1; M1=ceil(6.6*pi/th)+1;%M1=70 %n=[0:M1-1]; wt3=(wt1+wt2)/2;%wt3=11pi/90 hd=wide(wt3,M1); w=(hamming(M1))'; h=hd.*w; L=length(h); z1=zeros(M+L-1,N); z2=zeros(M+L-1,N); I=zeros(M,N); Q=zeros(M,N); K=ceil(L/2); for i=1:N for j=1:10 z1(1+(j-1)*M/10:j*M/10+L-1,i)=conv(h,y1(1+(j-1)*M/10:j*M/10,i)')'; z2(1+(j-1)*M/10:j*M/10+L-1,i)=conv(h,y2(1+(j-1)*M/10:j*M/10,i)')'; I(1+(j-1)*M/10:j*M/10,i)=z1(1+(j-1)*M/10+K:j*M/10+K,i); Q(1+(j-1)*M/10:j*M/10,i)=z2(1+(j-1)*M/10+K:j*M/10+K,i); end end w1=I(:)'; w2=Q(:)'; subplot(4,1,1); j=0:length(w1)-1; plot(j,w1),title('I路信号'); subplot(4,1,2); plot(j,w2),title('Q路信号'); subplot(4,1,3); w=w1+i*w2; plotpu(w),title('频谱'); subplot(4,1,4); plotpu(h),title('幅频特性'); function [x,y]=AD(z1,z2) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%模数转换 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% m1=min(min(z1)); m2=min(min(z2)); Vmax=6; N=12; x=Vmax/(2^N)*floor((z1-m1)*2^N/Vmax); y=Vmax/(2^N)*floor((z2-m2)*2^N/Vmax); [M,N]=size(z1); for i=1:N x(1:M,i)=Vmax/(2^N)*floor((z1(1:M,i)-m1)*2^N/Vmax); y(1:M,i)=Vmax/(2^N)*floor((z2(1:M,i)-m2)*2^N/Vmax); end w1=x(:)'; w2=y(:)'; subplot(2,1,1); plot(w1); subplot(2,1,2); plot(w2); function [I1,Q1]=maichongyasuo(I,Q,h) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%对正交两路信号进行脉冲压缩 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [M,N]=size(I); K=M/10; f0=3*10^7; fs=3*f0; Tr=600/f0; mt2=floor(3*Tr*fs/10); mt4=floor(3*Tr*fs/5); i1=I(mt2+1:mt2+K,1:N); q1=Q(mt2+1:mt2+K,1:N); for i=1:N S1(1:K,i)=fft(i1(1:K,i)); S2(1:K,i)=fft(q1(1:K,i)); w_ham=(hamming(K)); h1=h.*w_ham'; H1=fft(h); s(1:K,i)=(S1(1:K,i)+j*S2(1:K,i)).*H1'; S3(1:K,i)=ifft(s(1:K,i)); end i2=I(mt4+1:mt4+K,1:N); q2=Q(mt4+1:mt4+K,1:N); for i=1:N P1(1:K,i)=fft(i2(1:K,i)); P2(1:K,i)=fft(q2(1:K,i)); w_ham=(hamming(K)); h1=h.*w_ham'; H1=fft(h); p(1:K,i)=(P1(1:K,i)+j*P2(1:K,i)).*H1'; P3(1:K,i)=ifft(p(1:K,i)); end I1=[I(1:mt2,1:N)',real(S3)',I(mt2+K+1:mt4,1:N)',real(P3)',I(mt4+K+1:M,1:N)']'; Q1=[I(1:mt2,1:N)',imag(S3)',I(mt2+K+1:mt4,1:N)',imag(P3)',I(mt4+K+1:M,1:N)']'; w1=I1(:)'; w2=Q1(:)'; w3=sqrt(w1.^2+w2.^2); i=0:N*M-1; subplot(3,1,1); plot(i,w1),title('I路信号'); subplot(3,1,2); plot(i,w2),title('Q路信号'); subplot(3,1,3); plot(i,w3),title('两路合成一路'); function [s1,s2]=MTI(x1,x2) %做动目标检测 %%%%%%%%%%%%%%%%%%%%%%%%%%% [M,N]=size(x1); s1=zeros(M,N); for i=1:M for j=1:N if(i==1) s1(i,j)=x1(i,j); else s1(i,j)=x1(i,j)-x1(i-1,j); end end end for i=1:M for j=1:N if(i==1) s2(i,j)=x2(i,j); else s2(i,j)=x2(i,j)-x2(i-1,j); end end end w1=s1(:)'; w2=s2(:)'; i=0:length(w1)-1; subplot(2,1,1); plot(i,w1),title('MTI'); subplot(2,1,2); plot(i,w2); function y=qumo(x1,x2) %对两路信号进行取模运算合成一路信号 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [M,N]=size(x1); for i=1:N y(1:M,i)=sqrt(x1(1:M,i).^2+x2(1:M,i).^2); end w=y(:)'; i=0:length(w)-1; plot(i,abs(w)); function y=jilei(x) %做积累 %%%%%%%%%%%%%%%%%%%%%%%%%%% [M,N]=size(x); N1=N/10; y=zeros(M,N1); for i=1:N1 for j=1:10 y(1:M,i)=y(1:M,i)+x(1:M,i*j); end end w=y(:)'; i=0:length(w)-1; plot(i,w); function y=CFAR(x) %作恒虚警处理 %%%%%%%%%%%%%%%%%%%%%%% n=length(x); for i=1:n if(x(i)>=10^3) y(i)=x(i); else y(i)=0; end end w=y(:); i=1:length(w); plot(i,w); function y=u(n)%发射信号包络函数 %%%%%%%%%%%%%%%% if (n>=0) y=1; else y=0; end function h=wide(w,M) a=(M-1)/2; n=[0:1:(M-1)]; m=n-a+eps; h=sin(w*m)./(pi*m); function leidaxitong %总的雷达系统仿真 figure(1); [h,s1,g,f0,fs,f1]=huibo;%,title('目标回波'); figure(2); y1=gaofang(s1,f0,fs,f1);%,title('高放'); figure(3); y2=hunpin(y1,f0,fs,f1);%,title('混频'); figure(4); y3=zhongfang(y2,f0,fs,f1);%,title('中放'); figure(5); [I,Q]=xiangganjianbo(y3,fs,f0,f1);%,title('相干检波'); [x4,y4]=AD(I,Q); figure(6); [I1,Q1]=maichongyasuo(x4,y4,h);%,title('脉冲压缩'); figure(7); [y5,y6]=MTI(I1,Q1);%,title('MTI'); figure(8); y7=qumo(y5,y6),title('取模'); figure(9); y8=jilei(y7),title('脉冲积累'); figure(10); y9=CFAR(y8),title('恒虚警');