matlab--program 1.纽马克算法

1.纽马克算法，振动分析方法，计算多自由度系统振动响应，位移，速度，加速度 2.用于计算多自由度系统固有频率、固有振型、模态质量、模态刚度矩阵、正则模态向量、受迫响应振动 3.里兹法是一种缩减系统自由

clc
clear
M=[1234];       
m=diag（M）;          %质量对角线矩阵
k=[ 800   -800   0      0    ;
    -800  2400   -1600  0    ;
    0     -1600   4800 -3200 ;
    0     0      -3200 8000  ];    %计算刚度矩阵
c=0.05*m+0.02*k;%c=2*m*0.05*sqrt（k/m） ;            
f0= 100;
t1=5;                %力的作用时间                 
nt=2000;              %分2000步完成               
dt=0.01;              %时间步长               
alfa=0.25;             %γ=0.25           
beta=0.5;              %β=0.5            

a0=1/alfa/dt/dt;         %              
a1=beta/alfa/dt;         % 
a2=1/alfa/dt;            % 
a3=1/2/alfa-1;            % 
a4=beta/alfa-1;          % 
a5=dt/2*（beta/alfa-2）;     %  
a6=dt*（1-beta）;         % 
a7=dt*beta;            %  
d=zeros（4nt）;           %初位移为0                 
v=zeros（4

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-03-24 15:35  matlab程序\
     目录           0  2016-03-24 14:57  matlab程序\Newmark数值分析法\
     文件        2783  2016-03-24 11:10  matlab程序\Newmark数值分析法\Newmark_suanfa.m
     目录           0  2016-03-24 14:57  matlab程序\庞加莱截面图法\
     文件         124  2016-03-24 08:59  matlab程序\庞加莱截面图法\zjzdfun.m
     文件         732  2016-03-24 13:56  matlab程序\庞加莱截面图法\zjzd_solution.m
     目录           0  2016-03-24 14:57  matlab程序\线加速度法求系统振动响应\
     文件        1072  2016-03-24 14:23  matlab程序\线加速度法求系统振动响应\M635_1.m
     目录           0  2016-03-24 15:35  matlab程序\计算多自由度系统模态质量模态刚度矩阵\
     文件        1195  2016-03-24 15:34  matlab程序\计算多自由度系统模态质量模态刚度矩阵\jisuanduoziyouduxitong.m
     目录           0  2016-03-24 14:56  matlab程序\里兹法求多自由度系统前两阶固有频率及振型\
     文件         951  2016-03-24 14:47  matlab程序\里兹法求多自由度系统前两阶固有频率及振型\lizifaqiuguyoupinlv.m
     目录           0  2016-03-24 14:55  matlab程序\雨流计数法\
     文件        9216  2013-12-09 16:56  matlab程序\雨流计数法\rainflow.dll
     文件         926  2013-12-09 16:56  matlab程序\雨流计数法\rainflow.m
     文件         896  2016-01-18 10:32  matlab程序\雨流计数法\rfdemo2.m
     文件        1941  2013-12-09 16:56  matlab程序\雨流计数法\rfhist.m
     文件        3151  2013-12-09 16:56  matlab程序\雨流计数法\rfmatrix.m
     文件        2944  2013-12-09 16:56  matlab程序\雨流计数法\sig2ext.m
     文件          19  2016-03-24 14:55  matlab程序\雨流计数法\运算说明.txt
     目录           0  2016-03-24 14:58  matlab程序\频谱图法\
     文件         292  2016-01-22 19:39  matlab程序\频谱图法\pinputu.m