大体看了一下,
我认为你可以用方阵的变换把b提出来,去解不含有b的方阵,这样用eig函数很容易就可以实现了
这里提取b时,你最好能把方阵转为行列式,这里可能成立。
这里好像在说大话,不过我还是要说,关于建模编程有很多不是一两个小时能做好的,甚至有很多是解不出来的,我对题目了解的不多,所以给你去编程是没有意义的,还不如留给你自己去做,锻炼一下。
祝9.11能一路顺风。
我想建一个群,关于9.11时的技术交流,感兴趣的话给我留言
最后说一下你的编程
for
i=3:62
l(i,i-1)=m;
end
%这个好像超出范围了,以后小心
还有你们为什么不用m文件,我感觉用m文件做程序有很多优点的,建议要经常使用(我一般做一个程序都是用m文件de
)
a= [1.0000 0.5000 4.0000 3.0000 3.0000
2.0000 1.0000 7.0000 5.0000 5.0000
0.2500 0.1429 1.0000 0.5000 0.3333
0.3333 0.2000 2.0000 1.0000 1.0000
0.3333 0.2000 3.0000 1.0000 1.0000];
>> [v,d]=eig(a)
其中v对应特征向量阵。d对应特征值对角阵,