- 博客(26)
- 收藏
- 关注
原创 均匀量化函数
function [F] = quant_bit(bit)lon =pi;F = [0:2^bit-1]*(2*lon/2^bit)-lon+2*lon/2^(bit+1);end
2022-05-19 19:00:40 111
原创 LT码在BI-AWGN信道下的译码算法
clc;clear; K = 500; m = K; u = 0; degree = 40; %迭代次數 sigema =sqrt(0.5); % %(0.5/r)*10^(-0.1*(Eb_N0));% Eb_N0 = 1 ; %单位:dB % n = 10;% r = m/n;% P = ...
2021-01-04 17:32:26 426
原创 LT(7.21,Sun)于圖書館
clc;clear;m = 5000;K = 10000;n = 10000;c = 0.01;o = 0.5;x_1 = rand(1,m)<0.5; %信源碼字Pgenerate = rand(1,K); %隨機生成K個數count = 1;% for d = 1:K% index_d = 2...
2021-01-04 17:31:54 61
原创 LT 02
clc;clear;m = 300;K = 10000;n = 400;c = 0.01;o = 0.5;x = ones(1,m);x_1 = rand(1,m)<0.5; %信源碼字Pgenerate = rand(1,n); %隨機生成n個數count = 1;% for d = 1:K% index_d = 2...
2021-01-04 17:30:52 95
原创 LT_BP_performance
clc;clear;m = 500;K = 500;overhead = [550 600 700 800];i_ii =10000;for iii = 1:4 n = overhead(iii); for i_i = 1:i_ii if mod(i_i,100)==0 i_i end c = 0.0...
2021-01-04 17:29:29 163
原创 噴泉碼生成矩陣草稿_01
clc;clear;% k = ;x = rand(1,k)<0.5; %隨機生成碼元序列x%***************初始度分佈函數開始***************InitD = 2:K;if find(d==InitD) P = 1/d/(d-1);elseif d==1 P = 1/K;endPP(d) = P;%***************初始度分...
2021-01-04 17:29:12 56
原创 噴泉碼小程序
clc;clear;K = 10000;c = 0.2;o = 0.05;S = fix(c*sqrt(K)*log(K/o));for d = 1:50 %*******************************P(d)開始************************************ index_d = 2:K; if find(index...
2021-01-04 17:28:55 104
原创 稀疏矩陣壓縮存儲與還原
#1.模塊clc;clear;H = [1 1 1 1 0 1 0 1 1 0 1 0 1 0 1 1 0 1 1 1 0 1 1 1]C = [1 2 3 4 5 6];[m,n] = size(H);for i = 1:m k = 1; for j = 1:n if H(i,j)~=0 index_co...
2021-01-04 17:28:39 102
原创 項目二:LDPC碼高斯信道譯碼
clc;clear;degree_sum = 100;for sum_loop = 1:degree_sumload PEGH_3_6_1008.mat;[m,n] = size(H);C = zeros(1,n);for i = 1:nif C(i)==1C1(i) = -1;elseC1(i) = 1;endendC1 ;SNR = 2;b2 = 10^((-0....
2021-01-04 17:27:16 139
原创 項目一:LDPC在碼刪除信道下譯碼
clc;clear;for loop = 1:1000 load PEGH_3_6_1024.mat [m,n] = size(H); % H_simple = zhu_rank(H); %H化为行最简阵 %***********************高斯消元开始****************************...
2021-01-04 17:26:55 255 2
原创 噴泉碼LT譯碼程序
clc;clear;G = [ 1 1 0 1 0 1 1 1 0 1 1 0][m,n] = size(G);y = [0 1 1 1];x = ones(1,m)while(1) %***********判斷準備開始1***************** for i = 1:n front_judge = length(find(G(...
2021-01-04 17:26:30 118
原创 項目三:LT短信源数据包数目情况下,譯碼性能
clc;clear;m = 500;K = 500;overhead = [550 600 650 700 750 800];i_ii =10000;%0;for iii = 1:6 n = overhead(iii); for i_i = 1:i_ii if mod(i_i,100)==0 i_i end ...
2021-01-04 17:26:08 70
原创 fountain_BP_Gauss_decode
%*****************************************************************% 编码——BPSK——AWGN——高斯白(0,cigema平方)% mcisj :从编码比特ci到信息比特sj的LLR% msicj :从信息比特si到编码比特cj的LLR%%****************************************...
2021-01-04 17:24:47 152
原创 喷泉码raptor待续
% %**********生成矩陣G開始************* Pgenerate = rand(1,n); label = zeros(1,n); for j = 1:length(Pgenerate) for i = 1:length(sum_1) if Pgenerate(j)&l...
2021-01-04 17:21:15 816
原创 Ergodic_Capacity_CDF
clear all;close all;clc SNR_dB=10; SNR_linear=10.^(SNR_dB/10); N_iter=5000; sq2=sqrt(0.5); grps=['b:';'b-']; for Icase=1:2 if Icase==1 nT=2; nR=2; else nT=4;...
2021-01-04 17:20:51 383
原创 无线信道测量第一个程序
rnum = 12;freqpoint = 1601;rgetpoint = 100;for i_file = 1:rnum %接收位置点数 path_file = 'D:\research\SSM\1\'; file1 = dir(fullfile(path_fi...
2021-01-04 17:20:27 260
原创 SAGE算法1.0
clcclear%*****************理论对比**************************************************% load('coml_11.mat')% H = coml; %*******************实测数值***********************************************load('C:\...
2021-01-04 17:20:11 1338
原创 LT_BP_Gauss
clc;clear;K = 5000;m = K;u = 0;degree = 100; %迭代次數sigema =sqrt(0.5); % %(0.5/r)*10^(-0.1*(Eb_N0));% Eb_N0 = 1 ; %单位:dB % n = 10;% r = m/n;%sum_1= [0.007969,0.501539,0.6677...
2021-01-04 17:10:54 71
原创 LT_柱状图
clc;clear;K = 10000;%%%%%%%%%%%%%m = K;c = 0.01;o = 0.5;x = ones(1,m);S = c*sqrt(K)*log(K/o); for cou_num = 1:100 %½yÓ‹ %%%%%%%%% x_1 = rand(1,m)<0.5 ; %ÐÅÔ´´...
2019-08-23 21:43:13 110
原创 求矩阵秩
function [FF]=zhu_rank(H)[m,n] = size(H);a = H;for j = 1:m for i = j:m-1 for k = i+1:m if a(i,j)==0 && a(k,j)==1 ...
2019-08-13 09:22:09 842
原创 冒泡法排序
clear;clc%************冒泡法排序開始**************A = [3 1 4 7];n = length(A)for i = 1:n-1 %定值 for j = i+1:n %查找最小值 if A(i)>A(j) %交換 t = A(i); ...
2019-07-25 21:06:30 148
原创 選擇法排序
clear;clc%************選擇法排序開始**************A = [6 7 9 0 4 9 3 3 9 5];n = length(A)for i = 1:n %定值 for j = i:n %查找最小值 if A(i)>A(j) %交換 t = A(i);...
2019-07-25 21:04:09 112
原创 判斷素數
clear;clc%************判斷一個數是否是最小公倍數開始**************input =115if input <= 1 disp('請輸入大於2的整數');elseif input==2 disp('素數:yes');else n = input-1; for i = 2:n record(i-1) = ...
2019-07-25 20:48:10 78
原创 求最大公約數
clc%************求兩個整數最大公約數,最小公倍數開始**************m = 188;n = 188;while(1) r = mod(m,n); if r == 0 n disp('是最大公約數'); break; else m = n; n = r; ...
2019-07-25 20:22:32 86
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人