Matlab求解非线性超定方程组恰定方程组欠定方程组Word格式文档下载.docx

1、x0=,;fsolve（fun,x0,optimset（fsolve）注：.为续行符 m 文件必须以 function 为文件头，调用符为 ；文件名必须与定义的函数名相同； fsolve （）主要求解复杂非线性方程和方程组，求解过程是一个逼近过程。Matlab 求解线性方程组AX=B或 XA=B在MATLAB中，求解线性方程组时，主要采用前面章节介绍的除法运算符 “/和”。如：X=AB表示求矩阵方程AX= B的解;X= B/A表示矩阵方程XA=B的解。对方程组X= AB，要求A和B用相同的行数，X和B有相同的列数，它的行数等于矩阵A的列数，方程X= B/A同理。如果矩阵A不是方阵，其维数是 m

2、xn，则有:m= n恰定方程，求解精确解;mn 超定方程，寻求最小二乘解；mm。则方程组没有精确解，此时称方程组为超定方程组。 线性超定方程组经常遇到的问题是数据的曲线 拟合。对于超定方程，在 MATLAB中，禾I用左除命令（x=Ab）来寻求它的最小二乘解；还可以用广义逆来求，即 x=pinv（A）所得的解不一定满足Ax=b, x只是最小二乘意义上的解。 左除的方法是建立在奇异值分解基础之上， 由此获得的解 最可靠；广义逆法是建立在对原超定方程直接进行 householder 变换的基础上， 其算法可靠性稍逊与奇异值求解，但速度较快；求解超定方程组A=2 -1 3;3 1 -5;4 -1 1;

3、1 3 -13A=2-1 331 -54-1 1 1 3 -13 b= 3 0 3 -6； rank（A） ans=x1=Ab x1= x2=pinv（A）*b x2=A*x1-b ans=可见x1并不是方程Ax=b的精确解，用x2=pinv（A）*b所得的解与x1相同。三欠定方程组欠定方程组未知量个数多于方程个数，但理论上有无穷个解。 MATLAB将寻求一个基本解，其中最多只能有 m 个非零元素。特解由列主元 qr 分解求得。解欠定方程组A= 1-2 1 1;1 -2 1 -1;1 -2 1 51 -2 11 -2 1 b=1 -1 5 x1=AbWarning:Rank deficient

4、,rank=2 tol= x1= x2=pinv（A）*b x2=四方程组的非负最小二乘解 在某些条件下， 所求的线性方程组的解出现负数是没有意义的。 虽然方程组可以 得到精确解，但却不能取负值解。 在这种情况下， 其非负最小二乘解比方程的精 确解更有意义。在MATLAB中，求非负最小二乘解常用函数nnIs,其调用格式为:（1）X=nnls（A,b返回方程Ax=b的最小二乘解，方程的求解过程被限制在x的条件下；（2）X=nnls（A,b,TOl指定误差TOL来求解，TOL的默认值为TOL=max（size（A）*norm（A,1）*eps矩阵的1范数越大，求解的误差越大;（3）X,W=nnls

5、（A,b）当 x（i）=0 时，w（i）0 时，w（i）0,同时返回一个双 向量 w 。例 9 】求方程组的非负最小二乘解A=;b= ;X,W=nnls（A,b）X=W= x1=Ab x1=A*X-b ans=关于采用 matlab 进行指定非线性方程拟合的问题 （1）探1。优化工具箱的利用描述函数LSQLIN 有约束线性最小二乘优化LSQNONNEG 非负约束线性最小二乘优化问题阅优化工具箱文档中的相应利用这两个函数的例子。d. 非线性曲线拟合利用MATLAB的内建函数FMINBND 只解决单变量固定区域的最小值问题FMINSEARCH多变量无约束非线性最小化问题（ Nelder-Mead方

6、法）。面给出一个小例子展示一下如何利用 FMINSEARCH1 首先生成数据 t=0:.1:10; t=t（:）; Data=40*exp*t）+rand（size（t）; % 将数据加上随机噪声2 .写一个m文件，以曲线参数作为输入，以拟合误差作为输出fun cti on sse=myfit （p arams,l npu t,Actural_Out put）A=p arams（1）;lamda=p arams（2）;Fitted_Curve=A.*ex p（-lamda*l nput）;Error_Vector=Fitted_Curve-Actural_Out put;%当曲线拟合的时候，一

7、个典型的质量评价标准就是误差平方和 sse=sum（Error_Vector.A2）;% 当然，也可以将 sse 写作：sse=Error_Vector（:）*Error_Vector（:3.调用 FMINSEARCH Strart in g=ra nd（1,2）; op ti ons=op timset（Dis play,ite门; Estimates=fim in search（myfit,Strart ing,op ti on s,t,Data）; plot（t,Data,* hold on plot（t,Estimates（1）*ex p（-Estimates（2）*t）,rEstim

