专业论坛经典问题

信道 专栏收录该内容
4 篇文章 0 订阅

【转】专业论坛经典问题

  (2010-06-30 16:08:24)
标签: 

it

分类: 无线通信

恩,回答了我很多迷惑不解的问题~

 

 【 在 dypx (第一片雪) 的大作中提到: 】
: 我的目的是估计时域的信道冲击响应,然后低通滤波出去噪声
: 理论上时域的信道冲击响应长度以后的值应该很小,为噪声产生
: 我的做法是利用长训练序列估计出频域响应[-26,...,-1,1,...26]等
: ...................
我想你估计出来的信道频谱序列写的有问题了。再做IFFT之前首先要保证你的信道频谱序列是正确的。而

频谱序列的排列顺序应该是[1 2 ... 26 -26 ... -2 -1].才可以直接做IFFT.
我的愚见,欢迎拍砖.

 

 

发信站: BBS 水木清华站 (Wed Nov 29 18:17:04 2006), 站内

function mod_symbols=modulate_16QAM(interleaved_bits)

global ofdm_para;

users=ofdm_para.Users;                        

n_mod_syms=size(interleaved_bits,2)/4;         

mod_syms=[j*1/sqrt(5) -1/sqrt(5)+j*2/sqrt(5) 1/sqrt(5)+j*2/sqrt(5) j*3/sqrt(
5) ...

            1/sqrt(5)  2/sqrt(5)-j*1/sqrt(5) 2/sqrt(5)+j*1/sqrt(5)   3/sqrt(
5) ...

           -1/sqrt(5) -2/sqrt(5)+j*1/sqrt(5) -2/sqrt(5)-j*1/sqrt(5) -3/sqrt(
5) ...

         -j*1/sqrt(5)  1/sqrt(5)-j*2/sqrt(5) -1/sqrt(5)-j*2/sqrt(5) -j*3/sqr
t(5)];   

 

mod_symbols=zeros(users,n_mod_syms);

for k=1:users

    cor_bits=reshape(interleaved_bits(k,:),4,n_mod_syms);     

    mod_symbols(k,:)=mod_syms([8 4 2 1]*cor_bits+1);    这一句是什么意思

                                                      

end                                       

这是一个16QAM调制的程序 它怎么实现如下调制对应关系的:

0000 —— j*1/sqrt(5)

0001 ——  -1/sqrt(5)+j*2/sqrt(5)

0010 ——   1/sqrt(5)+j*2/sqrt(5)

0011 —— j*3/sqrt(5)

0100 ——   1/sqrt(5)

0101 ——   2/sqrt(5)-j*1/sqrt(5)

0110 ——   2/sqrt(5)+j*1/sqrt(5)

0111 ——   3/sqrt(5)

1000 ——  -1/sqrt(5)

1001 ——  -2/sqrt(5)+j*1/sqrt(5)

1010 ——  -2/sqrt(5)-j*1/sqrt(5)

1011 ——  -3/sqrt(5)

1100 —— -j*1/sqrt(5)

1101 ——   1/sqrt(5)-j*2/sqrt(5)

1110 ——  -1/sqrt(5)-j*2/sqrt(5)

1111 —— -j*3/sqrt(5)

发信人: wbear (熊熊), 信区: Commun      
标  题: Re: 请高手说一下软判决和硬判决的一些区别,谢谢!
发信站: BBS 水木清华站 (Mon Sep 16 13:50:05 2002), 站内信件

硬判决是根据汉明距离进行度量,软判决是根据欧式距离进行度量。

研学论坛:

有一个简单的A COST-207信道,有5径.各径的功率延迟系数为[a(1),a(2),a(3),a(4),a(5)],各径是独立的

复高斯随机过程.我用MATLAB这样编,可以吗
u=randn;
n=randn;
h=[];
for k=1:5
h1=a(k)*(u+j*n);
h=[h,h1];
end

