|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
$ ^/ [/ U; B1 X
function [chan_SR,chan_RD,dis_SR,dis_RD] = channel_settingmx(n,exp_path)
1 b; {3 E% }0 y$ j( }/ ]3 U%%距离和路径衰减6 r! u, f; {- d3 a. Y8 l7 @
s=[0 0];6 I% J* Q% M0 G5 ^2 ]+ S$ n* _" _) B
d=[10 0];
|) a& [/ \9 N3 {# v, ~3 k9 ur=[n 0];
+ E6 f z n4 j$ gdis_SR=norm(r-s);
1 a% Z& v; ^( a9 Odis_RD=norm(d-r);
6 j& \6 d" G9 U4 g2 V+ x7 Wpath_SR = 1/(dis_SR^exp_path);( J! f. w$ h6 r( o
path_RD = 1/(dis_RD^exp_path);
& q7 d( R" U8 f( w" Z5 N( e8 r2 e1 Q! I
%%定义信道增益
; y0 U+ q; m% d z( o; o3 @- |( tchan_SR_fast = (1 + sqrt(-1))/sqrt(2);
1 c, E% i! o( \+ WchanPW_SR_fast = abs(chan_SR_fast)^2;
+ X5 K8 @/ `# a' Mchan_SR = path_SR*chanPW_SR_fast;# F5 |/ V0 t* o
! c- p1 z) o$ u, s/ ]6 S
chan_RD_fast = (1 + sqrt(-1))/sqrt(2);
" G o+ l3 e8 g( z$ Q/ ^: h& Z0 UchanPW_RD_fast = abs(chan_RD_fast)^2;
" w) j/ T2 V8 x- Dchan_RD = path_RD*chanPW_RD_fast;7 C( l2 n: V t: p( d: r9 B
end S* a3 ~+ L: q3 g/ k, | L* |
1 W7 f$ S$ z" C; h- E* fw=20e6;%%频带宽度
9 U5 h2 m) \ J$ D; z5 ^5 L8 E' c! jep=3;%%路径衰减系数- m; G& d' n( m2 b5 N; b
N=4.0e-21;%%高斯白噪声功率
- {) r5 O, N# u* G/ bPr=100;
9 l) C$ u) j, j/ l# IPs=100;# p+ A1 m2 @ ^+ |8 c' v6 W6 k$ |
! M- g/ T0 t( v) M- vn=1:9
) \8 e8 y: ~9 y" @+ a: M: J: V[hsr,hrd]=channel_settingmx(n,3)
+ c( U4 @7 t9 ?' m( \; g; pfpl=(Ps*Pr*(hsr^2)*(hrd^2))/((Ps*(hsr^2)+Pr*(hrd^2)+N)*N);
" J# z, O7 k) C+ B% `8 IRsr=(1/2)*w*log2(1+((Ps*hsr^2)/N));. P8 E- t. R. y, q$ l$ Q, G) s. b7 S
Rrd=(1/2)*w*log2(1+((Ps*hrd^2)/N)+fpl);( I' d! H6 a/ Y% S
Rsd=(Rsr+Rrd)/2;/ Y6 P: p) I9 M% ?
Rsd1=Rsd$ F! o; x7 R* h3 w! U" Z' Q
plot(n,Rsd1)3 `0 m$ Y" Z$ _8 H; z) [
3 p4 u2 h& g7 z8 J0 R9 J9 u4 ?+ ] y+ [' Y S; r# G2 K
各位大佬 本人对MATLAB的矩阵思维不是很理解,这个程序报了矩阵维度错误,请问该怎样修改? |
|