協(xié)方差矩陣 在北太天元中可以使用內(nèi)置函數(shù) cov 計(jì)算， 相關(guān)系數(shù)矩陣再北太天元中暫時(shí)缺少，計(jì)劃補(bǔ)充 corr 函數(shù)來計(jì)算， 不過，利用相關(guān)系數(shù)和協(xié)方差的關(guān)系，下面的北太天元代碼用了兩行代碼也可以替代 corr 函數(shù) 的功能。

%北太天元代碼 展示 如何計(jì)算協(xié)方差矩陣和相關(guān)系數(shù)矩陣

% 創(chuàng)建一個(gè)包含身高和體重的數(shù)據(jù)集

X = [160, 170, 180, 190, 200]; ?%身高

Y = [ 50, 60, 70, 80, 90]; ? ? ? ? ?%體重

data = [ ?X' ? ?, ? ?Y' ];???

% 計(jì)算協(xié)方差矩陣

cov_matrix = cov(data);

% 輸出協(xié)方差矩陣

disp(cov_matrix);

% 標(biāo)準(zhǔn)差

std_matrix = std(data);

% 計(jì)算相關(guān)系數(shù)矩陣

corr_matrix = cov_matrix ./ ?( std_matrix' * std_matrix);

% 上面這一行的兩行代碼可以代替北太天元暫時(shí)沒有實(shí)現(xiàn)的

% ?的函數(shù) corr ?( MATLAB的函數(shù)corr 是計(jì)算相關(guān)系數(shù)矩陣的)

% 輸出相關(guān)系數(shù)矩陣

disp(corr_matrix);

plot(data(:,1), data(:,2), 'r-*', 'LineWidth' , 3);

xlabel('身高X'); ylabel('體重Y');

title("Y=X-110")

p.s.

昨天寫專欄的時(shí)候考慮不周，今天補(bǔ)上這些內(nèi)容，但偶感風(fēng)寒導(dǎo)致補(bǔ)得很潦草，后面再改進(jìn)。