[問題]解方程式的方法(用solve??)已回收
UQ=0.03 ; %表未來訂單成長率:這邊可能要再調整成一些變數的組合權重
Uq=0.02 ; %表現在產能成長率:同上要調整
SQ=0.2 ; %表未來訂單波動率
Sq=0.1 ; %表現在產能波動率
lo=0.5 ; %表相關係數
P=0.3 ; %表示單位價格
vc=0.2 ; %表示變動成本
q0=1 ; %現在產能在時間0的情況下
rf=0.01 ; %無風險利率
beta=0.3 ; %CAPM下的beta
rm=0.04 ; %市場報酬
R=rf+beta*(rm-rf) ; %一家公司的要求獲利率
D=UQ-Uq-0.5*(SQ^2-Sq^2) ;
V=SQ^2-2*lo*SQ*Sq+Sq^2 ; %V已經表示我上面寫的v^2
a=(-D+sqrt(D^2+2*R*V))/V;
b=(-D-sqrt(D^2+2*R*V))/V;
Profit=(P-vc)*q0 ; %假定起初的報酬
g=Profit/(R-0.5*V-D) ; %組合變數讓式子看起來較簡單
K0=0.1 ; %為撤廠的成本
K1=0.8 ; %為建廠的成本
%以上是變數的數值解,也可以改。
syms uT lT;
F=
(g*((b*(lT^b*uT^a-uT^(a+b)*lT)-a*(uT^(1+b)*lT^a-lT*uT^(a+b)))/(a*b*(uT^b*lT^a-lT^b*uT^a))-uT))+K1
;
G=
(g*((b*(lT^(a+b)*uT-uT^b*lT^(1+a))-a*(uT*lT^(a+b)-lT^(1+b)*uT^a))/(a*b*(uT^b*lT^a-lT^b*uT^a))-lT))+K0
;
[llT,uuT]=solve(F,G);
Warning: Warning, solutions may have been lost
Warning: Explicit solution could not be found.
> In solve at 140
In sym.solve at 49
%主要是這個方程式解解不出來
%希望大大幫忙看一下,問題出在哪。謝謝。
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.42.197.191
※ 編輯: hagry 來自: 114.42.197.191 (02/24 12:23)