clc;
clear?all;
m=250;n=[3005001000];%校驗矩陣的n是碼長
frame=500;
Eb_N0=[0:0.5:5];

for?ii=1:3
R（ii）=（n（ii）-m）/n（ii）
H（ii）=getH（mn（ii））???????%load?H;?H-->G
[G（ii）valid]=H2G（H（ii））;?
while?valid==0???????????????????%valid作為校驗矩陣是否為滿秩的標志若為非滿秩（valid=0）則返回重新利用1A隨機構造
H（ii）=getH（mn（ii））;?????????????????????%構造校驗矩陣H
[G（ii）valid]=H2G（H（ii））;????????????????%將校驗矩陣H轉化為生成矩陣G
end
end




????for?i=1:length（Eb_N0）
?????EN=10^（Eb_N0（i）/10）
????sigma=sqrt（1./（2*EN*R））;
????ber（i）=0;

??????for?num=1:frame??
????????num
????????x?=?（sign（randn（1size（G1）））+1）/2;?%?random?bits
????????y?=?mod（x*G2）;?????????????????????%?coding?
????????bpskmod?=bpsk（y）;??????????????????????????%BPSK?modulation
????????
????????z=bpskmod?+?sigma*randn（1length（bpskmod））;???%?AWGN?transmission
????????f1=1./（1+exp（-2*z/（sigma^2）））;?????????%?likelihoods
????????f0=1-f1;
????????[z_hat?success?k]?=?ldpc_decode（zf0f1H）;
????????x_hat?=?z_hat（size（G2）+1-size（G1）:size（G2））;%譯碼
????????x_hat?=?x_hat‘;???????????????????????????%轉置，保證x_hat與x的矩陣維數一致?
????????%erro=find（x_hat~=x）;?????????????????????%尋找錯誤信息位
????????err（num）=length（find（x_hat~=x））??????????%求出錯誤信息位位數
????????
????????end???%for?num
????????ber（i）=sum（err）/（length（x）*frame）????????%計算比特誤碼率BER
????????
????end?%for?i
???

???????
semilogy（Eb_N0ber‘b-o‘）;
xlabel（‘Eb/No（dB）‘‘fontweight‘‘bold‘）;
ylabel（‘BER‘‘fontweight‘‘bold‘）;
%title（‘ldpc在awgn下性能仿真‘‘fontsize‘12‘fontweight‘‘bold‘‘fontname‘‘黑體‘）
%c=legend（‘r=3/4L=512‘1）;
hold?on
grid?on;
%hold?off;
?%for?j

????????
????????
????????
????????
????????
????????
????????
????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1731??2011-05-22?19:52??LDPC_程序-50\1.asv

?????文件???????1767??2011-05-22?21:21??LDPC_程序-50\1.m

?????文件????????347??2011-05-20?00:22??LDPC_程序-50\bpsk.m

?????文件??????64037??2011-05-21?21:32??LDPC_程序-50\duibi.jpg

?????文件???????3302??2011-05-22?20:05??LDPC_程序-50\frame.asv

?????文件??????12010??2011-05-20?21:04??LDPC_程序-50\G.mat

?????文件????????333??2010-05-09?11:31??LDPC_程序-50\getG.asv

?????文件???????4255??2011-05-22?15:54??LDPC_程序-50\getH.asv

?????文件???????4219??2010-05-08?10:14??LDPC_程序-50\getH.m

?????文件???????8440??2011-05-22?22:09??LDPC_程序-50\H.mat

?????文件????????847??2010-05-12?19:53??LDPC_程序-50\H2G.m

?????文件???????3523??2011-05-20?20:46??LDPC_程序-50\ldpc_decode.asv

?????文件???????3523??2011-05-21?22:03??LDPC_程序-50\ldpc_decode.m

?????文件???????1457??2011-05-22?21:57??LDPC_程序-50\ldpc_decodeber.asv

?????文件???????1465??2011-05-22?22:01??LDPC_程序-50\ldpc_decodeber.m

?????文件??????????0??2011-05-21?17:38??LDPC_程序-50\Untitled.asv

?????文件??????63282??2011-05-21?17:13??LDPC_程序-50\untitled.jpg

?????文件??????65922??2011-05-20?22:30??LDPC_程序-50\圖.jpg

?????文件?????236278??2010-05-09?08:59??LDPC_程序-50\圖r=0.5\252_504_1000.bmp

?????文件?????252518??2010-05-09?11:27??LDPC_程序-50\圖r=0.5\252_504_10000（補充）.bmp

?????文件?????236278??2010-05-09?11:00??LDPC_程序-50\圖r=0.5\252_504_10000.bmp

?????文件?????236278??2010-05-11?12:48??LDPC_程序-50\圖r=0.5\256_512_20000（補充）.bmp

?????文件???????6784??2010-05-15?11:11??LDPC_程序-50\圖r=0.5\256_512_3000.fig

?????文件?????236278??2010-05-13?12:15??LDPC_程序-50\圖r=0.5\256_512_5000.bmp

?????文件???????8963??2010-05-14?18:24??LDPC_程序-50\圖r=0.5\256_512_5000.fig

?????文件?????252518??2010-05-09?15:26??LDPC_程序-50\圖r=0.5\96_192_10000.bmp

????..A.SH.?????27136??2010-05-13?15:03??LDPC_程序-50\圖r=0.5\Thumbs.db

?????文件???????6776??2010-05-15?16:49??LDPC_程序-50\圖r=0.5\untitled.fig

?????文件???????6740??2010-05-16?18:37??LDPC_程序-50\圖r=0.667\1.fig

?????文件???????6790??2010-05-14?20:52??LDPC_程序-50\圖r=0.667\2.fig

............此處省略38個文件信息