• 基于PSO粒子群优化的汽车刹车稳定性数据matlab仿真与分析

    欢迎订阅《FPGA学习入门100例教程》、《MATLAB学习入门100例教程》

    目录

    一、理论基础

    二、核心程序

    三、测试结果

    一、理论基础

    刹车类别A

    温度1

    测试1

    测试2

    .......

    测试n

    温度3

    测试1

    测试2

    .......

    测试n

    .......

    .......

    温度N

    测试1

    测试2

    .......

    测试n

            每次测试的温度间隔，根据实际操作的可行性，来选择。这个无所谓，从里面看到温度从20多度到60度，如果条件允许，可以再将温度扩展一下。即测试的温度范围再大点。然后，对于同一个温度，测试次数的确定，如果你每次测试的数值差别较大，那么测试次数最好多点，如果每次测试的差别不大，那么测试20次就够了。

          粒子群优化(PSO)算法是通过模拟鸟群觅食过程中的迁徙和群聚行为而提出的一种基于群体智能的全局随机搜索算法。PSO是将群体(swarm)中的个体看作是在D维搜索空间中没有质量和体积的粒子(particle)，每个粒子以一定的速度在解空间运动，并向自身历史最佳位置pbest和邻域历史最佳位置聚集，实现对候选解的进化。

           PSO从这种模型中得到启示并用于解决优化问题。PSO中，每个优化问题的潜在解都是搜索空间中的一只鸟，称之为粒子。所有的粒子都有一个由被优化的函数决定的适值( fitness value)，每个粒子还有一个速度决定它们飞翔的方向和距离。然后粒子们就追随当前的最优粒子在解空间中搜索。

           PSO初始化为一群随机粒子，然后通过迭代找到最优解。在每一次迭代中，粒子通过跟踪两个极值来更新自己；第一个就是粒子本身所找到的最优解，这个解称为个体极值；另一个极值是整个种群目前找到的最优解，这个极值是全局极值。另外也可以不用整个种群而只是用其中一部分作为粒子的邻居，那么在所有邻居中的极值就是局部极值。

    二、核心程序

    1. ................................................................
    2. parameter;
    3. ishow = (get(handles.checkbox1,'Value'));
    4. if choice == 1
    5. L2    = func_anti_noise(L1,T1);
    6. else
    7. L2    = func_anti_noise(L1,T2);    
    8. end
    9. hold off;
    10. axes(handles.axes4);
    11. color_table;
    12. for i = 1:len1
    13.     str = colors{i};
    14.     plot(L2(i,:),str);
    15.     hold on;
    16. end
    17. if ishow == 1
    18.     figure;
    19.     color_table;
    20.     for i = 1:len1
    21.         str = colors{i};
    22.         plot(L2(i,:),str);
    23.         hold on;
    24.     end
    25.     title('Leaks anti noise');
    26. end 
    27.  
    28.  
    29. % --- Executes on button press in pushbutton4.
    30. function pushbutton4_Callback(hObject, eventdata, handles)
    31. % hObject    handle to pushbutton4 (see GCBO)
    32. % eventdata  reserved - to be defined in a future version of MATLAB
    33. % handles    structure with handles and user data (see GUIDATA)
    34. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%多项式拟合%%%%%%%%%%%
    35. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%多项式拟合%%%%%%%%%%%
    36. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%多项式拟合%%%%%%%%%%%
    37. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%多项式拟合%%%%%%%%%%%
    38. parameter;
    39. ishow            =        (get(handles.checkbox1,'Value'));
    40. KK               = str2num(get(handles.edit1,'String')); 
    41.  
    42. if choice == 1
    43. [y1,x1,y0,x0,P1fit] =  func_fit_means(T1,L2,KK);
    44. else
    45. [y1,x1,y0,x0,P1fit] =  func_fit_means(T2,L2,KK);    
    46. end
    48. hold off;
    49. axes(handles.axes5);
    50. plot(x0,y0,'o');
    51. hold on;
    52. plot(x1,y1,'r','LineWidth',1);
    53. hold on;
    54.  
    55. if ishow == 1
    56.     figure;
    57.     plot(x0,y0,'o');
    58.     hold on;
    59.     plot(x1,y1,'r','LineWidth',1);
    60.     hold on;
    61.     title('多项式拟合曲线');
    62.     
    63. end 
    65. %计算拟合误差
    66. err1 = func_error(y0,x0,y1,x1);
    67. set(handles.edit2,'string',num2str(err1)); 
    68.  
    69.  
    70.  
    71.  
    72. % --- Executes on button press in pushbutton5.
    73. function pushbutton5_Callback(hObject, eventdata, handles)
    74. % hObject    handle to pushbutton5 (see GCBO)
    75. % eventdata  reserved - to be defined in a future version of MATLAB
    76. % handles    structure with handles and user data (see GUIDATA)
    77. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%最大值多项式拟合%%%%%%%%%%%
    78. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%最大值多项式拟合%%%%%%%%%%%
    79. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%最大值多项式拟合%%%%%%%%%%%
    80. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%最大值多项式拟合%%%%%%%%%%%
    81. parameter;
    82. ishow            =        (get(handles.checkbox1,'Value'));
    83. KK               = str2num(get(handles.edit1,'String')); 
    84. if choice == 1
    85. [y1max,x1max,y0max,x0max,P1max] =  func_fit_max(T1,L2,KK);
    86. else
    87. [y1max,x1max,y0max,x0max,P1max] =  func_fit_max(T2,L2,KK);    
    88. end
    89. hold off;
    90. axes(handles.axes5);
    91. plot(x0,y0,'o');
    92. hold on;
    93. plot(x1,y1,'r','LineWidth',1);
    94. hold on;
    95. plot(x1max,y1max,'g','LineWidth',1);
    96. hold on; 
    97. if ishow == 1
    98.     figure;
    99.     plot(x0,y0,'o');
    100.     hold on;
    101.     plot(x1,y1,'r','LineWidth',1);
    102.     hold on;
    103.     plot(x1max,y1max,'g','LineWidth',1);
    104.     hold on; 
    105.     title('多项式拟合曲线');
    106. end 
    113. function pushbutton6_Callback(hObject, eventdata, handles)
    114. % hObject    handle to pushbutton6 (see GCBO)
    115. % eventdata  reserved - to be defined in a future version of MATLAB
    116. % handles    structure with handles and user data (see GUIDATA)
    117. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%最小值多项式拟合%%%%%%%%%%%
    118. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%最小值多项式拟合%%%%%%%%%%%
    119. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%最小值多项式拟合%%%%%%%%%%%
    120. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%最小值多项式拟合%%%%%%%%%%%
    121. parameter;
    122. ishow            =        (get(handles.checkbox1,'Value'));
    123. KK               = str2num(get(handles.edit1,'String')); 
    124. if choice == 1
    125. [y1min,x1min,y0min,x0min,P1min] =  func_fit_min(T1,L2,KK);
    126. else
    127. [y1min,x1min,y0min,x0min,P1min] =  func_fit_min(T2,L2,KK);    
    128. end
    129. hold off;
    130. axes(handles.axes5);
    131. plot(x0,y0,'o');
    132. hold on;
    133. plot(x1,y1,'r','LineWidth',1);
    134. hold on;
    135. plot(x1max,y1max,'g','LineWidth',1);
    136. hold on; 
    137. plot(x1min,y1min,'m','LineWidth',1);
    138. hold on; 
    139. if ishow == 1
    140.     figure;
    141.     plot(x0,y0,'o');
    142.     hold on;
    143.     plot(x1,y1,'r','LineWidth',1);
    144.     hold on;
    145.     plot(x1max,y1max,'g','LineWidth',1);
    146.     hold on; 
    147.     plot(x1min,y1min,'m','LineWidth',1);
    148.     hold on; 
    149.     title('多项式拟合曲线');
    150.     
    151. end 
    154. % --- Executes on button press in pushbutton7.
    155. function pushbutton7_Callback(hObject, eventdata, handles)
    156. % hObject    handle to pushbutton7 (see GCBO)
    157. % eventdata  reserved - to be defined in a future version of MATLAB
    158. % handles    structure with handles and user data (see GUIDATA)
    159. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%PSO优化拟合%%%%%%%%%%%
    160. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%PSO优化拟合%%%%%%%%%%%
    161. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%PSO优化拟合%%%%%%%%%%%
    162. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%PSO优化拟合%%%%%%%%%%%
    163. parameter;
    164. ishow            =        (get(handles.checkbox1,'Value'));
    165. KK               = str2num(get(handles.edit1,'String')); 
    168. %进行PSO优化
    169. [V_score2,PP2] = func_pso(P1fit,x1,y0,x0);
    170. ypso           = polyval(PP2,x1);
    171.  
    172. hold off;
    173. axes(handles.axes6);
    174. plot(V_score2,'Linewidth',2);
    175.  
    176.  
    177.  
    178. hold off;
    179. axes(handles.axes7);
    180. plot(x0,y0,'o');
    181. hold on;
    182. plot(x1,y1,'r','LineWidth',1);
    183. hold on;
    184. plot(x1max,y1max,'g','LineWidth',1);
    185. hold on; 
    186. plot(x1min,y1min,'m','LineWidth',1);
    187. hold on; 
    188. plot(x1,ypso,'k--','LineWidth',2);
    189. hold on; 
    190.  
    191. if ishow == 1
    192.     figure;
    193.     subplot(121);
    194.     plot(V_score2,'Linewidth',2);
    195.     xlabel('迭代次数');
    196.     ylabel('拟合误差优化过程');
    197.     subplot(122);
    198.     plot(x0,y0,'o');
    199.     hold on;
    200.     plot(x1,y1,'r','LineWidth',1);
    201.     hold on;
    202.     plot(x1max,y1max,'g','LineWidth',1);
    203.     hold on; 
    204.     plot(x1min,y1min,'m','LineWidth',1);
    205.     hold on; 
    206.     plot(x1,ypso,'k--','LineWidth',2);
    207.     hold on; 
    208.     title('多项式拟合曲线');
    209.     legend('实际测量值','均值拟合曲线','最大测量值拟合曲线','最小测量值拟合曲线','优化后拟合曲线');
    210. end 
    211.  
    212. err2 = func_error(y0,x0,ypso,x1);
    213. set(handles.edit3,'string',num2str(err2)); 
    214.  
    215.  
    216. % --- Executes on button press in pushbutton8.
    217. function pushbutton8_Callback(hObject, eventdata, handles)
    218. % hObject    handle to pushbutton8 (see GCBO)
    219. % eventdata  reserved - to be defined in a future version of MATLAB
    220. % handles    structure with handles and user data (see GUIDATA)
    221. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%统计量计算%%%%%%%%%%%
    222. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%统计量计算%%%%%%%%%%%
    223. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%统计量计算%%%%%%%%%%%
    224. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%统计量计算%%%%%%%%%%%
    225. parameter;
    227. %计算均值和方差，将这两个计算结果存放到四个变量的后两位
    228. for i = 1:size(T1,1)
    229.     %计算均值
    230.     E1(i) = mean(T1(i,:));
    231.     E2(i) = mean(T2(i,:));
    232.     E3(i) = mean(P1(i,:));
    233.     E4(i) = mean(L2(i,:));
    234.     %计算标准差
    235.     D1(i) = std(T1(i,:));
    236.     D2(i) = std(T2(i,:));
    237.     D3(i) = std(P1(i,:));
    238.     D4(i) = std(L2(i,:));    
    239.     %计算方差
    240.     V1(i) = var(T1(i,:));
    241.     V2(i) = var(T2(i,:));
    242.     V3(i) = var(P1(i,:));
    243.     V4(i) = var(L2(i,:)); 
    244. end
    245.  
    246. figure(10);
    247. plot(E1,E4,'b-o');
    248. hold on;
    249. plot(E1,D4,'r-o');
    250. hold on;
    251. plot(E1,V4,'k-o');
    252. legend('均值','标准差','方差');
    253.  
    254.  
    255. fid = fopen('statistics.txt','wt');
    256. fprintf(fid,'均值       标准差     方差\n');
    257. for i = 1:length(E4)
    258.     fprintf(fid,'%f   ',E4(i));
    259.     fprintf(fid,'%f   ',D4(i));
    260.     fprintf(fid,'%f   ',V4(i));
    261.     fprintf(fid,'\n');
    262. end
    263. fclose(fid);
    264.  
    265.  
    266. % --- Executes on button press in pushbutton9.
    267. function pushbutton9_Callback(hObject, eventdata, handles)
    268. % hObject    handle to pushbutton9 (see GCBO)
    269. % eventdata  reserved - to be defined in a future version of MATLAB
    270. % handles    structure with handles and user data (see GUIDATA)
    271. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%拟合统计分析%%%%%%%%%%%
    272. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%拟合统计分析%%%%%%%%%%%
    273. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%拟合统计分析%%%%%%%%%%%
    274. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%拟合统计分析%%%%%%%%%%%
    275. parameter;
    276. ishow            =        (get(handles.checkbox1,'Value'));
    278. hold off;
    279. axes(handles.axes8);
    280. plot(x0,y0,'o');
    281. hold on;
    282. plot(x1,ypso,'k--','LineWidth',2);
    283. hold on; 
    284. %找到最接近0的温度
    285. Best_tmp = 0;
    286. Ind      = 0;
    287. for i = 2:length(ypso)-1
    288.     if ypso(i-1) < 0 & ypso(i)>=0   
    289.        Best_tmp = x1(i);
    290.        Ind      = i;
    291.     end   
    292. end
    293. plot(x1(Ind),-4:0.01:8,'r','LineWidth',1);
    294. hold on; 
    295. plot(x1,ypso(Ind),'r','LineWidth',1);
    296.  
    297.  
    298. %根据统计量分析
    299. %找到均值最小的值
    300. Tmpes          = abs(E4);
    301. [V,I]          = sort(Tmpes);
    302. Best_Team_mean = I(1);
    303.  
    304. %找到均值最小的值
    305. Tmpes          = abs(D4);
    306. [V,I]          = sort(Tmpes);
    307. Best_Team_std  = I(1);
    308.  
    309. if  ishow == 1
    310.     figure(11);
    311.     subplot(121);
    312.     plot(x0,y0,'o');
    313.     hold on;
    314.     plot(x1,ypso,'k--','LineWidth',2);
    315.     hold on; 
    316.     plot(x1(Ind),-4:0.01:8,'r','LineWidth',1);
    317.     hold on; 
    318.     plot(x1,ypso(Ind),'r','LineWidth',1); 
    319.     axis square;
    320.     
    321.     subplot(122)
    322.     plot(E1(Best_Team_mean),E4(Best_Team_mean),'--rs','MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',10);
    323.     hold on;
    324.     plot(E1,E4,'b-o'); 
    325.     hold on;
    326.     plot(E1(Best_Team_std),D4(Best_Team_std),'--rs','MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',10);
    327.     hold on;
    328.     plot(E1,D4,'k-o');       
    329.     hold on;
    330.     axis square;
    331.     legend('均值','标准差');
    332. end
    333.  
    334. if Best_Team_mean == Best_Team_std 
    335.    msgbox(['最稳定温度为：',num2str(Best_tmp),' 对应的测试组为：',num2str(Best_Team_mean)]);
    336. else
    337.    msgbox(['最稳定温度为：',num2str(Best_tmp),' 对应的测试组为：',num2str(Best_Team_mean),' 和 ',num2str(Best_Team_std)]); 
    338. end

    三、测试结果

     1：显示EXCEL中的数据，将EXCEL导入MATLAB，然后显示出其中的数据：

     

    2：3：显示两个温度和泄漏之间对应关系图，这个是根据样本直接显示的

    4：显示多个LEAK的样本值：

    5：去掉LEAK样本中的一些异常采样点：

     

    6：多项式拟合效果显示：

    7：粒子群优化后的误差收敛曲线

    8：粒子群优化后的多项式拟合效果：

    9：最后的稳定性分析结果：

    左图是根据优化拟合结果，得到实际中最有可能为稳定区域的温度值，如上面的对话框为31.43.

    但是这个在实际中可能没法测试到。

    所以上图中右图就是表示的这个情况，通过均值和方差最小化，得到测试组13和14两个测试温度区域为最佳的区域。

    通过分析，

    在GUI上的数据显示区域，会自动对最佳的数据测试组显示the best的字样。

    A28-13

  • 相关阅读:
    【Rust 笔记】14-集合（上）
    ubuntu非源码安装 openCV
    Monaco Editor教程（十六）：缩略图minimap的配置详解
    【C语言】数据的存储_数据类型：浮点型存储
    HTTP协议及Servlet类笔记
    软件加密系统Themida应用程序保护指南（九）：通过命令行进行保护
    Spring Boot ： ORM 框架 JPA 与连接池 Hikari
    职场人一起进阶吧（内附技术人进阶路径、Java开源项目完整推荐等）
    【文件系统】Linux文件系统的基本存储机制
    Android修行手册 - Toolbar搜索联动以及各种监听和小技巧
  • 原文地址：https://blog.csdn.net/ccsss22/article/details/127895045