% 我今天在龍芯+統(tǒng)信UOS的筆記本上安裝了

%北太天元， 然后寫了一個(gè)線性回歸的腳本。

% 第一種方法：使用了北太天元的內(nèi)置函數(shù) %polyfit

% 第二種方法：老老實(shí)實(shí)的使用了高中課本

%給 的公式

% 線性回歸分析

% 抽取了11個(gè)城市居民家庭的關(guān)于收入和食品支出的樣本，

% 通過線性回歸得到一個(gè)一元一次函數(shù) f(x) = b*x+a

% 可以用來預(yù)測(cè)某個(gè)收入對(duì)應(yīng)的食品支出

% 通過將點(diǎn)X0帶入到 f(x) 中得到 f(X0)

%

clear;

收入 = [82 93 105 130 144 150 160 180 270 350 600]; % 可以更改

支出 = [75 85 92 105 120 130 145 156 200 240 300]; % 可以更改

p = polyfit(收入,支出,1);

plot(收入,支出,'o')

x0 = [min(收入):max(收入)];

模擬支出 = p(1).*x0 + p(2);

hold on

plot(x0,模擬支出)

hold off

xlabel("收入")

ylabel("食品支出")

title("收入與食品支出的回歸分析")

["函數(shù)方程為：f(x) = ",p(1)," *x ",p(2)]

%按照高中課本上給出的公式計(jì)算 haty = a + hatb * x

x= 收入;

y =支出;

barx = mean(x);

bary = mean(y);

xx = x - barx;

yy = y - bary;

hatb = sum( xx .* yy) / sum( xx .* xx )

hata = bary - hatb * barx

"你會(huì)觀察到 hatb == p(1) , hata == p(2)"

p