KevinTq的专栏

%这是主程序,命名为Main1206.mclc;clear all;close all;%本程序研究的对象是瞬时点源污染扩散问题subplot(1,2,1);daspect([1 1 1]);xmin=-500;dx=10;xmax=500;ymin=-100;dy=5;ymax=100;zmin=-10;dz=1;zmax=10;[cx,cy,ct]=meshgrid(xmin:dx

clear syms a bc=[a,b]';A=[89677,99215,109655,120333,135823,159878,182321,209407,246619,300670];B=cumsum(A); %原始数据累加n=length(B);for i=1:(n-1) C(i)=(B(i)+B(i+1))/2; %生成累加矩阵

clearclca = 0.95k = [5;10;13;4;3;11;13;10;8;16;7;4];k = -k; % 模拟退火算法是求解最小值，故取负数d = [2;5;18;3;2;5;10;4;11;7;14;6];restriction = 46;num = 12;sol_new = ones(1,num); % 生成初始解E_current = in

%3.1_cwstd.m,用总和标准化法标准化矩阵 function std=cwstd(vector) cwsum=sum(vector,1); %对列求和 [a,b]=size(vector); %矩阵大小,a 为行数,b 为列数 for i=1:a for j=1:b std(i,j)= vector(i

function [xm,fv]=PSO_Kades(fitness,N,c1,c2,w,M,D)%待优化的函数：fitness%粒子数目：N%学习因子：c1%学习因子：c2%惯性权重：W%最大迭代次数：M%取得最小值时的自变量值：xm%目标函数的最小值：fv%初始化初始位置format long;for i=1:N for j=1:D x(i,j)=

%这是一个用GA算法来求函数极值的例子；clcclear allT=100;optionsOrigin=gaoptimset('Generation',T/2);[x,fval,reason,output,finnal_pop]=ga(@ch14_2f,2,optionsOrigin);%options1=gaoptimset('Generations',T/2,'Init

x1=[5.5 2.5 8 3 3 2.9 8 9 4 6.5 5.5 5 6 5 3.5 8 6 4 7.5 7]';x2=[31 55 67 50 38 71 30 56 42 73 60 44 50 39 55 70 40 50 62 59]';x3=[10 8 12 7 8 12 12 5 8 5 11 12 6 10 10 6 11 11 9 9]';x4=[8 6 9 16 15