- 积分
- 0
- 注册时间
- 2003-5-13
- 仿真币
-
- 最后登录
- 1970-1-1
|
!*******transient analysis*******
fini
/clear,nostart
/filname,frf_test
pi=acos(-1)
/prep7
et,1,21
et,2,14
r,1,1,1,1
r,2,2*pi*2*pi,0
n,1,
n,2,1
type,1
real,1
e,2
type,2
real,2
e,1,2
nsel,s,,,1,2
d,all,uy
d,all,uz
d,1,all
alls,
/solu
antype,modal
modopt,lanb,2
solve
fini
/solu
antype,trans
trnopt,msup,1
outres,nsol,all
autots,off
timint,on
!frequency range: 0.1 ~ 10 Hz
!so omage's range: 2*pi*0.1 ~ 2*pi*10
!define 400 frequency points
!middle values is 2*pi*0.1 /4
mid_val=2*pi*0.1 /4
time,mid_val
deltim,mid_val
kbc,0
lswrite,1
*do,i,2,400
time,i*mid_val
f,2,fx,sin(i*mid_val)
lswrite,i
*enddo
lssolve,1,400
!****matlab code*****
time=[];
ux=[];
m=(1/0.15708/400)1/0.15708/400):(1/0.15708);
n=400;
w1=window(@hamming,n);
jiatime=w1.*time;
jiaux=w1.*ux;
transftime=abs(fft(jiatime));
transfux=abs(fft(jiaux));
transf=transfux./transftime;
plot(m(1:200),transf(1:200));
loglog(m(1:200),transf(1:200)) |
|