科創(chuàng)數(shù)模樹模城市微小區(qū)地圖對(duì)應(yīng)數(shù)據(jù)txt

line_number=1;pointnumber=index_list=i=whiletline=ifstr2num(tline(1))==0pointnumber=0;data.X(line_number)=str2num(tline(1:21));line_number=line_number+1;pointnumber=pointnumber+1;ift=index_list=forindex_list=holdon;t=t+clearallcloseall;tp= tp=tp+cot((pi-a2+a1)/2/n).*F(k*L*2*(cos((2*pi*n*N1-a2+a1)/2)).^2);tp=tp-cot((pi+a2+a1)/2/n).*F(k*L*2*(cos((2*pi*n*N2-a2-a1)/2)).^2);tp=tp-cot((pi-a2-a1)/2/n).*F(k*L*2*(cos((2*pi*n*N1-a2-a1)/2)).^2);symstx=xall(i);ifelseifx>10elseif(x>=0.001)&&(x<=10);1111111111111-111111%[500;200]開始,(7次以內(nèi)clearallN=z_end=[300;350];%接收點(diǎn)fori=1:N%對(duì)于每個(gè)角度tht=2*pi/N*i;path{i}=N_rflzeros(1,2);%while f_seek_intersec(z_bgn,tht,rec_seg_y,'row'); f_seek_intersec(z_bgn,tht,rec_seg_x,'clm');ifdistance_y<distance_xintersec=intersec_y;typ_rfl=typ_rfl_y;

intersec=intersec_x;typ_rfl=typ_rfl_x;clearilne_min_yintersec_ytheta_intersec_ydistance_ytyp_rfl_yilne_min_xintersec_xtheta_intersec_xdistance_xtyp_rfl_xiftyp_rflN_rfl(typ_rfl)=N_rfl(typ_rfl)+

iff_ifExceedTimes(N_rfl)z_bgn=

tht=clearintersecp(:,1)=p(:,2)=ifj<7+1end;clearj;%hold end;clearifunction[ilne_min,intersec,theta_intersec,distance,typ_rfl]=f_seek_intersec(z,theta,segs,%z為起點(diǎn),thetasegstyprow'時(shí)為橫線段集合'clm'%縱線段時(shí),90度計(jì)算,再做相應(yīng)變換即可,因此首先只需考慮橫theta0=ifstrcmp(typ,'clm')z=flipud(z);%theta_intersec=-distance=1e10;fori=ydfr=segs(i,1)-ifsegs(i,4)*sin(theta)<0&&sin(theta)*ydfr>%外法方向與射線方向相反,intersec=z(1)-

intersec=z(1)+if(intersec-segs(i,2))*(intersec-segs(i,3))<=0%若交點(diǎn)段端點(diǎn)中間%?考慮繞射時(shí),須改disnow=(intersec-z(1))^2+ydfr^2;ifdisnow<distanceilne_min=intersec_min= distance=disnow;

if1

typ_rfl1;%typ_rfl2;%clearend;cleartheta_intersec=pi-theta0;%clear

theta_intersec=-theta0;%clearilne_min=-ifsin(theta)>ybiankuang=

ybiankuang=ifstrcmp(typ,'clm')intersec=z(1)-

intersec=z(1)+

typ_rfl=0;intersec=intersec=rec_seg_y=load('SegY.txt');rec_seg_x=load('SegX.txt');rec_seg_x=sortrows(rec_seg_x,1);if_E=0;ifN(1)>7+1||N(2)>if_E=%%close%p_bgnz_end之間,n_mostrfl%p_bgn=%z_end=%z_end=%n_mostrfl=%N=Np_bgn射出的均勻射線的條數(shù)[rec_seg_y,rec_seg_x]=f_rec_seg();eps_frl=1;fl_rec_Nrfl=fopen('Rec_Nrfl.txt','w');if_success=zeros(N,1);Nrfl_success=[];fori1:N%z_bgn=p_bgn;tht0=2*pi/N*i;tht=tht0;N_rfl=0;%反射次數(shù)while1 f_seek_intersec(z_bgn,tht,rec_seg_y,'row',z_end,eps_frl); f_seek_intersec(z_bgn,tht,rec_seg_x,'clm',z_end,eps_frl);ifdistance_y<ilne_min=[1;ilne_min_y];%注意:ilne_min

intersec=intersec_y;typ_rfl=typ_rfl_y;ilne_min=[2;ilne_min_x];intersec=intersec_x;typ_rfl=typ_rfl_x;clearilne_min_yintersec_ytheta_intersec_ydistance_ytyp_rfl_yilne_min_xintersec_xtheta_intersec_xdistance_xtyp_rfl_xswitchtyp_rflcase1%此時(shí)已經(jīng)到達(dá)終點(diǎn),if_success(i)=1;case0caseN_rflN_rfl1;case2%衍射

z_bgn=intersec;clearintersecforj=1:size(path{i},2)-1p(:,2)=ifj<n_mostrflend;clearj;hold%%forj=1:size(path{i},2)-1p(:,2)=ifj<n_mostrflend;clear hold Nrfl_success(end+1)=N_rfl; %pauseend;clearifunction[ilne_min,intersec,theta_intersec,distance,typ_rfl]=f_seek_intersec(p_bgn,theta,segs,typ,p_end,eps_ys)%p_bgn為起點(diǎn),thetasegs的最近交點(diǎn)以及反射線的角typrow'時(shí)為橫線段集合'clm'%縱線段時(shí),90度計(jì)算,再做相應(yīng)變換即可,theta0theta;z=p_bgn;ifstrcmp(typ,'clm')z=flipud(z);%theta_intersec=-distance=1e10;no_seg_ys=fori=ydfr=segs(i,1)-ifsegs(i,4)*sin(theta)<0&&sin(theta)*ydfr>%外法方向與射線方向相反,intersec=z(1)-

intersec=z(1)+if(intersec-segs(i,2))*(intersec-segs(i,3))<0%若交點(diǎn)段端點(diǎn)中disnow=(intersec-z(1))^2+%記錄更近的交點(diǎn)所ilne_min=i;distance=disnow;%如果交點(diǎn)太靠端點(diǎn),不考慮為反射現(xiàn)象,但需單獨(dú)記錄該線段的編號(hào),因?yàn)閕fabs(intersec-segs(i,2))<eps_ys||abs(intersec-segs(i,3))<eps_ysno_seg_ys=[no_seg_ys,i]; if typ_rfl1;% typ_rfl2;% clear

theta_intersec=pi-theta0;%cleartheta_intersec_mintheta_intersec=-theta0;%cleartheta_intersec_minilne_min=-ifsin(theta)>ybiankuang=

ybiankuang=ifstrcmp(typ,'clm')intersec=z(1)-

intersec=z(1)+

typ_rfl=0;typ_rfl=1;%判斷最近交點(diǎn)處是否為衍射現(xiàn)象,若是,typ_rfl=2;intersec=%%p_endp_bgnintersec之間,p_endscl=(p_end-