[繪圖] 畫圖可能遇到的問題?

看板MATLAB作者 (rffvb)時間7年前 (2017/02/19 23:23), 編輯推噓1(101)
留言2則, 1人參與, 最新討論串1/1
各位大大你們好 遇到個問題不知道有沒有大大 有沒有遇過類似經驗 -------------------------------------- 算 wishart matrices 的特徵值分佈 http://i.imgur.com/Tdn1Qmi.jpg
http://i.imgur.com/rrWbJx4.jpg
以下的圖都是用上面的 equation 所畫出來的 Code就蠻多行的,假如有大大好心能幫我看看可以站內信,或許是哪裡出錯了 … 我先假設 code (Matlab)是對的 這是我在PAPER上看到的圖 http://i.imgur.com/DPsw9rt.jpg
這是我畫出來的圖 http://i.imgur.com/70b5YZ2.jpg
看起來很類似 總共會有三個特徵值的分布 和一條 平均線 但現在我把q和p的大小改了 q=6,p=12 圖還是平順的 http://i.imgur.com/LyxenU5.jpg
q=6,p=18 綠色的線已經有點壞掉了 http://i.imgur.com/CZgVCta.jpg
q=6,p=48 3,4,5張圖已經壞的蠻大的 http://i.imgur.com/VTvXDUP.jpg
照理說公式算出來的圖應該是要平順的 但P的值越大 壞的越多 因為公式裡蠻多複雜的計算 還有階乘的運算( EX :48 ! ,就不知道幾位數了) 是不是要想辦法避開這麼多位數的值? 假如因為很多位數乘上很多位數 不知道在matlab 上運算會有會有 overflow or underflow 的問題 不知道這是不是造成圖畫錯的原因 ? 假如有大大有類似經驗 希望可以分享 或者可以從哪些地方下手 謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 223.140.89.42 ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1487517839.A.56F.html

02/20 17:09, , 1F
48!還好,簡單乘除應該不會出問題,不然你就用gammaln先取
02/20 17:09, 1F

02/20 17:10, , 2F
log吧,我會猜問題在det,但也要看你矩陣的本質。
02/20 17:10, 2F
文章代碼(AID): #1OgRYFLl (MATLAB)
文章代碼(AID): #1OgRYFLl (MATLAB)