【求助】块衰落信道建模
请各位做个块衰落信道模型的兄弟指点一下:仿真信道时y=ax+n中,对块衰落系数在一个block(或者帧

上)是恒定的,但各个块之间是怎么变化的我就不清楚了,我用方差为1的complex gaussian noise建模

,但这样对系数很小时,等价于信噪比非常低,这样这个块的误码率就非常高了。不知道是不是a系数在

块之间变化采用瑞丽建模时,方差有问题,还是模型不正确,如果哪位仿真过,麻烦给个回复,谢谢!

块衰落就是一块1个乘性因子到符号的幅度上,你这样做没有什么不对.

块衰就是相当与每个块产生一个随即的a,归一化就是方差为1的complex gaussian noise.

你说的情况是通信中的深度衰落的问题,当然误码率就非常高了.

块衰的仿真方法和快衰的仿真方法是不一样的.对发送符号来说就不是各态历经的,需要的是仿真outage

概率.要仿真大量的a才行.

【求助】Jakes仿真器可以仿真高斯型的多普勒功率谱吗?
一般来说, 信道仿真器的实现可以用三种方法, 一个是正旋波叠加的方法, 即SOS方法, 第二中常用的方

法是无线通信原理与应用一书中提到的基于FFT变换的方法, 就是根据Dopper谱的形状和表达式在频率域

实现, 这是个很好的方法,其缺点是不适合长时的连续的信道产生, 第三种方法是IIR滤波加插值的方法,

这个方法也很好, 很适合长时间的连续信道产生.
如果要仿高斯型多普勒功率谱, 我觉得第二,三种方法都可以, 第一种我不确定.
其中, 用第二种方法很简单, 如要用第三种方法, 如何设计IIR滤波器可能有点麻烦.

【求助】如何由Doppler shift和Doppler Spread求得maximum doppler frequence
在做DRM信道仿真,根据标准提供的信道,已知的参数是delay, path gain, Doppler Shift和Doppler

Spread。
参考"Mobile Radio Channels Modeling in Matlab"by Nikolay KOSTOV里提供的程序来实现每个可分辩

径的Rayleigh包络。
function r = ray_doppler(fm, M, dt, N)
T = N*dt-dt; t = 0:dt:T;
c = sqrt(2/M); w = 2*pi*fm; x = 0; y = 0;
for n = 1:M
alpha = (2*pi*n-pi+(2*pi*rand-pi))/(4*M);
ph1 = 2*pi*rand - pi;
ph2 = 2*pi*rand - pi;
x = x + c*cos(w*t*cos(alpha) + ph1);
y = y + c*cos(w*t*sin(alpha) + ph2);
end
r = sqrt(x.^2 + y.^2)/sqrt(2);
标准里Doppler Shift这个参数应该是PSD的中心频率,Doppler Spread这个参数是PDS的宽度...程序的参

数fm是maximum doppler frequence

那具体这个Doppler Spread的含义是什么,标准只是说the width of the PDS is quantified by a

number and is refered to as the Doppler spread

那如果要利用这个程序,也就是要从Doppler shift和Doppler spread来求出fm,因为概念不是很理解,

所以很疑惑该如何求得fm,特别是当Doppler shift不为0的时候。

多普勒shift可以理解为群偏移,而多普勒Spread可以理解为偏移量的摆动范围。因此最大多普勒频移为

Dsh+Dsp,最小为Dsh-Dsp。
二者的具体涵义我以前发帖说过的,有兴趣可以参考国外讲信道的书籍。不要看国内的

恩,谢谢xiexiang57 的提醒。
Jakes PSD Dsp=fmax/sqrt(2)
Gaussian PSD Dsp=fc/sqrt(2*ln2)
对于fc=sqrt(ln2)*fmax的情况,两者是一致的, Dsp=fmax/sqrt(2)

这个程序仿真多径(当然M=1就是单径咯)
fm是多普勒频移
N是采样的点数
dt采样间隔
M是多径的条数

