matlab09-线性方程式、统计

一：线性方程组的解

1. 高斯消去法
   增广矩阵通过行变换化为最简阶梯矩阵，得出等价方程组
   Ax=b
   R=rref([A b])

>> A=[1 2 1;2 6 1;1 1 4];
>> b=[2;7;3];
>> R=rref([A b])

R =

     1     0     0    -3
     0     1     0     2
     0     0     1     1

>> A\b   %左除

ans =

   -3.0000
    2.0000
    1.0000

2. Ax=b x=inv(A)b (det(A)~=0)

>> inv(A)*b

ans =

   -3.0000
    2.0000
    1.0000

x=0:0.5:40;
y=0:0.5:40;
[X Y]=meshgrid(x,y);
z1=-X-Y;
z2=Y-X;
z3=X/3;
colormap(hsv);
hold on
surf(X,Y,z1)
surf(X,Y,z2)
surf(X,Y,z3)
axis square;
 

二：统计

function	功能描述
mean()	平均数
median()	中位数
mode()	众数
prctile()	百分位点，median()相当于50%
max	最大值
min	最小值
var()	方差
std()	标准差

练习：

>> load stockreturns;
>> x4=stocks(:,4);
>> length(x4)

ans =

   100

>> mean(x4)

ans =

  -5.8728e-04

>> max(x4)

ans =

    6.4692

>> min(x4)

ans =

   -5.8764

>> mode(x4)

ans =

   -5.8764

>> median(x4)

ans =

    0.0617

>> prctile(x4,50)

ans =

    0.0617

>> prctile(x4,70)

ans =

    1.2359

>> prctile(x4,90)

ans =

    2.9449

>> var(x4)

ans =

    5.5649

>> std(x4)

ans =

    2.3590

2. 频数分布直方图：已知频数fq， bar(fq)

>> freq=[1 0 1 0 4 0 1 0 3 1 0 0 1 1];
>> length(freq)

ans =

    14

>> subplot(1,3,1);
>> bar(freq);
>> xlim([0 15]);
>> subplot(1,3,2);
>> area(freq);
>> xlim([0 15]);
>> subplot(1,3,3);
>> stem(freq);
>> xlim([0 15]);

练习：告诉原数据，画出频数分布图

求出频数分布：自定义函数freq()

function fq=freq(data)
ii=max(data)
i1=length(data)
fq=zeros(1,ii)
for jj=1:i1
    fq(data(jj))=fq(data(jj))+1
end
return  

执行代码：

>> x=[1 3 5 5 5 5 7 9 9 9 10 13 14];
>> fq=freq(x);
>> fq

fq =

     1     0     1     0     4     0     1     0     3     1     0     0     1     1

>> subplot(1,3,1)
>> bar(fq);
>> xlim([0 15]);
>> subplot(1,3,2)
>> area(fq);
>> xlim([0 15]);
>> subplot(1,3,3)

3. boxplot()

练习

>> load stockreturns;
>> boxplot(stocks);