跪求用Riccati传递矩阵法计算磁悬浮转子临界转速的Matlab的程序，越详细越好。万分感谢！

%求解转子系统前三个临界转速和主振型的传递矩阵法

%等截面轴参数

l1=0.12d=0.04A=pi*d*d/4

%轮盘参数

D=0.5h=0.025%盘轴材料参数（忽略轴的质量）

a=1u=0.3rou=7800E=2.0e11G=E/(2*(1+u))I=pi*(d^4)/64K1=2.0e7

v1=6*E*I/(a*G*A*l1*l1)

mi=rou*pi*D^2/4*h%轮盘的集质量

Jp=mi*D^2/8 Jd=Jp/2Ji=Jp-Jd%参数的数组形式

L=[l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 l1 0 0]M=[0 mi mi mi mi mi mi 0 0 0 0 0 mi mi 0]K=[K1 0 0 0 0 0 0 K1 0 0 0 K1 0 0 0]

v=[v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 v1 0 0]J=[0 Ji Ji Ji Ji Ji Ji 0 0 0 0 0 Ji Ji 0] k=0

Tit=['第一阶频率贺没漏的振型和弯矩图''第二阶频率的振型和弯矩图''第三阶频率的振型和弯矩图']

for w=0:0.01:4000for i=1:15

T(:,:,i)=[1+(L(i)^3)*(1-v(i))*(M(i)*w^2-K(i))/(6*E*I) L(i)+L(i)^2*J(i)*w^2/(2*E*I)

L(i)^2/(2*E*I) L(i)^3*(1-v(i))/(6*E*I) (L(i)^2)*(M(i)*w^2-K(i))/(2*E*I) 1+L(i)*J(i)*w^2/(E*I) L(i)/(E*I) L(i)^2/(2*E*I)

L(i)*(M(i)*w^2-K(i)) J(i)*w^2 1 L(i) M(i)*w^2-K(i) 0 0 1]end H=T(:,:,1)

for i2=2:15

H=T(:,:,i2)*Hend

F=H(3,1)*H(4,2)-H(3,2)*H(4,1)if F*(-1)^k <0 %求解临界转速察贺k=k+1 wi(k)=w w=wi(k)

ni(k)=wi(k)*30/piend end

for i1=1:3w=wi(i1)

for j=1:14

T(:,:,j)=[1+(L(j)^3)*(1-v(j))*(M(j)*w^2-K(j))/(6*E*I)

L(j)+L(j)^2*J(j)*w^2/(2*E*I)

L(j)^2/(2*E*I) L(j)^3*(1-v(j))/(6*E*I)

(L(j)^2)*(M(j)*w^2-K(j))/(2*E*I) 1+L(j)*J(j)*w^2/(E*I) L(j)/(E*I) L(j)^2/(2*E*I) L(j)*(M(j)*w^2-K(j)) J(j)*w^2 1 L(j) M(j)*w^2-K(j) 0 0 1]end

H=T(:,:,1)for j=2:15

H=T(:,:,j)*H end

b=-H(4,1)/H(4,2) X(:,1)=([1 b 0 0]') for n=2:16

X(:,n)=T(:,:,n-1)*X(:,n-1) %相邻两质点右边的传递关系end for j1=1:15y(j1)=X(1,j1) z(j1)=X(3,j1) x(j1)=(j1-1)*l1 end

y(16)=X(1,16)

x(16)=1.56z(16)=X(3,16)

y=y/max(abs(y))%归一禅烂化 z=z/max(abs(z))subplot(3,1,i1)

plot(x,y,'b-',x,z,'r:') title(Tit(i1,:))

xlabel('轴长'),ylabel('不平衡值') axis([0,1.56,-1.2,1.2]) grid on z

end

超引力磁悬浮蓝牙音箱主要使用说明：

开机启动：长按开机键2秒，听到飞碟音效后启动完成；

蓝牙匹配：开机后手机搜索“Super Gravity”，完成链接；

悬浮定位：双手平托飞碟，缓慢放置底托正中位置，待距离至2厘米处观察底托四角灯带，全部亮起后可放开；

音频播放：打开音乐或视频软件，即可开始播放；老侍余

底座通电：标配充电器，谈斗悬浮阶段由底座为磁悬浮模块供侍滚电；

飞碟充电：低电量提示音后，将飞碟放稳后连接充电。

海尔中央空调变频离心机组是磁悬浮

其特点是：节能高效

机组采用磁悬浮压缩机技术、变频控制技术、毁迟无油润滑等先进技术，产品能效比有了很大的提高。李余镇机组部分负荷最高能效比达11.3，综合能效比达到9.55（按照ARI标准测试）

静音无振动

完全无摩擦运行，运行噪音低，结构振动接近0，无须昂贵的减震配件

绿色环保

机组采用环保冷媒R134a，对臭氧层耗损值（ODP）为0。属正压型冷媒，避免了系统混入空气的危险

高效无摩擦

使用先进的磁轴承技术，实现机组的无油运行。完全避免常规压缩机轴承的高摩擦损失、润滑油的管理与控制，确保了卓越的能效及哪粗可靠、长期的无磨损运行

自由式容量调节

在冷凝温度下降或热负荷下降的情况下，降低压缩机的转速，从而可实现额定负荷10%-100%的宽负荷范围内自由调节冷媒的输出，优化压缩机的能耗

抗喘振

压缩机控制模块中提供了压缩机安全运行的控制曲线，通过实时监测压缩机的运行状态，计算判断后对转速进行及时调整，确保压缩机始终在安全区域内。

最小电网冲击

启动电流低，单台压缩机启动电流只有2A。由于采用磁悬浮轴承，单台压缩机开机时只需要2A的电流使转轴悬浮起来，启动扭矩很小，实现对电网的冲击最小

寿命长

压缩机由航空等级的铝制铸件及高强度的热塑电子外壳制造而成，可以长期高效运行

来源于http://www.haiersyrb.com/index.aspx

---