chrip的博客

原创 时变、状态空间系统

close all; clc; clear; %% %%系统%% m1=1;m2=2; c1=0.01;c2=0.02; k1=1000;k2=2000; M=[m1,0;0,m2]; for Theta=1:10 tau=5Theta c1=0.01+0.01tau/600 c2=0.01+0.01tau/600 for i=1:5 C{(Theta-1)5+i}=[c1,0;0,c2]; end end K=[k1,-k1;-k1,k1+k2]; n=size(M,1); %% %%生成状态空间系统%%