作者cmsuteki (WWJD)
看板comm_and_RF
标题Re: [问题] 如何用MATLAB模拟都普勒?
时间Tue Jul 14 21:48:32 2009
※ 引述《qeagle (请以分手为前提和我交往)》之铭言：
: : 我想模拟一个信号经过都普勒的变化,不知道大家如何用MATLAB或SIMULINK来模拟?
: : 直接乘一个估计的SINE波吗?或是乘一个EXP(JFT)?
: : 请高手赐教
: : --
: : ◆ From: 114.43.0.160
: : 推 lovewa:可以使用Jacke's Model... 140.113.144.28 07/08 01:03
: : 推 onlywig:可以用Jake's Model... 163.18.104.154 07/08 09:35
: : → lovewa:打错，拍谢～ 140.113.144.28 07/08 10:46
: : 推 onlywig:love 你没说我还没发现 单纯只是附和XD 163.18.104.154 07/08 13:54
: : → qeagle:请问那是什麽?我在matlab中没有找到... 114.43.9.7 07/08 22:54
: : → onicalz:google! 140.113.236.140 07/09 01:14
: : 推 manboy:matlab: help : Rayleighchan 140.112.17.232 07/09 23:04
: : → manboy:' DopplerSpectrum ' 你有很多可以选~ 140.112.17.232 07/09 23:05
: : → manboy:如果猛一点 你想用3gpp或是itu通道 140.112.17.232 07/09 23:05
: : → manboy:你可以查 help : stdchan 140.112.17.232 07/09 23:06
: : → manboy:会有很多standard chan.可以给你挑 版本08a 140.112.17.232 07/09 23:06
: 谢谢,想再请教一下,jacke's model好像是以多个反射的机率模型
: 但我想先从单一个反射的doppler来做起
Jakes' model
比较像真正的无线通讯环境，而不是单一个反射的。
: 想请教一下,单一个反射信号若有doppler在里面,应该是一个频移的结果
: 有人跟我说直接乘sine或cosine就会有频移,但这样做的频移是双边的
: 理论上应该是单边才合理,我也不想再後面再加滤波器滤掉一边(很难做那麽靠近)
恩
: 又有人跟我说用I/Q一个乘SINE一个乘COSINE就会有单边
: 但我不知道怎麽做
: 也有人跟我说直接乘EXP(-J*FT),就会只有单边频移
听起来这两个是同一件事。
只是一个是equivalent baseband.
假设passband signal为m(t)，equivalent baseband是m*(t)
假设原本的讯号的equivalent baseband(也就是complex envelope)为Mi+jMq
(Mi+jMq)*exp(-jFT)={Mi(cosFT)+Mq(sinFT)}+j{-Mi(sinFT)+Mq(cosFT)}
={real part}+j{imaginary part}
=I+jQ
(受限於bbs，我尽力了)
m(t)=Re{m*(t)exp{wct}}=I[cos(wct)]-Q[sin(wct)]
↑就照这个式子做吧！
wc:载波角频率=2*pi*fc
: 这个真的只有单边频移,但它的负频率不对称了,似乎有点令人疑惑
: 且乘上一个复数信号,原来的实数信号变成包含了虚数
: 我不知道这样怎麽做後续的处理
: 请各位指点一二,谢谢
有虚数就是equivalent baseband。
你把bandpass的讯号直接整个右移，左右当然不对称，
我再补充一点关於I-Q的东西，希望对你做的东西有帮助：
一开始的passband讯号可以把它写成m(t)=Mi(t)cos(wct)-Mq(t)sin(wct)
要把Mi和Mq分别出来，只要分别乘cos(wct)和sin(wct)再积分(或LPF)就行了
解出I和Q以後，再去乘exp(jFT)，至於电路上怎麽实作呢，把式子展开你就知道了
然後乘完以後，得到新的I和Q，
再套一个常见的公式 m(t)=Re{m*(t)exp(j*wct)}，就可以得到bandpass signal了
至於这个要怎麽实作在电路，一样，把式子展开就行了。
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.236.143
※ 编辑: cmsuteki 来自: 140.113.199.122 (07/14 22:19)
※ 编辑: cmsuteki 来自: 140.113.199.122 (07/14 22:25)
1F：推 qeagle:请问你的m(t)化成I/Q两项,是把THETA角抽出 114.43.16.178 07/14 23:21
2F：→ qeagle:弄成三角函数吗?又假设我是载波上有个方波 114.43.16.178 07/14 23:22
是的。
3F：→ qeagle:信号存在,也就是只有有时域信号,没有公式 114.43.16.178 07/14 23:22
4F：→ qeagle:有办法直接乘,而不用推导来得到答案吗? 114.43.16.178 07/14 23:23
看我上面写的吧！乘载波+LPF就是了。
5F：→ onicalz:楼上有修过数位通讯吗 @@?140.113.236.140 07/14 23:24
6F：→ onicalz:你的问题可以参考一下 David Tse 的140.113.236.140 07/14 23:25
7F：→ onicalz:Fundamentals of wireless comm.140.113.236.140 07/14 23:26
8F：→ onicalz:2.2.2 Baseband equivalent model140.113.236.140 07/14 23:27
9F：→ onicalz:有推出你想要的 channel model140.113.236.140 07/14 23:27
推楼上
10F：→ ninzaburo:m(t)化成I/Q两项，是带通讯号表示式 219.70.173.60 07/15 00:19
11F：→ ninzaburo:其中m*(t)是m(t)的complex-envolope 219.70.173.60 07/15 00:20
12F：推 ninzaburo:我记得Haykin书的後半附录有讲 219.70.173.60 07/15 00:24
13F：→ cmsuteki:我再把I-Q的部份做个补充吧140.113.198.184 07/15 00:44
※ 编辑: cmsuteki 来自: 140.113.198.184 (07/15 01:16)
※ 编辑: cmsuteki 来自: 140.113.198.184 (07/15 01:20)
14F：推 onlywig:oni大是吴老师的学生吗@@? 112.104.40.54 07/15 20:12
15F：→ onicalz:不是耶 @@" 是黄老师 QQ"140.113.236.140 07/15 22:33
16F：推 onlywig:还不赖喔^^ 羡慕 163.18.104.154 07/16 09:23