close all;

SNR = 3:3:18;
figure;
subplot(1,2,1)
% Channel = 1
BER1=[0.328788690476190,0.219756944444444,0.131839285714286,0.0653154761904762,0.0350436507936508,0.0189315476190476];
semilogy(SNR,BER1,'+-k','linewidth',2);
hold on;

% Channel = 2
BER2=[0.350351190476190,0.215301587301587,0.0998961640211640,0.0385562169312169,0.0112863756613757,0.00340013227513228];
PER2=[0.986000000000000,0.879333333333333,0.596000000000000,0.291666666666667,0.111333333333333,0.0293333333333333];
semilogy(SNR,BER2,'^-b','linewidth',2);


% Channel = 4
BER4=[0.368680555555556,0.218835317460317,0.0813068783068783,0.0177870370370370,0.00172156084656085,0.000355357142857143];
PER4=[0.998666666666667,0.947333333333333,0.671666666666667,0.270333333333333,0.0563333333333333,0.00933333333333333];
semilogy(SNR,BER4,'s-g','linewidth',2);

% Channel = 8
BER8 = [0.384902777777778,0.230017857142857,0.0675138888888889,0.00725793650793651,0.000438492063492064,1.66948853615520e-05];          %無窮bit量化
semilogy(SNR,BER8,'d-r','linewidth',2);
grid on;

xlabel('SNR/dB');
ylabel('BER');
title('\fontsize{12}多徑信道，Unquant軟判決');
legend('\fontsize{12}Flat Fading','\fontsize{12}L = 2','\fontsize{12}L = 4','\fontsize{12}L = 8','Location','SouthWest');
hold off;


subplot(1,2,2);
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');
axis([0 65 0 2]);
legend('\fontsize{12}L = 2','\fontsize{12}L = 8','\fontsize{12}L = 16');
title('\fontsize{12}Channel State Information');
xlabel('\fontsize{12}子載波');
ylabel('\fontsize{12}信道幅度H_{f}');
hold off;