|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
, W9 G% Y( n7 L; E7 zfunction [chan_SR,chan_RD,dis_SR,dis_RD] = channel_settingmx(n,exp_path)' {% `- n! J6 b2 ]
%%距离和路径衰减
( ]& F4 Q7 y: q5 v+ y$ Us=[0 0];; q" N& W$ v$ C" e
d=[10 0];& O, D; K, n; e2 X. U; \1 ^, _# u
r=[n 0];
% H% O8 H$ ~5 ]6 w9 [/ }dis_SR=norm(r-s);
) [0 R5 S5 u3 ^0 J1 o& o1 `dis_RD=norm(d-r);
/ l! _/ Q: C% w% }, dpath_SR = 1/(dis_SR^exp_path);
. P3 I8 w& |/ hpath_RD = 1/(dis_RD^exp_path);0 V/ ]4 W. ]4 A1 `4 l: y
) t6 L6 x$ `1 A7 ^" |
%%定义信道增益
+ m3 Y1 l7 e+ c& k+ J% `chan_SR_fast = (1 + sqrt(-1))/sqrt(2);2 e* E' K. f2 [8 ~5 N4 I9 W
chanPW_SR_fast = abs(chan_SR_fast)^2;" s% X! x) `4 j6 H# ~
chan_SR = path_SR*chanPW_SR_fast;1 u: D: e# @8 K$ {7 j
( Y0 f4 M' _& }) t/ {: `; Z
chan_RD_fast = (1 + sqrt(-1))/sqrt(2);& ~4 b& R' E4 |3 e' ~" N% g
chanPW_RD_fast = abs(chan_RD_fast)^2;5 Z7 x8 w# x% U$ K5 a9 [4 a& y
chan_RD = path_RD*chanPW_RD_fast;+ a1 Q: \7 [( m) L
end
2 `& ?/ L' S: a# P: a' s! Z7 u! j* a f5 C
w=20e6;%%频带宽度
& Q; ]# ^" _# B/ W3 Zep=3;%%路径衰减系数- P) N4 y) q% N# p" L* s* r" y7 K; F
N=4.0e-21;%%高斯白噪声功率
6 a: o6 E) _; e3 N/ LPr=100;3 `1 Y3 P1 D0 W" L4 z: i
Ps=100;
+ x6 H& f9 j. J1 u* t0 ^1 Q: I- k% k# Z/ u+ x- W, l
n=1:9 Y& r9 f1 e4 G' W3 a. J1 l& h
[hsr,hrd]=channel_settingmx(n,3)
( Y& X/ y5 [& w9 tfpl=(Ps*Pr*(hsr^2)*(hrd^2))/((Ps*(hsr^2)+Pr*(hrd^2)+N)*N);
: F7 M8 u2 z8 G2 y; J+ X: VRsr=(1/2)*w*log2(1+((Ps*hsr^2)/N));& T; v+ r' z% h7 D. B2 ^9 p
Rrd=(1/2)*w*log2(1+((Ps*hrd^2)/N)+fpl);
6 T9 o# ]+ f# F6 w) T# p3 dRsd=(Rsr+Rrd)/2;
0 s3 H& j0 N9 h+ S3 G4 ~7 SRsd1=Rsd
8 |' y8 g/ @* gplot(n,Rsd1)6 W% w/ L# m. V4 G% Q1 t, m1 {
' J- I5 X& ` {$ p2 w
. x/ @& u) E' t3 Z& U各位大佬 本人对MATLAB的矩阵思维不是很理解,这个程序报了矩阵维度错误,请问该怎样修改? |
|