clear?all;%%%%?一列情況下的LS?%%%
clc;??
carrier_count=128;%
bits_per_symbol=2;%用4進制調制
symbols_per_carrier=1;%%
I=10;%%
pilot=[1:I:carrier_count]%%%%%%%%%%%%%%%%%
if?max（pilot）~=carrier_count
???pilot_index=[pilotcarrier_count];
???else
???pilot_index=pilot;
end
M=?size（pilot_index2）%%%%%%%%%%%%%%%%%%%%

N=carrier_count;%總載波數?128
N_number=carrier_count*symbols_per_carrier*bits_per_symbol;%總二進制數
fd=0.01;%假設多普勒頻移為
Ts=1;%采樣周期
CP=11;%保護時間
L=10;
SNR_number=7;%SNR數量
err_rate=ones（1SNR_number）;
mse=ones（1SNR_number）;

h=zeros（L1）;
W_h=fft（eye（NL））/sqrt（N）;

%--------------------生成信號----------------------------------------------
????X_msg=randint（1N4）;
%--------------------QPSK調制----------------------------------------------
????X=pskmod（X_msg4）;????
%%%%%%%??串并變換?%%%%%%%%%%%
????X1=X.‘;
????X2=ifft（X1）;
%--------------------循環前綴----------------------------------------------
????for?i=1:carrier_count;
???????x（i+CP1）=X2（i1）;?
????end
????for?i=1:CP;
????????x（i1）=X2（i+carrier_count-CP1）;
????end??????????????????????????????????
%%%%%%%??并串變換?%%%%%%%%%%%
????X3=x.‘;
for?n=1:50;????
%-------------------?生成信道----------------------------------------------
????K=6;%稀疏度為4
????delay=[1?3?5?7?8?9];%各徑延遲
???%?a=[0?-0.12?-0.78?-1?-1.25?-2];
???%?a=[0?-0.87?-2.17?-3.04?-3.91?-4.34];?
???%a=[0?-1.87?-4.17?-7.04?-9.91?-10.34];?
??a=[0?-2.87?-7.17?-12.04?-17.91?-19.34];?
????RLchannel=rayleighchan（Tsfddelaya）;
????RLchan.ResetBeforeFiltering=0;?%必須設置為0這樣多徑信道對符號的拖尾就可以算到下一個符號上?
?????RLchannel.StorePathGains=1;
????