作者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