图表处理+试验+遗传算法(第3,4讲).ppt

图表处理方法 产品资料输入 变压器设计 电磁计算（铁芯和绕组数据等） 结构设计 方案输出 function y=pt(x) bb= [1.65 1.66 1.67 1.68 1.69 1.7 1.71 1.72 1.73 ]; ptb=[1.209 1.232 1.256 1.280 1.305 1.330 1.356 1.382 1.408 ]; for kk= 2:9 if xbb(kk) bb1=bb(kk-1); bb2=bb(kk); pt1=ptb(kk-1); pt2=ptb(kk); pt=pt1+(x-bb1)*(pt2-pt1)/(bb2-bb1); break; end end bb=[1.65 1.66 1.67 1.68 1.69 1.7 1.71 1.72 1.73 ]; ptb=[1.209 1.232 1.256 1.280 1.305 1.330 1.356 1.382 1.408]; a=polyfit(bb,ptb,3); bb=[1.65 1.66 1.67 1.68 1.69 1.7 1.71 1.72 1.73] ptb=[1.209 1.232 1.256 1.280 1.305 1.330 1.356 1.382 1.408]; x=[1.65:0.001:1.73]; a =[ -11.7845 61.9986 -106.0747 60.3786]; y=a(1)*x.^3+a(2)*x.^2++a(3)*x+a(4); plot(bb,ptb,o,x,y,-) 实验一、约束遍数法优化计算 一、实验目的 通过编程了解和掌握约束遍数法优化计算的基本原理和方法。 三、实验内容 铁心直径φ100mm，由三级硅钢片叠制而成，求截面最大时的各级片宽和叠厚尺寸（如有多种方案，应全部保留）。利用约束遍数法编程解决上述工程问题。 要求：铁心片宽为5的倍数， 末级片宽≥30mm。 实验一、约束遍数法优化计算 for A1=30:5:95 B1 D1计算 for A2=30:5:95 if A2A1 B2 D2计算 for A3=30:5:95 if A3A2 B3 D3计算 s= if ssmax smax=s 实验一、约束遍数法优化计算 clc;clear all; d=100; smax=0; for a1= 30:5: 95 b1=round((d^2-a1^2)^.5); d1=b1/2; for a2=30:5:95 if a2a1 b2=round((d^2-a2^2)^.5); d2=b2/2-d1; for a3=30:5:95 if a3a2 b3=round((d^2-a3^2)^.5); d3=(b3-b2)/2; s=2*(a1*d1+a2*d2+a3*d3); if ssmax smax=s; jg=[a1,b1,a2,b2,a3,b3,smax]; end end end end end end jg 遗传算法原理 二进制编码 二进制 0 1 0 1 实码 基数+二进制值*级差 90，95，100，105 编码 0000 0001 0010 0011 染色体与基因 单个二进制码 各个变量的编码连接 初始种群 随机 产生 遗传算法原理 评价 适应度=1/目标函数（最小值） 选择 计算个体适应度 总适应度计算 总选择概率=1 个体选择概率段 产生一个随机数 看该随机数落在哪个概率