24小时热门版块排行榜

返回列表

xuetb

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 39
帖子: 3
在线: 4.5小时
虫号: 1007926
注册: 2010-04-28
性别: GG
专业: 交通工程

[求助] matlab最小二乘问题求助怎么定义多变量

最小二乘问题：
非线性最小二乘问题：minf^2 其中， f=[x1+10*x2;sqrt(5)*(x3-x4);(x2-2*x3)^2;sqrt(10)*(x1-x4)^2];x0=[3,-1,0,1]
程序如下，可以求单变量，怎么定义x=[x1,x2,x3,x4],求解呢？请高手指教
function fv = Funval( f,varec,varval )
%UNTITLED Summary of this function goes here
% Detailed explanation goes here
tic;
var = findsym(f);
varc = findsym(varec);
s1=length(var);
s2=length(varc);
m=floor((s1-1)/3+1);
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,varec,varval);
end
end
x=[x1,x2,x3,x4];
function [x,minf]= minGN(f,x0,var,eps)%目标函数;f;初始点:x0，自变量向量;var,
  %精度：eps，目标函数最小值时的x以及目标函数最小值
  format long;
  if nargin==4
   eps=1.0e-5;
  end
  S=transpose(f)*f;
  k=length(f);
  n=length(x0);
  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);
   Ax=Funval(A,var,x0);
   gSx=transpose(Ax)*Fx;
   dx=-transpose(Ax)*Fx;
   x0=x0+dx;
   tol=norm(dx);
  end
  x=x0;
  minf=Funval(S,var,x);
  format short;
  toc;
  Iterations;
end

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

C语言变量问题，求助大侠！！已经有17人回复
matlab非线性参数拟合问题已经有7人回复
【求助】关于多相流mixture模型中材料定义问题已经有12人回复
【求助】Matlab非线性最小二乘拟合活度系数模型（Willson、NRTL、UNIQUAC）已经有17人回复
【求助】【已解决】Gaussview建分子内坐标使用变量出错已经有15人回复
【求助】MATLAB多个mat中数据画图问题已经有7人回复
【求助】matlab 二次规划的优化的问题已经有4人回复
【求助】最小二乘问题（急）已经有5人回复

1楼2011-12-05 15:04:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

xuetb

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 39
帖子: 3
在线: 4.5小时
虫号: 1007926
注册: 2010-04-28
性别: GG
专业: 交通工程

补充：偶在matlab中输入： syms x=[x1 x2 x3 x4];f=[x1+10*x2;sqrt(5)*(x3-x4);(x2-2*x3)^2;sqrt(10)*(x1-x4)^2];x0=[3,-1,0,1];[x,minf]=minGN(f,x0,[x],1.0e-5),S,A,k
错误： Error using ==> syms at 61
Not a valid variable name.