%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