|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
<P>function A=f(x)<BR>A=100*(x(2)-x(1)^2)^2+(1-x(2))^2;<BR>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<BR>function f_grad = grad(f,x)</P>
<P>% 求符号函数f在x处的梯度<BR>% 注意自变量必须为x1,x2,.....<BR>% x是列向量<BR> <BR>dim=size(x);<BR>for j=1:dim(1,1)<BR>for j=1:3<BR> grad_f(1,j) = diff(f,x(j,1)); % 符号偏导<BR>end<BR>f_grad = grad_f.'; % 将符号数字转化为实际数值,列向量<BR><BR>%%%%%%%%%%%%%%%%%%%%%%%%%%5<BR>上面是两个m文件<BR>下面运行<BR>sym x<BR>x=[1 1 2]' %列向量<BR>grad(f,x)<BR>结果报错<BR>麻烦个位大侠帮改下<BR>谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢谢<BR><BR></P> |
|