[閒聊] 有關高統作業及matlab

[focuward (泳者鬥惡龍)]

大家好

我今天花了一些時間把高統作業的程式碼寫了

希望大家可以一起討論一下，若要討論，請隨時發站內信給我，謝謝。


1.normdistribution 的 contour (100圈版)

clear all
mx=input('Input mean of x = ');
my=input('Input mean of y = ');
sx=input('Input standard variance of x = ');
sy=input('Input standard variance of y = ');
rho=input('Input correlation coefficient rho = ');
n=input('Input bound of x and y = ');
x=linspace(-n,n);
y=x;
[x,y]=meshgrid(x,y);
z=exp((-0.5/(1-rho^2))*(((x-mx)/sx).^2-2*rho*((x-mx)/sx).*((y-my)/sy)+((y-my)/sy).^2));

contour(x,y,z,100)//括號最後一項可控制圈數

2.p112 x marginal

clear all
l=input('Input lunda ');
n=input('Input upperbound of x ');
x=linspace(0,n);


z=l*exp(-l*x);

plot(x,z)

3.p112 y margianl

clear all
l=input('Input lunda ');
n=input('Input upperbound of y ');
y=linspace(0,n);

z=(y*l*l).*exp(-l*y);//注意前面刮弧後有一個頓號，這是為了確保兩刮弧是同秩的
                     //linespace
plot(y,z)

4.p112 joint pdf
這題我遇到一些困難，附上我的兩個版本，第一個版本是不考慮y>=x的情況

clear all
xlabel('x')
ylabel('y')
l=input('Input lunda ');
n=input('Input upperbound of   ');

x=linspace(0,n);
y=linspace(0,n);

    [x,y]=meshgrid(x,y);

    z=l*l*exp(-l*y);


meshz(x,y,z)

xlabel('X')
ylabel('Y')

第二個版本，我利用描點的方式：把x,y點放入由0到range n，以 d 為每行列間距的
矩陣裡，之所以用離散的表達是希望把每一點拿出來比大小，之後我再把 y 矩陣中的
每一個點換成 x 矩陣中的對應點減去 d，如此保證 y>=x，然後再畫出來就好了。
(其中d代表每行列間距，要訂為小於一，但如果訂得太小電腦會跑不動，我自己測試的時
候，訂為0.03可以得到近乎平滑的表面)

缺點：本程式的精確性不夠，尤其是在邊界值，因為我所用的演算法是把y 矩陣中的
每一個點換成 x 矩陣中的對應點減去 d，所以在 y=0 還有 x=n(boundary) 會有奇怪
的圖形，希望大家跟我討論，一起修正。



clear all

l=input('Input lunda ');
n=input('Input upperbound of y  ');
d=input('Input desity ');

[x,y]=meshgrid([0:d:n]);
for i=1:fix(1/d*(n));
    for j=1:fix(1/d*(n));

    y(i,j)=x(i,j)-d;

end
end

z=l*l*exp(-l*y);

meshz(x,y,z)
xlabel('X')
ylabel('Y')

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.112.7.59
推 Avle:推 泳者鬥餓龍~ 下週六 一起來龍門游吧~
推 yu0101:XD