matlab符号计算求解两自由度振子传递函数.doc

alljoyland

syms m1
m2
c
k1
k2
Y
w
F
den
num

M =[ m1 0; 0 m2]

C =[ c -c ; -c c]

K = [ k1,-k1;-k1,k1+k2]

U0 = -w^2*M + i*w*C +K

%拉普拉斯变换
然后把s 换为jw

U = inv(U0)*F

U = inv(U0)*F

detU=det(U)

numden(det(U))

[num,den]=numden(det(U))

% den = -den

collect(den,w)

U1 = U/F

det(U1)

temp_num1=[0;k2*Y]

%Cramer 法则,解代数方程一样

num1= [temp_num1,U1(:,2)]

%det(num1)/det(U1)

det(num1)/det(U0)

%缺少一个比例常数

num1=[temp_num1,U0(:,2)]

det(num1)/det(U0)

%这个才是对的了

%因为U1已经求过逆了

% because the U1 has been the inverse of the 2dof

% so you can you use the

F =[0;k2*Y]

U1*F

%另一种方法
就是矩阵的
求逆和乘法而已, 可以这样子求解啦

% the inv(U1)*F is just the numerator not the denominator

%其中非常值得注意的是,我们的这种谐波激励求解方法也默认结果是谐波

%的,所以其中
最后的解
是
取
虚部的 ,因为假定激励是sin 不是cos的

%否则应该取实部