使用遗传算法优化的RBF神经网络优化算法代码（matlab）
今天咱们聊点硬核的——如何用遗传算法给RBF神经网络调参。这玩意儿就像给机器人做基因改造，让它的"脑子"自己进化出最优结构。直接上实战，手把手教你用Matlab整活！
先看RBF网络的骨架。核心就三部分：输入层、隐层径向基函数、输出层加权和。关键参数是隐层中心点、宽度和输出权重。传统做法用k-means找中心，但容易陷入局部最优，这时候就该进化算法出场了。
% 初始化遗传算法参数 pop_size = 50; % 种群规模 max_gen = 100; % 进化代数 pc = 0.8; % 交叉概率 pm = 0.05; % 变异概率 gene_len = 30; % 每个中心点参数长度这个配置里藏着玄机：种群规模太小容易早熟，太大算到地老天荒。我试过把pm设到0.1，结果种群直接变异成外星生物，完全跑偏。建议新手先保持默认值，后面再微调。
适应度函数是进化的指挥棒，这里用预测误差的倒数：
function fitness = calc_fitness(pop, X, y) fitness = zeros(size(pop,1),1); for i=1:size(pop,1) centers = reshape(pop(i,:), [], 3); % 解码染色体 net = newrb(X, y, 0.01, 0.1, centers(:,1:2)', centers(:,3)); pred = sim(net, X); fitness(i) = 1/(1 + mse(y - pred)); % 防止除零 end end注意那个reshape操作，这里把染色体编码成N个中心点（包含二维坐标和宽度）。newrb函数里的0.01是目标误差，0.1是隐层节点间隔，这两个参数控制着网络复杂度。
进化过程最带劲的是选择操作，这里用轮盘赌：
% 轮盘赌选择 cum_fit = cumsum(fitness)/sum(fitness); new_pop = zeros(size(pop)); for i=1:pop_size r = rand(); new_pop(i,:) = pop(find(cum_fit >= r, 1), :); end这个实现有个坑——cumsum别忘记归一化！有次我忘了除总和，结果选择的都是垃圾个体，整个种群直接崩盘。
使用遗传算法优化的RBF神经网络优化算法代码（matlab）
交叉变异才是重头戏，来看这个两点交叉：
% 两点交叉 for i=1:2:pop_size if rand() < pc pts = sort(randi(gene_len,1,2)); temp = pop(i, pts(1):pts(2)); pop(i, pts(1):pts(2)) = pop(i+1, pts(1):pts(2)); pop(i+1, pts(1):pts(2)) = temp; end end随机选两个切点交换基因片段，比单点交叉探索空间更大。不过要注意索引越界问题，Matlab的randi函数是闭区间，别手抖写成0了。
最后测试效果时，拿sinc函数开刀：
% 生成训练数据 X = -5:0.5:5; y = sinc(X) + 0.1*randn(size(X)); % 传统RBF std_net = newrb(X, y, 0.01, 0.5); % 遗传算法优化RBF ga_net = train_ga_rbf(X, y); % 可视化对比 test_X = -5:0.1:5; plot(X, y, 'ro', test_X, sim(std_net, test_X), 'b--', test_X, sim(ga_net, test_X), 'g-');跑出来的结果特别有意思：传统方法（蓝色虚线）在边缘区域明显摆烂，遗传算法优化的绿色曲线稳稳贴合真实数据。这说明进化后的网络不仅记忆训练数据，还学会了泛化规律。
调参时发现个反直觉的现象——有时候增加变异概率反而收敛更快。因为RBF参数空间存在多个局部最优，适当"捣乱"能跳出陷阱。不过别过火，超过0.1的话，适应度曲线就跟心电图似的狂抖。
这种算法融合的思路其实可以玩出花：把进化策略改成粒子群，或者用神经网络指导遗传算法的参数自适应调整。不过那就是另一个故事了，下回咱们再接着唠。
