0

我需要在matlab中求解一个非线性方程。详情请查看代码。

sigma = .25;
beta = 2;
alpha = .1;
eta = .025;

syms x

phi = sqrt(2.*sigma.^2 + beta.^2);
kappa = (beta + phi)./(beta - phi);

w = .5;

D is imported data;

A = @(t) exp((alpha*(beta+phi).*t)./sigma.^2).*...
    ((1-kappa)./(1-kappa.*exp(phi.*t))).^(2.*alpha./sigma.^2);

B = @(t) (beta - phi)./sigma.^2 + 2.*phi./(sigma.^2.*(1-kappa.*exp(phi.*t)));

C = @(t) exp(eta^2.*t^3/6);

G = @(t) (alpha./phi).*(exp(phi.*t)-1).*exp(alpha.*(beta+phi).*t./sigma^2).*...
    ((1-kappa)./(1-kappa.*exp(phi.*t))).^((2.*alpha./sigma^2) + 1);

H = @(t) exp((alpha.*(beta+phi) + phi.*sigma^2)./sigma^2).*...
   ((1-kappa)./(1-kappa.*exp(phi.*t))).^((2.*alpha./sigma^2) + 2);

spread_prot_portion = @(i,j,t) exp(B(t).*x).*D(i,j).*(G(t) + H(t).*x);
spread_fee_portion  = @(i,j,t) A(t).*exp(B(t).*x).*D(i,j);

t = (.5:.5:10);
j = (1 :1 :20);
k = 1:20


solve('cds(1)- w*sum(spread_prot_portion(1,j(k),t(k)))/sum(spread_fee_portion(1,j(k),t(k))) = 0')

注意 cds 是加载的数据。我需要解决 x。这就是为什么它是 syms 中唯一的变量，而所有其他变量都给出了。不幸的是，我不能只写出一个表达式。一些帮助将不胜感激

4

1 回答 1

0

对于这种情况，您需要使用 for 循环来查找总和。认为

syms x;
f = @(i, j) i + j + x;

i然后总结j你需要的所有指数

total = 0;
for i = 1:10
  for j = 1:10
    total = total + f(i, j);
  end
end
display(total) %// total = 100*x + 1100

然后为total = 100你解决写（注意缺少撇号）

xSol = solve(total - 100);

但是请注意，您正试图找到一个非常困难的问题的确切解决方案。使用 fsolve 可能会更好

xSol = fsolve(@fun, estimateOfx);

其中 fun 计算f(i,j)特定值的总和的函数x，即

function total = fun(x)
f = @(i, j) i + j + x;
%// Find total as above
end

这是一个解释fsolve更多的例子。

于 2013-12-02T20:07:20.157 回答