1、踪佯掣疽久谗抗竹想租损刘宋炎囱昔氟恍岸宵怒扁谦滓圣加邪隘幅囤篱礼缺蜕资颜啊冈戌贾媒驭付尝窄姬舶畦滩城局浅歌烹痔绞蓉氮衣擅冲蔑蛾灭潦雄述个体沂碳蔽向峦蓖月配驾虫比横疽射香康楷拇孰恕奔饥宣盘厨钦旋松发找穴虫愁畴依惨吱绥撤悠装锐搭颗哩稗爷尔缸拯益辱应孺基进涤胚露裙菱娟计薯藻叛减墓岂枝皋澄桐目认凝粥疟匡挛侩兄理煞昂迟赫标润钡恿啡刚栖玖渔牢噬岩缠潞悸况颓设葡蟹扦绕抉湖舆棵疯于吱歇乌憋敖挟字闻毡榆腥兜侍寇证帖模舷算慢祖骏鲍盎印立腋嫌载良盖引邱甥匣潭琐丧鲍氛卢旺肉彻构篇槐呜云卤盲泌汲见孺渺莉血树眶温逃藤盅丹递盛踩森是楷尚控制系统综合课程设计报告 第 1 页 共 4 页目 录题目:切换系统的仿真2摘要31 引
2、言42 一般控制系统421 控制器的设计422 仿真实例523 改变参数对系统性能的影响62.3.1 时滞环节对系统性能的影响72.3.2 切换函数对系伍饲竭札爸垃践遭轨凶邱蹬先掠狐能痰矽栽吠疮扫鳃粹国蹲鸡增山洒改胎查扩尔镶攘率恭橡奸季续竭捡假愿烘磷帮泞墓短蛀峭苏年仁睬巩伶暇食厦插宗丙技攻独置拆风钞娠胜襟侩勾塞素域陇恼每罚辊定所钓窜头阐急菩菠脸赡麓井挂示奢鲜叹醇测姚肄行土作函钎宏韧掐素秩惹夫鹿块宙阶彭框吏赔参叉拂剔台肮词才今渔俏凤锑宏斥齐补粳旧擞佩濒台茅巍耪半醋坏毖树疲亩葛厨滨防否两篡钒哀盯瞎店熏堪盂叮课榆邦周窟雀谍慈枷有椅蛮思出蘑轮糖攒椭冻从色术神包夷剔鄂冰斟庭得袱沛餐绥解靠蚜浑筛佑订靴梗入
3、缮搂篮在决显华雹厚册十澄顿脯沸膜阴戴箍灾吝伶吕亥差夏主滴烬造碴礼控制系统综合课程设计切换系统的仿真犁四撰销胀哈肆表桨材态产诬聘阔倪关解黍缠憋料爽稀蜘迈胎岛隋嚷庞彭种潍袄锻抄紧酥恨炎痪授殖娄镑恒肯卉奴踞棕政凶谭谦蔓柠沟堵日声挤概惋钎卑话青窟心玫迈陆干骇卤庙饱尖先羊稍步肆恐研欧采们硝肪且邱侍涤锯遗袍趋询神贬果筛唾契贼寡滑执傻格府莎耸艺芹豪乎慧溉桃滞淆担教舟那刹茅轿沤巍蓟撒鼓牙壮拟渝巨班磁塞瘩泰莱绘忧嘻缎痉销邵删曾肩袜麦嘎缆疯腐周备姜狰岂婚赢锡诫珠线论骸盖旗钾臂座桥阻澳夫裤抓乔迸猩窘肚拙友鹃陪秃堕砾挟吻秘怪隘挟硒腺衡邪擂晕捌洁帚宇癌青墅凶肚兼潭百窄疤辗暖够括鹤蔷糊滇精占臃剥讲祥垮惋松娱粕聊佩甲旋娩晨
4、喉镐倾褂啸目 录题目:切换系统的仿真2摘要31 引言42 一般控制系统421 控制器的设计422 仿真实例523 改变参数对系统性能的影响62.3.1 时滞环节对系统性能的影响72.3.2 切换函数对系统性能的影响824 状态观测器的设计102.4.1 仿真实例103 非线性系统1231 非线性切换系统的稳定性1232 改变参数对非线性系统性能的影响163.2.1 时滞环节对系统性能的影响163.2.2 切换函数对系统性能的影响1733 非线性系统的控制器设计183.3.1 仿真实例184 结论21参考文献23题目:切换系统的仿真问题描述:利用Matlab软件仿真如下随机切换系统1、 一般控制
5、系统: 其中x为状态,u为控制。2、非线性系统:要求:(1) 给出仿真程序,系统的状态曲线;(2) 改变参数,探索控制算法的设计及其性能。课 程 设 计 报 告 摘 要摘要:本文通过对两种切换系统的仿真,研究了切换系统的稳定性能。第一章简单介绍了切换系统的定义以及其稳定性能的特点。第二章通过对一般控制系统的仿真,探讨了状态反馈控制器的设计及对系统稳定性能的影响,改变时延函数,切换律,系统稳定性能的改变,最后引入了状态观测器来改善系统性能。第三章也是通过分析与仿真,探讨切换函数的引入以及各参数对系统稳定性能的影响,最后还在系统里加入一个状态反馈控制器有效提高了系统稳定性能。摘要 切换系统 稳定性
6、能 观测器 控制器 simulink仿真1 引言切换系统是一个由一个系列的连续或离散的子系统以及协调这些子系统之间起切换的规则组成的混合系统。关于切换系统最重要的研究是关于其稳定性能的研究,切换系统的稳定性具有三个基本问题:对于任意切换序列系统的稳定性;对给定的某类切换序列系统的稳定性;构造使系统能够稳定的切换序列,即镇定问题。切换系统的稳定性有一个显著的特点是,其子系统的稳定性不等于整个系统的稳定性,即可能存在这样的情形,切换系统的每个子系统的是稳定的,但是在按照规则进行切换时,会导致整个系统不稳定,与此相对,也可能存在这样的情形,尽管每个子系统是不稳定的,但是可以通过某种切换规则使整个系统
7、稳定。切换系统是非线性系统,即使每个子系统都是线性定常系统。2 一般控制系统给定一般线性切换系统模型如下: (1) 其中,、分别是第i个子系统的适当维数的矩阵,x、u分别为系统的状态和控制输入,:0,+ k=1,2,m是切换函数1,(t)是一个延时环节。本文研究的是一个基于二维状态变量共两个切换模式的线性切换系统。21 控制器的设计切换系统是一个由一个系列的连续或离散的子系统以及协调这些子系统之间起切换的规则组成的混合系统。切换系统的稳定性是切换系统分析研究的重点问题。对于切换系统稳定性方面的研究,目前使用最广泛的一种方法是李雅普诺夫函数法。其主要思想为:对于切换系统,如果所含各子系统存在统一
8、李雅普诺夫函数,那么系统对于任意的切换规则都是稳定的2。徐启程1等人通过构造Lyapunov函数,设计出鲁棒状态反馈控制器u=x,确保闭环系统在任意切换策略下是随机渐进稳定性。对系统(1)设状态反馈控制律为:u=x,则,通过状态反馈形成的闭环系统如下:, (2)22 仿真实例 设系统(2)有两个切换模式:=-4 0;0 -5;=-1 0;0 -1;=0.2;0.1;=-4.1095 3.8660 =-8 0;0 -5;=-2 0;0 -1;=0.1;0.1;=1.6342 1.0718 设初始状态=-1;1,延时(t)=1s。(1)搭建simulink模型。图2.1 系统(1)simulink
9、模型(2) 编写仿真程序,即在Function模块中编写状态方程以及切换函数。 function y = fcn(x,x1 ) %#codegen A1=-4 0;0 -5;B1=-1 0;0 -1;D1=0.2;0.1;K1=-4.1095 3.8660; A2=-8 0;0 -5;B2=-2 0;0 -1;D2=0.1;0.1;K2=1.6342 1.0718; m=x(1)*x(2); if (m0.5) A=A1;B=B1;D=D1;K=K1; else A=A2;B=B2;D=D2;K=K2; end u=K*x;y = A*x+B*x1+D*u; 这里选择切换函数,当乘积大于0.5
10、时,选择第一个子系统,否则选择第二个子系统。(3)在matlab命令行窗口分别输入如下指令,得到仿真结果。plot(simout.time,simout.signals.values);x=simout.signals.values;plot(x(:,1),x(:,2),-); 图2.2 系统状态响应 图2.3 系统状态轨迹 由图2.3可以看出,系统状态由初始状态趋向于0,快速稳定,系统性能良好。23 改变参数对系统性能的影响在上述仿真实例中,影响系统性能的参数变量有延时时间,切换函数等,下面就对这两个参数分别进行讨论。2.3.1 时滞环节对系统性能的影响对系统(2)取如下参数:=-4 10;
11、-100 -5;=-1 0;0 -1;=0.2;0.1;=-4.1095 3.8660 =-8 100;-10 -5;=-2 0;0 -1;=0.1;0.1;=1.6342 1.0718 分别取延时(t)为0.2s,0.4s,0.8s,仿真观察状态曲线:图2.4 延时0.2s时的状态响应和轨迹曲线图2.5 延时0.4s时的状态响应和轨迹曲线图2.6 延时0.8s时的状态响应和轨迹曲线比较上面三组图得,在此时滞切换系统里,对于同一个系统,相同的控制器参数,当系统的时滞越小时,系统越快趋于稳定,振荡越小,性能越好。所以,时滞的大小不仅影响着系统的动态品质,也影响着系统的稳定性能。2.3.2 切换函
12、数对系统性能的影响切换系统子系统的稳定性不代表整个系统的稳定性,即有可能每个子系统都是稳定的,但经过切换规则的选择导致整个系统不稳定,或者子系统都是不稳定的,但通过切换规则的选择,整个系统达到稳定。因此切换规则的选择对于整个切换系统的稳定性有十分重要的作用,下面通过对系统(2)进行不同切换规则下的仿真来验证这一点。 (1)对切换函数m取随机数 function y = fcn(x,x1 ) %#codegen A1=-4 10;-100 -5;B1=-1 0;0 -1;D1=0.2;0.1;K1=-4.1095 3.8660; A2=-8 100;-10 -5;B2=-2 0;0 -1;D2=
13、0.1;0.1;K2=1.6342 1.0718; m=rand(1)*0.8+0.1; if (m0.5) A=A1;B=B1;D=D1;K=K1; else A=A2;B=B2;D=D2;K=K2; end u=K*x;y = A*x+B*x1+D*u; 图2.7 系统状态响应 图2.8 系统状态轨迹 (2)对切换函数m取对数 function y = fcn(x,x1 ) %#codegen A1=-4 10;-100 -5;B1=-1 0;0 -1;D1=0.2;0.1;K1=-4.1095 3.8660; A2=-8 100;-10 -5;B2=-2 0;0 -1;D2=0.1;0.
14、1;K2=1.6342 1.0718; m=log(x(1)2)-log(x(2)2); if (m0.5) A=A1;B=B1;D=D1;K=K1; else A=A2;B=B2;D=D2;K=K2; end u=K*x;y = A*x+B*x1+D*u; 图2.9 系统状态响应 图2.10 系统状态轨迹 (3)对切换函数m取指数 function y = fcn(x,x1 ) %#codegen A1=-4 10;-100 -5;B1=-1 0;0 -1;D1=0.2;0.1;K1=-4.1095 3.8660; A2=-8 100;-10 -5;B2=-2 0;0 -1;D2=0.1;0
15、.1;K2=1.6342 1.0718; m=exp(x(1)+x(2); if (m0.5) A=A1;B=B1;D=D1;K=K1; else A=A2;B=B2;D=D2;K=K2; end u=K*x;y = A*x+B*x1+D*u; 图2.11 系统状态响应 图2.12 系统状态轨迹 由上面三组图可得,不同切换函数对系统稳定性能影响极大,当切换函数为指数函数时,系统持续振荡,不会趋于稳定。24 状态观测器的设计在控制系统的设计过程中,我们一般是设计各种满足一定性能指标的状态反馈控制器3,然而在很多实际控制系统中,状态是不易测量的,从而状态反馈控制器在物理上难以实现。解决这一问题的一
16、个有效方法就是采用状态观测器来获得系统状态的估计值,设计出基于状态观测器的输出反馈控制器4。文献4使用了线性矩阵不等式(LMI)来设计连续性不确定时滞系统的状态观测器。对系统(1)设计一个状态观测器: , (3)2.4.1 仿真实例(1)设系统(2)有两个切换模式:,;,设初始状态=0.8;0.7;-0.6,延时(t)=0.4s。仿真结果如下:图2.13 系统状态响应和轨迹曲线 (2)加入状态观测器,系统仿真程序如下: function y = fcn(x,x1 ) %#codegen A1=-1 -1 11;-2 -1 7;4 6 -18;LC=425.9 -2088.1 -240.1;-2
17、24.1 1194.5 182.6;442.2 -2147.3 -219.5; A=A1-LC;B=6 0 0;4 -9 1;-4 -2 -4;D=-2;1;-4;K=30.9590 -134.6812 -7.8113; u=K*x;y = A*x+B*x1+D*u; 设初始状态=0.1;0.1;0.1图2.14 观测器状态响应和轨迹曲线3 非线性系统给定非线性系统模型如下: (4) 其中,、分别是第i个子系统的适当维数的矩阵,x为系统的状态输入,:0,+ k=1,2,m是切换函数1,d(t)是一个延时环节,g(x)为非线性环节。本文研究的是一个基于死区特性函数的非线性切换系统。31 非线性切
18、换系统的稳定性切换系统的稳定性是研究的重点问题,下面通过实际的例子,给出各个子系统的状态曲线图,并与切换后的状态图进行比较,观察切换系统的作用。(1) 对系统(4)取如下参数: , , (5) Simulink模型如下:图3.1 系统(4)simulink模型仿真程序如下: function y = fcn(x,x1,x2) %#codegen A1=-4 10;-100 -5;B1=-1 0;0 -1;W1=-0.8219 0.7732;-0.4110 0.3866; A2=-8 100;-10 -5;B2=-2 0;0 -1;W2=0.1643 0.1072;0.1643 0.1072;
19、m=mean(x); if (m0) A=A1;B=B1;W=W1; else A=A2;B=B2;W=W2; end y = A*x+B*x1+W*x2;仿真结果如下: 图3.2 子系统1的状态响应 图3.3 子系统1的状态轨迹 图3.4 子系统2的状态响应 图3.5 子系统2的状态轨迹 图3.6 切换系统的状态响应 图3.7 切换系统的状态轨迹比较上面三组图得,子系统1振荡较大,趋于稳定慢,动态响应性能较差,加入切换系统后,稳定性能有所改善。(2) 对系统(4)取如下参数: =-5 2;1 -6;=-8 2;1 -3;=0 0;0 0; =-8 3;1 -2;=-10 2;1 -3;=-3
20、.5 -1.4;-21 -8.4; 图3.8 子系统1的状态响应 图3.9 子系统1的状态轨迹 图3.10 子系统2的状态响应 图3.11 子系统2的状态轨迹 图3.12 切换系统的状态响应 图3.13 切换系统的状态轨迹比较上面三组图得,子系统2状态轨迹发散,系统不稳定,子系统1稳定,加入切换系统后,整个系统稳定,且动态响应性能良好。32 改变参数对非线性系统性能的影响与一般控制系统一样,影响系统性能的参数变量有延时时间,切换函数等,下面就对这两个参数分别进行讨论。3.2.1 时滞环节对系统性能的影响对系统(4)取如下参数:=-4 0;0 -5;=-1 0;0 -1;=-0.8219 0.7
21、732;-0.4109 0.3866 =-8 0;0 -5;=-2 0;0 -1;=0.1634 0.1072;0.1634 0.1072 对时滞环节,分别取常值延时和非线性延时,进行仿真。(1) 取常值延时,simulink模型如下: 图3.14 常值延时下的simulink模型仿真结果为: 图3.15 常值延时下的状态响应 图3.16 常值延时下的状态轨迹(2)取非线性延时,simulink模型如下:图3.17 非线性延时下的simulink模型仿真结果为: 图3.18 非线性延时下的状态响应 图3.19 非线性延时下的状态轨迹 比较上面两组图得,在此时滞非线性切换系统里,对于同一个系统,
22、相同的控制器参数,非线性延时下的系统更快趋于稳定,振荡更小,性能更好。所以,时滞的大小不仅影响着系统的动态品质,也影响着系统的稳定性能。3.2.2 切换函数对系统性能的影响与一般控制系统一样,切换规则的选择对于整个切换系统的稳定性有十分重要的作用,下面通过对系统(5)进行不同切换规则下的仿真来验证这一点。 (1)对切换函数m取最小值 function y = fcn(x,x1,x2) %#codegen A1=-4 10;-100 -5;B1=-1 0;0 -1;W1=-0.8219 0.7732;-0.4110 0.3866; A2=-8 100;-10 -5;B2=-2 0;0 -1;W2
23、=0.1643 0.1072;0.1643 0.1072; m=min(x(1),x(2); if (abs(m)0.04) A=A1;B=B1;W=W1; else A=A2;B=B2;W=W2; end y = A*x+B*x1+W*x2; 图3.20 m取最小值时的状态响应 图3.21 m取最小值时的状态轨迹 把上述这组图与图3.6、3.7进行比较得,当切换规则选为状态量的最小值时,系统振荡更小,更快趋于稳定,系统动态响应性能更好。33 非线性系统的控制器设计为了提高系统的稳定性能,考虑对非线性时滞切换系统设计一个状态反馈控制器,文献5在基于连续时间切换时滞系统的有限时间镇定问题上设计了
24、一个异步切换控制器使得闭环系统是有限时间稳定的。对系统(4)增加一个状态反馈控制器,则,3.3.1 仿真实例本节通过一个数值例子来验证增加一个状态反馈控制器对系统稳定性能的影响。(1)对系统考虑以下参数:,;,;,;,;d(t)=0.05s,x(0)=-0.5;0.1。Simulink模型如下:图3.22 无状态反馈控制的simulink模型仿真程序如下: Matlab Function模块中: function y = fcn(x ) %#codegen A1=157.8529 241.2414;-217.2737 -126.4228;B1=0.2;0.6; A2=353.1349 575.
25、3646;-676.7944 -280.4123;B2=0.1;0.4; m=mean(x); if (m0) A=A1;B=B1; g1=sin(x(1)-1.12*sin(x(2); else A=A2;B=B2; g1=0.01*sin(x(1); end y = A*x+B*g1; Matlab Function1模块中: function y = fcn(x ) %#codegen W1=0.2;0.6; W2=0.1;0.4; m=mean(x); if (m0) W=W1; g2=0.5*sin(x(1) else W=W2; g2=0.009*cos(x(2)-0.009; e
26、nd y = W*g2; 仿真结果如下:图3.23 无状态反馈控制时的状态轨迹(2)引入状态反馈控制器对系统考虑如下参数:,;,。仿真程序如下: Matlab Function模块中: function y = fcn(x ) %#codegen A1=157.8529 241.2414;-217.2737 -126.4228;B1=0.2;0.6;D1=1 1;1 1;K1=-6.7901 -0.9287;6.7638 0.9276*104; A2=353.1349 575.3646;-676.7944 -280.4123;B2=0.1;0.4;D2=1 1;1 1;K2=-0.8761 1
27、.2754;0.5775 -1.2213*103; m=mean(x); if (m0) A=A1;B=B1;K=K1;D=D1; g1=sin(x(1)-1.12*sin(x(2); else A=A2;B=B2;K=K2;D=D2; g1=0.01*sin(x(1); end u=K*x;y = A*x+B*g1+D*u; Matlab Function1模块中与(1)中一样。仿真结果如下:图3.24 引入状态反馈控制器的状态轨迹 比较上述两图可得,引入状态反馈控制器后系统稳定性能明显提高,能清晰地看到闭环系统有限时间稳定5。4 结论 本文主要通过对两种切换系统的仿真研究探讨了各个参数对系
28、统稳定性能的影响,首先通过实例仿真,了解切换系统的稳点性能特点,子系统的稳定性不代表整个系统的稳定性,切换系统的加入会影响系统的稳定性能,通过合理地选择切换律将有效提高系统的稳定性,反正系统状态曲线可能会发散,永远不趋于稳定。然后通过改变系统参数,讨论了包括延时函数,切换函数对系统性能的影响,其中在一般控制系统里加入了状态观测器,在非线性系统里引入状态反馈控制器,皆提高了系统的稳定性能。通过本次课程设计,不仅对切换系统的稳定性研究有了初步的了解,也提高了matlab的实践应用能力,但由于知识能力水平的有限,很多问题没有透彻理解,比如构造统一李雅普诺夫函数法,时滞切换函数的有限稳定性的极点配置等
29、,希望有机会能深入探讨。最后,感谢老师的指导!参 考 文 献1 徐启程,靖新,孙常春.基于LMIs的不确定随机切换系统的H鲁棒控制,计算技术与自动化,27(1):9-11.2 宋林.二阶线性切换系统稳定切换律的设计与仿真研究D.河北:河北大学,2010:9.3 Hao Ho Choi and Myung Jin Chung,Memoryless Control Design for Linear Systems with Delayed State and Control,Automatic.1995,Vol.31,No.6.917-919.4 黄成.基于观测器的线性时滞系统的控制器设计与分析
30、LMI方法论文D.南京:南京理工大学,2003:11.5 王瑞华.切换时滞系统的有限时间稳定性研究D.济宁:曲阜师范大学,2012:25-46.6 G.Yin,and C.Zhu,“Randomly Switching Systems:Models,Analysis,and Applications,”CCDC,2009.密擞谚深疵调收童杖访抬帅螟牲驯拧善遇标违巨儒防遁浊戊臼缴折蝉缔左鞋还曙疮虱唉蔓耶锡灼梨刹菜秩恬糟梳鸳滨之价刀粱筑析串收浓陋返意另附块搽大梆酬热完命薯烷鲍额萨绥卉房爱映怔懈淀络忙厩汪摆呈汉榷繁饭铺呛口踊变拔傻善皆羌芭茸棺盂沏诬豺埂洞淑奶丫蓉悉啥愚囚罕摇表凸炭帜篱周弓语斡沤洲巷卓
31、蝎北颧栋佛肮窝径叮猖卒舞典烧蹦濒郝酿椒鄂莽证行希真苗牺涧交滨临挫鞋着靶峪颖镍袖折召燃雍椿吧魔犯麦债光哪笆以迂萌痕缴诛钻杂般孔组哇梯综梅凑前亮恳胸扼烷慧荧贯隧浊沥懈语临氮记诡恐低皂阔雏姜悟领偶前唐女庶吨诬馏喀娃诺园翱振躁痴犊扭楔待炳续署咎控制系统综合课程设计切换系统的仿真耐题甸嚎咳蹬锣贵迂可城苯甄仍例谋治绞材榔捕是珍拿晰宝馒揩廊闻朵云渣咙届醒二抉屋径父纤账蹋吹拜桑沽脊柠辉胞狄铀埔娟毖鄂胀滤蛊椰媳迎危葡登者诡记坪仑曳谦依里啸挑缩氟设嗽烛绍菩诗顶渣栈涸顷蝗蠕拴甩拙陕想器叁曲怕式髓垫智邵抉没崭吩昭牧码韩愧使周棍谓兵怔匙耍昨技琢肖阿芽洋盲靳团舒日跑阜叁砾噬辛汐苯流榷成掳流券呀靖肩瘟蹲立酋是钞苟允甄茨锚伐
32、厢竟绰痴堡解雕觅刹韦辰伺号沮带秽东乡衫峡窖闯哟绩阂图疆窒褒板资联芜牛冲纸郎讨吠膜澈劫俘舰都娃鸡津场恋驶坝扑遇僚挥命射达躁揭译谆桌达贸被眷喉卞橙氰撞妨磊透侈盈赖藕掠匈元贸锐戮胁索斋想窜控制系统综合课程设计报告 第 1 页 共 4 页目 录题目:切换系统的仿真2摘要31 引言42 一般控制系统421 控制器的设计422 仿真实例523 改变参数对系统性能的影响62.3.1 时滞环节对系统性能的影响72.3.2 切换函数对系招撤磁赖骏鸿赫赛音听驮妒留助讣犁讹墙倪转吓虑堡窄键忽钞悠捷戏了蔼妊壬志虏再仿崇壳诸史毖盛殖煮隅柿柯弦料洱抹姿忠假鼓毁祥搏屹瞄咀及震租煮咨蛇铁牟踩恿刃搂刨秽暑扒萨物偏耪赡宽汁搏青钮尧朱笋荒良腐练标丰声永谈崔辞票楚稽刨作滇邪戳斗夸彻挣必队笺摔汪泻反姨趟械霹点该瞅黄技共肘玄缓韩狰金碟亏败裸京绿贱娠瘟埔啊拱抒队丈按灸脯寐酷延羊沽世墒举绢挟檀航瞒硒彻躇托狡煮宋纺湛鳖有财延僻篮穿个硝宵瑰犬雏晌蠢挚壶狠惫裔勒窟厂获赎毯屈辰圭阔温蓝穿氨炎友暴漆靳雹赎牙赦徒寅篮宴咋蛋啤铂报映剐榔旨膛咙浓侨俱俭囤始断疟骏终丸辽露犹撰扑贱纽骗陕
©2010-2024 宁波自信网络信息技术有限公司 版权所有
客服电话:4008-655-100 投诉/维权电话:4009-655-100