clear?all
close?all
clc?
SNR_db?=[101520];
SNR?=10.^（SNR_db/10）;
color?=?[‘k‘;‘r‘;‘b‘;‘k‘;‘b‘;‘r‘];
shape1?=?[‘-s‘;‘-d‘;‘-p‘;‘-h‘];
shape2?=?[‘-*‘;‘-o‘;‘-<‘;‘-+‘];
%wucha?=?[00.010.040.090.160.250.360.490.640.811]

wucha?=?[1e-41e-31e-21e-11]
wuchashu=[12345]
y=1e2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%等收發天線組合，不同信噪比
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%發送天線數目等于接收天線數目dt=dr
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

dt=[1?2?3?4];?%發射天線數
dr=[4?4?4?4];?%接受天線數
for（i=1:length（SNR））
??for（k=1:length（wucha））
????for（x=1:y）
????[Cf（x）Cg（x）U（x）]=xindaorongliangbuwanquan（wucha（k）SNR（i）44）;?
????end
????Cwf（k）=sum（Cf）/y;?%注水算法信道容量
????Cdg（k）=sum（Cg）/y;?%等功率算法信道容量
????u（k）=sum（U）/y;????%返回不同信噪比時的注水線u
????Cwd1（k）=Cwf（k）-Cdg（k）;%注水-等功率增益差
????Cwd2（k）=Cdg（k）/Cwf（k）;%等功率/注水增益
??end
????f1=figure（1）;
????hold?on
????grid?on
????plot（wuchashuCwfshape2（i:）‘color‘‘k‘）????????????????%注水曲線
%????plot（SNR_dbCdgshape2（i:）‘color‘‘k‘）???%等功率曲線
?%???legend（‘注水3*3‘‘等功率3*3‘‘注水4*4‘‘等功率4*4‘‘注水5*5‘‘等功率5*5‘‘注水6*6‘‘等功率6*6‘‘Location‘‘northwest‘）
????????legend（‘SNR=10dB‘‘SNR=15dB‘‘SNR=20dB‘‘0.16‘‘Location‘‘northwest‘）
????xlabel（‘SNR/dB‘）;
????ylabel（‘Capacity/bits/s/Hz‘）;
????f2=figure（2）
????hold?on
????grid?on
???????plot（wuchashuCdg‘color‘‘k‘）
??
????
?%???plot（SNR_dbCwd1shape2（i:）‘color‘‘k‘）?????????????????%注水-等功率增益差曲線
%????legend（‘注水-等功率3*3‘‘注水-等功率4*4‘‘注水-等功率5*5‘‘注水-等功率6*6‘‘Location‘‘northeast‘）??%注水-等功率增益差曲線
legend（‘4發4收‘‘5發5收‘‘6發6收‘‘Location‘‘northeast‘）??%注水-等功率增益差曲線
????xlabel（‘SNR/dB‘）;
????ylabel（‘Capacity/bits/s/Hz‘）;
????f3=figure（3）
????hold?on
????grid?on
%????plot（SNR_dbCwd2shape2（i:）‘color‘‘k‘）?????????????????%注水-等功率增益比曲線
%????legend（‘等功率/注水3*3‘‘等功率/注水4*4‘‘等功率/注水5*5‘‘等功率/注水6*6‘‘Location‘‘northwest‘）??%注水-等功率增益比曲線
????????legend（‘EPA/WF4發4收‘‘EPA/WF5發5收‘‘EPA/WF6發6收‘‘Location‘‘northwest‘）??%注水-等功率增益比曲線
????xlabel（‘SNR/dB‘）;
????ylabel（‘Ratio/無單位‘）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-10-21?17:46??MIMO\
?????文件????????2428??2011-06-02?12:34??MIMO\buwuanquanxindao.m
?????文件????????6008??2011-05-28?23:11??MIMO\mimo?yuan.m
?????文件????????5905??2011-10-16?19:44??MIMO\mimo.m
?????文件????????3793??2011-05-31?20:07??MIMO\mimotianxianguding.m
?????文件????????2367??2011-05-31?18:39??MIMO\mimoxinzaobi.m
?????文件????????1667??2011-06-01?20:30??MIMO\sanzhongxitong.m
?????文件?????????784??2011-05-24?10:18??MIMO\xindaorongliang.m
?????文件?????????229??2011-05-25?19:33??MIMO\xindaorongliang1.m
?????文件?????????137??2011-05-26?00:49??MIMO\xindaorongliang2.m
?????文件?????????812??2011-06-02?12:11??MIMO\xindaorongliangbuwanquan.m