24小时热门版块排行榜    

查看: 1496  |  回复: 3

xuetb

铁虫 (初入文坛)

[求助] 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuetb

铁虫 (初入文坛)

补充:偶在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.
2楼2011-12-05 15:54:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhfzh

木虫 (正式写手)

【答案】应助回帖


余泽成(金币+1): 谢谢参与应助!欢迎常来程序语言版! 2011-12-08 14:06:03
你可以用strcat、  eval函数通过循环来实现多个变量的定义
3楼2011-12-05 17:11:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuetb

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by zhfzh at 2011-12-05 17:11:33:
你可以用strcat、  eval函数通过循环来实现多个变量的定义

能帮助编写下么? 如何实现?谢谢啦
4楼2011-12-06 00:54:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xuetb 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +7 大米饭! 2026-03-15 7/350 2026-03-16 10:25 by 了了了了。。
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
信息提示
请填处理意见