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