04b-0-机械振动2015121.pptVIP

强迫振动-3 2015.12.11 韩清凯 * clc clear % syms a b XX=(a+i*b)*(a-i*b) expand(XX) % syms m c k f0 w syms A phi H1=1/(k-m*w*w+i*c*w) H_below=((k-m*w*w+i*c*w)*(k-m*w*w-i*c*w)) H_up=(k-m*w*w-i*c*w) H2=H_up/H_below realH2=(k-m*w*w)/H_below imagH2=-c*w/H_below A=sqrt(realH2^2+imagH2^2) tanPhi=imagH2/realH2 响应 振幅放大因子 相位差角 syms phi s ksi syms k H=((1-s^2)-i*(2*ksi*s))/((1-s^2)^2+(2*ksi*s)^2) realH=(1-s^2)/((1-s^2)^2+(2*ksi*s)^2) imagH=-(2*ksi*s)/((1-s^2)^2+(2*ksi*s)^2) phi=atan(imagH/realH) %注意公式中的角度为负号 beta=1/sqrt(((1-s^2)^2+(2*ksi*s)^2)) a=beta*cos(phi) b=beta*sin(phi) a+i*b % =