MATLAB experiment＜3.1＞

最新推荐文章于 2024-04-01 10:17:12 发布

地表最温柔阿福

最新推荐文章于 2024-04-01 10:17:12 发布

阅读量667

点赞数 1

分类专栏： matlab 文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/weixin_60787500/article/details/127433958

版权

matlab 专栏收录该内容

8 篇文章 3 订阅

订阅专栏

5. (经济预测) 在某经济年度内，各经济部门的投入产出表如下表3.5（单位：亿元）

假设某经济年度工业，农业及第三产业的最后需求均为17亿元，预测该经济年度工业，农业及第三产业的产出（提示：对于一个特定的经济系统而言，直接消耗矩阵和Leontief矩阵可视作不变）。

B=[6,2,1;2.25,1,0.2;3,0.2,1.8];x=[25 5 20]'

>> C=B/diag(x) %第1列除以x(1),以此类推

C =

    0.2400    0.4000    0.0500

    0.0900    0.2000    0.0100

    0.1200    0.0400    0.0900

>> A=eye(3,3)-C

A =

    0.7600   -0.4000   -0.0500

   -0.0900    0.8000   -0.0100

   -0.1200   -0.0400    0.9100

>> D=[17 17 17]';x=A\D  %尽管D变了，但这里假设C,A不变

x =

   37.5696

   25.7862

   24.7690

__________

6. 求下列矩阵的行列式、逆、特征值和特征向量

___________________________________________________________________

1,行列式

h=det(A)

2，逆

n=inv(A)

3,特征值eigenvalue和特征向量eigenvector

[v,d]=eig(a)

//A*V = V*D

V特征向量矩阵，D特征对角阵

————————————————

>> a=[4 1 -1;3 2 -6;1 -5 3];det(a),inv(a),[v,d]=eig(a)


————————————————————————————————————————————————————————-------
 

%Exercise 6(2)

>> a=[1 1 -1;0 2 -1;-1 2 0];det(a),inv(a),[v,d]=eig(a)

 

%Exercise 6(3)

>> A=[5 7 6 5;7 10 8 7;6 8 10 9;5 7 9 10]


>> det(A),inv(A), [v,d]=eig(A)



%Exercise 6(4)(以n=5为例，n=50或500类似) 

%关键是矩阵的定义，当n很大时不能直接写出。

%方法一（三个for）

n=5;

for i=1:n, a(i,i)=5;end

for i=1:(n-1),a(i,i+1)=6;end

for i=1:(n-1),a(i+1,i)=1;end

a

%方法二（一个for）

n=5;a=zeros(n,n);

a(1,1:2)=[5 6];

for i=2:(n-1),a(i,[i-1,i,i+1])=[1 5 6];end

a(n,[n-1 n])=[1 5];

a

%方法三（不用for）

n=5;a=diag(5*ones(n,1));

b=diag(6*ones(n-1,1));

c=diag(ones(n-1,1));

a=a+[zeros(n-1,1),b;zeros(1,n)]+[zeros(1,n);c,zeros(n-1,1)]

%下列计算

>> det(a)


>> inv(a)


>> [v,d]=eig(a)

判断第6题各小题是否可以相似对角化，如果是，求出对角矩阵和对应的相似变换矩阵

inv(v)*a*v %验算

因为A*V = V*D

所以inv(v)*a*v=inv(v)*v*a

Exercise 7(1)

>> a=[4 1 -1;3 2 -6;1 -5 3];[v,d]=eig(a)

v =

    0.0185   -0.9009   -0.3066

   -0.7693   -0.1240   -0.7248

   -0.6386   -0.4158    0.6170

d =

   -3.0527         0         0

         0    3.6760         0

         0         0    8.3766

>> det(v)

ans =

   -0.9255    %v行列式正常, 特征向量线性相关，可对角化

>> inv(v)*a*v  %验算

ans =

   -3.0527    0.0000   -0.0000

    0.0000    3.6760   -0.0000

   -0.0000   -0.0000    8.3766

8. 判断第6题各小题是否为正定矩阵。

正定矩阵

对称且特征值全部大于零

(只有3满足)

9. 求下列向量组的秩和它的一个最大线性无关组，并将其余向量用该最大无关组线性表示。

>> a=[4 -3 1 3;2 -1 3 5;1 -1 -1 -1;3 -2 3 4;7 -6 -7 0]

>> rank(a)

ans =

     3

>> rank(a(1:3,:))

ans =

     2

>> rank(a([1 2 4],:)) %1,2,4行为最大无关组

ans =

     3

>> b=a([1 2 4],:)';、
c=a([3 5],:)';                    %其他向量

>> b\c                            %线性表示的系数

ans =

    0.5000    5.0000

   -0.5000    1.0000

         0   -5.0000

10. (二次型标准化) 用正交变换化下列二次型为标准形

二次型标准化

标准型：X^T*A*X

1，写出二次型的矩阵；


a =

     1    -2     2
    -2    -2     4
     2     4    -2

2，求出矩阵的特征值和特征向量；

[v,d]=eig(a)

v =

    0.3333    0.9339   -0.1293

    0.6667   -0.3304   -0.6681

   -0.6667    0.1365   -0.7327

d =

   -7.0000         0         0

         0    2.0000         0

         0         0    2.0000

3，
%正交变换为y=vx, 标准型-7*y1^2+2*y2^2+2*y3^2

(电路网) 图3.1是连接三个电压已知终端的电路网，求a, b, c点的电压。

%设经过6个电阻的电流分别为i1, ..., i6. 列方程组如下

%20-2i1=a; 5-3i2=c; a-3i3=c; a-4i4=b; c-5i5=b; b-3i6=0;

%i1=i3+i4;i5=i2+i3;i6=i4+i5;

%计算如下

>> A=[1 0 0 2 0 0 0 0 0;0 0 1 0 3 0 0 0 0;1 0 -1 0 0 -3 0 0 0; 1 -1 0 0 0 0 -4 0 0;

0 -1 1 0 0 0 0 -5 0;0 1 0 0 0 0 0 0 -3; 0 0 0 1 0 -1 -1 0 0;0 0 0 0 -1 -1 0 1 0;

0 0 0 0 0 0 -1 -1 1]; 

>>b=[20 5 0 0 0 0 0 0 0]'; A\b

ans =

   13.3453

    6.4401

    8.5420

    3.3274

   -1.1807

    1.6011

    1.7263

    0.4204

    2.1467