Question

我需要解决一个二维方程组。我有一个列向量SpreadCLPbp，其中N=length(SpreadCLPbp)返回500。

我想解决以下等式：

0 = sum_{i=1}^{N}{(-1/mu2 +exp(mu2*delta)/mu2)*(SpreadCLPbp(k)-(-mu1/mu2+(mu1/mu2+SpreadCLPbp(k))*exp(mu2*delta)))/(mu1*sigma/(2*mu2^2)-sigma^2/mu2*(mu1/mu2+SpreadCLPbp(k))*exp(mu2*delta)+sigma^2/mu2*(mu1/(2*mu2)+SpreadCLPbp(k))*exp(2*mu2*delta))}

0 = sum_{k=1}^{N}{(mu1/mu2+delta*exp(mu2*delta)*(mu1/mu2+SpreadCLPbp(k))-exp(mu2*delta)*mu1/mu2)*(SpreadCLPbp(k)-(-mu1/mu2+(mu1/mu2+SpreadCLPbp(k))*exp(mu2*delta)))/(mu1*sigma/(2*mu2^2)-sigma^2/mu2*(mu1/mu2+SpreadCLPbp(k))*exp(mu2*delta)+sigma^2/mu2*(mu1/(2*mu2)+SpreadCLPbp(k))*exp(2*mu2*delta))}

我尝试了以下

syms k mu1 mu2;

S1 = symsum((-1/mu2 +exp(mu2*delta)/mu2)*(SpreadCLPbp(k)-(-mu1/mu2+(mu1/mu2+SpreadCLPbp(k))*exp(mu2*delta)))/(mu1*sigma/(2*mu2^2)-sigma^2/mu2*(mu1/mu2+SpreadCLPbp(k))*exp(mu2*delta)+sigma^2/mu2*(mu1/(2*mu2)+SpreadCLPbp(k))*exp(2*mu2*delta)),1,N);

S2 = symsum((mu1/mu2+delta*exp(mu2*delta)*(mu1/mu2+SpreadCLPbp(k))-exp(mu2*delta)*mu1/mu2)*(SpreadCLPbp(k)-(-mu1/mu2+(mu1/mu2+SpreadCLPbp(k))*exp(mu2*delta)))/(mu1*sigma/(2*mu2^2)-sigma^2/mu2*(mu1/mu2+SpreadCLPbp(k))*exp(mu2*delta)+sigma^2/mu2*(mu1/(2*mu2)+SpreadCLPbp(k))*exp(2*mu2*delta)),1,N);

S = solve([S1 == 0, S2 == 0],[mu1, mu2])

但是我收到以下错误：

使用sym / subsindex时出错（第732行）
      索引或函数定义无效。定义功能时，请确保执行   参数是符号变量，函数体是SYM
  表达。建立索引时，输入必须是数字，逻辑或＆＃39;：＆＃39;。

analisisdatoscir出错（第42行）
  S1 = symsum（（ - 1 / mu2 + exp（mu2 * delta）/ mu2）*（SpreadCLPbp（k） - （ - mu1 / mu2 +（mu1 / mu2 + SpreadCLPbp（k））* exp（mu2 * delta））） /（MU1 *西格玛/（2 * MU2 ^ 2）^ -sigma 2 / MU2 *（MU1 / MU2 + SpreadCLPbp（K））* EXP（MU2 *增量）+Σ-^ 2 / MU2 *（MU1 /（2 * MU2）+ SpreadCLPbp（K））* EXP（2 * * MU2三角洲）），1，N）;

如果有人能帮助我，我将非常感激。

如何用内部的symsum求解方程

0 个答案: