有關(guān)信道編碼中turbo碼的原理與實(shí)現(xiàn)-資料下載頁
2025-06-29 02:35本頁面
【正文】 a(k1)))。 %計(jì)算輸入為1時的gamma的值； gamma(last_state(state2,2))= (rec_s(2*k3)+rec_s(2*k2)*last_out(state2,4)).... +L_a(k1)log(1+exp(L_a(k1)))。 %根據(jù)gamma值和前一時刻的alpha值計(jì)算前向遞推計(jì)算當(dāng)前時刻的alpha值 if(sum(exp(gamma+Alpha(k1,:)))1e300) Alpha(k,state2)=Infty。 else Alpha(k,state2) = log( sum( exp( gamma+Alpha(k1,:) ) ) )。 end end %alpha值歸一化，對數(shù)形式為減去最大值 tempmax(k) = max(Alpha(k,:))。 Alpha(k,:) = Alpha(k,:) tempmax(k)。end % 后向遞推計(jì)算beta值for k = L_total1:1:1 for state1 = 1:nstates gamma = Infty*ones(1,nstates)。 %輸入為0時的gamma值 gamma(next_state(state1,1)) = (rec_s(2*k+1)+rec_s(2*k+2)*next_out(state1,2)).... log(1+exp(L_a(k+1)))。 %輸入為1時的gamma值 gamma(next_state(state1,2)) = (rec_s(2*k+1)+rec_s(2*k+2)*next_out(state1,4)).... +L_a(k+1)log(1+exp(L_a(k+1)))。 %根據(jù)gamma值和前一時刻的beta值計(jì)算前向遞推計(jì)算當(dāng)前時刻的beta值 if(sum(exp(gamma+Beta(k+1,:)))1e300) Beta(k,state1)=Infty。 else Beta(k,state1) = log(sum(exp(gamma+Beta(k+1,:))))。 end end %beta值歸一化，對數(shù)形式為減去最大值 Beta(k,:) = Beta(k,:) tempmax(k+1)。end% 計(jì)算對數(shù)似然比形式的軟輸出for k = 1:L_total for state2 = 1:nstates gamma0 = (rec_s(2*k1)+rec_s(2*k)*last_out(state2,2)).... log(1+exp(L_a(k)))。 gamma1 = (rec_s(2*k1)+rec_s(2*k)*last_out(state2,4))... +L_a(k)log(1+exp(L_a(k)))。 temp0(state2) = exp(gamma0 + Alpha(k,last_state(state2,1)) + Beta(k,state2))。 temp1(state2) = exp(gamma1 + Alpha(k,last_state(state2,2)) + Beta(k,state2))。 end L_all(k) = log(sum(temp1)) log(sum(temp0))。end
電大資料相關(guān)推薦
鄂ICP備17016276號-1