Matlab.zip_dungluong_luong_mo在matlab

hold on 'Do rong pho' w=15;'m' 'Goc den so voi truc pho' u=20;'do' 'Tan so' f=880;'MHz' 'Ban kinh' r=0.1:0.1:8;'km' 'Do cao anten di dong' hm=1.5;'m' 'Khoang cach giua cac toa nha doc theo duong truyen' b=25;'m' 'Do cao tram BTS' hb=30;'m' 'Do cao toa nha' hr=15;'m' 'Ton hao khong gian tu do' lf=32.4+20*log10(r)+20*log10(f);'dB' if (0<=u & u<55) l0=-9.646; else l0=2.5+0.075*(90-55); end if (hb>=hr) lbsh=-18*log10(11)+(hb-hr); else lbsh=0; end if (hb>hr) ka=54; else if (hb<=hr & r>=0.5) ka=54-0.8*hb; else ka=54-1.6*(hb-hr)*r; end end if (hb<hr) kd=18; else kd=18-15*((hb-hr)./(hr-hm)); end 'Nhieu xa mai nha-pho va ton hao tan xa' lts=l0+20*log10(hr-hm)+10*log10(f)-10*log10(w)-16.7; 'Trung tam thanh pho' kf=4+1.5*((f./925)-1); 'Thanh pho trung binh va vung ngoai o' kf1=4+0.7*((f./925)-1); 'Doi voi trung tam thanh pho' 'Ton hao cac vat che chan' lms=lbsh+ka+kd*log10(r)+kf*log10(f)-9*log10(b); 'Ton hao cho phep' lcp=lf+lts+lms; subplot(2,1,1) grid on hold on plot(r,lcp,'g*-') xlabel('Ban kinh [km]') title('Suy hao thanh pho lon theo mo hinh Walfsch-Ikegami ') ylabel('Suy hao [dB]') 'Doi voi thanh pho trung binh va vung ngoai o' 'Ton hao cac vat che chan' lms1=lbsh+ka+kd*log10(r)+kf1*log10(f)-9*log10(b); 'Ton hao cho phep' lcp1=lf+lts+lms1; subplot(2,1,2) grid on hold on plot(r,lcp1,'k*-') xlabel('Ban kinh [km]') title('Suy hao thanh pho trung binh va ngoai o theo mo hinh Walfsch-Ikegami ') ylabel('Suy hao [dB]') hold off