[问答]

萌新求助，求单纯型法MATLAB程序的编程

931 matlab

问答对人有帮助，内容完整，我也想知道答案 0 萌新求助，求单纯型法matlab程序的编程 0
2021-11-22 07:43:08　　评论淘帖0 邀请回答您可以邀请以下用户，快速回答问题 × heks 该类别下有 50 个回答。邀请回答 hgimtk 该类别下有 44 个回答。邀请回答新星之火12138 该类别下有 42 个回答。邀请回答 wang21cj 该类别下有 41 个回答。邀请回答 chm5 该类别下有 40 个回答。邀请回答 hjfjsdgfjdsf 该类别下有 37 个回答。邀请回答 werywer 该类别下有 35 个回答。邀请回答 h1654155957.9520 该类别下有 35 个回答。邀请回答 fdjslkjd 该类别下有 35 个回答。邀请回答 uvysdfydad 该类别下有 35 个回答。邀请回答维生素B2 该类别下有 34 个回答。邀请回答凤毛麟角该类别下有 34 个回答。邀请回答 jenny042 该类别下有 34 个回答。邀请回答 dfzvzs 该类别下有 34 个回答。邀请回答刘洋该类别下有 33 个回答。邀请回答 meihuacg 该类别下有 33 个回答。邀请回答 hfgdzc 该类别下有 33 个回答。邀请回答 ggfvxv 该类别下有 33 个回答。邀请回答河南顺之航该类别下有 32 个回答。邀请回答 wenminglang 该类别下有 32 个回答。邀请回答举报回头太晚相关推荐 • 萌新求助，求50Hz陷波器的Matlab程序 1148 • 怎样去编写一种有约束复合型法MATLAB主程序呢 1494 • 萌新求助，求大佬详细介绍一下stm32S型加减速算法 610 • 萌新求助，求一个matlab例程 552 • 萌新求助，求云模型及发生器matlab代码 803 • 什么是单纯形法？单纯形法有何作用 1602 • 萌新求助，求使用cubeMX生成IO的配置程序 652 • 萌新求助，求大佬分享一级倒立摆MATLAB仿真程序 1302 • 萌新求助，求基于51单片机心形流水灯C语言源程序 1687 • 萌新求助，求继电器控制照明设备的程序代码 1843 1个回答

答案对人有帮助，有参考价值 0 写这篇博文是因为一些相关的程序编程用到了个人改动后的函数作为子函数运行，在个人其他博文介绍时会涉及到，所以附程序如下，分享出来方便大家交流。 clear all;clc%% 初始参数设定，求f最小值syms x1 x2 xf=x1^2+2x2^2-4x1-2x1x2;% f 函数符号表达式x0=[0 0];% 给定初始点h=4; % 初始单纯形构造因子alpha=1; % alpha 反射因子gama=2; % gama 扩张因子beta=0.5;% beta 收缩因子eps=0.5;% eps 精度值%% 构造初始单纯形x0=x0';X(:,1)=x0;Num_var_f=length(symvar(f));%二维函数for i = 1:Num_var_f %二维函数，需要三个顶点坐标 x = zeros(Num_var_f,1);%预分配内存 x(i) = x0(i) + h; X(:,i+1)=x;% 将单纯形存储endn = 0; %迭代循环次数m = 1; %精度差值k = 0; %数据记录g = 0; %操作记录 %% 计算函数值 for i = 1:Num_var_f+1 F_val(i)=double(subs(f,symvar(f),X(:,i)')); end [F_sort,F_index]=sort(F_val);%从小到大排列函数 % 最好点 f_best_0=F_sort(1); X_best_0=X(:,F_index(1)); % 最差点 f_bad_n=F_sort(end); X_bad_n=X(:,F_index(end)); % 次差点 f_nextbad_1n=F_sort(end-1); X_nextbad_1n=X(:,F_index(end-1)); % 计算形心 Xc_n1=1/Num_var_fsum(X(:,F_index(1:end-1)),2); f_Xc_n1=double(subs(f,symvar(f),Xc_n1')); %% 循环迭代语句while m > eps %% 反射 flag = 0; X_reflect_n2=Xc_n1+alpha(Xc_n1-X_bad_n); f_reflect_n2=double(subs(f,symvar(f),X_reflect_n2')); %% 反射成功 if f_reflect_n2= f_nextbad X_reflect_shrink_n4=Xc_n1+beta(X_reflect_n2-Xc_n1);% 收缩点 f_reflect_shrink_n4=double(subs(f,symvar(f),X_reflect_shrink_n4')); X(:,F_index(end)) = X_reflect_shrink_n4; g=2;%收缩点 end end %% 计算迭代完后函数序列，并计算精度 for i = 1:Num_var_f+1 F_val(i)=double(subs(f,symvar(f),X(:,i)')); end [F_sort,F_index]=sort(F_val);%从小到大排列函数 % 最好点 f_best_0=F_sort(1); X_best_0=X(:,F_index(1)); % 最差点 f_bad_n=F_sort(end); X_bad_n=X(:,F_index(end)); % 次差点 f_nextbad_1n=F_sort(end-1); X_nextbad_1n=X(:,F_index(end-1)); % 计算形心 Xc_n1=1/Num_var_fsum(X(:,F_index(1:end-1)),2); f_Xc_n1=double(subs(f,symvar(f),Xc_n1')); % 计算精度与迭代次数 m=abs(double(subs(f,symvar(f),X(:,F_index(end))'))-f_Xc_n1); n=n+1; % 运算数据记录 k=k+1; A(k)=g; B(:,k+1)=X(:,F_index(end))'; M(k)=m;end %% 结果输出minx=X_best_0;A=A';B=B';M=M';ABMnminxm

2021-11-22 09:27:27 评论举报高彬