数学试验七.pdfVIP

数学实验七 陈海涛 核31 2013011784 实验目的： 掌握matlab优化工具箱的使用方法，对不同算法进行分析和比较。 练习使用无约束优化的方法建立模型，求解模型。 第五题： 如果单纯在平面结构上讨论，于是我们可以取第一个点为坐标原点，对于第i个 原子，我们用（x,y ）来表示它的具体位置，于是当目标函数 i i 2 2 2 ((x x ) (y y ) d )  i j i j i,j i,j 取得最小值时候，我们可以认为这时候得到的x,y就是所要求的坐标 于是，我们在matlab中输入 然后，我们再定义一个函数： function f = question5(x,d,c) f=0; for i=1:52 f(i)=f+((x(1,d(i,1))-x(1,d(i,2))).^2+(x(2,d(i,1))-x(2,d(i,2))).^2-d(i, 3)).^2; end end c=zeros(1,25); x=lsqnonlin(@question5,x0,d,c) 警告:Lengthoflowerboundsislength(x);ignorin extrabounds. Incheckboundsat27 Inlsqnsetupat62 Inlsqnonlinat181 警告:Lengthofupperboundsislength(x);fillin inmissin upperboundswith+Inf. Incheckboundsat47 Inlsqnsetupat62 Inlsqnonlinat181 Exitin duetoinfeasibility: 25lowerboundsexceedthecorrespondin upperbounds. 在这里出现警告，并且输出的x为一个2×25的全为1的矩阵，肯定是在哪个环节出现了问 题 [x,norm,res]=lsqnonlin(@question5,x0,[],[],[],d,c) 警告:Trust-region-reflective algorithm requires at least as many equations as variables; usin Levenberg-Marquardtalgorithminstead. Inlsqncommonat55 Inlsqnonlinat236 Initialpointisalocalminimum. Optimizationcompletedbecausethesizeofthegradientattheinitialpoint islessthanthedefaultvalueofthefunctiontolerance. stoppin criteriadetails x= Columns1through23 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 Columns24through25 1 1 1 1 norm= 1.9688e+03 res= 44.3710 这里误差平方和和误差向量都是非常大的，可能是我的初始值取得不对 我们取 X=[0.1:0.1:2.5;2.5:-0.1:0.1];之后得到这个结果 于是，我们在图中画出这些点 2.5