可参考文涛,基于Matlab语言的多自由度振动系统的固有频率及主振型计算分析,2007 对于无阻尼系统 [VEC,VAL]=eig(inv(A)*K) 对于有阻尼系统,参考振动论坛计算程序 输入M,D,K function [v,w,zeta]=vbr_sf(m,d,k) %vbr_sf vbr_sf(m,d,k) % [v,w,zeta]=vbr4(m,d,k) % function vbr_sf finds the mode shapes and natural…
function [ y ] = moving_average( x, win_size ) y1=filter(ones(1,win_size/2+1)/win_size,1,x); y2=filter(ones(1,win_size/2+1)/win_size,1,fliplr(x)); y=y1+fliplr(y2)-(1/win_size)*x; end L = win_size/2; yi = (1/win_size)∑(j=i-L:i+L)xj 本系列文章允许转载,转载请保留全文…