看不懂这个瑞利程序:(
瑞利信道是一个可辨径(多个不可辨径的合成),能把1:M的加和看作是多径的表达吗?

这个仿真程序是出自论文"Mobile Radio Channels Modeling in Matlab"by Nikolay KOSTOV
程序中的M指的是正弦波的数目,应该不小于8,M越大约精确
而不是指径数

对的 我理解的有问题,:)
这个程序的输出是瑞丽衰落因子,要仿真多径 还需要时延等参数
当瑞丽衰落因子的仿真方法有两种:
(1)正弦波叠加的方法 (上述程序所示)
(2)Jakes多普勒谱滤波的方法
对于第一种方法 详细的可以参照 patzold的mobile radio fading一书 以及一些相关英文文献
第二种方法 T.S.Rappapot的无线通信原理于应用中有提到,书中也有提到整个多径仿真的框图


这两个公式从何而来,请指示。我怎么看见如此漠生,能不能提供一下参考文献。
公式在Mobile Fading Channels page 39...


是用在描述DRM信道属性时的文字。。

1.描述为什么选择WSSUS信道模型时说的:
The justification for the stationary approach with different parameter sets is, that results

on real channels lead to BER curves between best and worst cases found in the simulation.
这句话是什么意思呢?

2.描述经过信道的信号的PDS和方差的:
The variance is a measure for the average signal power which is received via this path and

is defined by the relative attenuation ρk and the PDS determines the average speed of

variation in time. The width of the PDS is quantified by a number and is referred to as the

Doppler spread Dsp of that path - listed in table B.1.

the PDS determines the average speed of variation in time. 这一句又是怎么理解。

第二句话比较好办,其中 “the PDS determines the average speed of variation in time” 意即

"PSD 决定 (了) 随时间变化的平均速度"。

第一句话完全按字面翻译的话比较别扭。其大意是说,配以不同参数的平稳模型的合理性 (或正当性) 在

于由从真实信道中得到的结果所画出来的 BER 曲线(正好) 位于从 (采用这种平稳模型所进行的) 模拟中

得到的最好的和最坏的结果的中间。意思就是说,从真实信道中得到的 BER 曲线与从这种模型中得到的

结果很吻合,所以该模型是合理的。

xixi,先谢谢未昏...

可还是一直很糊涂,为什么PSD决定了 随时间变化的平均速度。
还有一个描述是,Doppler Spread是谱展宽的一个测量值,它是移动无限信道的时间变化率的一种度量。