8、ates将是一个包含了对原数据集进行估计的参数值的向量。附图见后面:FMINSEARCH通常能够用来解决不连续情况，特别是如果他们不出现在解的附近的时候。它得到的通常也是局部解。FMINSEARCH只能够最小化实数值（也就是说，解的域必须只能包括实数，函数的输出只能够为实数值）。 当感兴趣的是复数变量的域的时候，他们必须被分 割为实部与虚部。的FIGURE窗口：最基本的拟合界面与数据统计工具MATLAB通过基本的拟合界面也支持基本曲线拟合。利用这个界面，你可以快速地在简单易用的环境中实现许多基本的曲线拟合。这个界面可以实现以下功能：a 通过比样条插值（ spline interpolant ）

9、、 hermite 插值、或者是高达 10 阶的多项式插值 实现数据的拟合；b对给定数据同时实现多样插值的绘制；c绘制残差图；d检查拟合结果的残差的数值；e通过内插值或者外推插值评价一个拟合结果；f对拟合结果和残差的模进行图形绘制；g将拟合结果保存入 MATLAB工作空间。开发你的拟合应用的时候，你可以通过基本拟合（ Basic Fitt ing）界面，也可以通过命令行函数，也可以同时作用。你可以通过基本拟合界面只能够实现2D 数据的拟合。然而，如果你用 subplot 绘制多个数据集，只要有至少一个数据集是2D的，那么就可以用基本拟合界面。可以通过如下步骤激活基本拟合界面：1 绘制数据；2

10、从 figure 窗口的 Tools 菜单条下面选择 Basic Fitting菜单项；有关Basic Fitting界面的更多信息，请查阅 MATLAB帮助文档的相应部分。对于HP，IBM以及SGI平台，（）以及的基本拟合界面不受支持。统计界面激活：1 制数据；2 从 figure 窗口的 Tools 菜单条下面选择 Data Statistics 菜单项；关于采用 matlab 进行指定非线性方程拟合的问题 （2）一。优化工具箱函数LSQNONLIN 解决非线性最小二乘法问题，包括非线性数据拟合问题LSQCURVEFIT解决非线性数据拟合问题面给出利用这两个函数的例子：LSQNONLIN禾

11、U用这个函数最小化连续函数只能够找到句柄解。下面的例子说明利用LSQNONLIN函数用下面的函数进行拟合:f = A + B exp（C*x）+D*exp（E*x）对数据集x与y进行拟合，其中y是在给定x的情况下的期望输出（可以是方程给出数组,也可以是单独数据组成的数组 ）。为了解决这个问题，先建立下面的名为的函数，它利用数据x与y,将他们作为优化输入参数传递给 LSQNONLIN利用给定的数据x 计算 f 的值，再与原始数据 y 进行比较。经验值与实际计算出的值之间的差异作为输出值返回。LSQNOLIN函数就是最小化这些差的平方和。function diff = fit_simp（x,X,Y

12、）%此函数被LSQNONLIN调用% x 是包含等式系数的向量% X 与 Y 是作为操作数传递给 lsnonlinA = x（1）;B = x（2）;C = x（3）;D = x（4）;E = x（5）;diff = A + B.*exp（C.*X）+D.*exp（E.*X）-Y;面的脚本是利用上面定义的函数的一个例子： 定义你打算拟合的数据集合 X=0:.01:.5; Y=.*exp.*X）+.*exp.*X）+.*rand（size（X）; 初始化方程系数 X0=1 1 1 1 1;% 设置用中等模式（ memdium-scale ）算法 options=optimset（Largesca

13、leoff%通过调用LSQNONLIN重现计算新的系数 x=lsqnonlin（fit_simp,X0,options,X,Y）;%调用LSQNONLIN结果输出表明拟合是成功的Optimization terminated successfully:Gradient in the search direction less than tolFunGradient less than 10*（tolFun+tolX）% 绘制原始数据与新的计算的数据 Y_new=x（1）+x（2）.*exp（x（3）.*X）+x（4）.*exp（x（5）.*X）; plot（X,Y,+r,X,Y_new,b注意

