[未答复] 方程组求解

[复制链接]
bobo85 发表于 2021-6-1 17:24:34
10 财富积分

本帖最后由 bobo85 于 2021-6-1 17:26 编辑

ee=0;
le=1;
Mu=[3.5 3.5];
Sigma=[0.3 0.3];
LB=[0 5];
UB=[5 10];
g=@(x,t)  x(:,1).^2.*x(:,2)-5*x(:,1).*t(:,1)+(x(:,2)+1).*(t(:,1).^2)+x(:,1).*t(:,1).*t(:,2)-25;
gt{1}=@(x,t) -5.*x(:,1)+2.*(x(:,2)+1).*t(:,1)+x(:,1).*t(:,2);
gt{2}=@(x,t) x(:,1).*t(:,1);

b0=@(t) g(Mu,t);
a{1}=@(x,t) 2*x(:,1).*x(:,2)-5*t(:,1)+t(:,1).*t(:,2);
a{2}=@(x,t) x(:,1).^2+t(:,1).^2;
b{1}=@(t) a{1}(Mu,t).*Sigma(1);
b{2}=@(t) a{2}(Mu,t).*Sigma(2);

Db0{1}=@(t) gt{1}(Mu,t);
Db0{2}=@(t) gt{2}(Mu,t);
Da{1,1}=@(x,t) -5+t(:,2);
Da{1,2}=@(x,t) t(:,1);
Da{2,1}=@(x,t) 2.*t(:,1);
Da{2,2}=@(x,t) 0;
Db{1,1}=@(t) Da{1,1}(Mu,t).*Sigma(1);
Db{1,2}=@(t) Da{1,2}(Mu,t).*Sigma(1);
Db{2,1}=@(t) Da{2,1}(Mu,t).*Sigma(2);
Db{2,2}=@(t) Da{2,2}(Mu,t).*Sigma(2);

Dbb{1}=@(t1,t2) 0;
bb{1}=@(t1,t2) 0;
Dbb{2}=@(t1,t2) 0;
bb{2}=@(t1,t2) 0;
for i=1:length(LB)
    for j=1:length(b)
        Dbb{i}=@(t1,t2) Dbb(t1,t2)+Db{j,i}(t1).*b{j}(t2);
        bb{i}=@(t1,t2) bb(t1,t2)+b{j}(t1).*b{j}(t2);
    end
end
f{1}=@(t) Db0{1}(t)-(ee+b0{1}(t)).*Dbb{1}(t,t)./bb{1}(t,t);
f{2}=@(t) Db0{2}(t)-(ee+b0{2}(t)).*Dbb{2}(t,t)./bb{2}(t,t);

想要联立f{1}和f{2}求t,其中t是二维的,该怎么做呢?求大佬指教!

您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
热门教程
站长推荐
快速回复 返回顶部 返回列表