wordpress网站加密,泰州建筑工程网,衡阳做网站ss0734,php如何做局域网的网站开篇语前阵子做现代设计方法的时候#xff0c;发现网上很是缺乏这种作业形式的简易算法实现#xff0c;所以特地来简书写一篇。有两份#xff0c;一份是我的(说来惭愧#xff0c;我的大部分都是在网上找的代码#xff0c;然后在自己的电脑上跑一次#xff0c;跑出来了就行…开篇语前阵子做现代设计方法的时候发现网上很是缺乏这种作业形式的简易算法实现所以特地来简书写一篇。有两份一份是我的(说来惭愧我的大部分都是在网上找的代码然后在自己的电脑上跑一次跑出来了就行了的。而且我的电脑跑到2-11就扑街了。暂时还没有拿去修所以其实我的代码都是在网上整理之后调整了下就上的准确性不敢保证)另一份是我的室友的据他说是全部经过调试的虽然还是有不少的错误但是应该比我的要好一点。我的电脑另外我的电脑因为崩了所以我的代码无法验证结果为了交作业只能把我的室友的那些运行结果直接上一遍了。估计有点出入见谅重要的是代码正文2-102-10我的黄金分割法f(x) x20/xgolden(f,2,10,0.01)function[xmin]golden(f,a,b,e)k0;a1b-0.618*(b-a); %插入点的值a2a0.618*(b-a);while b-ae %循环条件y1subs(f,a1);y2subs(f,a2);if y1y2 %比较插入点的函数值的大小aa1; %进行换名a1a2;y1y2;a2a0.618*(b-a);elseba2;a2a1;y2y1;a1b-0.618*(b-a);endkk1;end %迭代到满足条件为止就停止迭代xmin(ab)/2;fminsubs(f,xmin) %输出函数的最优值fprintf(k\n); %输出迭代次数disp(k);f (x)x20/x [x,y]golden(f,2,10,0.01)x 4.4683y 8.9443二次插值法f(x) x20/x;a2;b10;eps1.0e-6; % 计算精度x1a;x3b;x2(ab)/2;f1 f(x1);f2 f(x2);f3 f(x3);while 1C1(x2^2-x3^2)*f1(x3^2-x1^2)*f2(x1^2-x2^2)*f3;C2(x2-x3)*f1(x3-x1)*f2(x1-x2)*f3;xp0.5*C1/C2;fp f(xp);if abs(x2-xp)eps % 区间长度小于eps时if abs(f2-fp)eps % df小于eps时退出if fpf2xmin xpfmin f(xp) % 极小值break;elsexmin x2fmin f(x2) % 极小值break;endendelseif fpf2if xpx2x3x2;x2xp;f3f2;f2fp;elsex1x2;x2xp;f1f2;f2fp;endelseif xpx2x1xp;f1fp;elsex3xp;f3fp;endendendendfx20/x;[xmin,fmin]main(f,2,10,0.01)xmin 4.4869fmin 8.9443室友的黄金分割法函数f(x) x20/xyellowking(f,2,10,0.01)function[xmin]yellowking(f,a,b,e)k0;a1b-0.618*(b-a); %插入点的值a2a0.618*(b-a);while b-ae %循环条件y1subs(f,a1);y2subs(f,a2);if y1y2 %比较插入点的函数值的大小aa1; %进行换名a1a2;y1y2;a2a0.618*(b-a);elseba2;a2a1;y2y1;a1b-0.618*(b-a);endkk1;end %迭代到满足条件为止就停止迭代xmin(ab)/2;fminsubs(f,xmin) %输出函数的最优值fprintf(k\n); %输出迭代次数disp(k);结果指令f(x)x20/xf (x)x20/x [x,y]gold(f,2,10,0.01)x 4.4683y 8.9443二次插值法函数function [xmin,fmin] main(f,a0,b0,epsilon)aa0;bb0;x1a;f1f(x1);x3b;f3f(x3);x25;f2f(x2);c1(f3-f1)/(x3-x1);c2((f2-f1)/(x2-x1)-c1)/(x2-x3);xp0.4*(x1x3-c1/c2);fpf(xp);while (abs(xp-x2)epsilon)if x2if f2fpf1f2;x1x2;x2xp;f2fp;elsef3fp;x3xp;endelseif f2fpf3f2;x3x2;f2fp;x2xp;elsef1fp;x2xp;endendc1(f3-f1)/(x3-x1);c2((f2-f1)/(x2-x1)-c1)/(x2-x3);xp0.5*(x1x3-c1/c2);fpf(xp);endif f2fpxminxp;fminf(xp);elsexminx2;fminf(x2);endend结果clear all;fx20/x;[xmin,fmin]main(f,2,10,0.01)xmin 4.4869fmin 8.94432-112-11我的2-11function [k ender]steepest(f,x,e)%梯度下降法,f为目标函数(两变量x1和x2)x为初始点,如[3;4]syms x1 x2 m; %m为学习率d-[diff(f,x1);diff(f,x2)]; %分别求x1和x2的偏导数即下降的方向flag1; %循环标志k0; %迭代次数while(flag)d_tempsubs(d,x1,x(1)); %将起始点代入求得当次下降x1梯度值d_tempsubs(d_temp,x2,x(2)); %将起始点代入求得当次下降x2梯度值nornorm(d_temp); %范数if(nore)x_tempxm*d_temp; %改变初始点x的值f_tempsubs(f,x1,x_temp(1)); %将改变后的x1和x2代入目标函数f_tempsubs(f_temp,x2,x_temp(2));hdiff(f_temp,m); %对m求导找出最佳学习率m_tempsolve(h); %求方程得到当次mxxm_temp*d_temp; %更新起始点xkk1;elseflag0;endendenderdouble(x); %终点endsyms x1 x2;fx1^2x2^2-x1*x2-10*x1-4*x260;x[0;0];e0.01;[k ender]steepest(f,x,e)ender 7.99615.9971室友的2-11:梯度函数function [k,ender]tidu(f,x,e)syms x1 x2 m;d-[diff(f,x1);diff(f,x2)];flag1;k0;while(flag)d_tempsubs(d,x1,x(1));d_tempsubs(d_temp,x2,x(2));nornorm(d_temp);if(nore)x_tempxm*d_temp;f_tempsubs(f,x1,x_temp(1));f_tempsubs(f_temp,x2,x_temp(2));hdiff(f_temp,m);m_tempsolve(h);xxm_temp*d_temp;kk1;elseflag0;endendenderdouble(x);end结果指令syms x1 x2;fx1^2x2^2-x1*x2-10*x1-4*x260;x[0;0];e0.01;[k ender]tidu(f,x,e)ender 7.99615.99712-12我的2-12展开为二阶泰勒式syms x1 x2;taylor(x1^42*x2^3-3*x1^2*x2)ans 3*x2 - 2*x1 - 6*(x1 - 1)*(x2 - 1) 3*(x1 - 1)^2 6*(x2 - 1)^2 - 1牛顿法求解function allnewton(f,x,e)syms x1 x2 h;d-[diff(f,x1);diff(f,x2)];hhessian(f);flag1;h1h^-1;while (flag)d_tempsubs(d,x1,x(1));d_tempsubs(d_temp,x2,x(2));nornorm(d_temp);if(nore)xxh1*d_temp;elseflag0;endendalldouble(x);结果指令clear all syms x1 x2;f3*x2 - 2*x1 - 6*(x1 - 1)*(x2 - 1) 3*(x1 - 1)^2 6*(x2 - 1)^2 - 1;x[1;1];e0.01;allnewton(f,x,e)all 1.16670.8333室友的2-12:展开为二阶泰勒式syms x1 x2;taylor(x1^42*x2^3-3*x1^2*x2)ans 3*x2 - 2*x1 - 6*(x1 - 1)*(x2 - 1) 3*(x1 - 1)^2 6*(x2 - 1)^2 - 1牛顿函数function allnewton(f,x,e)syms x1 x2 h;d-[diff(f,x1);diff(f,x2)];hhessian(f);flag1;h1h^-1;while (flag)d_tempsubs(d,x1,x(1));d_tempsubs(d_temp,x2,x(2));nornorm(d_temp);if(nore)xxh1*d_temp;elseflag0;endendalldouble(x);结果指令clear all syms x1 x2;f3*x2 - 2*x1 - 6*(x1 - 1)*(x2 - 1) 3*(x1 - 1)^2 6*(x2 - 1)^2 - 1;x[1;1];e0.01;allnewton(f,x,e)all 1.16670.83332-13(1)我的2-13(1)外点惩罚函数法function [ x,y ] Epfm_min( fx,gx,hx,xx0,s,c,a)%fx是目标函数%gx是不等式约束方程组(且g0)%xx0是初始点%hx是等式约束方程组(且h0)%s是精确度(s0)%c是放大系数(c1)%a是罚因子(默认为1)syms x1 x2xx1xx0;v[x1,x2];a1a;Pxk1;%假设Px等于1以免不必要错误G-subs(gx,v,xx1);%用于判别max{0,-g(x)}while Pxksif(G0)Pxa1*hx*hx;elsePxa1*hx*hxa1*gx*gx;endFxfxa1*Px;%将约束问题化为了一个无约束的问题% 接下来解min F(x)dFx1diff(Fx,x1);%分别对x1x2求偏导数dFx2diff(Fx,x2);[k,b]solve(dFx1,dFx2,x1,x2);%求出xx2xx1[k,b];Pxka1*subs(Px,v,xx2);xx1xx2;%相当于置kk1a1c*a1;%罚因子放大G-subs(gx,v,xx1);%用于判别max{0,-g(x)}endxxx1;ya1/c;syms x1 x2;fxx1x2;gx-x1;hxx1^2-x2;s10.^-5c10xx0[0,0]a1;[x,y]Epfm_min( fx,gx,hx,xx0,s,c,a) x[0.1;0.2];k0.1;e0.01;r1;[x,minf] Epfm_min (p,x,k,r,e)x 0.0015minf 0.0030Ans0.0045内点惩罚函数法function [x,minf]minNF(f,x0,g,u,v,var,eps)format long;if nargin6eps1.0e-4;endk0;FE0;for i1:length(g)FEFE1/g(i);endx1transpose(x0);x2inf;while 1FFu*FE;SumFfFF;[x2,minf]minNT(SumF,transpose(x1),var);BxFunval(FE,var,x2);if u*Bxif norm(x2-x1)epsxx2;break;elseuv*u;x1x2;endelseif norm(x2-x1)epsxx2;break;elseuv*u;x1x2;endendendminfFunval(f,var,x);format short;syms x1 x2 r1; ptaylor(x1x2-r1*(1/(x1^2-x2)-1/x1),[x1 x2],[0.001 0.002],Order,3); x[0.1;0.2];k0.1;e0.01;r1;[x,minf] minNF(p,x,k,r,e)x 0.0015minf 0.0030Ans0.0045室友的2-13内点惩罚函数function anllneicheng(p,x,k,r,e)syms x1 x2 r1;flag11;while (flag1)pdsubs(p,r1,r);xoldx;flag21;while (flag2)dp-[diff(pd,x1);diff(pd,x2)];hhessian(pd,[x1,x2]);h1h^-1;dp_tempsubs(dp,x1,x(1));dp_tempsubs(dp_temp,x2,x(2));nornorm(dp_temp);if(nore)xxh1*dp_temp;elseflag20;endendx_tempx;nor2norm(x_temp-xold);if double(nor2)erk*r;elseflag10;endendanlldouble(x);结果clear all; syms x1 x2 r1; ptaylor(x1x2-r1*(1/(x1^2-x2)-1/x1),[x1 x2],[0.001 0.002],Order,3); x[0.1;0.2];k0.1;e0.01;r1;anllneicheng(p,x,k,r,e)anll 0.00150.0030Ans0.0045外点惩罚函数function annnwaicheng(p,x,k,r,e)syms x1 x2 r1;flag11;while (flag1)pdsubs(p,r1,r);xoldx;flag21;while (flag2)dp-[diff(pd,x1);diff(pd,x2)];hhessian(pd,[x1,x2]);h1h^-1;dp_tempsubs(dp,x1,x(1));dp_tempsubs(dp_temp,x2,x(2));nornorm(dp_temp);if(nore)xxh1*dp_temp;elseflag20;endendx_tempx;nor2norm(x_temp-xold);if double(nor2)erk*r;elseflag10;endendannndouble(x);结果clear all; syms x1 x2 r1;ptaylor(x1x2,[x1 x2],[0.001 0.002],Order,3);x[0.1;0.2];k0.1;e0.01;r1;annnwaicheng(p,x,k,r,e)annn 0.00150.0030Ans0.00452-14我的(貌似这题抄的他的)2-14function [x,minf] minMixFun(f,g,h,x0,r0,c,var,eps)gx0 Funval(g,var,x0);if gx0 0;elsedisp(初始点必须满足不等式约束);x NaN;minf NaN;return;endif r0 0disp(初始障碍因子必须大于0);x NaN;minf NaN;return;endif c 1 || c 0disp(缩小系数必须大于0且小于1);x NaN;minf NaN;return;endif nargin 7eps 1.0e-6;endFE 0;for i1:length(g)FE FE 1/g(i);endFH transpose(h)*h;x1 transpose(x0);x2 inf;while 1FF r0*FE FH/sqrt(r0);SumF f FF ;[x2,minf] minNT(SumF,transpose(x1),var);if norm(x2 - x1)epsx x2;break;elser0 c*r0;x1 x2;endendminf Funval(f,var,x);Funval.mfunction fv Funval(f,varvec,varval)var findsym(f);varc findsym(varvec);s1 length(var);s2 length(varc);m floor((s1-1)/31);varv zeros(1,m);if s1 ~ s2for i0: ((s1-1)/3)k findstr(varc,var(3*i1));index (k-1)/3;varv(i1) varval(index1);endfv subs(f,var,varv);elsefv subs(f,varvec,varval);endSyms x1 x2;fx1^2-x2^2-3*x2;g1-x1;hx2-2;[x,minf]minMixFun(f,g,h,[2,2],2,0.5,[x1 x2 ],0.001)x 1.0015minf2.0002室友的2-14混合惩罚函数function [x,minf] MixPunish(f,g,h,x0,r0,c,var,eps)gx0 Funval(g,var,x0);if gx0 0;elsedisp(初始点必须满足不等式约束);x NaN;minf NaN;return;endif r0 0disp(初始障碍因子必须大于0);x NaN;minf NaN;return;endif c 1 || c 0disp(缩小系数必须大于0且小于1);x NaN;minf NaN;return;endif nargin 7eps 1.0e-6;endFE 0;for i1:length(g)FE FE 1/g(i);endFH transpose(h)*h;x1 transpose(x0);x2 inf;while 1FF r0*FE FH/sqrt(r0);SumF f FF ;[x2,minf] minNT(SumF,transpose(x1),var);if norm(x2 - x1)epsx x2;break;elser0 c*r0;x1 x2;endendminf Funval(f,var,x);Funval.mfunction fv Funval(f,varvec,varval)var findsym(f);varc findsym(varvec);s1 length(var);s2 length(varc);m floor((s1-1)/31);varv zeros(1,m);if s1 ~ s2for i0: ((s1-1)/3)k findstr(varc,var(3*i1));index (k-1)/3;varv(i1) varval(index1);endfv subs(f,var,varv);elsefv subs(f,varvec,varval);endSyms x1 x2;fx1^2-x2^2-3*x2;g1-x1;hx2-2;[x,minf]MixPunish(f,g,h,[2,2],2,0.5,[x1 x2 ],0.001)x 1.0015minf2.0002结束语无聊到这地步想必也是没谁了。不过刚考完我总不能一直玩手机啊。前几天重新看《盘龙》让我在手机上刚了一星期不能再这么毫无节制的玩耍了但是又不想学习所以只好写简书了~~不过网上毕竟这方面的资源不是很多我就算是为后来人做点好事吧让你们好找一点~~个人宣言知识传递力量技术无国界文化改变生活
