函数问题

cpupygu

8 l+ E9 P4 U4 h4 }9 Iclear;
t=1.43;
t0=0.562;%初始参数
, p" x& z2 G. v/ v1 \zeta_m=sqrt(2)/2;%ζm
A=(4*t+t0)/(2*t*t0);
B=t0^2*t;
C=t0+t;%简化参数
syms x1;
syms x2;
& v1 O0 s, m' R  X; ]/ \+ rw0=A/(zeta_m*x1);
a=A*(1-1/x1);
kp=0.5*(a*zeta_m+w0)*a*w0*B-1;
ki=a^2*w0^2*B*0.25;
kd=0.25*(a^2+4*a*zeta_m*w0+w0^2)*B-C;%三个PID参数
% A, |0 a4 h7 y% o& mKiw=(kd*(1i*x2)^2+kp*(1i*x2)+ki)/(1i*x2);
# @" _% D7 W; q/ ]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);
0 W' a* o9 `4 o3 D4 A" @5 r' T5 G5 L; @Siw=1/(1+Kiw*Giw);
f=norm(Siw,inf);
7 u, r' l0 ~0 o. G9 Q( T- G) @- ox=[7;1];
5 R. i+ e9 b1 E" Y- Ze=10^(-2);
[k, ender]=steepest(f,x,e);
  |6 r+ d; k  u" _: q: b; t  }( t* efunction [k, ender]=steepest(f,x,e)
! }; c  i2 Z; Q8 w% L$ q; K- U%梯度下降法,f为目标函数（两变量x1和x2），x为初始点,如[3;4]
5 ~; \+ r. m& x+ E& \syms x1 x2 m; %m为学习率
d=-[diff(f,x1);diff(f,x2)];  %分别求x1和x2的偏导数，即下降的方向
flag=1;  %循环标志
& z3 K5 N' D0 F0 ak=0; %迭代次数
while(flag)
6 {2 C$ X& Z9 R( c2 V    d_temp=subs(d,x1,x(1));      %将起始点代入，求得当次下降x1梯度值
; T  c) [/ d/ V! j  p8 Z  ]    d_temp=subs(d_temp,x2,x(2)); %将起始点代入，求得当次下降x2梯度值
    nor=norm(d_temp); %范数
    if(nor>=e)
) u! k. A2 V* S        x_temp=x+m*d_temp;            %改变初始点x的值
. V1 C; ~% Y2 R9 n, |/ G        f_temp=subs(f,x1,x_temp(1));  %将改变后的x1和x2代入目标函数
        f_temp=subs(f_temp,x2,x_temp(2));
        h=diff(f_temp,m);  %对m求导，找出最佳学习率
2 D1 }/ D+ M2 ?' \+ B& ~$ U% u        m_temp=vpasolve(h);   %求方程，得到当次m
        x=x+m_temp*d_temp; %更新起始点x
0 D3 l3 k- K6 T( |        k=k+1;
    else
* ?/ B7 u% _6 S6 p9 x& C        flag=0;
$ Z# I7 k2 \; t) F+ j) s    end
- N, |% P2 p, I- Cend
ender=double(x);  %终点
0 s7 ]# T" q1 p0 ~0 @end
% l& j% [/ q5 m1 p; G
( @7 K  `& ^" `# |: I# W

9 o% M7 y% U5 c运行后提示
beiyong
错误使用 symengine
The dimensions do not match.

/ B0 B* j5 M- u7 u出错 sym/privBinaryOp (line 946)
+ `; h, |/ a+ e& _            Csym = mupadmex(op,args{1}.s, args{2}.s, varargin{:});
7 W, ?% x* L' g
出错  *  (line 273)
        X = privBinaryOp(A, B, 'symobj::mtimes');

$ }/ ~5 g7 S; _0 `; S6 h; X出错 beiyong>steepest (line 38)
        x=x+m_temp*d_temp; %������ʼ��x
: G& z7 V2 c6 H% S
出错 beiyong (line 21)
[k, ender]=steepest(f,x,e);
7 s" B9 `  `! k3 r4 S/ T求助，感谢！

求助

cichishia

帮你顶一下

IBB-EUT

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

zzz.dan

来学习学习

nkkopd

楼上正解，楼主可以试试