关于fminbnd的不解

cnzjtzszf

Q1=A(:,:,1)*(( transpose(A(:,:,1))*A(:,:,1))\(transpose(A(:,:,1))*Q(:,:,1)))
Q_2=(Q(:,:,1)-Q1).*(Q(:,:,1)-Q1);
f=@(tc1)sum(Q_2);
[tc1,f]=fminbnd(f,90,10800)
想通过以上程序求得变量tc1在区间(90,10800)内sum(Q_2)的最小值,可是出现以下问题,

??? Function 'le' is not defined for values of class 'sym'.
Error in ==> le at 18
  [varargout{1:nargout}] = builtin('le', varargin{:});
Error in ==> fminbnd at 277
    if fu <= fx
但是如果先计算出sum(Q_2)然后粘贴至f方程中又可以完成求值过程,请问是否因为存在变量的原因呢?如果想一气呵成应该怎样做呢?

请大家帮帮我(Q_2是含有tc1字符变量的变量)

minivet

“le”是什么？是已经赋值的一个数值变量么？

cnzjtzszf

le 好象是小于等于的逻辑符号,我在想怎么样才能将符号变量Q_S转化成fminbnd所需要的函数fun呢,谁来帮帮我,急啊

lengyunfeng

在Matlab中，函数（不管你是M函数或者inline函数，或者你这里的anonymous function）指的是能够创建局部独立空间的程序语句或者文件，在函数中定义的所有变量名对函数而言均是局部的，既不会影响其它程序、函数以及命令窗口工作空间的任何变量也不受其影响。你可以用下面语句做个试验：
a=8;
f=@(a)a^2+2*a+1;
f(2)
f(a)
看到了吧，在f定义之前的量对f没有任何影响影响。所以你在f=@(tc1)sum(Q_2);中定义的tc1并非就是你Q_2中的tc1，其具体区别在哪，你可以看看rocwoods大侠在本版中推销的那本书（:lol）。
想要一次性解决问题的话，建议你用M函数来构建，不用匿名函数，那玩意有点绕～～

cnzjtzszf

谢谢,我试试看

lin2009

1# cnzjtzszf

描述问题要完整。
你的程序中的A＝？
没看到在Q_2表达式中没看到tc1的影子。

le是Matlab内置函数
<=   Less than or equal.


Q1=A(:,:,1)*(( transpose(A(:,:,1))*A(:,:,1))\(transpose(A(:,:,1))*Q(:,:,1)))
Q_2=(Q(:,:,1)-Q1).*(Q(:,:,1)-Q1);
f=@(tc1)sum(Q_2);
[tc1,f]=fminbnd(f,90,10800)

cnzjtzszf

按五楼的意思是变量传递出现的问题,不知道能否设置全局变量来解决问题呢

lengyunfeng

我在前面所说的只是一个思路，你要是想用全局变量来解决这一问题的话，建议你按6#的意思先把代码给全再说哈～～

7# cnzjtzszf

cnzjtzszf

sorry,本人是初学者.
A(:,:,1)是一个含有字符参数tc1的4×9矩阵如下:
A(:,:,1)=[3 6 8 tc1 5 7 6 4 2;
             6 4 5 1 3 8 2 4 tc1;
             6 4 8 3 0 tc1 8 3 5;
             5 8 7 2 6 tc1 6 3 4;]
Q(:,:,1)=[5;7;3;4]