EDA365电子论坛网
标题:
如何求解如下优化问题
[打印本页]
作者:
smileqq
时间:
2020-12-2 14:55
标题:
如何求解如下优化问题
) z' \% o0 W# K. r! _) T9 ^
f=min(abs(0.5*ay_tp1)+abs(0.5*ay_tp2));
8 D; ?7 O: o+ t: J) ~# _
s.t.
( W Y7 _6 }- O5 q) W/ \
0<p<10,-10<q<0,q^2-4*p<=0;2*pi/(4*p-q^2)^(1/2)-tf=0;
3 S* X; z% F: P6 ~6 n
%%其中：
4 C/ r! X, ]* X- J# P6 p
w=3.75;
) N, J0 {! E* l0 @7 m0 U5 {
tf=5;
% X G8 H0 \- B5 s: R7 x
tp1=0;
9 a$ X2 t1 v# v- k' b. M9 d+ n( G% B
tp2=2*(pi-atan(q*(4*p-q^2)^(1/2)/(q^2-4*p)))/((4*p-q^2)^(1/2));
, f1 M9 o, G% w1 V. m1 O6 F
A1=(4*p/(4*p-q^2))^(1/2);
! C! v, n$ Z$ F/ L U
A2=((4*p-q^2)/2)^(1/2);
5 M8 y3 q( P# ]; J$ |' ~ w$ l
A3=q/((4*p-q^2)^(1/2));
+ }8 S" {! I' s6 [( X/ L
, `2 f4 H/ j/ J/ i$ \& n7 |
ay_tp1=(w*p*exp(q*tp1/2)*A1*sin(((A2/(2^(1/2)))*tp1)+1/A3))/(1-exp(q*tf/2)*A1*sin(A2*tf-atan(A3)));
, Q: D, z+ g+ F$ N: z
ay_tp2=(w*p*exp(q*tp2/2)*A1*sin(((A2/(2^(1/2)))*tp2)+1/A3))/(1-exp(q*tf/2)*A1*sin(A2*tf-atan(A3)));
2 ~5 ?" _8 u! V0 @; D
使用fmincon函数求解时显示 相邻两次迭代点的变化小于预先给定的容忍度。
* ~2 `" v6 W, I9 x% o3 R
求助大佬。感谢
) }& a% H7 I. a: N# x9 S
作者:
SsaaM7
时间:
2020-12-2 15:56
路过，帮顶一下
作者:
ExxNEN
时间:
2020-12-2 15:56
看不懂啊，有心无力哦
作者:
haidaowang
时间:
2020-12-2 15:57
下面结果参考下：
+ }. s0 @# ] C9 ]; h2 D
Objective Function (Min.): 1.28973902645915E-5
1 z$ ?& ]4 j# w# w4 t2 \/ x
p: 0.394784176495481
+ y: t. ~# x( W7 g
q: -4.3001505065871E-5
- Z' d$ c% l% O
& R, `. T, L" V* ~
Constrained Functions:
8 B2 D0 U( B# G/ M
1: q^2-4*p-0 = -1.5791367041328
- [! h+ U& K# a; d0 w
2: 2*pi/(4*p-q^2)^(1/2)-5-0 = 6.57038867757365E-11
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2