线性回归方程式与线性系统

最新推荐文章于 2023-03-19 19:34:22 发布

关关雎鸠儿

最新推荐文章于 2023-03-19 19:34:22 发布

阅读量311

点赞数

分类专栏： MATLAB

本文链接：https://blog.csdn.net/weixin_42260102/article/details/100849636

版权

MATLAB 专栏收录该内容

13 篇文章 1 订阅

订阅专栏

Gaussian Elimination:rref()
在这里插入图片描述
LU Factorization:lu()

Matrix Left Division:\or mldivide()

>> A=[1 2 1;2 6 1;1 1 4];
>> b=[2;7;3];
>> x=A\b

在这里插入图片描述
Cramer’s(Inverse)Method:
Solving Equations Using Cramer’s Method:

>> A=[1 2 1;2 6 1;1 1 4];
>> b=[2;7;3];
>> x=inv(A)*b

在这里插入图片描述
Singular:
inv(A)未必存在

>> A=[1 2 3 4;2 4 6 8;9 8 7 6;1 3 2 8];
>> inv(A)

在这里插入图片描述

>> x=0:40;
>> y=0:40;
>> [X,Y]=meshgrid(x,y);
>> Z1=-X-Y;
>> surf(X,Y,Z1);
>> hold on
>> Z2=-X+Y;
>> surf(X,Y,Z2);
>> Z3=-1/3*X;
>> surf(X,Y,Z3);
>> zlim([-20,20]);

在这里插入图片描述
Problem with Cramer’s Method:
Functions to Check Matrix Condition:
k值越大越不健康

Linear System:
Eigenvalues and Eigenvectors:
Interpretation of Figenvalues and Eigenvectors:
Solving Eigenvalues and Eigenvecors:
eig()

Matirix Exponential:expm()

A=[0 -6 -1;6 2 -16;-5 20 -10];
x0=[1 1 1]';
X=[];
for t=0:.01:1
    X=[X expm(t*A)*x0];
end
plot3(X(1,:),X(2,:),X(3,:),'-o');
xlabel('x_1');
ylabel('x_2');
zalbel('x_3');
grid on;
axis tight square;