| 查看: 587 | 回复: 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不久，希望大家帮忙下.新人没有多少金币。。 |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有242人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
