Matlab'da Yapısal Bir Sistemin Dürtü Yanıtını Elde Etme

Kütle(Mass), Rijitlik(Stiffness) ve Sönümleme(Damping) matrisleri bilinen bir dinamik sistemin dürtü yanıtını(Impulse Response) çizen bir Matlab fonksiyonu. Yapısal(Structural) sistemler üstünde çalışırken bu kod parçasını sıklıkla kullanırım. Tek yapmanız gereken yarattığınız sistemin boyutunu ve yapısal özelliklerini sırayla Matlab formatında girmek.

clear all
N = input('model order :  '); 
M = input('mass matrix :  ') ;
C = input('damping matrix :  ') ;
K = input('stiffness matrix :  ') ;

A1=inv(M)*C;
A0=inv(M)*K;
B0=inv(M);
I=eye(N,N);
Z=zeros(N,N);

A=[Z I;-A0 -A1]; 
B=[Z;B0];   
C1=[I Z];   % continous time system
D=0;

sys=ss(A,B,C1,D);

dt=0.01; %sampling time 
duration=20; %duration of the signal
t=0:dt:duration;
y=impulse(sys,t); %impulse response