振动理论与应用 matlab 程序

振动理论及应用（第5版）——国际著名力学图书影印版系列，matlab 程序

%This program produces the mass and stiffness matrix for 
%a bar where each element is the same length but may
%vary in its properties. The bar is pinned at the left end. 
clear
ll=input（‘ What is the length of the bar?  ‘）
n=input（‘ How many elements would you like? ‘）
for i=1:n
 l（i）=ll/n;
end
mass=input（‘Enter  the mass of the bar as a vector  ‘）
for i=1:n
m（i）=mass（i）*l（i）;
end
E=input（‘Enter the elastic modulus as a vector ‘）
a=input（‘Enter the cross-sectional area as a vector  ‘）
for i=1:n
k（i）=E（i）*a（i）/l（i）;
end
mm=zeros（n+1）;
kk=zeros（n+1）;
mm（11）=2*m（1）;
mm（12）=m（1）;
for i=2:n
 mm（ii-1）=m（i-1）;
 mm（ii）=2*m（i-1）+2*m（i）;
 mm（ii+1）=m（i）;
end
mm（n+1n）=m（n）;
mm（n+1n+1）=2*m（n）;
mm=（1./6.）*mm;
kk（11）=k（1）;
kk（12）=-k（1）;
for i=2:n
kk（ii-1）=-k（i-1）;
kk（ii）=k（i-1）+k（i）;
kk（ii+

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    .......      5606  1997-10-27 00:00  thoery of vibration\Archive.sit

    .......       975  1997-10-27 00:00  thoery of vibration\bar.m

    .......      1304  1997-10-27 00:00  thoery of vibration\beam.m

    .......       917  1997-10-27 00:00  thoery of vibration\choljac.m

    .......       143  1997-10-27 00:00  thoery of vibration\f.m

    .......       142  1997-10-27 00:00  thoery of vibration\force1.m

    .......       143  1997-10-27 00:00  thoery of vibration\force2.m

    .......       143  1997-10-27 00:00  thoery of vibration\force3.m

    .......       621  1997-10-27 00:00  thoery of vibration\iterate.m

    .......       765  1997-10-27 00:00  thoery of vibration\mykl.m

    .......       741  1997-10-27 00:00  thoery of vibration\polyn.m

    .......      1373  1997-10-27 00:00  thoery of vibration\runga.m

    .......       398  1997-10-27 00:00  thoery of vibration\tor.m

     目录          0  2011-01-27 17:05  thoery of vibration

----------- ---------  ---------- -----  ----

                13271                    14