fly.m

来自「四面体有限元程序,单元刚度的计算及总刚度矩征集成」· M 代码 · 共 9 行

%%program fly.m
x=[4 3 7 -5 1 -3 -4 2]';
Y=fft(x)
X=ifft(Y)
Y(1)=[];
n=fix(length(Y)/2);
freq=(1:n)/length(Y);
power=abs(Y(1:n).^2/(length(Y).^2));
power=power'