matlab 用eig和diag解微分方程
来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/06/28 09:09:10
matlab 用eig和diag解微分方程
![](http://img.wesiedu.com/upload/1/4b/14beafdad31fbc465fda813bd0670e6c.jpg)
![](http://img.wesiedu.com/upload/1/4b/14beafdad31fbc465fda813bd0670e6c.jpg)
A=[-5 3; 100 -301];
[V D]=eig(A);
lambda=diag(D);
C=V\[50;100];
t=linspace(0,1,1000)';
y1=C(1)*V(1,1)*exp(lambda(1)*t)+C(2)*V(1,2)*exp(lambda(2)*t);
y2=C(1)*V(2,1)*exp(lambda(1)*t)+C(2)*V(2,2)*exp(lambda(2)*t);
plot(t,[y1 y2]);
legend('y1','y2');
再问: 如果不看hint可以知道要用什么矩阵吗
再答: 矩阵就是方程的系数啊
再问: 我是说为什么v1等于这个矩阵的行向量。算了基本明白了
[V D]=eig(A);
lambda=diag(D);
C=V\[50;100];
t=linspace(0,1,1000)';
y1=C(1)*V(1,1)*exp(lambda(1)*t)+C(2)*V(1,2)*exp(lambda(2)*t);
y2=C(1)*V(2,1)*exp(lambda(1)*t)+C(2)*V(2,2)*exp(lambda(2)*t);
plot(t,[y1 y2]);
legend('y1','y2');
![](http://img.wesiedu.com/upload/4/2c/42cfc1ca8d951b3501b98d5ccf5b3e80.jpg)
再问: 如果不看hint可以知道要用什么矩阵吗
再答: 矩阵就是方程的系数啊
再问: 我是说为什么v1等于这个矩阵的行向量。算了基本明白了