matlab优化求解的问题

Uifhjvv

关于matlab优化求解的问题请教大家；程序报错截图如下，我分析报错的原因是迭代点（s,d1）改变时，（expect1，expect2）未改变，具体问题描述如下：
# T) m6 K9 R( Y$ n$ c! M1.expect1，expect2是s,d1的函数，当fsolve函数按照优化算法找到新的点（s',d1'）时，我的问题里时需要expect1，expect2也变化，这样fsolve的目标函数就也变化了；
2.根据上述1的描述，这个优化问题是否应该将（expect1，expect2）也纳入优化求解目标中去？这样该如何实现？
clc;
/ {' E1 D1 [# t% K* p: s& Tclear all;
%%%  赋值外生给定的参数  %%%
wn=0.3;
wt=1-wn;
2 `4 D8 g* q( I7 P8 n/ sr=0.05;
beta=1/(1+r);
epislon=0.05;
phi=0.3;
y1=1;
7 @  o. J  g  _" K7 gyn2=1;
. ^0 v; ~$ Z0 m: r4 Z" |y2mean=1;
2 T4 d3 d* T6 e- {% Q. _y3=1;
$ t- I- d- U2 h  T, g, U5 m# ?theta=0.9;
. C$ m, i) S+ S! G6 S

% 赋值(y2,s,d1)
% y6 p* E0 S+ S1 D. J/ Xy2=0.95:0.01:1.05;
( b$ Q3 x: c7 Fs=0.1; %给s,d1的任意初值
& W0 `( I7 E9 x* r5 m9 p3 M3 wd1=0.1;
0 S# |6 Z' c; M2 z7 yx0 =[s,d1];


2 r- ^9 Y  @/ L, P( ]* Y. M- W) k%在这个for循环里得到均值expect1,expect2，改变(s,d1),均值expect1,expect2也改变


3 q! R7 m( U7 H  a' Mfor i=1:1:length(y2);
s=0.01;
d1=0.1;

2 {1 J7 E2 [4 @8 P& L
/ N# l* L% a- M) W9 ]0 }%%% 判断是否紧约束 %%%
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);
9 o( J* }" P3 ~7 x3 sif check(i)>standard(i)
    ct2(i)=((1-s)*y2(i)-d1+phi*(1-s)*y2(i))/(1-phi*wn/wt);
, v% ]* B$ N, H* S- F( z    d2(i)=phi*(1+r)*(wn/wt*ct2(i)+(1-s)*y2(i));
' {' i  {1 h8 x    ct3(i)=(1-s)*y3-d2(i);
' @% ]) E4 h; z  u& ?    mu(i)=wt/((1-s)*y2(i)-d1+d2(i)/(1+r))-beta*wt*(1+r)/((1-s)*y3-d2(i));
) X" V, _( a7 h- O   
3 ^$ _( [! M6 k7 Velse
1 H+ |, p3 c* L# z; c' I5 ~    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);
    ct3(i)=(1-s)*y3-d2(i);


    mu(i)=0;
end
e1(i)=wt/ct2(i);%每给定一个y2(i)，计算得到一个e1(i)，e2(i)
% R4 b: k4 t2 a. re2(i)=wt/ct2(i)*y2(i)+mu(i)*phi*y2(i)+beta*wt/ct3(i)*y3;

- ^5 E, M) W, I( N/ \: y
4 A* K" |! y0 `* G%%%求均值%%%
expect1=mean(e1(1,：));
expect2=mean(e2(1,：));
end

2 R# ^, M1 K* a5 ?: t" C% i1 `
* w+ s: D/ Q& G% I5 Y; H%我的问题中，expect1，expect2实质上是s,d1的函数
%在用fsolve优化求解的时候，每迭代一个新的（s,d1),都应该运行一次上述for这个循环，得到新的expect1、expect2再进行优化求解，
%也即（s,d1)变化时，expect1、expect2也要变化，如此fsolve的目标函数也变化
%这样一个问题没有编程思路，请求大家给予指导，如果我的问题没有表述清楚，十分欢迎追问，谢谢。
' ?- x! ]4 v7 F) q

[x,fval,exitflag,output,]=fsolve(@(x)foc(x,expect1,expect2),x0);
& P8 X, y2 T) z! p4 g/ c$ C

% y' S4 m/ s, h) Ifunction f=foc(x,expect1,expect2)
$ ~1 X) N8 k" ]' qwn=0.3;
% v) J8 i4 _# E+ @& x2 X2 S* |6 G. fwt=1-wn;
r=0.05;
0 T0 @; d$ i, @4 v& h7 Ybeta=1/(1+r);
) _  g4 F% [% b+ ~; g/ @epislon=0.05;
. V8 W" J8 X6 ~- j( H9 Uphi=0.3;
y1=1;
yn2=1;
" t3 U" r, Z) T1 ^5 \4 my2mean=1;
2 U. S" b# c7 J0 e+ Q  k% Qy3=1;
theta=0.9;
s=0.01;
d1=0.1;
4 u3 n) N8 R5 cct1=s*(1-theta)*y1+d1/(1+r);% ct1也与s,d1相关；
  {+ d0 Q! d; A/ j

f(1)=wt/ct1-beta*(1+r)*expect1;
f(2)=wt/ct1*y1*(1-theta)-beta*expect2;
end
5 I+ S3 |: X. U0 ?  x6 [

pTDbn25

帮你顶一下

zzz.dan

你这个写的啥，你foc的函数的变量f1和f2里面x的变量呢，你这个方程写的有很大问题

nkkopd