clear?all;
Numusers?=?1;
Nc?=?16;????%擴頻因子
ISI_Length?=?1;??%碼元長度
EbN0db?=?2;??????%信噪比
Tlen?=?10;???????%數據長度
?
%我們假設出的每徑功率因子
power_unitary_factor1?=?sqrt（?6/9?）;
power_unitary_factor2?=?sqrt（?2/9?）;
power_unitary_factor3?=?sqrt（?1/9?）;
s_initial?=?randsrc（?1?Tlen?）;?????%產生數據源1或-1
figure（1）
stairs（s_initial‘Color‘‘r‘）;
ylim（[-1.51.5]）;
title（‘信號源序列‘‘Color‘‘b‘）
xlabel（‘數據長度‘）
ylabel（‘幅度‘）
%根據擴頻因子等于16產生walsh矩陣
wal2?=?[?1?1;?1?-1?];
wal4?=?[wal2?wal2;?wal2?wal2*（-1）];
wal8?=?[wal4?wal4;?wal4?wal4*（-1）];
wal16?=?[wal8?wal8;?wal8?wal8*（-1）];
%擴頻
s_spread?=?zeros（?Numusers?Tlen*Nc?）;
ray1?=?zeros（?Numusers?2*Tlen*Nc?）;
ray2?=?zeros（?Numusers?2*Tlen*Nc?）;
ray3?=?zeros（?Numusers?2*Tlen*Nc?）;

for?i?=?1:Numusers
????x0?=?s_initial（?i:?）.‘*wal16（?8:?）;
????x1?=?x0.‘;
????s_spread（?i:?）?=?（?x1（:）?）.‘;
end
%將擴頻后的輸出重復然后進行延遲
ray1（?1:2:2*Tlen*Nc?-?1?）?=?s_spread（?1:Tlen*Nc?）;
ray1（?2:2:2*Tlen*Nc?）?=?ray1（?1:2:2*Tlen*Nc?-?1?）;
?
%產生第二第三徑信號
ray2（?ISI_Length?+?1:2*Tlen*Nc?）?=?ray1（?1:2*Tlen*Nc?-?ISI_Length?）;
ray2（?2*ISI_Length?+?1:2*Tlen*Nc?）?=?ray1（?1:2*Tlen*Nc?-?2*ISI_Length?）;
?
for?nEN?=?1:length（?EbN0db?）
????en?=?10^（?EbN0db（nEN）/10?）;?????%信噪比轉換
????sigma?=?sqrt（?32/（2*en）?）;
????%收到信號
????demp?=?power_unitary_factor1*ray1+...
???????????power_unitary_factor2*ray2+...
???????????power_unitary_factor3*ray3+...
???????????（?rand（?12*Tlen*Nc?）+randn（?12*Tlen*Nc?）*i?）*sigma;
????dt?=?reshape（?demp32Tlen?）‘;
????%將walsh碼重復
????wal16_d（1:2:31）?=?wal16（81:16）;
????wal16_d（2:2:32）?=?wal16（81:16）;
????%解擴后第一徑輸出
????rdata1?=?dt*wal16_d（1:）.‘;
????%將walsh碼延遲
????wal16_delay1（12:32）?=?wal16_d（11:31）;
????%解擴后第二徑輸出
????rdata2?=?dt*wal16_delay1（1:）.‘;
????%將walsh碼延遲
????wal16_delay2（13:32）?=?wal16_d（11:30）;
????wal16_delay2（11:2）?=?wal16_d（131:32）;
????%解擴后第三徑輸出
????rdata3?=?dt*wal16_delay2（1:）.‘;
????
????p1?=?rdata1‘*rdata1;
????p2?=?rdata2‘*rdata2;
????p3?=?rdata3‘*rdata3;
????p?=?p1?+?p2?+?p3;
????u1?=?p1/p;
????u2?=?p2/p;
????u3?=?p3/p;
????
????%最大比合并
????rd_m1?=?real（?rdata1*u1+rdata2*u2+rdata3*u3）;
????%等增益合并
????rd_m2?=?（real（rdata1+rdata2+rdata3））/3;
????%選擇式合并
????u?=?[u1u2u3];
????maxu?=?max（u）;
????if（maxu==u1）
????????rd_m3?=?real（rdata1）;
????else
????????if（maxu==u2）
????????????rd_m3?=?real（rdata2）;
????????else?rd_m3?=?real（rdata3）;
????????end
????end
?
????%輸出數據
????r_Data1?=?sign（rd_m1）‘;
????r_Data2?=?sign（rd_m2）‘;
????r_Data3?=?sign（rd_m3）‘;
????
end
figure（2）
subplot（311）
stairs（r_Data1‘Color‘‘r‘）;
ylim（[-1.51.5]）;
title（‘最大比合并方式‘‘Color‘‘b‘）
xlabel（‘碼元長度‘）
ylabel（‘幅度‘）
subplot（312）
stairs（r_Data2‘Color‘‘r‘）;
ylim（[-1.51.5]）;
title（‘等增益合并‘‘Color‘‘b‘）
xlabel（‘碼元長度‘）
ylabel（‘幅度‘）
subplot（313）
stairs（r_Data3‘Color‘‘r‘）;
ylim（[-1.51.5]）;
title（‘選擇式合并‘‘Color‘‘b‘）
xlabel（‘碼元長度‘）
ylabel（‘幅度‘）
figure（3）
subplot（311）
stairs（s_initial‘Color‘‘r‘）;
ylim（[-1.51.5]）;
hold?on
stairs（r_

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

?????文件???????3777??2018-05-06?10:21??新建文件夾?（4）\fashe.m

?????文件???????3183??2017-05-28?23:56??新建文件夾?（4）\jieshou.m

?????目錄??????????0??2018-05-06?10:26??新建文件夾?（4）

-----------?---------??----------?-----??----

?????????????????6960????????????????????3