搜索
1.已知系统质量矩阵M和刚度矩阵K，试用矩阵迭代幂法编程上机计算其固有频率与对应正则振型（要求独立完成，并打印出源程序与结果清单）。30010-3-2M=010K=-38-40020-47matlab编程如下：M=[300;010;002];%质量K=[10-30;-38-4;-2-47];%刚度U0=[1;2;3];%定义一个值D=K\M;fori=1:3;lamda=max(D*U0);u1=1/lamda*D*U0;while(norm((u1-U0),2)1e-4)U0=u1;lamda=max(D*u1);u1=1/lamda*D*u1;endw(i)=1/sqrt(lamda);%公式1disp(w(i));%打印出固有频率M1=u1'*M*u1;phai(:,i)=1/sqrt(M1)*u1;%公式2disp(phai(:,i));%打印出固有频率对应的阵型D=D-lamda*phai(:,i)*phai(:,i)'*M;end结果如下1.2057%第一阶0.2401%对应的阵型0.45120.55831.9265%第二阶0.5665%对应的阵型0.1850-0.03901.0985e+004%第三阶0.5665%对应的阵型0.1850-0.03902.已知系统质量矩阵M和刚度矩阵K，试用矩阵迭代幂法编程上机计算其固有频率与对应正则振型（要求独立完成，并打印出源程序与结果清单）。102-1M=K=02-12编程如下：M=[10;02];%质量K=[2-1;-12];%刚度U0=[1;2];%定义一个值D=K\M;fori=1:2;lamda=max(D*U0);u1=1/lamda*D*U0;while(norm((u1-U0),2)1e-4)U0=u1;lamda=max(D*u1);u1=1/lamda*D*u1;endw(i)=1/sqrt(lamda);%公式1disp(w(i));%打印出固有频率M1=u1'*M*u1;phai(:,i)=1/sqrt(M1)*u1;%公式2disp(phai(:,i));%打印出固有频率对应的阵型D=D-lamda*phai(:,i)*phai(:,i)'*M;end结果如下：0.7962%第一阶固有频率0.4597%第一阶对应的阵型0.62801.5382%第二阶固有频率0.8881%第二阶对应的阵型-0.32513.已知系统质量矩阵M和刚度矩阵K，试用矩阵迭代幂法编程上机计算其固有频率与对应正则振型（要求独立完成，并打印出源程序与结果清单）。100003-100001000-12-100M=00100K=0-12-100001000-12-100001000-13matlab编程如下：M=[10000;01000;00100;00010;00001];%质量K=[3-1000;-12-100;0-12-10;00-12-1;000-13];%刚度U0=[1;2;3;4;5];%定义一个值D=K\M;fori=1:5;lamda=max(D*U0);u1=1/lamda*D*U0;while(norm((u1-U0),2)1e-4)U0=u1;lamda=max(D*u1);u1=1/lamda*D*u1;endw(i)=1/sqrt(lamda);%公式1disp(w(i));%打印出固有频率M1=u1'*M*u1;phai(:,i)=1/sqrt(M1)*u1;%公式2disp(phai(:,i));%打印出固有频率对应的阵型D=D-lamda*phai(:,i)*phai(:,i)'*M;end结果如下：0.61800.19540.51170.63250.51170.19541.17560.37170.60150.0001-0.6015-0.37181.6180-0.5116-0.19550.6325-0.1953-0.51171.90210.6014-0.3715-0.00030.3720-0.60162.0000-0.44750.4474-0.44720.4470-0.4469