close all;
clc;

% % 平衰落 
% % 有編碼
% SNR = 3:3:18;
% BER_1=[0.329420634920635,0.215067460317460,0.127912698412698,0.0730079365079365,0.0422420634920635,0.0209880952380952];
% PER_1=[0.926000000000000,0.732000000000000,0.498000000000000,0.290000000000000,0.166000000000000,0.0940000000000000];
% 
% figure;
% semilogy(SNR,BER_1,'*-g','linewidth',2);
% hold on;
% BER_3 = [0.352567460317460,0.284365079365079,0.219750000000000,0.152115079365079,0.104226190476190,0.0674801587301587,0.0466230158730159,0.0297103174603175];
% semilogy(4:2:18,BER_3,'^-.g','linewidth',2);
% 
% semilogy(SNR,PER_1,'+-b','linewidth',2);
% grid on;
% xlabel('\fontsize{14}SNR/dB');
% ylabel('\fontsize{14}BER');
% 
% % without coding
% BER_2=[0.259912109375000,0.197408203125000,0.137386718750000,0.0886406250000000,0.0558398437500000,0.0302968750000000];
% PER_2=[1,1,0.994000000000000,0.962000000000000,0.788000000000000,0.546000000000000];
% semilogy(SNR,BER_2,'d-r','linewidth',2);
% semilogy(SNR,PER_2,'s-k','linewidth',2);
% 
% legend('\fontsize{14}BER with coding(Unquant)','\fontsize{14}BER with coding(Hard)','\fontsize{14}PER with coding','\fontsize{14}BER without coding','\fontsize{14}PER without coding','Location','SouthWest');
% 

% AWGN
% simulation in Eb/n0
figure;

BER_0 = [0.140981635066842,0.118997407465924,0.0977418537374869,0.0774530602925490,0.0586237372834044,0.0418927600464623,0.0278713278451503,0.0169667343687604,0.00924721374147442,0.00439033608735211,0.00175415061789273,0.000564706106481744,0.000138658688812619,2.42337854663159e-05,2.76320800168778e-06];
semilogy(0:14,BER_0,'-b','linewidth',2);

grid on;
hold on;

BER_3=[0.350591269841270,0.216067460317460,0.0655555555555556,0.00673412698412698,0.000202380952380952,1.12433862433862e-06];
semilogy(1:2:11,BER_3,'*-r','linewidth',2);

BER_2=[0.377821428571429,0.233333333333333,0.0452698412698413,0.00201984126984127,4.04383975812547e-05];
semilogy(-1:2:7,BER_2,'^-g','linewidth',2);

xlabel('\fontsize{14}Eb/n0(dB)');
ylabel('\fontsize{14}BER');
legend('\fontsize{12}Not coding','\fontsize{12}Hard Decision,Rc=0.5','\fontsize{12}Unquant Decision,Rc=0.5','Location','SouthWest');

