波条干扰‘.rar_干扰_干扰效果_波条干扰_雷达 干扰_雷达干扰

function[echo]=Chaff(azimuthRadar,elevationRadar,vRadar,xRadar,yRadar,zRadar,xChaff,yChaff,zChaff,tChaff,N) % clc % clear all % ============== parameters ==================== global Pt G c f0 waveLength vChaff vMuzzle % xRadar yRadar zRadar azimuthRadar elevationRadar vRadar % xChaff yChaff zChaff elevation azimuth rChaff deltaAzimuth deltaElevation % xVector yVector zVector rChaffRadar thetaRadarChaff thetaRadarChaffMax % height % iMax jMax iMid jMid thetaNormal % the taChaff rrChaffRadar fai Rcs PtEcho fd xMissile=xRadar; yMissile=yRadar; zMissile=zRadar; % xChaff=120; % yChaff=160; % zChaff=200; % r=10; % azimuthRadar=0*pi/180; % elevationRadar=100*pi/180; % Pt=1e3; % Gdb=20; % G=10^(Gdb/10); % vRadar=1e3; % c=3e8; % f0=60e6; % waveLength=c/f0; yita=0.8;%鸟窝效应造成的损失 N=yita*N; thegma=0.17*waveLength*waveLength; %单根箔条平均反射系数 elevation=linspace(-pi/2,pi/2,101); azimuth=linspace(-pi/2,pi/2,101); deltaAzimuth=pi/101; deltaElevation=pi/101; [rChaff]=rContrail(vMuzzle,vChaff,tChaff); [height]=undulate(rChaff,rChaff/10); % deltaRcs=4*pi*r^2/(101*101); % x=r.*cos(elevation).*cos(azimuth); % y=r.*cos(elevation).*sin(azimuth); % z=r.*sin(elevation); % height=zeros(100,100); xVector=-xMissile+xChaff; yVector=-yMissile+yChaff; zVector=-zMissile+zChaff; xRadarUnity=cos(elevationRadar)*cos(azimuthRadar); yRadarUnity=cos(elevationRadar)*sin(azimuthRadar); zRadarUnity=sin(elevationRadar); %天线方向上考虑参数 rChaffRadar=sqrt(xVector^2+yVector^2+zVector^2); %雷达和箔条中心的距离 costhetaRadarChaffMid=(xVector*xRadarUnity+yVector*yRadarUnity+zVector*zRadarUnity)/rChaffRadar^2;%雷达天线与箔条中心夹角 thetaRadarChaffMid=acos(costhetaRadarChaffMid); % thetaRadarChaffMid=0; thetaNormal=0; echo=0; % rChaff=100; thetaRadarChaffMax=asin(rChaff/rChaffRadar);% 箔条占有的雷达宽度（弧度） iMax=round(2*(pi/2-thetaRadarChaffMax)/deltaElevation); %箔条占有雷达宽度俯仰点数 jMax=iMax; iMid=fix(iMax/2); jMid=fix(jMax/2); for i=1:iMax for j=1:jMax % ii=mod(i+M,2*M)+1; % jj=mod(j+N,2*N)+1; % x(i,j)=(r+height(ii,jj)).*cos(elevation(j)).*cos(azimuth(i)); % y(i,j)=(r+height(ii,jj)).*cos(elevation(j)).*sin(azimuth(i)); % z(i,j)=(r+height(ii,jj)).*sin(elevation(j)); % thetaChaff(i,j)=atan((tan(deltaAzimuth*i))^2+(tan(deltaElevation*j))^2); rrChaff(i,j)=rChaff+height(i,j);%每块的起伏半径 vectorAzimuth=[cos(azimuth(i)),sin(azimuth(i)),0]; %横向量 vectorElevation=[cos(elevation(j)),0,sin(elevation(j))];%纵向量 thetaChaff(i,j)=acos(dot(vectorAzimuth,vectorElevation));%每块与箔条中心的夹角 rrChaffRadar(i,j)=sqrt(rrChaff(i,j)^2+rChaffRadar^2-2*rrChaff(i,j)*rChaffRadar*cos(thetaChaff(i,j)));%回波路径 fai(i,j)=acos((rrChaffRadar(i,j)^2+rChaffRadar^2-rrChaff(i,j)^2)/(2*rrChaffRadar(i,j)*rChaffRadar));%每个部分与雷达夹角 deltaRcsReal(i,j)=rChaff^2*deltaAzimuth*(sin(elevation(j+1))-sin(elevation(j))); deltaRcs(i,j)=deltaRcsReal(i,j)*(1-exp(-N*thegma/deltaRcsReal(i,j))); %箔条干扰的单位Rcs Rcs(i,j)=deltaRcs(i,j)*cos(fai(i,j)+thetaChaff(i,j))*cos(thetaNormal); PtEcho(i,j)=Pt*G^2*Rcs(i,j)/(4*pi*rrChaffRadar(i,j)^2)^2; % fd(i,j)=vRadar*cos(thetaRadarChaff-fai(i,j)); fd(i,j)=2*vRadar*cos(thetaRadarChaffMid-azimuth(i))*cos(elevation(j))/waveLength; [deltaEchoChaff]=echoChirp(rrChaffRadar(i,j),PtEcho(i,j),fd(i,j)); echo=echo+deltaEchoChaff; end end plot(echo)