【正文】 Delay: 平均延遲%****************************定義終端狀態(tài)常數(shù)以及仿真結(jié)束參數(shù)****************************STANDBY =0。 %等待TRANSMIT =1。 %傳輸COLLISION =2。 %碰撞TOTAL=5000。 %成功傳輸多少數(shù)據(jù)包后仿真結(jié)束%******************************** 定義信道參數(shù) *************************************brate=6e6。 %比特速率Srate=。 %符號速率Plen =500。 %包長(符號數(shù))Ttime=Plen/Srate。 %每個數(shù)據(jù)包的傳輸時間Dtime=。 %歸一化傳播延遲delay=Dtime*Ttime。 %實際延遲alfa =3。 %路徑損耗指數(shù)sigma=6。 %陰影衰落標(biāo)準(zhǔn)差 [dB]%**********************************定義接入點信息***********************************r =100。 %服務(wù)區(qū)域半徑 [m]bxy =[0,0,5]。 %接入點位置坐標(biāo) (x,y,z)[m]t =10。 %接入點進行正確信號解調(diào)所需要的最低信號功率 [dBm]%************************************定義終端信息 ***********************************Mnum =100。 %終端數(shù)目m =30。 %終端在服務(wù)區(qū)域邊緣時，接入點接收到的信號功率[dBm]mpow =10^(m/10)*sqrt(r^2+bxy(3)^2)^alfa。 %終端的發(fā)射信號功率h=0。 %終端高度mxy=[randsrc(2,Mnum,[r:r])。randsrc(1,Mnum,[0:h])]。 %隨機生成終端坐標(biāo)while 1 d=sqrt(sum(mxy(1:2,:).^2))。 %判斷終端與接入點的水平距離是否超過r [tmp,indx]=find(dr)。 if length(indx)==0 break end mxy(:,indx)=[randsrc(2,length(indx),[r:r])。mxy(3,indx)]。 %超過r重新生成位置坐標(biāo)enddistance=sqrt(sum(((ones(Mnum,1)*bxy).39。mxy).^2))。 %終端距離接入點的距離mrnd=randn(1,Mnum)。 %每個終端的陰影衰落G=[::1,2:10,20:20:40]。 %理論業(yè)務(wù)量for indx=1:length(G) %*************************************** 初始化相關(guān)參數(shù) ******************************** Tint =Ttime/log(1G(indx)/Mnum)。 %數(shù)據(jù)包產(chǎn)生間隔的期望值 Rint =Tint。 %數(shù)據(jù)包重傳間隔的期望值 Spnum=0。 %成功傳輸?shù)陌鼈€數(shù) Splen=0。 %成功傳輸?shù)姆柕膫€數(shù) Tplen=0。 %成功傳輸?shù)姆枖?shù) Wtime=0。 %傳輸延遲時間(s) mgtime =Tint*log(1rand(1,Mnum))。 %初始數(shù)據(jù)包產(chǎn)生時刻 Mstime =zeros(1,Mnum)inf。 %數(shù)據(jù)包傳輸時刻 mtime =mgtime。 %終端狀態(tài)改變傳輸時刻 Mstate =zeros(1,Mnum)。 %終端狀態(tài) Mplen(1:Mnum)=Plen。 %每個終端傳輸?shù)臄?shù)據(jù)包長度大小 now_time =min(mtime)。 %*************************************** 仿真循環(huán) ************************************** while 1 idx=find(mtime==now_time amp。 Mstate==TRANSMIT)。 %成功傳輸數(shù)據(jù)包的終端ID if length(idx)0 Spnum =Spnum+1。 Splen =Splen+Mplen(idx)。 Wtime =Wtime+now_timemgtime(idx)。 Mstate(idx)=STANDBY。 mgtime(idx)=now_timeTint*log(1rand)。 %下一個數(shù)據(jù)包產(chǎn)生時刻 mtime(idx) =mgtime(idx)。 %下一次狀態(tài)改變時刻 end idx=find(mtime==now_time amp。 Mstate==COLLISION)。 %數(shù)據(jù)包傳輸失敗的終端ID if length(idx)0 Mstate(idx)=STANDBY。 mtime(idx)=now_timeRint*log(1rand(1,length(idx)))。 %重傳等待時間 end idx=find(mtime==now_time amp。 Mstate==STANDBY)。 %開始傳輸數(shù)據(jù)包的終端ID if length(idx)0 Tplen=Tplen+sum(Mplen(idx))。 for ii=1:length(idx) jj=idx(ii)。 %載波監(jiān)聽 idx1=find((Mstime+delay)=now_time amp。 now_time=(Mstime+delay+Ttime))。 if length(idx1)==0 %信道空閑 Mstate(jj)=TRANSMIT。 %傳送數(shù)據(jù)包 Mstime(jj)=now_time。 %數(shù)據(jù)包開始傳輸時間 mtime(jj) =now_time+Mplen(jj)/Srate。 %數(shù)據(jù)包傳輸結(jié)束時間 else %信道忙 mtime(jj)=now_timeRint*log(1rand)。 %等待一段時間再進行監(jiān)聽 end end end if Spnum=TOTAL %如果成功傳輸?shù)臄?shù)據(jù)包達(dá)到設(shè)定條件，仿真結(jié)束 break end % 有數(shù)據(jù)包傳輸或發(fā)生碰撞的終端ID idx=find(Mstate==TRANSMIT|Mstate==COLLISION)。 if capture==0 %不考慮捕獲效應(yīng) if length(idx)1 Mstate(idx)=COLLISION。 %同時傳輸數(shù)據(jù)包的終端大于1，發(fā)生碰撞 end else %考慮捕獲效應(yīng) if length(idx)1 dxy =distance(idx)。 %比較發(fā)生碰撞的終端的距離 %計算接入點收到的各個終端信號功率,其中考慮了陰影衰落的影響 pow =mpow*dxy.^alfa.*10.^(sigma/10*mrnd(idx))。 [maxp no]=max(pow)。 if Mstate(idx(no))==TRANSMIT if length(idx)==1 =10*log10(maxp)。 else =10*log10(maxp/(sum(pow)maxp+1))。 end Mstate(idx)=COLLISION。 if =t %接收到的信號功率大于捕獲門限 Mstate(idx(no))=TRANSMIT。 %傳輸成功 end else Mstate(idx)=COLLISION。 end end end now_time=min(mtime)。 %更新時刻 end Traffic(indx)=Tplen/Srate/now_time。 %統(tǒng)計實際產(chǎn)生的業(yè)務(wù)量 S(i