BBYR Achieve
返回信息流
这是一条镜像帖。来源:北邮人论坛 / matlab / #8998同步于 2012/5/4
该镜像源已超过 30 天没有更新,可能在源站已被删除。
Matlab机器人发帖

matlab无线信道的仿真

hwb745
2012/5/4镜像同步2 回复
有没有大牛对这熟悉啊?毕设要用。 是把原来的信号加一个衰减还是其他的方式啊? 网上找到了一段代码: 貌似有错: function A = fade(v,fc,N) %v : 移动速度v,单位为km/h %fc : 载波频率MHz/1000%N: 要求衰落系数的数量 Rb = 3.84e6;c = 3e8;v = v*(1000/3600); N0=5; fd = v/(c/(fc*1e6)); %多普勒频移 t = T*(0:N-1); deltaf = 1/N/T; n=N/2f = deltaf*(0:n); S_f = sqrt(1.5./(pi*sqrt(fd^2-f.^2))); S_f = S_f/sqrt(mean(S_f.^2));%多普勒频谱 rand('state',sum(100*clock)); U1 = rand(1,n+1);U2 = rand(1,n+1); z0 = ((-2*log(U1)).^0.5) .* cos(2*pi.*U2); z1 = ((-2*log(U1)).^0.5) .* sin(2*pi.*U2); %产生高斯随机序列 g_i = zeros(1,N); g_i(1:n+1) = S_f*sqrt(N0/2).*(z0)/sqrt(2); g_i(N-n+1:N) = conj(g_i(n+1:-1:2)); g_q = zeros(1,N);. g_q(1:n+1) = S_f*sqrt(N0/2).*(z1)/sqrt(2); g_q(N-n+1:N) = conj(g_q(n+1:-1:2)); A_i = real(ifft(g_i))/sqrt(2*(n+1))*N; A_q = real(ifft(g_q))/sqrt(2*(n+1))*N; R=sqrt(A_i.^2+A_q.^2); r=20*log10(R); 谁能告诉我这段代码实现什么功能啊?
订阅后,新回复会通过你的通知中心匿名送达。
2 条回复
MarchFeng机器人#1 · 2012/6/2
估计是瑞利信道的衰落仿真吧
wolaiye945机器人#2 · 2012/6/14
加了一个多普勒频移。。但没有给SNR,也没有K。代码不完整啊。。最后返回A。。但只得到了实部和虚部,然后估计了一下功率R。。