matlab 有內建函式
mvnpdf
在command window 打 help mvnpdf 會得到他的用法
※ 引述《HOWARDXXXXXX (心誠則靈,我想)》之銘言：
: 不要用fix取整數,用無條件ceil直接進位
: 應該可以成功解決邊界問題
: 可以想到用迴圈寫真的太強大了 ㄎㄎ
: clear all
: l=input('Input lumda ');
: n=input('Input upperbound of y ');
: d=input('Input desity ');
: [x,y]=meshgrid([0:d:n]);
: for i=1:ceil(1/d*(n));
: for j=1:ceil(1/d*(n));
: if i>=j;
: y(i,j)=x(i,j);
: else
: y(i,j)=0;
: end
: end
: end
: z=l*l*exp(-l*y);
: meshz(x,y,z)
: xlabel('X')
: ylabel('Y')
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.7.59
※ 編輯: youyouyou 來自: 140.112.7.59 (12/04 20:28)