早教吧作业答案频道 -->其他-->
Matlab问题求解.我写了个程序e=0.001;c=[5;4;3;2;1]Q=[5 4 3 2 1;4 4 3 2 1;3 3 3 2 1;2 2 2 2 1;1 1 1 1 1];x=[x1;x2;x3;x4;x5];f = @(x) (c.'*x + 1/2 * x.' * Q * x);df = @(x) (c + Q *x);x1=0;x2=0;x3=0;x4=0;x5=0;flag=1;k=0;Q=inv(Q);d=-Q*df(x);while
题目详情
Matlab问题求解.
我写了个程序
e=0.001;
c=[5;4;3;2;1]
Q=[5 4 3 2 1;4 4 3 2 1;3 3 3 2 1;2 2 2 2 1;1 1 1 1 1];
x=[x1;x2;x3;x4;x5];
f = @(x) (c.'*x + 1/2 * x.' * Q * x);
df = @(x) (c + Q *x);
x1=0;x2=0;x3=0;x4=0;x5=0;
flag=1;k=0;
Q=inv(Q);
d=-Q*df(x);
while flag
k=k+1;
if (norm(d) norm at 28
[varargout{1:nargout}] = builtin('norm',varargin{:});
然后再次用原程序运行(无syms的)就可以正常算出来了
我的Matlab版本是7.0.1的 请问这是什么问题能怎么解决
我写了个程序
e=0.001;
c=[5;4;3;2;1]
Q=[5 4 3 2 1;4 4 3 2 1;3 3 3 2 1;2 2 2 2 1;1 1 1 1 1];
x=[x1;x2;x3;x4;x5];
f = @(x) (c.'*x + 1/2 * x.' * Q * x);
df = @(x) (c + Q *x);
x1=0;x2=0;x3=0;x4=0;x5=0;
flag=1;k=0;
Q=inv(Q);
d=-Q*df(x);
while flag
k=k+1;
if (norm(d) norm at 28
[varargout{1:nargout}] = builtin('norm',varargin{:});
然后再次用原程序运行(无syms的)就可以正常算出来了
我的Matlab版本是7.0.1的 请问这是什么问题能怎么解决
▼优质解答
答案和解析
分析
首先,恕我冒昧,这程序真的是你写的?
如果是,这么简单的问题不应该看不出来啊.
1、之所以【运行时会出现?Undefined function or variable 'x1'.】,是因为程序中x1(还包括x2 x3 x4 x5)在使用的时候确实未定义.用到x1的最靠前的代码是:
x=[x1;x2;x3;x4;x5];
你可以看看,前面有x1定义吗?
2、之所以【加入一个syms x1 x2 x3 x4 x5;再运行,会提示?Function 'norm' is not defined for values of class 'sym'.】,是因为求范数的函数norm对sym类型没有定义,而你用syms命令定义x1-x5是sym类型,进而x、再进而d都是sym类型,对d求范数当然会出错.
3、之所以【然后再次用原程序运行(无syms的)就可以正常算出来了】,是因为你前面加了syms命令后,尽管执行到norm的时候出错,但在此之前,已经执行了下面这一行:
x1=0;x2=0;x3=0;x4=0;x5=0;
这样你再次运行的时候,x1-x5都已经定义,当然就不会出现第一个错了;而此时x和d也不再是sym类型,求范数没问题,也就不会出现第二个错了.
修改
经过上面的分析,相信你应该知道怎么修改了,对,没错,就是把这一句
x1=0;x2=0;x3=0;x4=0;x5=0;
移到
x=[x1;x2;x3;x4;x5];
的前面就可以了.
首先,恕我冒昧,这程序真的是你写的?
如果是,这么简单的问题不应该看不出来啊.
1、之所以【运行时会出现?Undefined function or variable 'x1'.】,是因为程序中x1(还包括x2 x3 x4 x5)在使用的时候确实未定义.用到x1的最靠前的代码是:
x=[x1;x2;x3;x4;x5];
你可以看看,前面有x1定义吗?
2、之所以【加入一个syms x1 x2 x3 x4 x5;再运行,会提示?Function 'norm' is not defined for values of class 'sym'.】,是因为求范数的函数norm对sym类型没有定义,而你用syms命令定义x1-x5是sym类型,进而x、再进而d都是sym类型,对d求范数当然会出错.
3、之所以【然后再次用原程序运行(无syms的)就可以正常算出来了】,是因为你前面加了syms命令后,尽管执行到norm的时候出错,但在此之前,已经执行了下面这一行:
x1=0;x2=0;x3=0;x4=0;x5=0;
这样你再次运行的时候,x1-x5都已经定义,当然就不会出现第一个错了;而此时x和d也不再是sym类型,求范数没问题,也就不会出现第二个错了.
修改
经过上面的分析,相信你应该知道怎么修改了,对,没错,就是把这一句
x1=0;x2=0;x3=0;x4=0;x5=0;
移到
x=[x1;x2;x3;x4;x5];
的前面就可以了.
看了 Matlab问题求解.我写了...的网友还看了以下:
如果x^2-px+q=(x+a)(x+b)那么p等于() 2020-04-27 …
已知x^2+px+q=(x+a)(x+b),当a 2020-05-22 …
设函数-x^2+7x-12>0的解集为A.(1)求集合A(2)设p:x属于A,q:x>a.且q是p 2020-05-23 …
√(a-6)+(b+4)^2=0,求|x+a|+|x+b|的最小值 2020-06-03 …
已知集合M=(x/x^2+px+2=0)N=(x/x^2-x-q=0)已知M={x|x^2+px+ 2020-06-14 …
两个集合相交为空集,它的含义是什么?解释的清晰明白的我再给高分,不胜感谢……设集合P{x^-x-6 2020-06-15 …
已知集合P={X|X^2-X-6<0},Q={X|a-1<x<a+1},当a为何值时,P交Q=空集 2020-07-30 …
p是q的什么条件(1)p:未位数是2的正整数,q:可以被2整除的整数(2)p:角A与角B是对顶角, 2020-07-30 …
已知a>1,设命题P:a(x-2)+1>0,命题Q(x-1)^2>a(x-2)+1.求使得P,Q都 2020-08-01 …
高一复数题(会任意一道都行)1.设复数z1、z2、A(A≠0)满足关系式z1z2拔+A拔z1+Az2 2020-10-31 …