matlab优化求解的问题

Uifhjvv

关于matlab优化求解的问题请教大家；程序报错截图如下，我分析报错的原因是迭代点（s,d1）改变时，（expect1，expect2）未改变，具体问题描述如下：
1.expect1，expect2是s,d1的函数，当fsolve函数按照优化算法找到新的点（s',d1'）时，我的问题里时需要expect1，expect2也变化，这样fsolve的目标函数就也变化了；
; D) F: z; r5 [* m, Q2.根据上述1的描述，这个优化问题是否应该将（expect1，expect2）也纳入优化求解目标中去？这样该如何实现？
, n- ]2 v8 Y) i( f" k' b+ eclc;
& ^$ W1 d0 ]% r! x7 }& |clear all;
%%%  赋值外生给定的参数  %%%
wn=0.3;
wt=1-wn;
0 B/ L" k6 @% `5 `* K$ m  R/ l; ?r=0.05;
  p; t8 V3 N7 ?% K1 h$ Pbeta=1/(1+r);
epislon=0.05;
phi=0.3;
y1=1;
yn2=1;
  E; Q# Q9 K) v- P$ h7 G; wy2mean=1;
y3=1;
" V) h# B" L3 l6 C: Dtheta=0.9;
4 |' a! t& l% e2 a4 y

% 赋值(y2,s,d1)
y2=0.95:0.01:1.05;
/ [( o$ ~" D$ M8 c% j4 is=0.1; %给s,d1的任意初值
$ n! }) _4 U! L7 G9 J9 \7 _3 H7 ?( O- j  ld1=0.1;
7 a7 }8 c& U% Y" [, |% rx0 =[s,d1];

% f  x% P+ P3 i6 b* }$ ~
%在这个for循环里得到均值expect1,expect2，改变(s,d1),均值expect1,expect2也改变
1 ]: q" [( V6 s

. j$ y+ `3 A5 C3 ~# S( X. g* ifor i=1:1:length(y2);
s=0.01;
7 j- A/ h' @% {6 a& Q5 a4 {- ]# ed1=0.1;


! F& ]! m$ Q0 C" f* L%%% 判断是否紧约束 %%%
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)
    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));
    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));
# z$ y1 b' }7 w" |  g   
; n- M, d* J: Xelse
    ct2(i)=((1-s)*y2(i)-d1+(1-s)*y3/(1+r))/(1+beta);
/ y4 q0 e1 ?" {; S    d2(i)=((1-s)*y3-beta*(1+r)*(1-s)*y2(i)+beta*(1+r)*d1)/(1+beta);
( \$ o5 {$ e& S6 k* t3 H' ]8 |- L. Y    ct3(i)=(1-s)*y3-d2(i);

1 g, M5 i* z, C
    mu(i)=0;
; L( |" u0 n5 o7 nend
, W9 ?0 l# p% f# o2 \e1(i)=wt/ct2(i);%每给定一个y2(i)，计算得到一个e1(i)，e2(i)
e2(i)=wt/ct2(i)*y2(i)+mu(i)*phi*y2(i)+beta*wt/ct3(i)*y3;
# t% j4 a# ~  j* V, P7 ^

%%%求均值%%%
3 H* ]% H% }  nexpect1=mean(e1(1,：));
expect2=mean(e2(1,：));
* ?* {2 L* P2 ~end
$ P  H7 x. ?! u" X; w/ `; D
8 |2 x7 D: A/ Z1 m" m
%我的问题中，expect1，expect2实质上是s,d1的函数
%在用fsolve优化求解的时候，每迭代一个新的（s,d1),都应该运行一次上述for这个循环，得到新的expect1、expect2再进行优化求解，
%也即（s,d1)变化时，expect1、expect2也要变化，如此fsolve的目标函数也变化
+ S7 j/ Y4 Q6 j# J: P& E1 r$ R%这样一个问题没有编程思路，请求大家给予指导，如果我的问题没有表述清楚，十分欢迎追问，谢谢。
7 R5 Z) F1 P- L' u* N
$ ]$ X# A8 }8 G$ P8 _, d! ~
" ]& c% L# @6 e7 x; [% b[x,fval,exitflag,output,]=fsolve(@(x)foc(x,expect1,expect2),x0);


function f=foc(x,expect1,expect2)
wn=0.3;
% G' ]5 c- \6 A8 C8 R# Dwt=1-wn;
r=0.05;
beta=1/(1+r);
epislon=0.05;
phi=0.3;
y1=1;
yn2=1;
& X0 m  j4 M) B0 Wy2mean=1;
- e" t; K+ q2 V1 fy3=1;
* ^- z! h) I: m. _( e- u5 V" Vtheta=0.9;
9 a* c* [9 E  H' `! ^) R1 Us=0.01;
- u- P2 S! {( ^- O4 Gd1=0.1;
& g' s6 p0 S5 [. C0 }. bct1=s*(1-theta)*y1+d1/(1+r);% ct1也与s,d1相关；


f(1)=wt/ct1-beta*(1+r)*expect1;
f(2)=wt/ct1*y1*(1-theta)-beta*expect2;
end
: \. `( H* s( l( [  H
$ J$ o& x  S' A8 c# v9 E
5 K  y0 p; ^, L2 d4 K

pTDbn25

帮你顶一下

zzz.dan

1 E! S7 D; E/ e- w你这个写的啥，你foc的函数的变量f1和f2里面x的变量呢，你这个方程写的有很大问题

nkkopd