一:线性方程组的解
- 高斯消去法
增广矩阵通过行变换化为最简阶梯矩阵,得出等价方程组
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
- 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
- 频数分布直方图:已知频数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);