| 查看: 492 | 回复: 0 | ||
[求助]
用Gauss-Neuton算法解非线形最小二乘法问题
|
|
%xdata，ydata1是测得的数据 syms z1 z2 z3 z4 z5 xdata = S1o; f = ydata1 - exp(z1)./(exp(z1)+exp(-z1)).*exp(-(xdata./exp(z3)).^exp(z2))+(1-exp(z1)./(exp(z1)+exp(-z1))).*exp(-(xdata./exp(z5)).^exp(z4)); x0 = [1.7380;0.9555;-1.3863;2.0149;-2.5257]; var = findsym(f); eps = 1.0e-3; minGN %minGN是m文件里的一个函数，用于实现G-N算法，这个程序是自己网上找的，后面的注释是为了我自己方便理解自己加的，具体如下： function [x,minf] = minGN(f,x0,var,eps) format long; if nargin == 3 %如果没有设置eps，则eps=1.0e-6 eps = 1.0e-6; end S = transpose(f)*f; %trnspose是转置 k = length(f); % n = length(x0); %n为变量个数 x0 = transpose(x0); %将行向量转置为列向量 tol = 1; %初始容差 A = jacobian(f,var); while tol>eps Fx = zeros(k,1); for i=1:k Fx(i,1) = Funval(f(i),var,x0); end Sx = Funval(S,var,x0); %求得f转置与其本身值 Ax = Funval(A,var,x0); %雅克比矩阵的值 gSx = transpose(Ax)*Fx; dx = -transpose(Ax)*Ax\gSx; x0 = x0 + dx; tol = norm(dx); %tol=dx绝对值之和 end minf = Funval(S,var,x); format short; x = x0; %%Funval函数也是在网上找到的资源，后面的注释原本就有，具体如下： function fv = Funval(f,varvec,varval) var = findsym(f);%找出表达式 即f函数包含的变量，即s,t varc = findsym(varvec);%找出传递参数的变量，即[t s]中的t,s s1 = length(var);%函数的个数，例子中即s t 2个 s2 = length(varc);%变量个数，即[t s]两个 m =floor((s1-1)/3+1);%floor(3.6)=3,即最靠近左边的整数 varv = zeros(1,m); if s1 ~= s2 for i=0: ((s1-1)/3) k = findstr(varc,var(3*i+1)); index = (k-1)/3; varv(i+1) = varval(index+1); end fv = subs(f,var,varv); else fv = subs(f,varvec,varval); %如果原来函数变量个数和传递的参数中变量个数一致 这调用subs函数 即计算在给点点处的函数值 end 现在在命令窗口输入上述命令，提示： 错误使用 minGN (line 6) 输入参数的数目不足。 谁能告诉我这是什么原因呢，参数数目应该是够的？另外，麻烦大侠再解释一下 Funval函数里m =floor((s1-1)/3+1);这条语句在函数的作用，尤其是为什么取((s1-1)/3+1）？谢谢了！刚接触matlab不久，希望大家帮忙下.新人没有多少金币。。 |
回复此楼» 猜你喜欢
哈尔滨理工大学物理系招收物理学考研调剂
已经有7人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
物理学I论文润色/翻译怎么收费?
已经有144人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有31人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有68人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
10