【matlab】基于粒子群PSO算法实现限定自变量取值范围下的函数寻优

任务描述:对bp训练得到的模型进行寻优,得出多个自变量在不同的取值范围下的函数的最小值或者最大值

PSO.m
此函数可用于寻找最小值

clc
clear
close all
load('net.mat')
data = xlsread('data.xlsx');
P_train = data(1:16,1:3)';
T_train = data(1:16,4:5)';
[inputn, ps_input] = mapminmax(P_train,0,1);%归一化
[outputn, ps_target] = mapminmax(T_train,0,1);%归一化

%速度更新
c1 = 2;
c2 = 1.5;
%
maxgen = 50; %最大更新代数
sizepop = 20; %种群大小
%
Vmax = 1;  %规定速度的大小范围
Vmin = -1;
%种群初始化
popmax1 = 12;  %规定自变量x 的取值范围
popmin1 = 2;
for i = 1:sizepop
    pop(i,1) = (popmax1-popmin1)*rand(1,1)+popmin1;
    V(i,1) = 2*Vmax*rand(1,1)+Vmin;
end
popmax2 = 4;  %规定自变量x 的取值范围
popmin2 = 1;
for i = 1:sizepop
    pop(i,2) = (popmax2-popmin2)*rand(1,1)+popmin2;
    V(i,2) = 2*Vmax*rand(1,1)+Vmin;
end
popmax3 = 40;  %规定自变量x 的取值范围
popmin3 = 10;
for i = 1:sizepop
    pop(i,3) = (popmax3-popmin3)*rand(1,1)+popmin3;
    V(i,3) = 2*Vmax*rand(1,1)+Vmin;
end
for i = 1:sizepop
    fitness(i) = fun(pop(i,:),ps_input,ps_target,net);
end

[best_fitness,best_index] = min(fitness);
zbest = pop(best_index,:);%群体最优个体
gbest = pop;%个体最优个体
fitness_gbest = fitness;%个体最优个体函数值
fitness_zbest = best_fitness;%群体最优函数值
%
w_start = 0.9;
w_end = 0.4;
for i = 1:maxgen
    w = w_start - (w_start-w_end)*(i/maxgen)^2;
    for j = 1:sizepop
        %速度更新
        V(j,:) = V(j,:) + c1*rand*(gbest(j,:)-pop(j,:)) + ...
                          c2*rand*(zbest-pop(j,:));
        V(j,find(V(j,:)>Vmax)) = Vmax;
        V(j,find(V(j,:)<Vmin)) = Vmin;
        pop(j,:) = pop(j,:) + w*V(j,:);
    end
    for j=1:3
        %粒子更新
        if j==1
            pop(find(pop(:,j)>popmax1),j) = popmax1;
            pop(find(pop(:,j)<popmin1),j) = popmin1;            
        end
        if j==2
            pop(find(pop(:,j)>popmax2),j) = popmax2;
            pop(find(pop(:,j)<popmin2),j) = popmin2;            
        end
        if j==3
            pop(find(pop(:,j)>popmax3),j) = popmax3;
            pop(find(pop(:,j)<popmin3),j) = popmin3;            
        end
        
    end
    for j = 1:sizepop
        fitness(j) = fun(pop(j,:),ps_input,ps_target,net);
    end
    %个体与种群极值更新
    for j = 1:sizepop
        %个体
        if fitness(j) < fitness_gbest(j)
            gbest(j,:) = pop(j,:);
            fitness_gbest(j) = fitness(j);
        end
        %种群
        if fitness(j) < fitness_zbest
            zbest = pop(j,:);
            fitness_zbest = fitness(j);
        end
    end
    xx(i,:) = zbest;
    result(i) = fitness_zbest;
end

%% 结果显示
plot(result);
title(['最优函数值:',num2str(fitness_zbest)]);
zbest

fun.m
bp模型的输出是2维的,所以这里分别进行处理,对每个输出都进行了寻优

function y = fun(x,ps_input,ps_target,net)
input_n = mapminmax('apply', x', ps_input);
an=sim(net,input_n);
BPoutput = mapminmax('reverse', an, ps_target); %反归一化
% y=BPoutput(1,1);%求MRR时用这个
y=BPoutput(2,1);%求Ra时用这个
end


  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
### 回答1: 粒子群优化算法(Particle Swarm Optimization, PSO)是一种进化计算算法,常用于解决优化问题。PSO算法模拟了鸟群或鱼群等动物群体协作求解问题的行为,通过不断调整粒子的位置来搜索最优解。 要实现基于PSO算法的圆形障碍物避障问题,可以按以下步骤进行: 步骤1:初始化参数 有关粒子群优化算法的参数有很多,其中包括粒子数、搜索空间的范围、迭代次数、惯性权重等。根据具体场景和需求确定这些参数,并进行初始化。 步骤2:生成初始粒子群 根据搜索空间的范围,在该范围内随机生成一定数量的粒子,并将其初始位置和速度设定为随机值。 步骤3:计算适应度函数 根据具体的圆形障碍物避障问题,定义适应度函数。适应度函数的设计应考虑到避免与圆形障碍物碰撞并尽可能接近目标点。 步骤4:更新速度和位置 根据粒子的当前速度和位置以及全局最优和个体最优解,更新粒子的速度和位置。更新公式一般包含三个部分:惯性项、个体经验项和群体经验项,具体的公式可以根据实际情况进行设计。 步骤5:检查是否满足终止条件 判断是否满足终止条件,例如迭代次数达到预设值或达到了预期的适应度值。 步骤6:返回最优解 迭代完成后,返回全局最优解,即找到的避开圆形障碍物并接近目标点的粒子位置。 以上是基于粒子群优化算法实现圆形障碍物避障的一般步骤,具体实现时需要根据问题的具体要求进行参数调整和代码编写。请注意,由于字数限制,无法提供完整的MATLAB代码。 ### 回答2: 粒子群优化算法PSO)是一种仿生智能算法,由于其简单性和高效性在优化问题上得到广泛应用。要基于PSO算法实现圆形障碍物避障,可以按照以下步骤进行: 首先,在Matlab中定义问题,包括定义圆的位置、大小和数量,以及设定目标函数和目标位置。 然后,定义粒子群的初始状态,包括设定粒子的初始位置和速度。可以将每个粒子的位置表示为二维坐标(x,y),并将速度表示为每个维度上的变化量。 接下来,对于每个粒子,根据PSO算法公式更新其速度和位置。公式中包含了个体最优和全局最优位置的影响因素,可以通过适应度函数来评估粒子的表现,并记录最优位置。 在更新位置之后,需要检查粒子是否与障碍物相交。如果粒子与障碍物相交,则根据一定策略调整粒子的位置,使其避免障碍物。 最后,根据设定的停止准则,判断是否继续迭代。如果满足停止准则,则结束迭代,输出最优解;否则,继续进行下一轮迭代。 通过以上步骤,就可以实现基于粒子群优化算法的圆形障碍物避障。具体实现过程中需要根据实际情况进行细节调整,如权重和学习因子的选择,障碍物避障策略的确定等。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值