我大概觀察到問題出在哪一點了就是我在計算 Wzf=(H^(H)*H)^(-1)*H^(H) 的程式碼是這樣去計算的 for i=1:1:N/nTx B(:,:,i)=H(:,:,i)'; BB(:,:,i)=B(:,:,i)*H(:,:,i); BBB(:,:,i)=inv(BB(:,:,i)); BBBB(:,:,i)=BBB(:,:,i)*B(:,:,i); %2x2 xHat(:,:,i)=BBBB(:,:,i)*y(:,i); end 他在BBB的時候變的超級大 像 1.0e+07 * 4.7093 - 0.0000i -4.7097 - 0.0001i -4.7097 + 0.0001i 4.7102 - 0.0000i 而改用h的BBB則是長 4.3620 + 0.0000i 4.2843 - 1.9992i 4.2843 + 1.9992i 5.8087 + 0.0000i 簡單來說我用H配上noise的時候 他好像就算不出來(H^(H)*H)^(-1) 可是今天如果是用h配上noise的時候算出來的值就都是0.XXX之類的小很多的值 可是明明noise就長一樣 H跟h也長超像的-..-.... 救命啊 誰能救救我 Orz PS 以我這次模擬的例子來看 H(:,:,1)長 -0.4550 - 1.7656i -0.4550 - 1.7655i -0.3984 - 0.6744i -0.3982 - 0.6744i h(:,:,1)長 -0.0600 - 0.0846i 0.3701 + 0.3301i -0.5598 - 1.2732i 0.8337 + 0.7209i 明明就超像為啥Wzf算出來差這麼多...... --------------------------------觀察更新--------------------- 目前觀察到是因為用 rayleighchan函式產生出來的H陣列都太接近了 差別都在小數點後四位附近，所以在算反矩陣的時候都變超大的導致解不出正確的結果? 那該怎麼辦哩，不用rayleighchan就不能用都普勒效應了，而且改16QAM 64QAM好像都沒 特別大的變化。 持續觀察中... ※ 引述《xm3vul3h06 (蠶哥)》之銘言： : 是這樣的 小弟之前做 Zero Forcing的h通道都是這樣產生的 : h = 1/sqrt(2)*[randn(nRx,nTx,N/nTx) + sqrt(-1)*randn(nRx,nTx,N/nTx)]; : 然後雜訊部分則是這樣產生 : noise = 1/sqrt(2)*[randn(nRx,N/nTx) + sqrt(-1)*randn(nRx,N/nTx)]; : 這樣做起來很順利，可是今天想換成有都普勒效應的通道就把h改成用H : chan=rayleighchan(ts,fd); : for raychan=1:1:nRx : ray_signal(raychan,:)=filter(chan,sMod_test(raychan,:)); : end : look=ray_signal./sMod_test; : H=reshape(look,[nRx,nTx,N/nTx]); : 大概來說就是把經過filter的信號除上原始信號就可以得到rayleigh通道的H係數 : 有把noise關閉過，可以完全正確的解回來，可是+上noise就一直降不下來。 : 觀察過 H 跟 h 感覺也沒啥不一樣的地方就是個亂數的矩陣... : 我知道講得有點攏統啦 但是我實在找不到Bug了 : 有沒有高手可以提點一下可能會錯在哪邊 Orz -- 你不可不知道的使用簡體字的後果： #1InxrhHQ (Gossiping) → chinaviva:按你的邏輯　三人為众　繁體的眾字也能叫眾？ 01/04 → chinaviva:繁體哪個三個人？那個中間的也算人？ 01/04 → chinaviva:繁體的憂無心　簡體的忧有心　無心何來憂？ 01/04 → chinaviva:英文的愛沒有心　阿拉伯的愛也沒有心　比繁體字低級? 01/04 → chinaviva:愛要有心　那哀要不要有心啊？ 01/04 --

※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 111.249.43.101 ※ 編輯: xm3vul3h06 來自: 111.249.43.101 (02/14 01:43)