|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
2 L% _ c! @* r9 H, t8 Q
function [chan_SR,chan_RD,dis_SR,dis_RD] = channel_settingmx(n,exp_path)/ O0 o1 E' ?) C5 ]
%%距离和路径衰减9 }3 [- u/ a; R
s=[0 0];4 j, t+ l) [5 ^
d=[10 0];- o$ K: `3 k0 E" u
r=[n 0]; D9 w5 U" B0 \) ]* v$ a
dis_SR=norm(r-s);
3 U/ `8 e' }8 ddis_RD=norm(d-r);- ~/ P6 ^* T' E3 z' w' s5 a
path_SR = 1/(dis_SR^exp_path);4 _ h3 o h# ^3 P0 E6 f# |. A# J5 i
path_RD = 1/(dis_RD^exp_path);/ J% m z; G( a9 n
& Z7 ?% `& r0 B- M/ N
%%定义信道增益
; _8 F# Y7 j1 Y `: H' m% [- _chan_SR_fast = (1 + sqrt(-1))/sqrt(2); j0 o6 j; F! @! Y2 j
chanPW_SR_fast = abs(chan_SR_fast)^2;4 B9 B7 q7 I+ V6 M$ m
chan_SR = path_SR*chanPW_SR_fast;7 G1 ~" {9 D, N+ Z
& e+ G1 `8 M8 f6 Q, b/ t% c+ |) p
chan_RD_fast = (1 + sqrt(-1))/sqrt(2);1 ~, E' b# o; j( F' o. ]
chanPW_RD_fast = abs(chan_RD_fast)^2;7 _( T1 D! b. w( R7 a F
chan_RD = path_RD*chanPW_RD_fast;
7 s) r$ [( |6 Y* f+ L0 Mend
0 x. u9 p9 {: k( X
. n0 B% Z( p; I4 i, ^ e) hw=20e6;%%频带宽度3 v5 x: l. ~4 ]$ z4 z
ep=3;%%路径衰减系数- R5 W! w4 W8 d
N=4.0e-21;%%高斯白噪声功率
' ]% J$ Q% M# o- x. ]/ tPr=100;7 h4 z n3 y6 M
Ps=100;$ e q# I7 D8 i j5 _" {! g; B2 B
$ Z+ D( H# c5 `. o& cn=1:9' o7 q7 A1 L5 {/ T. G+ ]/ A$ Q
[hsr,hrd]=channel_settingmx(n,3). o7 H. v8 o& X3 J( Y# `% ~" r# [. B
fpl=(Ps*Pr*(hsr^2)*(hrd^2))/((Ps*(hsr^2)+Pr*(hrd^2)+N)*N);1 \" y: Z- W2 ~7 [3 t
Rsr=(1/2)*w*log2(1+((Ps*hsr^2)/N));! ?* L! l1 c6 h" L! X& T$ C
Rrd=(1/2)*w*log2(1+((Ps*hrd^2)/N)+fpl);
3 {' m: l3 u2 i2 {5 \) z9 G! JRsd=(Rsr+Rrd)/2;! c/ U5 f% N" e8 S, H' w" O
Rsd1=Rsd
% U: \4 H7 J m$ f# Hplot(n,Rsd1)
~! z% b( W3 p B: `* [
. f1 K5 q. ~( R1 `5 r1 w B, n' y5 Y6 l. i9 [7 w2 n( F
各位大佬 本人对MATLAB的矩阵思维不是很理解,这个程序报了矩阵维度错误,请问该怎样修改? |
|