% figure;
% channel_16=[0.863902726369044 - 0.0616236545751702i,1.06512737328433 - 0.593905603391635i,0.670119792016308 - 1.17774317983073i,...
%     -0.0477610843620582 - 1.23633668532691i,-0.512439451867941 - 0.827010762843347i,-0.625616195997298 - 0.368171670943697i,...
%     -0.581112725422812 + 0.0608189580893372i,-0.256846669219316 + 0.519671332966920i,0.480776633406690 + 0.620088954894494i,...
%     1.08348248325568 - 0.0469354669168924i,0.793946529209606 - 0.992767031757949i,-0.165029842349707 - 1.20532282453285i,...
%     -0.720333518375221 - 0.556945154178910i,-0.451248497547026 - 0.00470591793360009i,-0.121316975261849 - 0.131444200162212i,...
%     -0.364786916536001 - 0.345932509705114i,-0.705462491282661 - 0.00715763819009363i,-0.495436790134623 + 0.486572360664072i,...
%     -0.0681196238157121 + 0.478493514600307i,-0.0835596897956207 + 0.255171688633970i,-0.269845651130643 + 0.494373676878428i,...
%     0.121118371612717 + 0.924239253700895i,0.857605011944152 + 0.661500898751606i,0.894523541711494 - 0.194741434084709i,0.0947814840126928 - 0.477313066207145i,...
%     -0.347953444255392 + 0.259725504227689i,0.344005071125462 + 0.922604938988274i,1.31478623563640 + 0.438084469875479i,1.30492784579245 - 0.698124792018651i,...
%     0.396199689372622 - 1.21859890721692i,-0.325515472631402 - 0.856077035143908i,-0.414409111865024 - 0.410553036963963i,-0.424149399935030 - 0.327319173499613i,...
%     -0.727589031121530 - 0.160870778961843i,-0.876765264441552 + 0.457870370621013i,-0.382210649281151 + 1.09797482137475i,0.481538987127257 + 1.10326181856093i,...
%     0.992014413068645 + 0.445660012436776i,0.785142596717620 - 0.285325735542510i,0.137156684834479 - 0.503476945689422i,-0.340799603302686 - 0.103867575014349i,...
%     -0.217263531210172 + 0.458134557398054i,0.311776517042344 + 0.542519791124116i,0.512522888156559 + 0.0518057462760559i,-0.0154162446180450 - 0.276896278630157i,-0.625053642098652 + 0.270029468377727i,-0.302562962457561 + 1.26239481458011i,0.882638825685606 + 1.47207674271852i,1.66471400911845 + 0.511927378758850i,1.21719376428883 - 0.534260771926236i,0.265309833456126 - 0.538630381147646i,0.0369887324030674 + 0.201805920114371i,0.622846226337742 + 0.546816660320334i,1.09706980583487 + 0.128506097447887i,0.950918842992287 - 0.398187865186306i,0.576403701933507 - 0.517951344035582i,0.387178990546905 - 0.425983393952344i,0.277775634435312 - 0.380487288640205i,0.127705866077018 - 0.257531151817409i,0.141462810500491 - 0.00276037524390843i,0.382765565969342 + 0.102305330955160i,0.535165705380052 - 0.0483985170636607i,0.463849071618318 - 0.131964374907676i,0.513106650715634 - 0.00298323511918960i;];
% 
% channel_8=[1.26856028945971 - 0.0468173372575202i,1.08931030480661 - 0.513122260040985i,0.680399715493310 - 0.801148898471985i,0.190721986111827 - 0.815134784299697i,-0.207525939076146 - 0.567608440255340i,-0.388666767033586 - 0.168547195770737i,-0.320607064879906 + 0.225107521506543i,-0.0702744243281025 + 0.476398763130209i,0.232725575395658 + 0.522421543852763i,0.456256638764014 + 0.393061486497838i,0.522097992346710 + 0.188474300444481i,0.434225672205139 + 0.0278552758834841i,0.268300359782743 - 0.00585108434035564i,0.129761463352949 + 0.0998303150564655i,0.101818326471606 + 0.287924109780253i,0.208158242815568 + 0.465409978776939i,0.407048076937376 + 0.551142273226651i,0.617814018976620 + 0.514328922968052i,0.765002805304813 + 0.385114793779852i,0.817294163782667 + 0.233066193910057i,0.801286026448500 + 0.125427389606255i,0.783070546127633 + 0.0871669483295377i,0.826877608940645 + 0.0841323480364849i,0.952067897017070 + 0.0393261629349462i,1.11149328679248 - 0.124273488051722i,1.20461192952325 - 0.434864293658241i,1.12206136790449 - 0.841905636442379i,0.802440730316382 - 1.22256148149758i,0.274500489763766 - 1.42551241046334i,-0.337261032340175 - 1.33526841005483i,-0.850140862074246 - 0.929101120738848i,-1.09106376662543 - 0.299998770139855i,-0.967269835461710 + 0.367432001601180i,-0.509487807218740 + 0.861100186581451i,0.132315063912589 + 1.02073774774488i,0.740907361345681 + 0.799964049067732i,1.11148308933851 + 0.285050365720853i,1.12697166011866 - 0.337907716146302i,0.798427972892818 - 0.854135971163958i,0.254607692348898 - 1.10037316774176i,-0.313753431618960 - 1.02138852801123i,-0.730575642195801 - 0.680687031874010i,-0.899487302155500 - 0.223183548220484i,-0.830845496189020 + 0.191902592130053i,-0.621653446957339 + 0.457626431680960i,-0.400174370839840 + 0.556588403569984i,-0.262112170421599 + 0.554301924996666i,-0.227896381001468 + 0.554358736279060i,-0.239897555598405 + 0.638293711684735i,-0.199359966924434 + 0.819442865897970i,-0.0242304044069520 + 1.03291602915893i,0.300623570566247 + 1.16648346002909i,0.705692454227613 + 1.11771174890914i,1.06125925787722 + 0.850033453265324i,1.23380613694155 + 0.420430960635055i,1.14886273535274 - 0.0358768008119293i,0.831000615095164 - 0.358673119410408i,0.401337545269790 - 0.431122818737622i,0.0317477826018500 - 0.232867730358840i,-0.125460967137610 + 0.146561752411322i,-0.00217601252633193 + 0.541293056776283i,0.354946263738461 + 0.776241259385012i,0.801837073130448 + 0.739477284582257i,1.15544502542203 + 0.428892155206848i;];
% channel_2=[0.895896296566544 + 0.213343528723103i,0.911618847373584 + 0.181776499889498i,0.924171579906107 + 0.148820395459082i,0.933433604487804 + 0.114792600733158i,0.939315722758814 + 0.0800208219764219i,0.941761286704696 + 0.0448399304264887i,0.940746744207768 + 0.00958873730081571i,0.936281865866883 - 0.0253932691405343i,0.928409650901213 - 0.0597691930524213i,0.917205913044268 - 0.0932079754957822i,0.902778550416196 - 0.125387582716190i,0.885266506405906 - 0.155998107504948i,0.864838431570272 - 0.184744753774928i,0.841691059437084 - 0.211350675608032i,0.816047311853624 - 0.235559643432733i,0.788154152127421 - 0.257138511654910i,0.758280206634590 - 0.275879463977374i,0.726713177800985 - 0.291602014784414i,0.693757073370569 - 0.304154747316937i,0.659729278644645 - 0.313416771898634i,0.624957499887909 - 0.319298890169644i,0.589776608337976 - 0.321744454115526i,0.554525415212303 - 0.320729911618598i,0.519543408770953 - 0.316265033277713i,0.485167484859066 - 0.308392818312043i,0.451728702415705 - 0.297189080455098i,0.419549095195297 - 0.282761717827026i,0.388938570406539 - 0.265249673816736i,0.360191924136559 - 0.244821598981102i,0.333586002303455 - 0.221674226847914i,0.309377034478754 - 0.196030479264455i,0.287798166256577 - 0.168137319538251i,0.269057213934113 - 0.138263374045420i,0.253334663127073 - 0.106696345211815i,0.240781930594550 - 0.0737402407813992i,0.231519906012853 - 0.0397124460554753i,0.225637787741843 - 0.00494066729873904i,0.223192223795961 + 0.0302402242511942i,0.224206766292889 + 0.0654914173768672i,0.228671644633774 + 0.100473423818217i,0.236543859599444 + 0.134849347730104i,0.247747597456389 + 0.168288130173465i,0.262174960084461 + 0.200467737393873i,0.279687004094751 + 0.231078262182631i,0.300115078930385 + 0.259824908452610i,0.323262451063573 + 0.286430830285715i,0.348906198647032 + 0.310639798110416i,0.376799358373236 + 0.332218666332593i,0.406673303866067 + 0.350959618655057i,0.438240332699672 + 0.366682169462097i,0.471196437130088 + 0.379234901994620i,0.505224231856012 + 0.388496926576317i,0.539996010612748 + 0.394379044847327i,0.575176902162681 + 0.396824608793209i,0.610428095288354 + 0.395810066296281i,0.645410101729704 + 0.391345187955396i,0.679786025641591 + 0.383472972989726i,0.713224808084952 + 0.372269235132781i,0.745404415305359 + 0.357841872504709i,0.776014940094118 + 0.340329828494419i,0.804761586364097 + 0.319901753658785i,0.831367508197202 + 0.296754381525597i,0.855576476021903 + 0.271110633942138i,0.877155344244080 + 0.243217474215934i;];
% 
% plot(1:64,abs(channel_2),'--k');
% hold on;
% 
% plot(1:64,abs(channel_8),'-r');
% hold on;
% plot(1:64,abs(channel_16),'-.b');
% 
% 
% 
%