14、：LSQNONLIN只可以处理实数变量。在处理包括复数变量的实例的拟合的时候，数 据集应该被切分成实数与虚数部分。下面给出一个例子演示如何对复数参数进行最小二乘拟 合。为了拟合复数变量，你需要将复数分解为实数部分与虚数部分， 然后把他们传递到函数中去,这个函数被LEASTSQ乍为单个输入调用。首先，将复数分解为实部与虚部两个向量。其次,复数数据，并用你想拟合的复数方程计算。 将输出向量分解实部与虚部， 将这两部分连接为一个单一的输出向量传递回 LEASTSQ下面，给出一个例子演示如何根据两个复数指数拟合实数X与Y。建立方程:fun ction zero = fit2（x,X,Y） %根据输入x

15、重建复数输入 cmpx = x（1:4）+i.*x（5:8）;%利用复数计算函数 zeroco mp = cmp x（1）.*ex p（cmp x（2）.*X） + cmp x（3）.* exp（cmp x（4）.*X）-Y;%将结果转换成一个列向量% 其中第一部分是实部，第二部分是虚部numx = length（X）; % 实部长度zero=real（zerocomp）; % 实部zero（numx+1:2*numx）=imag（zerocomp）; % 虚部为了评价计算这个函数， 需要X与丫数据集。LSQNONLIN将根据它拟合出下面方程中的参 数 a，b，c 与 d：Y = a*exp（

16、b*X）+c*exp（d*X）;其中， a， b，c 与 d 是复数变量。5; Y=sin（X）; Y=Y+.1*rand（size（Y）; cmpx0=1 i 2 2*i; x0（1:4）=real（cmpx0）; x0（5:8）=imag（cmpx0）; x=leastsq（fit2,x0,X,Y）; cmpx=x（1: Y1=real（cmpx（1）.*exp（cmpx（2）.*X）+cmpx（3）.*exp（cmpx（4）.*X）; plot（X,Y1,+二。LSQCURVEFIT利用此函数可以在最小二乘意义上解决非线性曲线拟合（数据拟合）问 题。也就是说，给定输入数据xdata，以及

17、观测的输出数据 ydata，找到系数x，使得函数F（x,xdata）能够最好的拟合向量值。 LSQCURVEFI利用与LSQNONLIN相同的算法。它的目的在于专门为数据拟合问题提供一个接口。3维参数拟在拟合的时候，2维、3维或者N维参数拟合是没有什么差别的。下面给出一个 合的例子。待拟合函数是:z = a1*y.*x.A2+a2*si n（x）+a3*y.A3;建立的的函数如下:function F = myfun（a, data）;x = data（1,:y = data（2,:F= a（1）*y.*x.2+a（2）*si n（x）+a（3）*yA3;面的脚本展示了这么利用上面的函数： d

18、ata=xdata; ydata; a0= 10, 10, 10; % 初识揣测 a, resnorm = lsqcurvefit（myfun,a0,data,zdata）Maximum number of function evaluations exceeded;increase resnorm = +004 format long a option=optimset（MaxFunEvals,800）; a, resnorm = lsqcurvefit（myfun,a0, data, zdata, , , option）Relative function value changing by

19、 less than resnorm = +004统计工具箱函数函数名nlinfit （非线性回归） 采用 Gauss-Newton 法进行非线性最小二乘数据拟合Iscov （线性回归） 根据已知协方差矩阵进行最小二乘估计regress 多元线性回归regstats 回归诊断ridge 脊回归（ Ridge regress） rstool 多维响应表面可视化（ RSV） stepwise 交互式逐步回归solve 求解出所有的解。具体例子请参阅相应文档在 Matlab 中如果求解一个二元二次方程，现在只有两个方程可以用 但是当存在多余 2 个方程的时候（即为超定方程组时） ，该如何求解呢还能用

20、 Solve 函数吗比如下面的方程：a1*xA2+b1*x+c1*yA2+d1*y=e1 a2*xA2+b2*x+c2*yA2+d2*y=e2 a3+b3*x+c3*y=e3x,y 为未知数，其他的为已知数，求解 x,y.用fmin求e1A2+e2A2+e3A2的最小解你的意思是将方程转换成：a1*xA2+b1*x+c1*yA2+d1*y-e仁f1a2*xA2+b2*x+c2*yA2+d2*y-e2=f2a3+b3*x+c3*y-e3=f3再用 fminsearch 求解最小值 min（f1+f2+f3 ），得到最终的未知数的值是不是？很有道理。但是这个 fminsearch 是需要初始解才能得到最终的解的。假如不知道初始解该如何求解呢1.fminsearch求解最小值min（卄人2+2人2+3人2 ）2.2.一般非线性问题的数值解往往都是迭代法，因而需要初值3.3.初值可以试出来或者对原始方程作简化然后定性分析出一个解的大体范围来确定初值。4.4.非线性问题很少有一劳永逸的一统天下的解法5.5.正因为如此，才让我们有事做，有饭吃