|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
( K. a2 V4 h' Y' P* k
function [chan_SR,chan_RD,dis_SR,dis_RD] = channel_settingmx(n,exp_path)
# z4 ?1 `9 T$ q6 u' _%%距离和路径衰减
* Q9 D" E- \3 _; t5 B% M2 ]s=[0 0];; b. [) q9 f2 e1 l0 v8 ^. V
d=[10 0];. }, O+ Y$ s& b, |3 V6 @
r=[n 0];. D L% P% ^8 N$ x4 G+ }3 Y# G/ I
dis_SR=norm(r-s);7 b/ ], _8 }+ y: ?6 y
dis_RD=norm(d-r);
+ q, J6 h* @; i0 X* Mpath_SR = 1/(dis_SR^exp_path);
& D6 k2 |3 u1 ^: l2 A7 ~path_RD = 1/(dis_RD^exp_path);
/ _8 b) m' f# t
. X7 d- b# Y6 I" l%%定义信道增益, z7 q4 I6 ?4 k$ F3 [0 z0 C
chan_SR_fast = (1 + sqrt(-1))/sqrt(2);% D0 K' v7 f; @* t- s9 W
chanPW_SR_fast = abs(chan_SR_fast)^2;
1 V7 l' i1 N( B' }chan_SR = path_SR*chanPW_SR_fast;
0 A1 A3 U ] ] p2 g0 C @3 U$ t; E7 ?5 i
chan_RD_fast = (1 + sqrt(-1))/sqrt(2);% z1 ?1 k7 @! d6 B# c5 K
chanPW_RD_fast = abs(chan_RD_fast)^2;: N6 m( g6 n8 i, E9 ~- p! R3 q: u# |
chan_RD = path_RD*chanPW_RD_fast;7 E% o: @1 H H: i6 f8 c X! s
end
$ d6 `: x1 W. G, i! _
; G* G: ^) l+ [; Tw=20e6;%%频带宽度
; t. G l7 m6 f8 E# xep=3;%%路径衰减系数
/ v( l9 |9 n+ e1 HN=4.0e-21;%%高斯白噪声功率
& U8 R6 @ j' TPr=100;/ `4 t/ C" o: A/ I; ?* H7 v
Ps=100;
0 b8 T. s' d+ E/ S% [5 X* r) K/ ~) v, Y
n=1:9
$ c% U2 k; Z+ F* p* Z[hsr,hrd]=channel_settingmx(n,3)
( M$ B- t8 Q. V6 q+ ^fpl=(Ps*Pr*(hsr^2)*(hrd^2))/((Ps*(hsr^2)+Pr*(hrd^2)+N)*N);
1 E" I+ Z w' c, q ?- s1 y& MRsr=(1/2)*w*log2(1+((Ps*hsr^2)/N));
c+ X3 K# I# z; U" Z0 Z& gRrd=(1/2)*w*log2(1+((Ps*hrd^2)/N)+fpl);
6 p# b# W3 T. l$ tRsd=(Rsr+Rrd)/2;
; P. D. O+ S4 p1 W5 @ R6 PRsd1=Rsd
* U( T. P* L& L" P$ [, Uplot(n,Rsd1), W7 w6 X" f5 j+ d
5 u' Q9 K2 Y2 J4 F8 h* s* P4 Q8 ^7 E8 `# B# a
各位大佬 本人对MATLAB的矩阵思维不是很理解,这个程序报了矩阵维度错误,请问该怎样修改? |
|