函数问题

cpupygu

) x8 A$ \1 D' b$ E1 _- d) uclear;
' G  R- H2 K8 pt=1.43;
t0=0.562;%初始参数
' s* D9 r6 u5 d" U5 lzeta_m=sqrt(2)/2;%ζm
A=(4*t+t0)/(2*t*t0);
$ w8 [4 @# ?! y. h1 Y" s7 `" y- h% ^B=t0^2*t;
6 R+ |( e/ C* Q! V5 O" Z7 J) o% B; KC=t0+t;%简化参数
syms x1;
syms x2;
w0=A/(zeta_m*x1);
# o* [% j! d$ ^% _3 ia=A*(1-1/x1);
$ n) e8 e: b  A, u% K0 Lkp=0.5*(a*zeta_m+w0)*a*w0*B-1;
. U8 c5 A. B) _7 Xki=a^2*w0^2*B*0.25;
kd=0.25*(a^2+4*a*zeta_m*w0+w0^2)*B-C;%三个PID参数
2 E4 m9 v$ k5 f; QKiw=(kd*(1i*x2)^2+kp*(1i*x2)+ki)/(1i*x2);
Giw=((1i*x2)^2-10.68*(1i*x2)+37.99)/(1.43*(1i*x2)^3+16.27*(1i*x2)^2+65.01*(1i*x2)+37.99);
Siw=1/(1+Kiw*Giw);
f=norm(Siw,inf);
x=[7;1];
e=10^(-2);
9 j# q+ F* B% c; e! C, K[k, ender]=steepest(f,x,e);
function [k, ender]=steepest(f,x,e)
1 n- m. w5 E7 r/ A" W%梯度下降法,f为目标函数（两变量x1和x2），x为初始点,如[3;4]
syms x1 x2 m; %m为学习率
d=-[diff(f,x1);diff(f,x2)];  %分别求x1和x2的偏导数，即下降的方向
/ J6 }1 q1 o; Lflag=1;  %循环标志
% ~# q9 H' Y+ mk=0; %迭代次数
while(flag)
3 s$ _; ~& O8 `# H- Z  Q  [    d_temp=subs(d,x1,x(1));      %将起始点代入，求得当次下降x1梯度值
    d_temp=subs(d_temp,x2,x(2)); %将起始点代入，求得当次下降x2梯度值
3 W1 C5 N* I* L' j2 o" C: q    nor=norm(d_temp); %范数
! Y( F: c+ Q, q. J    if(nor>=e)
        x_temp=x+m*d_temp;            %改变初始点x的值
5 a4 @$ O7 L. j        f_temp=subs(f,x1,x_temp(1));  %将改变后的x1和x2代入目标函数
2 w& s8 L) W4 @: y+ O1 e) w2 q: X        f_temp=subs(f_temp,x2,x_temp(2));
        h=diff(f_temp,m);  %对m求导，找出最佳学习率
        m_temp=vpasolve(h);   %求方程，得到当次m
) K2 _' t: X: X' q& a2 w1 |        x=x+m_temp*d_temp; %更新起始点x
3 u) {4 h1 N% a: Z8 j! ~        k=k+1;
    else
        flag=0;
: J" d7 r' t* e( a5 J" K/ g' H    end
end
  M: T& m! B" \6 m6 M; Aender=double(x);  %终点
end
" T: Z9 K* i- h# A9 U4 z) q
; Z5 w7 L, |  W- `3 t

% K% D* z' h9 ]8 ?运行后提示
beiyong
错误使用 symengine
The dimensions do not match.

  ]1 y/ }% J  A: [' M2 {5 C( c+ \! q6 `出错 sym/privBinaryOp (line 946)
            Csym = mupadmex(op,args{1}.s, args{2}.s, varargin{:});
6 D8 q4 x/ M8 E' f
7 a6 [$ P2 s( g2 b* G( a2 j出错  *  (line 273)
        X = privBinaryOp(A, B, 'symobj::mtimes');
$ Y) K1 N# U$ C$ E7 \
出错 beiyong>steepest (line 38)
        x=x+m_temp*d_temp; %������ʼ��x

出错 beiyong (line 21)
1 G' w9 p0 e# m9 J, ][k, ender]=steepest(f,x,e);
# R7 g8 Y9 g8 ?. k求助，感谢！
/ ^  R0 Z8 h% o, f$ J5 z. e+ x, ]  W
: o$ N* P9 C5 S  G" \) v$ l求助

1 q- p, e# [' ~7 a2 |

nkkopd

楼上正解，楼主可以试试

zzz.dan

来学习学习

IBB-EUT

个人认为在自定义函数求解方程时m_temp=vpasolve(h);   %求方程，得到当次m存在问题，求解结果为空，
改用：m_temp=solve(h,m,'ReturnConditions',true);

cichishia

帮你顶一下