matlab优化求解的问题

Uifhjvv

关于matlab优化求解的问题请教大家；程序报错截图如下，我分析报错的原因是迭代点（s,d1）改变时，（expect1，expect2）未改变，具体问题描述如下：
2 M" h: C; Q) @, Z& ]7 N( C1.expect1，expect2是s,d1的函数，当fsolve函数按照优化算法找到新的点（s',d1'）时，我的问题里时需要expect1，expect2也变化，这样fsolve的目标函数就也变化了；
$ f( w! }  U! z6 |! Z2.根据上述1的描述，这个优化问题是否应该将（expect1，expect2）也纳入优化求解目标中去？这样该如何实现？
clc;
5 p2 O- _% l  r/ U' P& [1 c- k5 Jclear all;
%%%  赋值外生给定的参数  %%%
( n2 M% H( L- ^( D/ wwn=0.3;
* P( t/ i  Z  |  [: P  Wwt=1-wn;
r=0.05;
beta=1/(1+r);
epislon=0.05;
9 v4 s; R. m! J$ x6 `6 W, K! Aphi=0.3;
y1=1;
1 D! F* Q7 q; u. s: Zyn2=1;
y2mean=1;
y3=1;
theta=0.9;

9 ]$ e/ E3 R: n+ P+ o# Y
# Z4 W  [$ _& d7 W3 x% 赋值(y2,s,d1)
y2=0.95:0.01:1.05;
s=0.1; %给s,d1的任意初值
d1=0.1;
x0 =[s,d1];
( G0 Z2 c( L# W
4 _+ y! b7 [7 U% `' o
%在这个for循环里得到均值expect1,expect2，改变(s,d1),均值expect1,expect2也改变
/ @2 _$ {( U9 h6 s1 }( [
' o, y4 P, `! A; U  g; G
for i=1:1:length(y2);
s=0.01;
d1=0.1;
8 ~8 V$ p- ?6 c3 [& A

* L4 @) V5 Y9 `; C: m* p%%% 判断是否紧约束 %%%
standard(i)=y2(i)*((beta/(1+beta)+phi/wt/(1-phi*wn/wt))-y3/((1+beta)*(1+r)))/(beta/(1+beta)+phi*wn/wt/(1-phi*wn/wt));
check(i)=d1/(1-s);
if check(i)>standard(i)
1 {8 I+ W5 d9 }. L2 ?    ct2(i)=((1-s)*y2(i)-d1+phi*(1-s)*y2(i))/(1-phi*wn/wt);
    d2(i)=phi*(1+r)*(wn/wt*ct2(i)+(1-s)*y2(i));
$ {! a: Q/ H4 e) [4 B& N) d( G    ct3(i)=(1-s)*y3-d2(i);
    mu(i)=wt/((1-s)*y2(i)-d1+d2(i)/(1+r))-beta*wt*(1+r)/((1-s)*y3-d2(i));
   
3 T5 h" Z7 P* [/ P$ d" o% i7 @else
    ct2(i)=((1-s)*y2(i)-d1+(1-s)*y3/(1+r))/(1+beta);
    d2(i)=((1-s)*y3-beta*(1+r)*(1-s)*y2(i)+beta*(1+r)*d1)/(1+beta);
9 i9 A* h# B; t3 x    ct3(i)=(1-s)*y3-d2(i);


5 }3 [5 }* F* c; r: I    mu(i)=0;
/ o/ g* F1 G7 S0 I) rend
! J. T& ^/ s0 U3 t3 D% We1(i)=wt/ct2(i);%每给定一个y2(i)，计算得到一个e1(i)，e2(i)
+ f* o$ T: @/ F7 Se2(i)=wt/ct2(i)*y2(i)+mu(i)*phi*y2(i)+beta*wt/ct3(i)*y3;
( [8 m" M) j3 N
. M; X9 |1 N6 ?" t. \. `
+ ?% L+ N% s' }3 t. y2 w9 H  [%%%求均值%%%
expect1=mean(e1(1,：));
/ T# L2 C1 t% Z! a  Pexpect2=mean(e2(1,：));
! m& S8 d9 c. L3 ?2 K: b$ Zend
) O% }& M" p* N3 `4 y

%我的问题中，expect1，expect2实质上是s,d1的函数
%在用fsolve优化求解的时候，每迭代一个新的（s,d1),都应该运行一次上述for这个循环，得到新的expect1、expect2再进行优化求解，
8 Z5 G1 s( r. ]& q  X, Y! t%也即（s,d1)变化时，expect1、expect2也要变化，如此fsolve的目标函数也变化
# J. u% W# s' Y%这样一个问题没有编程思路，请求大家给予指导，如果我的问题没有表述清楚，十分欢迎追问，谢谢。
8 g' J! T8 L- s! x; d
' Y+ {& m! Q" t8 T) N, v8 ^
6 H( Y6 Q( ^  L4 |9 S[x,fval,exitflag,output,]=fsolve(@(x)foc(x,expect1,expect2),x0);
2 v3 J" N; X4 Y% Q2 W+ X6 m! Z
# I2 a* y! M$ y
function f=foc(x,expect1,expect2)
% v0 P& r  s6 w6 Rwn=0.3;
wt=1-wn;
2 e" w+ p% C& V$ ~9 }9 i$ hr=0.05;
6 P% m8 G) [" k) E: t3 b4 ^7 qbeta=1/(1+r);
epislon=0.05;
2 A/ G! Q& T" {! z& I5 Nphi=0.3;
1 I/ ~+ {, c, m0 F% C  a! Ay1=1;
yn2=1;
) h) C+ Q+ r5 f! xy2mean=1;
1 I3 C# j+ e5 T+ `! Z, t7 o1 fy3=1;
theta=0.9;
s=0.01;
; Y, {6 D8 O! E# e4 n7 T2 Q: k8 _( Zd1=0.1;
$ Y) t, f; y4 p8 Q1 ^. g0 j9 Wct1=s*(1-theta)*y1+d1/(1+r);% ct1也与s,d1相关；

/ q- L& k, `' ~3 x
5 E' B0 f% V+ T# c: a0 gf(1)=wt/ct1-beta*(1+r)*expect1;
: u; L5 }# D5 I( J" `' J+ T% Uf(2)=wt/ct1*y1*(1-theta)-beta*expect2;
end

  H  D0 v* ?5 G) u; Q9 j/ |! v8 `

pTDbn25

帮你顶一下

zzz.dan

0 E7 f: N( E+ p. [你这个写的啥，你foc的函数的变量f1和f2里面x的变量呢，你这个方程写的有很大问题

nkkopd