matlab求数值线性代数问题

yzddaniel

考虑一端固支的悬臂梁，长度为1米，截面工字，相关数据为高度0.05米，上、下宽都为0.034米，厚度分别为0.00225米、0.0038米和0.0038米。自由端有一质量为2千克的集中质量，材料的杨氏模量为2.0e11，密度为7800kg/m3，Possion率0.33(图例见下页)。将此工字梁划分为40个节点，用有限元方法算出其前10阶振动频率和模态，以及质量阵和刚度阵。具体数据见design.mat.
1.  redgongwithstiff.mat(K)为提取出的刚度矩阵
2.  redgongwithmass.mat(M)为提取出的质量矩阵
3.  redgongwithmode.mat(X)为提取出的前十阶模态矩阵
4.  gongwithfre.mat为提取出的前十阶频率矩阵
5. gongwitheig.mat(D)为利用频率转换得到的特征值矩阵
这些是给出的已知条件,机子没装MATLAB放不上来
1.分别计算质量矩阵M和刚度矩阵K的谱半径；
2.**给出一种计算关于质量矩阵和刚度矩阵的按模最小的非零广义特征值的方法(误差小于3%)；
3.试估计前五个特征值所对应的特征向量，(此时可利用已给的测量特征值数据)，并给出估计结果与测量数据之间的相对误差；
4***.一般说来，测量数据应满足KX=MXD，但题目中给出的数据并不严格满足这一条件,试给出一种方法降低误差，要求不能改变X,D，可以改变K,M.
第1个不用求了,主要是第2和第3个,第二个最好用QZ方法做,还有第3个误差用什么方法算,最好给下程序,或者阐述下具体做的思想

sigma665

没看明白你的问题是什么

yzddaniel

额,其实前面的无所畏的,可以直接看下面的四个问题,比如2.**给出一种计算关于质量矩阵和刚度矩阵的按模最小的非零广义特征值的方法
给个程序源代码什么的,我可以自己找台机子做结果
矩阵太大传不上来

[ 本帖最后由 yzddaniel 于 2008-6-13 19:31 编辑 ]

sogooda

:@L 太专业了，不懂。

yzddaniel

555555555555,太伤心,都没人能帮我解决吗,MATLAB的高手们千万不要吝啬啊