不理解这个时间变化率是个什么样的概念...=(

一个随时间变化的随机信号,其变化的快慢可以用它的自相关函数反映出来。比如说,如果自相关函数是

一个冲激函数 (delta),那么它的变化是无穷快的;反之如果自相关函数是一条平坦的直线,那么它基本

是不变的。而 PSD 是自相关函数的付氏变换,当然也反映 (或决定) 着信号随时间变化的快慢。PSD 的

非零部分的宽度 (也就是 Doppler Spread) 越宽,信号变化就越快。 所以 Doppler Spread 的确 "是移

动无线信道的时间变化率的一种度量。"

更简单地说,一个信号随时间变化得越快,它的高频成分就越多,变化得越慢,它的高频成分就越少。而

PSD 就反映了高频成分的多少,所以 PSD 反映 (或决定) 着信号随时间变化的快慢。

 

【求助】为什么射频信号发送时,只需要发送实部信号?
为什么射频信号发送时,只需要发送实部信号,就能包含原始复数基带信号的所有信息了?
想不明白...
请大家释疑...
谢了~~~
基带复信号xb=xi+j*xq
调制后射频信号xrf=real(xb*exp(j*w*t))=xi*cos(w*t)-xq*sin(w*t)
解调后y=xrf*exp(-j*w*t)=(xi+j*xq)/2+xi*cos(2w)/2+....
滤掉2w分量就得到xb

【求助】信道的设计
我要对一个信道进行仿真,这个多径信道路径数目L=16,他的系数是独立的复高斯随机变量,方差是这

样定义的(英文原文):variance defined by an
exponential power profile,i.e. sigma(这里的上标是平方2,下标是k,不好表达就这么形容了)=exp(-

k/10),for k=0,1,...,L-1

我的疑问就是power profile,汉语应该是功率延迟分布吧,到底他的物理意义怎么理解啊,实在很很困

惑的,
然后不明白的就是,这个信道到底要怎么根据给的方差设计呢?是不是共有16条路径,假设要测试3次,

对应第一径,产生一组随机数字,就是3×1的矩阵他们的方差是k=0时对应的方差
sigma=1,然后对于第二径,同样,产生一组随机数字,他们的方差是k=0时对应的方差sigma=exp(-1/10),
以此类推,产生完16个这样的合并为一个3×16的矩阵,对应其中每次一仿真,信道系数取16个,就是这

个大矩阵的每一行?
还有个疑问就是如果这一点话是不是意味每次信道都定了接收到16径
不知道我的意思表达清楚没有,希望大家帮帮我,困扰好久了,我不是爱偷懒的好孩子,希望大家的帮助

下让我走的更远

"对应其中每次一仿真,信道系数取16个,就是这个大矩阵的每一行?" 恩,对的, 每一列的统计特性应

该服从Rayleigh或Rician分布特性,每一行代表信道取样(channel snapshot)

对应其中每次一仿真,信道系数取16个,就是这个大矩阵的每一行,16个独立同分布的随机变量!

个人看法,首先power profile 就是不同时延路径的功率分布。你要做的就是如何让各个路径的复高斯过

程独立。还有对与各个高斯过程中的参数计算方法都会对psd有依赖,这个psd是怎么定义的。还有就是你

如何定义这16条路径,时延间隔时多少。这些都很模糊,建议你读读 matthias patzold的mobile fading

channel.这样也许对你会有些帮助

【求助】请问多谱勒频移,多谱勒扩展,频率色散和频率弥散之间的区别?? [精华]

http://bbs.matwav.com/post/view?bid=28&id=577193&sty=1&tpg=4&age=0

【求助】信道多谱勒扩展和多谱勒频移的区别

这个问题困扰了我很久,最近终于有些理解了(参 Mobile Fading Channels - Modelling Analysis &

simulation ,第三章)

首先反对上述“最大多普勒频移=多普勒扩展”,从概念上讲不是这样的。

多普勒频移是针对一条路径的信号的,多普勒扩展是多径叠加产生的综合效应。文章中有讲,多普勒扩展

是平均多普勒频移(The average Doppler
shift (Doppler spread) describes the average frequency shift (frequency spread) that a

carrier signal experiences during transmission.原话),扩展在不同的谱下定义不同。比如高斯谱

下,定义为谱线的-3dB宽度。

申请技术分!!!

对多普勒频移和扩展的产生描述原文如下:
The transmitted electromagnetic
waves mostly do not, at least in urban areas, arrive at the vehicle antenna
of the receiver over the direct path. On the other hand, due to reflections from
buildings, from the ground, and from other obstacles with vast surfaces, as well as
scatters from trees and other scatter-objects, a multitude of partial waves arrive
at the receiver antenna from different directions. This effect is known as multipath
propagation. Due to multipath propagation, the received partial waves increase or
weaken each other, depending on the phase relations of the waves. Consequently, the
received electromagnetic field strength and, thus, also the received signal are both
strongly fluctuating functions of the receiver’s position [Lor85] or, in case of a moving
receiver, strongly fluctuating functions of time. Besides, as a result of the Doppler
effect, the motion of the receiver leads to a frequency shift (Doppler shift)1 of the
partial waves hitting the antenna. Depending on the direction of arrival of these partial
waves, different Doppler shifts occur, so that for the sum of all scattered (and reflected)
components, we finally obtain a continuous spectrum of Doppler frequencies, which is
called the Doppler power spectral density.

 

【求助】怎么多普勒滤波器长度?--做过WSSUS信道滤波法仿真的来帮忙!
在WSSUS模型中,单径的多普勒频移是由高斯白噪声经过滤波器生成的。那么,这个滤波器的时域抽头数

怎么选择?其频率响应是给定了的高斯形状。

滤波以后生成的衰落参数,怎么和实际的输入信号采样速率进行匹配呢?直接使用插值么?小弟做了几天

了,已经晕了...

 
呵呵,不知楼主所用到的多谱勒功率谱服从什么分布?是高斯型还是jakes型,若是高斯型,我做过,建

议去网站:http://www.qsl.net/ae4jy/pathsim.htm
上面看看,这篇文献很详细地讲解了构造高斯型锐利衰落信道的过程。
如果是后者,MATLAB上面不是有个函数rayleighchan直接就可以拿来用么,如果你自己想要编写程序的话

可以把源程序调出来看啊,不过我认为最简单的方法就是直接用MATLAB中自带的函数,还有我想问的是,

你为什么要用这种方法来仿真WASSUS?毕竟滤波器的设计是有一定难度的,建议采样其他方法。


任务在手,避免不了啊。我做的就是高斯型的信道,理论资料一堆i,实现资料很少,最后卡在滤波器与

插值的问题上了。作出来的衰落参数相关性摆动很大,不能拟合贝塞尔函数。滤波器长度太大,并且上插

值倍数造成插值后图线不平滑,问题很多。

做滤波器确实不容易,我查了些中文讲滤波器的文章,都是用matlab直接设计,对这种特殊窄带响应的没

有多少描述。你有这方面的资料么?麻烦发送一点,非常感谢,我对滤波器设计了解太少了。

xiexiang57@uestc.edu.cn


任务在手,避免不了啊。我做的就是高斯型的信道,理论资料一堆i,实现资料很少,最后卡在滤波器与

插值的问题上了。作出来的衰落参数相关性摆动很大,不能拟合贝塞尔函数。滤波器长度太大,并且上插

值倍数造成插值后图线不平滑,问题很多。

【讨论】频率选择性与多径是否可分离?
在CDMA系统中,经常使用不可分离多径,可分离多径。有的文献说,频率选择性衰落等效为一个FIR滤波器,

具有可分离性。但一般对频率选择性信道中不使用RAKE接收,是因为具有较多不可分离多径,如何理解可分

离与不可分离多径呢?请大侠们讨论!!!!!!!


频率选择性信道等价于一模拟带通滤波器, 变换到数字滤波器,其冲激响应为无穷多个脉冲,难以用

Rake分离主径与旁径.一般采用均衡器解决,均衡器所均衡的是信道的频率特性,而不是相位特性.

多径信道的延迟特性明显,可用Rake或采用均衡器解决,均衡器所均衡的是信道的相位特性,而不是频率

特性.

Rake使用的前提是信道各径的延迟特性明显,个数有限.

可分离多径是指到达接收端的两条路径的信号的时延相差超过一个chip时间(cdma系统),因为信道的相

干带宽近似等于最大时延差的倒数,所以信道的相干带宽小于信号的带宽,所以带宽系统会引起频率选择

形衰落。

Rake是时间分集,它适用于时延超过1 chip的多径信号,对时延小于1个 chip的一般采用时域或者频域均

衡。


有的资料将时间分集和多径分集分别定义,而RAKE归类于多径分集。关于频率分集的说法也是从一些资料

看到的,一时找不到是那份资料了,大概是这样说的:当带宽为W的信号可以分辨多径分量,将若干独立

的衰落信号路径提供给接收机,其时间分辨率为1/W;于是,宽带信号的使用刚好可看作阶数为L=W/fc的

频率分集的另一种方法(这里用fc表示相干带宽)。


LAS-CDMA基本原理是在选取的码字集(比如WALSH)中,对码字按照一
定规律插0,从而
得到的码字性能有一个0相关窗口。在传统的直扩或者是跳频,自相关
和互相关经证明
是不可能同时达到最优的,但是LAS-CDMA经过特殊的构造,成为一个三
值序列,其自相关
和互相关可以同时达到最优+1或者-1,不过同时带来的问题是,经过
这样构造的码字长度
太长,16长的WALSH经过这样的构造后形成的LAS-CDMA码字有4000多,
可想而知,其在频谱?的带宽的改变,不过如果将用户平均传输速率考
虑在内,其频谱利用率依然超过目前所有
的CDMA协议。 我觉得其中精妙之处在于0相关窗口,在相关主瓣与旁瓣
之间存在一段0值,
这个窗口宽度超过了最大弥散时间,因此可以达到其没有ISI,最重要的
是没有MAI的优势。
 太长,16长的WALSH经过这样的构造后形成的LAS-CDMA码字有4000多
,可想而知,其在频破咨?: 的带宽的改变,不过如果将用户平均传输
速率考虑在内,其频谱利用率依然超过目前所有
: 的CDMA协议。 我觉得其中精妙之处在于0相关窗口,在相关主瓣与旁
瓣之间存在一段0值?问题在于主瓣是有一定的宽度的呀!也就是说,
多径带来的ISI怎么可能没有呢?除非这个
窗度很小,这是说自相关,如果你需要区分多径就要它们分得开在时间
上。对于互相关,
如果能够做到这个0区域在零点附近,那么如果各个用户之间的相对时
延不是很大的话,
那的确没有MAI。那这样的话就需要对用户到达基站的时间进行严格的
控制,但是每个用户
离基站各不相同,路径时延各不相同,产生的多径各不相同,要像做到
无MAI谈何容易。
: 这个窗口宽度超过了最大弥散时间,因此可以达到其没有ISI,最重要
的是没有MAI的优势?还有就是这个码跟CDMA怎么联系?它这么长应该
是用来做扩频码用得吧!是用它来替代
IS-95种m序列与WALSH码结合的扩频码。

 

【求助】jakes模型的适用范围

大家好!

我最近有个问题比较伤神。

就是jakes、watterson信道模型的适用条件各是什么?

它们各在什么频率段使用?

我看国外的很多信道模拟产品都是用的jakes模型,Agilent的信号源里也用jakes模型。

但是也发现有一些用的watterson模型,

那位大侠告知一下二者的区别?

如果有标准之类的文档,麻烦给我发一份。

laoma-mail123@163.com

 

我了解的好像jakes模型适用于一般的无线移动通信(传输距离不太远)

watterson信道适用于短波窄带信号(跟一般无线信道的区别貌似在于短波经过电离层放射传输,传输距离远)

 

watterson模型一般用于30MHZ以下的中短波信道,而jakes模型用于一般无线移动信道的仿真。

 

【求助】WSSUS信道单径增益的时间相关性怎么仿真?

 

我在进行WSSUS信道的仿真,使用时域抽头滤波模型,每个抽头的系数对应一径的小尺度影响。

 

s输入

-------z----------z-------z---------z---------z--------

|c0.....|c1......|c2...............|cn-1.......|cn <-每一个参数对应一个抽头

-----------------------求和-------------------------

得到 y

 

每一径的衰落ci(t)是使用 滤波法生成的,用具有多普勒功率谱的滤波器对白噪声进行滤波,产生ci(t)。

 

文献提出,Ci(t)的时间相关函数 r(dt)= Et(ci(t)xci*(t+dt))趋近于一类0阶贝塞尔函数,但是我仿真出来拟合程度很差。

滤波器功率谱使用标准多普勒频移谱线

用两个独立高斯噪声源生成Ci(t)的实部和虚部

 

那位做过相关研究?望不吝赐教。

 

主要参考文献

MODELLING AND SIMULATION OF A FADING CHANNEL

可以直接从google找到。

 

希望各位在做相关研究的研友交流一下信道仿真的经验。我的邮箱xiexiang_57@163.com,欢迎交流

 

多做几次平均呢

 

我以前用正弦波叠加的方法做过类似的 效果一般

滤波方法没试过

 

楼上那位,作平均也没有很大的作用,到底是怎么回事呢?从理论上来说,为什么会出现贝塞尔函数呢?

 

M. Patzold, “Mobile Fading Channels”, John-Wiley, New York, NY, 1999

这本书上有讲到

多普勒谱满足Jakes谱形状 则其时域相关为贝塞尔函数

 

楼主,因为信道衰落参数的相关性是重要的信道衡量指标,相关性差,信道产生的US影响就产生了畸变。

请问你用正弦波叠加是否是用的Jacks模型的生成方法?相关性不好的原因在哪里呢?是否有分析过?谢谢。

 

 

【求助】DVB-T的能量扩散里的一个小问题

在能量扩散里PRBS的周期为什么是1053字节?我理解的是,在后面的7个包里每个包的同步字节,是没有PRBS的输出,可周期怎么会是1053。望知道的大虾回答下,谢谢

 

首先是1503

randomize的时候是不会对第一个同步头(B8)做加绕和解绕的,所以1503=188×8-1

 

后面的同步头0x47经PRBS加扰后在须强行还原成0x47.

 

第一个188字节的第一个byte胃sync,同步字,直接取反为B8hex

 

【讨论】信道估计的mmse算法

 

这个困扰很久了

就是这个mmse算法里面涉及到信道的自相关矩阵,这要是已知信道才能求得

但是要是知道的话那还估计什么信道阿?

我想那可能是针对导频位置的信道说的,这样可以解释通。但是这个又涉及到一个问题就是公式里面要求信道自相关函数,那好比有个信道h=[1 2 3]'

按照公式应该信道自相关函数是个3*3的矩阵

用命令xcorr是不能得到的

用文献里看得E{h*h},可是这样的话对他求逆必然行列式=0,没法求!

看文献里都大部分引用mmse算法,似乎很基础很简单,从没有人给过详细方法步骤

 

请大家讨论

 

信道估计是为了恢复信号阿

最简单的来说Y=X*H+noise

得到接收符号Y知道了信道H的话就可以恢复出信号X阿

 

我个人的理解是,mmse算法假设知道信道冲击响应(注意这只是假设,和真实的信道响应有一定的差距,但也有一定的相似).

可以参考IEEE文献:Jan-Jaap van de Beek On channel estimation in OFDM systems

 

但是请教楼上的同学(暂且这样称呼)

那信道自相关矩阵该怎么求呢

今天下午和导师争论了一下午,他说回家好好看看我研究的论文再说

他给我一个第四版自适应滤波器里面的公式,可是我觉得还是没法求这个东西

请多多指点

 

我觉得首先要区分随机量(random)和确定量(deterministic)。如果信道是确定量(如有线信道),mmse算法等同于ml算法。如果是随机量,mmse算法使信道估计有更小的估计方差。

考虑基於帧(packet/frame)传送的无线通讯系统,Rayleigh信道信道系数可以看作是复高斯随机变量,但是这个随机特性对於每个帧而言是表现不出来的(假设信道缓慢时变)。也就是说,相对每个帧而言,信道系数可以看作是随机变量的样本(确定量)。

但把许多帧作为整体来看,随机特性又显现出来。例如在GSM系统中,对於rural area, typical urban 等典型环境,GSM标准中给出的power-delay profile,其实就给出各信道系数的自相关函数。这些物理信道系数是统计独立的。当然由於接收器滤波和采样,等价信道系数(包括成形滤波,物理信道,接收滤波和采样)一般不再统计独立,所以有了信道自相关矩阵的概念。

推荐: "A Comparison of Pilot-Aided Channel Estimation

Methods for OFDM Systems", Morelli and Mengali, IEEE Trans. Signal Processing, pp. 3065-, Dec. 2001

 

互相关是接收信号和本地导频信号的相关运算,自相关当然是本地导频信号的自相关,至于相关运算你可以自己写个代码,不用xcorr,这样心里有数

 

MMSE在信道估计时有一个重要的前提是良好的帧同步,否则会引入pilot的旋转,导致MMSE效果下降.这也是为什么MMSE文献上到处都是而实际中应用很少的主要原因.

在实际中根本无法获得准确的帧同步,(FFT窗口落在GI处),此时简单的低通滤波即可达到很好的信道估计.

 

当信道的ISI不是很严重时,低通滤波即可较好实现均衡。至于帧同步,虽然很多时候是困难了一点,但是并非不可实现。为了获得良好的性能,不论是盲信道估计还是盲均衡都是很有价值的。

实际应用的算法很多时候都是ad hoc的,能满足应用即可。所以做工程长了水平很难提高啊。

 

自相关矩阵的求法其实很简单,你可以将信道响应作一个自相关运算,也就是经xcorr 或者autocorr都可以,基长度为2*N-1(假设信道响应的长度为N的序列)那么可以发现其运算后的序列是在第N点对称的,然后取其后边的N位就是其自相关函数了(这里用xcorr函数),在构造自相矩阵的时候,其实我们可以发现,那是一个具有一定对称性质的toeplitz矩阵,因此你可以用matlab自带的toeplitz函数将自相关函数进行处理,而后得到的矩阵就是一个信道响应的自相关矩阵

 

经典的做法是利用自相关函数和功率谱的傅立叶变换关系,首先去确定功率谱,而根据jakes <Microwave mobile communications>,功率谱具有经典jake型。有时候近似为矩形功率谱,那么自相关函数(自相关矩阵)为sinc型。

 

要用时间平均代替集平均计算自相关矩阵,误差很大,尤其矩阵的边界点。

 

【求助】信道估计的创新

 

我毕业论文是OFDM信道估计,要求一定有创新思想

以一篇2005年IEEE文章为基础,加入自己的创新,这篇文章的估计思路我很喜欢而且也仿真出来了,主体思路为

1. 粗估计,用MLE

2.精确估计,设计个代价函数对信道求导,迭代,精度很高

 

我的改进思想其实很简单拉,就是吧第一步粗估计的MLE用MMSEE代替,因为正好找到一篇2001年论文对比这辆种算法性能,而且仿真结果不错,改进后提高1dB,适合低信噪比,星座图明显比MMSEE和2005年这篇论文效果要好,误码率也好,

但是我的导师把我一顿挫,说这没什么意义,有人做过你还作什么啊,而且显然MMSE比MLE好嘛,想想也是,感觉信道估计MMSEE就像老祖宗一样,挺难确很经典,如果用MMSEE感觉象倒退了,

我考虑答辩可能老师嫌我的改进复杂,我想用奇异值分解MMSEE,这样算法精度应该高也简单些了,可心里还是忐忑,不知道自己研究的是不是要被全盘否定,来着里问问高人和热心肠的大家,走过路过麻烦您给我点建议,或者MMSEE有什么更好改进的方法?等等了

谢谢大家啊谢谢拉

 

MMSE 理论上不算新东西,可是实现上还是很难的.

于是有很多sub-optimal的方法来简化MMSE,或者用adaptive的方法逐步逼近MMSE.

 

我觉得你有时间的话,可以试着找到一种性能比MMSE更好,实现上也比较可能的方案,这样看起来,很有创新性了.

 

在我看来MMSE的性能还不够好,可是考虑到实现复杂度,这个已经很高了.

 

【讨论】复高斯随机变量的pdf

令Z=X+jY,其中X、Y都是实高斯RV,并假设其为标准正态。显然Z为复高斯RV,则Z的概率密度函数pdf如何来确定,f(Z)=?

请高人指教。

 

f(z)=1/(pi*sigma^2)*exp(-|z-u|^2/sigma^2)

其中u=E(z)=Ux+j*Uy

 

关于复高斯随即变量的pdf, 好像没有这么简单。f(x,y)准确表示他的分布。至于苦参网友的结论,似乎只是x,y服从i.i.d.情况下的特殊情况(circularly symmetric Gaussian complex R.V.)。

 

推荐楼主参考David Tse的“Fundamentals of Wireless Communication”的附录A, 里面有详细的讨论。

  • 0
    点赞
  • 0
    评论
  • 1
    收藏
  • 一键三连
    一键三连
  • 扫一扫,分享海报

©️2021 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值